JP6852823B1 - Motor control device - Google Patents

Abstract

PROBLEM TO BE SOLVED: To improve the vibration damping effect of a motor. In a motor control device 100, a control switching determination unit 15 determines whether or not a control region of a motor M is in a voltage saturation region, and a voltage command value generator 14 determines a speed command value and a motor speed. The voltage command value of the motor is generated based on the above, and when the control switching determination unit 15 determines that the control region of the motor M is in the voltage saturation region, the total torque command value and the maximum voltage that can be output to the motor M. The voltage vector angle of the output voltage applied to the motor is obtained from the limit value of, and the voltage command value is generated based on this voltage vector angle. [Selection diagram] Fig. 2

Description

The present disclosure relates to a motor control device.

In the compressor used in the air conditioner, the load torque fluctuates periodically during one rotation of the rotor of the motor that drives the compressor. This periodic load torque fluctuation occurs due to the pressure change of the refrigerant gas between the suction, compression, and discharge strokes, and is sometimes referred to as the fluctuation of the rotation speed of the motor (hereinafter, simply referred to as "speed fluctuation"). ) Causes motor vibration. When a compressor in which such load torque fluctuations occur is used, "torque control (periodic disturbance suppression control)" is performed in order to suppress fluctuations in the rotational speed of the motor.

Usually, the vibration of the motor appears remarkably in the low rotation region (for example, the normal control region where the maximum torque / current control of the motor is performed). However, depending on the specifications of the motor and the load conditions, vibration occurs even in a high rotation region (for example, a voltage saturation region where weak magnetic flux control is performed), and the peak current of the motor increases due to the generation of vibration. Further, when the peak current of the motor increases, the protection function of the inverter may operate to prevent demagnetization of the motor, and the motor may stop.

Therefore, as torque control in the voltage saturation region where weak magnetic flux control is performed, the voltage vector angle indicating the phase of the output voltage fluctuates in synchronization with the speed fluctuation while limiting the output voltage below the DC voltage that can be output by the inverter. A technique for performing torque control in a voltage saturation region has been proposed.

JP-A-2017-158414 JP-A-2017-158415

However, when torque control is performed in the voltage saturation region by fluctuating the voltage vector angle of the output voltage with respect to the speed fluctuation, it has been necessary to take man-hours to adjust the phase of the voltage vector angle fluctuation by tuning. Therefore, even if the phase of the voltage vector angle fluctuation is tuned under certain conditions to obtain a vibration damping effect, the voltage for generating the optimum output torque if the specifications of the inverter or motor, load conditions, etc. change. The damping effect may not be fully exerted without the vector angle fluctuating.

Therefore, the present disclosure proposes a technique capable of improving the vibration damping effect of the motor.

The motor control device of the present disclosure includes a voltage command value generator and a control switching determination unit. The control switching determination unit determines whether or not the control region of the motor is in the voltage saturation region. The voltage command value generator generates the voltage command value of the motor from the torque command value based on the speed command value and the speed of the motor, and when the control switching determination unit determines that the control region is in the voltage saturation region. Generates a voltage command value based on the torque command value, the limit value of the maximum voltage that can be output to the motor, and the voltage vector angle of the output voltage applied to the motor.

According to the present disclosure, it is possible to improve the vibration damping effect of the motor.

FIG. 1A is a diagram provided for explaining an operation example of the motor control device according to the first embodiment of the present disclosure. FIG. 1B is a diagram provided for explaining an operation example of the motor control device according to the first embodiment of the present disclosure. FIG. 2 is a diagram showing a configuration example of the motor control device according to the first embodiment of the present disclosure. FIG. 3 is a diagram showing a configuration example of the control switching determination unit according to the first embodiment of the present disclosure. FIG. 4 is a diagram showing a configuration example of the correction torque generator according to the first embodiment of the present disclosure. FIG. 5 is a diagram showing a configuration example of the voltage saturation region voltage command value generator according to the first embodiment of the present disclosure. FIG. 6 is a diagram showing a configuration example of the output voltage limit command value generator according to the first embodiment of the present disclosure. FIG. 7A is a diagram provided for explaining an operation example of the current command value calculator according to the first embodiment of the present disclosure. FIG. 7B is a diagram provided for explaining an operation example of the current command value calculator according to the first embodiment of the present disclosure. FIG. 8 is a diagram provided for explaining an operation example of the voltage vector angle calculator according to the first embodiment of the present disclosure. FIG. 9 is a diagram provided for explaining an operation example of the MTPI voltage amplitude limiting processor according to the first embodiment of the present disclosure. FIG. 10 is a diagram showing an example of the output voltage waveform of the first embodiment of the present disclosure. FIG. 11 is a diagram showing a configuration example of the normal control region voltage command value generator according to the first embodiment of the present disclosure. FIG. 12 is a diagram showing a configuration example of the current error correction value generator according to the first embodiment of the present disclosure. FIG. 13 is a diagram showing a configuration example of the voltage saturation region voltage command value generator according to the second embodiment of the present disclosure.

Hereinafter, examples of the present disclosure will be described with reference to the drawings. In the following examples, the same configurations are designated by the same reference numerals.

In the present disclosure, a motor control device that controls the torque of a permanent magnet synchronous motor (PMSM (Permanent Magnet Synchronous Motor)) that drives a compressor having periodic load torque fluctuations by position sensorless vector control, for example, an air conditioner or a low temperature. A motor control device used as a storage device or the like will be described as an example. However, the disclosed technique is widely applicable to a motor control device that controls torque of a motor that drives a load having periodic load torque fluctuations.

[Example 1]
<Operation of motor control device>
1A and 1B are diagrams for explaining an operation example of the motor control device according to the first embodiment of the present disclosure.

The constant induced voltage ellipse shown in FIGS. 1A and 1B (a part of the ellipse is shown in FIGS. 1A and 1B) is a current vector locus at which the induced voltage Vo of the motor becomes equal, and it is determined that the estimated electrical angle angular velocity ωe becomes large. The diameter of the induced voltage ellipse becomes smaller. The constant induced voltage ellipse shown in FIG. 1A shows the current vector locus when the estimated electrical angle angular velocity ωe is constant. Further, the constant induced voltage ellipse shown in FIG. 1B shows a current vector locus when the electric angle estimated angular velocity ωe fluctuates due to a load torque fluctuation. In FIG. 1B, a constant-induced voltage ellipse at the maximum value of the estimated electric angular velocity ωe, a constant-induced voltage ellipse at the minimum value of the estimated electric angular velocity ωe, and a constant-induced voltage ellipse at the average value of the estimated electric angular velocity ωe are shown. Shown.

^{As shown in FIG. 1B, the motor control device of the present disclosure has a constant torque curve T *} (= To) that fluctuates by ± ΔT due to torque control when performing torque control in the voltage saturation region where the weakening magnetic flux control of the motor is performed. ^{The q-axis current command value Iq *} and the d-axis current command value Id ^* are calculated based on the intersection of * ± ΔT) and the fluctuating constant-induced voltage ellipse. For example, the motor control device of the present disclosure has a constant torque curve which is a locus of a current ^{in which the total torque command value T * obtained} by adding the fluctuation torque command value ΔT which is a correction torque to ^{the average torque command value To * is constant and an output. Based on the intersection of the induced voltage command value V0 *} for setting the voltage amplitude Va (the amplitude of the voltage command value) to the desired amplitude and the constant induced voltage ellipse which is the locus of the current at which the estimated electrical angle velocity ωe is constant. The shaft current command value Id ^* and the q-axis current command value Iq ^* are calculated. The constant induced voltage ellipse and constant torque curve are not uniquely determined based on motor parameters such as reactance, but change from moment to moment depending on the operating state of the motor.

<Structure of motor control device>
FIG. 2 is a diagram showing a configuration example of the motor control device according to the first embodiment of the present disclosure. In FIG. 2, the motor control device 100 includes subtractors 11 and 38, a speed controller 12, an adder 13, a voltage command value generator 14, a control switching determination unit 15, and d−q / u, v. , W converter 23, PWM (Pulse Width Modulation) modulator 24, and IPM (Intelligent Power Module) 25. The IPM 25 is connected to the motor M. PMSM is an example of the motor M.

Further, the motor control device 100 includes a shunt resistor 26, current sensors 27a and 27b, and a 3φ current calculator 28. The motor control device 100 may have either a shunt resistor 26 or current sensors 27a and 27b.

Further, the motor control device 100 includes a u, v, w / dq converter 29, an axis error calculator 30, a PLL (Phase Locked Loop) controller 31, a position estimator 32, and 1 / Pn processing. It has a device 33 and a correction torque generator 34.

The voltage command value generator 14 includes a normal control region voltage command value generator 14a, a voltage saturation region voltage command value generator 14b, a switch SW1, and a switch SW2. The switch SW1 has contacts 14c-1, 14c-2, 14c-3. The switch SW2 has contacts 14c-4, 14c-5, 14c-6.

The subtractor 11 is the current estimated angular velocity output from the 1 / Pn processor 33 from ^{the mechanical angular velocity command value ωm *} input to the motor control device 100 from the outside of the motor control device 100 (for example, a higher-level controller). The angular velocity error Δω is calculated by subtracting a certain mechanical angular velocity estimated angular velocity ωm, and the calculated angular velocity error Δω is output to the speed controller 12.

^{The speed controller 12 generates an average torque command value To *} so that the angular velocity error Δω input from the subtractor 11 approaches zero, and outputs the generated average torque command value To ^* to the adder 13.

The adder 13 calculates the ^{total torque command value T *} by adding the average torque command value To ^* output from the speed controller 12 and the variable torque command value ΔT output from the correction torque generator 34. , The calculated total torque command value T ^* is output to the voltage command value generator 14.

The voltage command value generator 14 has a d-axis voltage command value Vd ^* and a q-axis voltage command value Vq ^{* based on the total torque command value T *} output from the adder 13 in each of the normal control region and the voltage saturation region. Is generated, and the generated d-axis voltage command value Vd ^* and q-axis voltage command value Vq ^* are output. The voltage saturation region is a region in which the output voltage amplitude Va is saturated in the high rotation region of the motor M and the weakening magnetic flux control is performed. The normal control region is a region other than the voltage saturation region, and is a region in which the motor M is controlled by varying the output voltage. In the normal control region, maximum torque / current control and the like are performed.

When the control signal CONTROL_TYPE: A (normal control) is output from the control switching determination unit 15, the voltage command value generator 14 connects the contact 14c-1 and the contact 14c-3 of the switch SW1 and switches. By connecting the contact 14c-4 of SW2 and the contact 14c-6, the d-axis voltage command value Vd ^* and the q-axis voltage command value Vq ^* generated by the normal control region voltage command value generator 14a are dq / Output to u, v, w converter 23. On the other hand, when the control signal CONTROL_TYPE: B (voltage saturation control) is output from the control switching determination unit 15, the voltage command value generator 14 connects the contact 14c-2 and the contact 14c-3 of the switch SW1. At the same time, the contact 14c-5 and the contact 14c-6 of the switch SW2 are connected to generate the d-axis voltage command value Vd ^* and the q-axis voltage command value Vq ^{* generated by the voltage saturation region voltage command value generator 14b.} Output to the d−q / u, v, w converter 23.

In the control switching determination unit 15, the current control area of the motor M is the normal control area and the voltage saturation area based on the output voltage limit value Vdq_limit, the d-axis voltage command value Vd ^*, and the q-axis voltage command value Vq ^*. Which of the above is determined. Then, when the control switching determination unit 15 determines that the current control area of the motor M is the normal control area, the control switching determination unit 15 outputs the control signal CONTROL_TYPE: A (normal control) to the voltage command value generator 14 and outputs the control signal CONTROL_TYPE: A (normal control) to the motor M. When it is determined that the current control region of the above is the voltage saturation region, the control signal CONTROL_TYPE: B (voltage saturation control) is output to the voltage command value generator 14. The output voltage limit value Vdq_limit is a DC voltage Vdc supplied to the IPM25 from the outside of the IPM25 (for example, a power converter (not shown)) converted into a voltage value in the dq rotation coordinate axis system which is a control system.

The d−q / u, v, w converter 23 outputs the two-phase d-axis voltage command value Vd ^* and the q-axis voltage command value Vq ^* output from the voltage command value generator 14 from the position estimator 32. Based on the current rotor position, the electric angle phase (dq axis phase) θe, the three-phase U-phase output voltage command value Vu ^* , V-phase output voltage command value Vv ^*, and W-phase output voltage command value Vw ^* Convert to. Then, the d−q / u, v, w converter 23 outputs the U-phase output voltage command value Vu ^* , the V-phase output voltage command value Vv ^*, and the W-phase output voltage command value Vw ^* to the PWM modulator 24. ..

The PWM modulator 24 generates a 6-phase PWM signal based on the U-phase output voltage command value Vu ^* , the V-phase output voltage command value Vv ^* , the W-phase output voltage command value Vw ^{*, and the PWM carrier signal.} , The generated 6-phase PWM signal is output to the IPM25.

The IPM 25 applies the DC voltage Vdc supplied from the outside of the IPM 25 to each of the U phase, V phase, and W phase of the motor M by converting the DC voltage Vdc supplied from the outside of the IPM 25 based on the 6-phase PWM signal output from the PWM modulator 24. AC voltage is generated, and each AC voltage is applied to the U phase, V phase, and W phase of the motor 10.

When the bus current is detected by the 1-shunt method using the shunt resistance 26, the 3φ current calculator 28 uses the 6-phase PWM switching information output from the PWM modulator 24 and the detected bus current to determine the motor. The U-phase current value Iu, the V-phase current value Iv, and the W-phase current value Iw of M are calculated. Alternatively, when the U-phase current and the V-phase current are detected by the current sensors 27a and 27b, the 3φ current calculator 28 calculates the remaining W-phase current value Iw based on Kirchhoff's law of “Iu + Iv + Iw = 0”. .. The 3φ current calculator 28 outputs the calculated phase current values Iu, Iv, Iw of each phase to the u, v, w / dq converter 29.

The u, v, w / dq converter 29 is a three-phase U phase output from the 3φ current calculator 28 based on the electric angle phase θe indicating the current rotor position output from the position estimator 32. The current value Iu, the V-phase current value Iv, and the W-phase current value Iw are converted into the two-phase d-axis current Id and the q-axis current Iq. Then, the u, v, w / dq converter 29 outputs the d-axis current Id and the q-axis current Iq to the voltage command value generator 14 and the axis error calculator 30.

The axis error calculator 30 outputs the d-axis voltage command value Vd ^* and the q-axis voltage command value Vq ^* output from the voltage command value generator 14 and the u, v, w / d-q converter 29. The axis error Δθ (difference between the estimated rotation axis and the actual rotation axis) is calculated using the d-axis current Id and the q-axis current Iq. Then, the axis error calculator 30 outputs the calculated axis error Δθ to the PLL controller 31.

The PLL controller 31 calculates the electric angle estimated angular velocity ωe, which is the current estimated angular velocity, based on the axis error Δθ output from the axis error calculator 30, and uses the calculated electric angle estimated angular velocity ωe as the position estimators 32 and 1. Output to / Pn processor 33.

The position estimator 32 estimates the electric angle phase θe and the mechanical angle phase θm based on the electric angle estimated angular velocity ωe output from the PLL controller 31. Then, the position estimator 32 outputs the estimated electric angle phase θe to the d−q / u, v, w converter 23 and the u, v, w / d−q converter 29, and estimates the mechanical angle phase θm. Is output to the voltage command value generator 14 and the correction torque generator 34.

The 1 / Pn processor 33 calculates the machine angle estimated angular velocity ωm by dividing the electric angle estimated angular velocity ωe output from the PLL controller 31 by the pole logarithm Pn of the motor M, and calculates the machine angle estimated angular velocity ωm. Output to subtractors 11 and 38.

The subtractor 38 calculates the machine angle estimated angular velocity fluctuation Δωm by subtracting the machine angle ^{velocity command value ωm *} from the machine angle estimated angular velocity ωm output from the 1 / Pn processor 33, and calculates the machine angle estimated angular velocity fluctuation Δωm. Is output to the correction torque generator 34.

^{The correction torque generator 34 has a speed fluctuation allowable value | Δωm | *} , which is a speed fluctuation range in which the vibration of the motor M can be tolerated, a mechanical angle estimated angular velocity fluctuation Δωm output from the subtractor 38, and an output from the position estimator 32. on the basis of the mechanical angle phase .theta.m, periodic rate is the variation of mechanical angle estimated angular velocity change Derutaomegaemu speed variation tolerance | Δωm | ^* generates fluctuation torque command value ΔT for suppressing below, the resulting variations The torque command value ΔT is output to the adder 13. The permissible speed fluctuation value | Δωm | ^* is stored in the motor control device 100. Further, the mechanical angle estimated angular velocity fluctuation (velocity fluctuation) Δωm is different only in the positive and negative signs from the value of the above angular velocity error Δω.

<Structure of control switching judgment unit>
FIG. 3 is a diagram showing a configuration example of the control switching determination unit according to the first embodiment of the present disclosure. In FIG. 3, the control switching determination unit 15 has a voltage amplitude calculator 15a and a control switching determination device 15b, and determines whether the current control region of the motor is a normal control region or a voltage saturation region as follows. To do.

The voltage amplitude calculator 15a calculates the output voltage amplitude Va according to the equation (1) based on the d-axis voltage command value Vd ^* and the q-axis voltage command value Vq ^{* output from the voltage command value generator 14.}

The control switching determination device 15b compares the peak value of the output voltage amplitude Va calculated by the voltage amplitude calculator 15a with the output voltage limit value Vdq_limit.

When the peak value of the output voltage amplitude Va is less than the output voltage limit value Vdq_limit, the control switching determination device 15b determines that the current control area of the motor M is the normal control area, and sets the control signal CONTROL_TYPE: A to a voltage. Output to the command value generator 14.

On the other hand, when the peak value of the output voltage amplitude Va is equal to or greater than the output voltage limit value Vdq_limit, the control switching determination device 15b determines that the current control region of the motor M is the voltage saturation region, and the control signal CONTROL_TYPE: B is output to the voltage command value generator 14.

<Structure of correction torque generator>
FIG. 4 is a diagram showing a configuration example of the correction torque generator according to the first embodiment of the present disclosure. In FIG. 4, the correction torque generator 34 includes a speed fluctuation component separator 34a, a speed fluctuation amplitude calculator 34b, a subtractor 34c, a correction torque amplitude calculator 34d, a speed fluctuation phase corrector 34e, and an orthogonal component. It has a separator 34f and a correction torque demodulator 34g.

In the correction torque generator 34, the fluctuation torque command value (correction torque) ΔT so that the speed fluctuation amplitude | Δωm | is within the speed fluctuation allowable value | Δωm | ^{* within the range where the vibration of the motor M does not pose a problem in actual use.} (Corrected torque amplitude | ΔT |) and phase are adjusted for each mechanical angle period.

The velocity fluctuation component separator 34a sets the mechanical angle estimated angular velocity fluctuation Δωm to the two Fourier coefficients ωsin (sin component) and ωcos (sin component) which are the fundamental wave components of Δωm according to the equations (2.1) and (2.2). cos component) is separated. By calculating the Fourier coefficient of the fundamental wave component of the machine angle estimated angular velocity fluctuation Δωm for each machine angle period, the harmonic component of the machine angle estimated angular velocity fluctuation Δωm is eliminated and the fundamental wave component of the machine angle estimated angular velocity fluctuation Δωm is accurate. It can be extracted well. ωsin and ωcos are values that are updated every mechanical angle period.

The velocity fluctuation amplitude calculator 34b calculates the velocity fluctuation amplitude | Δωm | according to the equation (3) based on the Fourier coefficients ωsin and ωcos. Since ωsin and ωcos are values that are updated every machine angle period, the velocity fluctuation amplitude | Δωm | is also updated every machine angle period.

The subtractor 34c calculates the velocity fluctuation error | Δωm | err by subtracting the ^{velocity fluctuation tolerance | Δωm | *} from the velocity fluctuation amplitude | Δωm | output from the velocity fluctuation amplitude calculator 34b. The permissible speed fluctuation value | Δωm | ^* defines the permissible speed fluctuation amplitude | Δωm | in the range in which the vibration of the motor M can be tolerated.

The corrected torque amplitude calculator 34d adjusts the corrected torque amplitude | ΔT | for each mechanical angle period according to the error between the speed fluctuation amplitude | Δωm | and the speed fluctuation allowable value | Δωm | ^*. For example, the correction torque amplitude calculator 34d multiplies the speed fluctuation error | Δωm | err, which is the error between the speed fluctuation amplitude | Δωm | and the speed fluctuation tolerance | Δωm | ^{*, by the correction gain k according to the equation (4).} , The correction torque amplitude | ΔT | is calculated by adding the multiplication result and | ΔT | _old. | ΔT | _old in the equation (4) is the corrected torque amplitude | ΔT | in the previous mechanical angle period. By setting the correction gain k appropriately, the speed fluctuation | Δω | can be hunted at ^{the boundary of the speed fluctuation tolerance | Δωm | *} , and the speed fluctuation | Δω | can be the speed fluctuation tolerance | It becomes ^{larger than Δωm | * and the} occurrence of vibration can be suppressed.

The velocity fluctuation phase corrector 34e corrects the phase of the mechanical angle estimated angular velocity fluctuation Δωm acquired for each mechanical angle period. For example, the velocity fluctuation phase corrector 34e multiplies each of the Fourier coefficients ωsin and ωcos by the correction gain k according to the equations (5.1) and (5.2), and adds ωsin_i_old and ωcos_i_old to the respective multiplication results. To do. The ωsin_i_old in the equation (5.1) is ωsin_i in the previous mechanical angle period, and the ωcos_i_old in the equation (5.2) is ωcos_i in the previous mechanical angle period. Then, the speed fluctuation phase corrector 34e calculates the inverse tangent (Arctangent) of ωsin_i and ωcos_i as the speed fluctuation correction phase φωi according to the equation (5.3). This speed fluctuation correction phase φωi serves as a reference for the phase when torque control is performed, and the phase retarded by π / 2 with respect to this reference becomes the phase (correction torque phase) of the fluctuation torque command value ΔT.

Based on the corrected torque amplitude | ΔT | and the speed fluctuation correction phase φωi, the orthogonal component separator 34f and the sin component (ωsin_i) of the speed fluctuation correction phase φωi according to the equations (6.1) and (6.2). The cos component (ωcos_i) is calculated. This process also has a role of preventing divergence at the time of phase correction by the calculation of the equations (5.1) and (5.2).

The correction torque demodulator 34g sets the fluctuation torque command value ΔT according to the equations (7.1) and (7.2) based on the sin component (ωsin_i) and the cos component (ωcos_i) of the speed fluctuation correction phase φωi. calculate. By this processing, the speed fluctuation correction phase φωi is converted into a correction torque phase retarded by π / 2, and an instantaneous value of the fluctuation torque command value ΔT at the mechanical angle phase θm is generated.

The correction torque demodulator 34g may calculate the instantaneous value of the variable torque command value ΔT according to the equation (8) instead of the equation (7.1) and the equation (7.2).

Then, the adder 13 calculates the ^{total torque command value T *} by adding the variable torque command value ΔT to ^{the average torque command value To *} output from the speed controller 12 according to the equation (9).

<Configuration of voltage saturation region voltage command value generator>
FIG. 5 is a diagram showing a configuration example of the voltage saturation region voltage command value generator according to the first embodiment of the present disclosure. In FIG. 5, the voltage saturation region voltage command value generator 14b includes an output voltage limit command value generator 14b1, an induced voltage command value calculator 14b2, a current command value calculator 14b3, and a temporary voltage command value calculator 14b4. , A voltage vector angle calculator 14b5 and a voltage command value calculator 14b6.

FIG. 6 is a diagram showing a configuration example of the output voltage limit command value generator according to the first embodiment of the present disclosure. In FIG. 6, the output voltage limit command value generator 14b1 includes an MTPI current command value calculator 14b1-1, an MTPI voltage command value calculator 14b1-2, an MTPI voltage amplitude calculator 14b1-3, and an average output voltage generation. It has a device 14b1-4, an MTPI voltage fluctuation component extractor 14b1-5, an adder 14b1-6, 14b1-8, and an MTPI voltage amplitude limiting processor 14b1-7.

The output voltage limit command value generator 14b1 outputs based on the total torque command value T ^* , the electric angle estimated angular velocity ωe, the output voltage limit value Vdq_limit, the d-axis current Id, the q-axis current Iq, and the mechanical angle phase θm. Generates the voltage limit command value Va ^*. This output voltage limit command value Va ^* adjusts the fluctuation amplitude of the output voltage within the range up to the output voltage limit value Vdq_limit, and sets the fluctuation phase of the output voltage as the fluctuation phase of the output voltage in the normal control region (MTPI control region). It is a voltage to match.

In FIG. 6, the MTPI current command value calculator 14b1-1 is the MTPI which is the intersection of the constant torque curve and the MTPI curve (maximum torque / current control curve) which are the loci of the current at which the ^{total torque command value T * is constant.} The assumed d-axis current command value Id_mtpi ^* and the MTPI assumed q-axis current command value Iq_mtpi ^* are calculated. The intersection of the constant torque curve and the MTPI curve is calculated using, for example, the motor torque equation of the equation (10) and the equation (11) which is the d-axis current equation on the MTPI curve when the q-axis current is known. Will be done.

By eliminating the d-axis current Id from the equations (10) and (11), a quartic equation relating to the q-axis current Iq can be obtained as shown in the equation (12).

By using, for example, Newton's method for the quartic equation shown in equation (12) as the solution of the quartic equation shown in equation (12), the intersection of the constant torque curve ^{of the total torque command value T * and the MTPI curve.} It is possible to derive a solution corresponding to the MTPI assumed q-axis current command value Iq_mtpi ^{* in.}

^{The MTPI current command value calculator 14b1-1 is based on the MTPI assumed q-axis current command value Iq_mtpi *} , which is the solution of the equation (12), and according to the d-axis current equation of the equation (11), the MTPI assumed d-axis current command value Id_mtpi ^* Calculate.

The MTPI voltage command value calculator 14b1-2 uses equations (13.1) and (13) based on the MTPI assumed d-axis current command value Id_mtpi ^* , the MTPI assumed q-axis current command value Iq_mtpi ^{*, and the estimated electrical angle angular velocity ωe.} According to the PMSM voltage equation shown in .2), the MTPI assumed d-axis voltage Vd_mtpi ^* and the MTPI assumed q-axis voltage Vq_mtpi ^* are calculated. In addition, "p" in equation (13.1) and equation (13.2) is a differential operator.

In equations (13.1) and (13.2), the voltage drops “p ・ Ld ・ Id” and “p ・ Lq ・ Iq” (p-term voltage) due to the inductance caused by the current change due to torque control are It is being considered. “Ld” indicates the d-axis inductance of the motor M, and “Lq” indicates the q-axis inductance of the motor M.

Here, the p-term voltage is shown using the time derivative of the current change. However, if the amount of change in the detected current is used as a differential value as it is, the MTPI assumed d-axis voltage Vd_mtpi ^* and the MTPI assumed q-axis voltage Vq_mtpi ^* react sensitively to current noise. Therefore, the differential value is generated based on the fluctuation of the current fundamental wave, for example, as follows.

In order to explain the generation of the p-term voltage, first, the variable components ΔIda and ΔIqa of the d-axis current Id and the q-axis current Iq are defined as equations (14.1) and (14.2).

Here, when one periodic fluctuation occurs in one rotation of the mechanical angle, “Ida” and “φd” included in the equation (14.1) indicate the fluctuation amplitude and the initial phase of ΔIda, respectively, and the equation (14. “Iqa” and “φq” included in 2) indicate the fluctuation amplitude and initial phase of ΔIqa, respectively, and “θm” included in equations (14.1) and (14.2) are instantaneous mechanical angle phases. Indicates the value.

Therefore, the p-term voltage generated by the fluctuation of the current fundamental wave is expressed by Eqs. (15.1) and (15.2). That is, the phase of the d-axis current fluctuation and the q-axis current fluctuation is advanced by π / 2, and the d-axis current fluctuation and the q-axis current fluctuation whose phase is advanced by π / 2 are multiplied by the mechanical angle estimated angular velocity ωm. , A differential value (p-term voltage) can be generated.

^{The MTPI voltage amplitude calculator 14b1-3 calculates the MTPI assumed output voltage Va_mtpi *} according to the equation (16) based on the MTPI assumed d-axis voltage Vd_mtpi ^* and the MTPI assumed q-axis voltage Vq_mtpi ^* .

The average output voltage generator 14b1-4 traces the MTPI curve (maximum torque / current control curve) so that the average values of the d-axis current Id and the q-axis current Iq, which fluctuate with each rotation of the rotor of the motor M, trace the MTPI curve (maximum torque / current control curve). Outputs the average output voltage command value Va0 ^* adjusted to. For example, the average output voltage generator 14b1-4 calculates the d-axis current Idt on the MTPI curve from the current q-axis current Iq so that the error between the calculated d-axis current Idt and the current d-axis current Id is eliminated. The average output voltage command value Va0 ^* is adjusted by PI control or the like. The average output voltage generator 14b1-4 calculates the ^{average output voltage command value Va0 *} according to, for example, equations (17.1) and (17.2). The average output voltage generator 14b1-4, when the average output voltage command value Va0 ^* exceeds the output voltage limit value Vdq_limit is according to equation (18), the average output voltage command value Va0 ^* output voltage limit value Limited to Vdq_limit. By limiting the average output voltage command value Va0 ^* to the output voltage limit value Vdq_limit, the weak magnetic flux control is performed. “Ψa” indicates the interlinkage magnetic flux of the motor M.

The MTPI voltage fluctuation component extractor 14b1-5 calculates ^{the fluctuation amplitude | ΔVa_mtpi | of the MTPI assumed output voltage Va_mtpi *} and the instantaneous value ΔVa_mtpi, for example, as follows.

^{The MTPI voltage fluctuation component extractor 14b1-5 first sets the fundamental wave component of the MTPI assumed output voltage Va_mtpi *} according to the equations (19.1) and (19.2), the Fourier coefficient Va_mtpi_sin which is a sin component, and cos. It separates from the component Va_mtpi_cos. MTPi voltage variation component extractor 14b1-5, by calculating the Fourier coefficients of the fundamental wave component of MTPi assumed output voltage Va_mtpi ^* the mechanical angle for each cycle, MTPi assumed output voltage Va_mtpi ^* basic harmonic component is removed Wave components can be extracted.

Next, the MTPI voltage fluctuation component extractor 14b1-5 is based on the Fourier coefficients Va_mtpi_sin and Va_mtpi_cos calculated according to the equations (19.1) and (19.2), and the MTPI assumed output voltage Va_mtpi ^{* according to the equation (20).} Calculate the amplitude | ΔVa_mtpi | of the fundamental wave component of. Since the Fourier coefficients Va_mtpi_sin and Va_mtpi_cos are values that are updated every machine angle period, the amplitude | ΔVa_mtpi | is also updated every machine angle period.

Then, the MTPI voltage fluctuation component extractor 14b1-5 calculates the instantaneous value ΔVa_mtpi of the fundamental wave component of the ^{MTPI assumed output voltage Va_mtpi * according to the equation (21).}

^{The adder 14b1-6 adds the amplitude | ΔVa_mtpi | of the fundamental wave component of the average output voltage command value Va0 *} and the MTPI assumed output voltage Va_mtpi ^* according to the equation (22), thereby adding the MTPI assumed output voltage fluctuation peak value Va_mtpi_peak. Is calculated.

The MTPI voltage amplitude limiting processor 14b1-7 is adjusted so that the MTPI assumed output voltage fluctuation peak value Va_mtpi_peak, which is the addition result of the adder 14b1-6, is equal to or less than the output voltage limit value Vdq_limit. Is generated, and the generated variable output voltage limit command value ΔVa_limit_mtpi is output.

For example, the MTPI voltage amplitude limiting processor 14b1-7 compares the average output voltage command value Va0 ^* , the MTPI assumed output voltage fluctuation peak value Va_mtpi_peak, and the output voltage limit value Vdq_limit with the amplitude ratio scale of the output voltage fluctuation component. Is calculated, and the amplitude ratio scale is multiplied by the MTPI assumed output voltage fluctuation component ΔVa_mtpi to generate the variable output voltage limit command value ΔVa_limit_mtpi. By doing so, it is possible to generate a variable output voltage limit command value ΔVa_limit_mtpi whose phase is matched with the MTPI assumed output voltage fluctuation component ΔVa_mtpi.

For example, the MTPI voltage amplitude limiting processor 14b1-7 calculates the amplitude ratio scale of the output voltage fluctuation component according to the equations (23.1) to (23.3), and based on the calculated amplitude ratio scale, the equation ( According to 23.4), the variable output voltage limit command value ΔVa_limit_mtpi is generated.

The adder 14b1-8 calculates the ^{output voltage limit command value Va * by adding the average output voltage command value Va0 *} and the variable output voltage limit command value ΔVa_limit_mtpi according to the equation (24). The adder 14b1-8 outputs the calculated output voltage limit command value Va ^* to the induced voltage command value calculator 14b2 and the voltage command value calculator 14b6.

In FIG. 5, the induced voltage command value calculator 14b2 is a motor represented by the equations (25.1) and (25.2) based on the current d-axis current Id, q-axis current Iq, and electric angle estimated angular velocity ωe. according to the model equations to calculate the induced voltage command value Vo ^* based on the output voltage limit command value Va ^*. The details of the calculation of the induced voltage command value Vo ^{* are shown below.}

The PMSM voltage equation (d-axis voltage Vd, q-axis voltage Vq), the theoretical formula of the output voltage amplitude Va, and the theoretical formula of the induced voltage Vo of the motor M are shown in equations (25.1) to (27). Is done.

Further, from the equations (25.1) to (27) ^{, the equation for associating the output voltage limit command value Va *} with the induced voltage command value Vo ^* is as shown in the equation (28). Therefore, the induced voltage command value calculator 14b2 calculates the induced voltage command value Vo ^* according to the equation (28), and outputs the calculated induced voltage command value Vo ^* to the current command value calculator 14b3.

The current command value calculator 14b3 has ^{a constant torque curve which is a current locus in which the total torque command value T *} is constant, and a current locus in which the induced voltage command value Vo ^* and the electric angle estimated angular velocity ωe are constant. ^{The q-axis current command value Iq *} and the d-axis current command value Id ^* are calculated based on the intersection with the induced voltage ellipse (see FIG. 1B). The current command value calculator 14b3 outputs the calculated q-axis current command value Iq ^* and d-axis current command value Id ^* to the temporary voltage command value calculator 14b4.

The intersection of the constant torque curve and the constant induced voltage ellipse can be calculated using, for example, the motor torque equation shown in the equation (29) and the induced voltage equation shown in the equation (30).

By eliminating the d-axis current Id from the equations (29) and (30), a quartic equation relating to the q-axis current Iq can be obtained as in the equation (31). However, in the formula (31), "ΔL = Ld-Lq".

By using, for example, Newton's method for the quartic equation shown in equation (31) as the solution of the quartic equation shown in equation (31), the total torque command value T ^* is a constant current locus. A solution corresponding to ^{the q-axis current command value Iq *} at the intersection of the torque curve and the constant induced voltage ellipse, which is the locus of the current at which the induced voltage Vo and the estimated electrical angle velocity ωe are constant, can be derived. (See FIG. 1B).

Current command value calculator 14b3 after calculating a q-axis current command value Iq ^*, in accordance with the equation the induced voltage type shown in (30) deformed into d-axis current equation Equation (32), the q-axis current command value Iq ^* Based on this, the d-axis current command value Id ^* is calculated.

Here, in Eq. (32), whether positive or negative is adopted as the sign related to √ is the point M (−Ψa / Ld, 0) parallel to the Iq axis and the center of the constant induced voltage ellipse. The torque at the intersection of the straight line passing through (hereinafter sometimes referred to as "M point boundary line") and the constant induced voltage ellipse (hereinafter sometimes referred to as "M point boundary torque T_M") It can be calculated and determined by comparing the torque T_M on the boundary of the M point with the total torque command value T ^*.

The procedure for calculating the d-axis current command value Id ^* and the q-axis current command value Iq ^* is shown below. 7A and 7B are diagrams for explaining an operation example of the current command value calculator according to the first embodiment of the present disclosure.

The current command value calculator 14b3 first calculates the d-axis current Id_M on the M point according to the equation (33).

Next, the current command value calculator 14b3 calculates the q-axis current Iq_M at the point where the M point boundary line and the constant induced voltage ellipse intersect. Since the q-axis current Iq_M can be calculated by substituting the d-axis current Id_M on the M point into the equation (30), it is calculated according to the equation (34).

Therefore, the current command value calculator 14b3 calculates the torque T_M on the boundary of the M point according to the equation (35).

Then, the current command value calculator 14b3 uses the d-axis current command according to the equations (36.1) and (36.2) based on the magnitude relationship between ^{the total torque command value T * and the torque T_M on the M point boundary.} Determine the value Id ^{*. Equation (36.1) shows the d-axis current command value Id *} (see FIG. 7A) in the case of "total torque command value T ^* ≤ torque on the M point boundary T_M", and equation (36.2) shows. , The d-axis current command value Id ^* (see FIG. 7B) in the case of "total torque command value T ^{*> torque on M point boundary T_M" is shown.}

The current command value calculator 14b3 outputs the d-axis current command value Id ^* and the q-axis current command value Iq ^* calculated as described above to the temporary voltage command value calculator 14b4.

In FIG. 5, the temporary voltage command value calculator 14b4 has equations (37.1) and (37.2) based on the estimated electrical angle angular velocity ωe, the d-axis current command value Id ^*, and the q-axis current command value Iq ^*. ), The temporary d-axis voltage command value Vd_m and the temporary q-axis voltage command value Vq_m are calculated by feed forward. The temporary voltage command value calculator 14b4 outputs the calculated temporary d-axis voltage command value Vd_m and the temporary q-axis voltage command value Vq_m to the voltage vector angle calculator 14b5. By calculating the temporary voltage command value by feedforward, it is possible to prevent windup (saturation phenomenon) that occurs in the integrator such as PI control due to input saturation.

In equations (37.1) and (37.2), the voltage drops “p ・ Ld ・ Id” and “p ・ Lq ・ Iq” (p-term voltage) due to the inductance caused by the current change due to torque control are It is being considered.

The voltage vector angle calculator 14b5 calculates the voltage vector angle δ according to the equation (38) based on the temporary d-axis voltage command value Vd_m and the temporary q-axis voltage command value Vq_m. The voltage vector angle calculator 14b5 outputs the calculated voltage vector angle δ to the voltage command value calculator 14b6. That is, as shown in FIG. 8, the voltage vector angle calculator 14b5 calculates the angle formed by the output voltage vector having the amplitude Va calculated by the equation (1) as the voltage vector angle δ. By doing so, in the voltage saturation region where the output voltage amplitude is limited to the DC voltage (DC voltage) or less that can be output by the inverter, the voltage vector angle δ corresponding to the ^{total torque command value T * can be generated by calculation.} Therefore, since the vibration damping control can be performed without tuning the voltage vector angle fluctuation, the vibration damping effect of the motor M can be improved without tuning the voltage vector angle fluctuation. FIG. 8 is a diagram provided for explaining an operation example of the voltage vector angle calculator according to the first embodiment of the present disclosure.

The voltage command value calculator 14b6 performs coordinate conversion from polar coordinates to Cartesian coordinates according to equations (39.1) and (39.2) based on the voltage vector angle δ and the output voltage limit command value Va ^*. The d-axis voltage command value Vd ^* and the q-axis voltage command value Vq ^* are calculated accordingly.

<Operation of MTPI voltage amplitude limiting processor>
FIG. 9 is a diagram provided for explaining an operation example of the MTPI voltage amplitude limiting processor according to the first embodiment of the present disclosure.

For example, as shown in the case (a) of FIG. 9, ^{when the peak value Va_mtpi_peak of the MTPI assumed output voltage fluctuation component ΔVa_mtpi that fluctuates around the average output voltage command value Va0 *} is equal to or less than the output voltage limit value Vdq_limit, the formula is used. Since the condition of (23.2) is satisfied, the MTPI voltage amplitude limiting processor 14b1-7 sets the amplitude ratio scale of the output voltage fluctuation component to “1”. Then, the MTPI voltage amplitude limiting processor 14b1-7 sets “scale = 1” in the equation (23.4), and outputs the MTPI assumed output voltage fluctuation component ΔVa_mtpi as it is as the variable output voltage limiting command value ΔVa_limit_mtpi. As a result, the output voltage limit command value Va ^* coincides with the MTPI assumed output voltage fluctuation component ΔVa_mtpi.

Further, for example, as shown in the case (b) of FIG. 9, the peak value Va_mtpi_peak of the MTPI assumed output voltage fluctuation component ΔVa_mtpi that fluctuates around ^{the average output voltage command value Va0 * exceeds the output voltage limit value Vdq_limit and is average.} If the output voltage command value Va0 ^* does not exceed the output voltage limit value Vdq_limit, the condition of Eq. (23.3) is met. Therefore, the MTPI voltage amplitude limiting processor 14b1-7 uses the amplitude ratio scale of the output voltage fluctuation component. Is "(Vdq_limit-Va0 ^* ) / | ΔVa_mtpi |". Then, the MTPI voltage amplitude limiting processor 14b1-7 outputs the variable output voltage limiting command value ΔVa_limit_mtpi with ^{“scale = (Vdq_limit-Va0 *) / | ΔVa_mtpi |” in the equation (23.4). As a result, an output voltage limit command value Va *} is generated in which the phase matches the MTPI assumed output voltage fluctuation component ΔVa_mtpi and the peak value due to the fluctuation amplitude is equal to or less than the output voltage limit value Vdq_limit.

Further, for example, as shown in the case (c) of FIG. 9, when the average output voltage command value Va0 ^* of the MTPI assumed output voltage fluctuation component ΔVa_mtpi is equal to or more than the output voltage limit value Vdq_limit, the condition of the equation (23.1) is satisfied. Therefore, the MTPI voltage amplitude limiting processor 14b1-7 sets the amplitude ratio scale of the output voltage fluctuation component to “0”. Then, the MTPI voltage amplitude limiting processor 14b1-7 outputs “scale = 0” in the equation (23.4) and outputs the variable output voltage limiting command value ΔVa_limit_mtpi as “0”. As a result, the output voltage limit command value Va ^* matches the output voltage limit value Vdq_limit.

By controlling the output voltage limit command value Va ^* in this way, as in the example of the output voltage waveform shown in FIG. 10, the output is output even immediately after the control region of the motor M transitions from the normal control region to the voltage saturation region. While keeping the voltage amplitude Va below the output voltage limit value Vdq_limit, the output voltage amplitude Va can be matched between the normal control region and the voltage saturation region. Therefore, it is possible to reduce the switching shock at the time of transition from the normal control region to the voltage saturation region. Further, since the motor control device 100 has the voltage saturation region voltage command value generator 14b, the average output voltage command value Va0 ^{*, which is the center of fluctuation of the output voltage limit command value Va *, is the output voltage limit in the voltage saturation region.} It can also be used when it is limited by the value Vdq_limit.

<Configuration of normal control area voltage command value generator>
FIG. 11 is a diagram showing a configuration example of the normal control region voltage command value generator according to the first embodiment of the present disclosure. In FIG. 11, the normal control region voltage command value generator 14a includes a current command value calculator 14a1, adders 16, 17, 21, 22, subtractors 18, 19, voltage command value calculator 20, and IIR. It has (Infinite Impulse Response) filters 35a and 35b, a decoupling controller 36, and a current error correction value generator 37.

The current command value calculator ^{14a1 has a q-axis current command value Iq *} and a d-axis current command value Id based on the intersection of the constant torque curve and the MTPI curve, which are the loci of the current at which the ^{total torque command value T * is constant. *} Calculate.

Here, the intersection of the constant torque curve and the MTPI curve is, for example, an equation showing the relationship between the motor torque equation shown in the equation (29) and the d-axis current Id and the q-axis current Iq in the MTPI curve (Equation (17. It can be calculated using 1) and. On the right side of the equation (29), the first term represents the magnet torque, the second term represents the reluctance torque, and the magnet torque includes only the q-axis current Iq, and the reluctance torque. Includes both the q-axis current Iq and the d-axis current Id. Therefore, by appropriately controlling the q-axis current Iq and the d-axis current Id, it is possible to generate an appropriate torque in the motor M.

By eliminating the d-axis current Id from the equations (29) and (17.1), the equation (40), which is a quartic equation relating to the q-axis current Iq, can be obtained.

By using, for example, Newton's method for the quartic equation shown in equation (40) as the solution of the quartic equation shown in equation (40), the intersection of the constant torque curve ^{of the total torque command value T * and the MTPI curve.} A solution corresponding to the q-axis current command value Iq ^{* can be derived.} Further, the current command value calculator 14a1 calculates the d-axis current command value Id ^{* by substituting the q-axis current command value Iq *} calculated according to the equation (40) into the equation (17.1).

The adder 17 has a q-axis current command value Iq ^* output from the current command value calculator 14a1 and a q-axis current error correction value ΔIq output from the current error correction value generator 37 according to the equation (41.1). Is added to calculate the q-axis current correction command value Iq_FF ^*. The adder 16 has a d-axis current command value Id ^* output from the current command value calculator 14a1 and a d-axis current error correction value ΔId output from the current error correction value generator 37 according to the equation (41.2). Is added to calculate the d-axis current correction command value Id_FF ^*.

The subtractor 19 subtracts the q-axis current Iq output from the u, v, w / d-q converter 29 from the q-axis current correction command value Iq_FF ^{* output from the adder 17, thereby subtracting the q-axis current Iq.} The q-axis current error Iq_diff, which is the error between the current correction command value Iq_FF ^{* and the q-axis current Iq, is calculated.} The subtractor 18 subtracts the d-axis current Id output from the u, v, w / dq converter 29 from the d-axis current correction command value Id_FF ^{* output from the adder 16, thereby subtracting the d-axis current Id.} The d-axis current error Id_diff, which is the error between the current correction command value Id_FF ^{* and the d-axis current Id, is calculated.}

The voltage command value calculator 20 calculates the q-axis voltage command value Vqt before decoupling by performing PI (Proportional Integral) control on the q-axis current error Iq_diff (Iq_FF ^{* −Iq) according to the equation (42.1).} To do. Further, the voltage command value calculator 20 calculates the pre-interference d-axis voltage command value Vdt by performing PI control on ^{the d-axis current error Id_diff (Id_FF * −Id) according to the equation (42.2).} Note that kp_q in Eq. (42.1) and kp_d in Eq. (42.2) are proportional constants, and ki_q in Eq. (42.1) and ki_d in Eq. (42.2) are constants of integration.

The adder 22 adds the q-axis decoupling correction value Vqa represented by the equation (43.1) to the q-axis voltage command value Vqt before decoupling according to the equation (43.3), thereby q. Calculate the shaft voltage command value Vq ^*. The adder 21 adds the d-axis decoupling correction value Vda represented by the equation (43.2) to the d-axis voltage command value Vdt before decoupling according to the equation (43.4), thereby d. Calculate the shaft voltage command value Vd ^*. As a result, the q-axis voltage command value Vq ^* and the d-axis voltage command value Vd ^{* in} which the interference between the dq axes is canceled by feedforward are calculated.

The IIR filter (Infinite Impulse Response Filter) 35a removes the noise of the d-axis current Id output from the u, v, w / d-q converter 29, and outputs the d-axis response current Id_ir after removing the noise. The IIR35b removes the noise of the q-axis current Iq output from the u, v, w / d-q converter 29, and outputs the q-axis response current Iq_ir after removing the noise. The IIR filters 35a and 35b are examples of noise reduction filters.

^{The decoupling controller 36 is based on the electric angular velocity command value ωe *} input from the outside of the motor control device 100 (for example, the upper controller) and the q-axis response current Iq_ir, and the d-axis voltage before decoupling. A d-axis decoupling correction value Vda for correcting the command value Vdt is generated. Further, the decoupling controller 36 is a q-axis decoupling correction value Vqa for correcting the q-axis voltage command value Vqt before decoupling based on the electric angular velocity command value ^{ωe * and the d-axis response current Id_ir.} To generate. The d-axis non-interference correction value Vda and the q-axis non-interference correction value Vqa are correction values for canceling the interference term between the dq axes by feedforward. Here, in order to achieve stable control, it is desirable that the non-interference correction value is a DC value. Therefore, in generating the non-interference correction value, the electric angular velocity command value ωe ^* is used for the speed, and for the d-axis current Id and the q-axis current Iq, the d-axis response in which the variable component is removed by the IIR filter. The current Id_ir and the q-axis response current Iq_ir are used.

The current error correction value generator 37 outputs the d-axis current command value Id ^* and the q-axis current command value Iq ^* output from the current command value calculator 14a1 and the u, v, w / dq converter 29. The d-axis current error correction value ΔId and the q-axis current error correction value ΔIq are generated based on the d-axis current Id and the q-axis current Iq and the mechanical angle phase θm output from the position estimator 32.

The current error correction value generator 37 integrates fluctuation errors (phase error and amplitude error) that occur when the dq-axis current cannot follow the current command value due to the response delay of the current command value calculator 14a1 and the interference of the dq axis. Then, the inverted output of the integrated value is generated as a current error correction value (d-axis current error correction value ΔId and q-axis current error correction value ΔIq). Here, the d-axis current error correction value ΔId is ^{a feed-forward component for correcting the fluctuation error between the d-axis current command value I *} and the d-axis current Id, and the q-axis current error correction value ΔIq is the q-axis current error correction value ΔIq. This is a feed-forward component for correcting the fluctuation error between the current command value Iq ^{* and the q-axis current Iq.}

<Configuration of current error correction value generator>
FIG. 12 is a diagram showing a configuration example of the current error correction value generator according to the first embodiment of the present disclosure. In FIG. 12, the current error correction value generator 37 includes subtractors 37a and 37e, a q-axis current error component separator 37b, a q-axis current error accumulator 37c, and a q-axis current error correction value demodulator 37d. It has a d-axis current error component separator 37f, a d-axis current error accumulator 37g, and a d-axis current error correction value demodulator 37h.

The subtractor 37a calculates the q-axis current fluctuation error Iq_err, which is the error between the q-axis current Iq and the q-axis current command value Iq ^{*, according to the equation (44).}

The q-axis current error component separator 37b has two Fourier coefficients Iq_err_sin (sin component) and IQ_err_cos (cos), which are fundamental wave components of the q-axis current fluctuation error Iq_err, according to equations (45.1) and (45.2). Component) is calculated for each machine angle period.

The q-axis current error integrator 37c is provided with the sin component Iq_err_sin of the q-axis current fluctuation error Iq_err and the cos component Iq_err_cos of the q-axis current fluctuation error Iq_err according to the equations (46.1) and (46.2), respectively. The correction gain k is multiplied, and the respective multiplication results are added to Iq_err_sin_i_old and Iq_err_cos_i_old. Iq_err_sin_i in the equation (46.1) is the integrated value of IQ_err_sin up to the current machine angle period, and Iq_err_cos_i in the equation (46.2) is the integrated value of Iq_err_cos up to the current machine angle period. Further, IQ_err_sin_i_old in the equation (46.1) is IQ_err_sin_i up to the previous machine angle period, and Iq_err_cos_i_old in the equation (46.2) is Iq_err_cos_i up to the previous machine angle period.

The q-axis current error correction value demodulator 37d calculates the q-axis current error correction value ΔIq according to the equations (47.1) and (47.2). As a result, the phase of the q-axis current fluctuation error is inverted, and the instantaneous value of the q-axis current error correction value ΔIq at the mechanical angle phase θm is calculated.

The subtractor 37e calculates the d-axis current fluctuation error Id_err, which is the error between the d-axis current Id and the d-axis current command value Id ^{*, according to the equation (48).}

According to equations (49.1) and (49.2), the d-axis current error component separator 37f has two Fourier coefficients, Id_err_sin (sin component) and Id_err_cos (cos), which are fundamental wave components of the d-axis current fluctuation error Id_err. Component) is calculated for each mechanical angle period.

According to the equations (50.1) and (50.2), the d-axis current error integrator 37g has the sin component Id_err_sin of the d-axis current fluctuation error Id_err and the cos component Id_err_cos of the d-axis current fluctuation error Id_err. The correction gain k is multiplied, and the respective multiplication results are added to Id_err_sin_i_old and Id_err_cos_i_old. Id_err_sin_i in the equation (50.1) is an integrated value of Id_err_sin up to the current machine angle period, and Id_err_cos_i in the equation (50.2) is an integrated value of Id_err_cos up to the current machine angle period. Further, Id_err_sin_i_old in the equation (50.1) is Id_err_sin_i up to the previous machine angle period, and Id_err_cos_i_old in the equation (50.2) is Id_err_cos_i up to the previous machine angle period.

The d-axis current error correction value demodulator 37h calculates the d-axis current error correction value ΔId according to the equations (51.1) and (51.2). As a result, the phase of the d-axis current fluctuation error is inverted, and an instantaneous value of the d-axis current error correction value ΔId at the mechanical angle phase θm is generated.

The first embodiment of the present disclosure has been described above.

[Example 2]
Example 1 of the present disclosure is effective when the distortion of the induced voltage of the motor M is small. However, the actual induced voltage waveform has an induced voltage distortion due to the structure of the motor M, and when this induced voltage distortion is large, the harmonic of the current of the motor M becomes large. As a result, the followability of the current control of the motor M may deteriorate. Therefore, in the second embodiment, the followability of the current control of the motor M is enhanced to further improve the stability of the control. Therefore, in the second embodiment, the configuration of the voltage saturation region voltage command value generator 14b is partially different from that of the first embodiment.

<Configuration of voltage saturation region voltage command value generator>
FIG. 13 is a diagram showing a configuration example of the voltage saturation region voltage command value generator 14b according to the second embodiment of the present disclosure. In FIG. 13, the voltage saturation region voltage command value generator 14b includes an output voltage limit command value generator 14b1, an induced voltage command value calculator 14b2, a current command value calculator 14b3, and a temporary voltage command value calculator 14b4. The point that the voltage vector angle calculator 14b5 and the voltage command value calculator 14b6 are provided is the same as that of the first embodiment (FIG. 5). In the second embodiment, the voltage saturation region voltage command value generator 14b further includes a current proportional controller 14b7, adders 14b8, 14b9, and subtractors 14b10, 14b11.

In FIG. 13, the voltage saturation region voltage command value generator 14b outputs the d-axis voltage command value Vd ^** calculated by the voltage command value calculator 14b6 to the adder 14b8 and calculates it in the same manner as in the first embodiment. The q-axis voltage command value Vq ^** is output to the adder 14b9.

In the first embodiment, the d-axis voltage command value and the q-axis voltage command value output from the voltage command value calculator 14b6 are described as Vd ^* and Vq ^* , respectively, whereas in the second embodiment, the voltage command value is described. The d-axis voltage command value and the q-axis voltage command value output from the calculator 14b6 are expressed as Vd ^** and Vq ^** , respectively. That is, the d-axis voltage command value Vd ^** and the q-axis voltage command value V ^** in the second embodiment correspond to the d-axis voltage command value Vd ^* and the q-axis voltage command value Vq ^* in the first embodiment.

The current command value calculator 14b3 outputs the calculated q-axis current command value Iq ^* and d-axis current command value Id ^* to the temporary voltage command value calculator 14b4. Further, the current command value calculator 14b3 outputs the calculated q-axis current command value Iq ^* to the subtractor 14b11, and outputs the calculated d-axis current command value Id ^* to the subtractor 14b10.

The subtractor 14b10 calculates the d-axis current error Id_p by subtracting the d-axis current Id from the d-axis current command value Id ^{*, and outputs the calculated d-axis current error Id_p to the current proportional controller 14b7.} The subtractor 14b11 calculates the q-axis current error Iq_p by subtracting the q-axis current Iq from the q-axis current command value Iq ^{*, and outputs the calculated q-axis current error Iq_p to the current proportional controller 14b7.}

The current proportional controller 14b7 calculates the d-axis compensation voltage Vd_p by multiplying the d-axis current error Id_p by the proportionality constant kp_d, and outputs the calculated d-axis compensation voltage Vd_p to the adder 14b8. Further, the current proportional controller 14b7 calculates the q-axis compensation voltage Vq_p by multiplying the q-axis current error Iq_p by the proportionality constant kp_q, and outputs the calculated q-axis compensation voltage Vq_p to the adder 14b9. That is, the d-axis compensation voltage Vd_p is calculated as “kp_d ・ Id_p”, and the q-axis compensation voltage Vq_p is calculated as “kp_q ・ Iq_p”.

The adder 14b8 adds the d-axis compensation voltage Vd_p output from the current proportional controller 14b7 to the d-axis voltage command value Vd ^** output from the voltage command value calculator 14b6 to obtain the final d-axis. Calculate the voltage command value Vd ^*. Further, the adder 14b9 finally adds the q-axis compensation voltage Vq_p output from the current proportional controller 14b7 to the q-axis voltage command value Vq ^{** output from the voltage command value calculator 14b6.} Calculate the q-axis voltage command value Vq ^*.

By adding the current proportional controller 14b7 in this way, the current harmonics due to the induced voltage distortion can be suppressed, and the followability to the current command value can be further improved.

The second embodiment has been described above.

As described above, the motor control device (motor control device 100 of the first embodiment) of the present disclosure includes a voltage command value generator (voltage command value generator 14 of the first embodiment) and a control switching determination unit 15. .. The control switching determination unit 15 determines whether the control region of the motor (motor M of the first embodiment) is in the voltage saturation region or the normal control region. The voltage command value generator is based on the speed command value (mechanical angular velocity command value ωm ^{* in} Example 1) and the motor speed (machine angle estimated angular velocity ωm in Example 1), and the voltage command value of the motor (Example 1). The d-axis voltage command value Vd ^* and the q-axis voltage command value Vq ^* ) are generated. ^{Further, the voltage command value generator outputs the torque command value (total torque command value T * of the} first embodiment) to the motor when the control switching determination unit determines that the control region of the motor is in the voltage saturation region. The voltage vector angle of the output voltage applied to the motor (voltage vector angle δ of Example 1) is obtained from the limit value of the maximum possible voltage (output voltage limit command value Va ^{* of Example 1), and this voltage vector angle is obtained.} Generates a voltage command value based on.

By doing so, it is possible to optimize the output torque fluctuation with respect to the speed fluctuation while suppressing the fluctuation of the voltage amplitude, so that it is possible to improve the vibration damping effect of the motor at the time of torque control in the voltage saturation region. ..

Further, the voltage command value calculator is based on the intersection of the constant torque curve of the motor and the constant induced voltage ellipse of the motor, and the current command value of the motor (q-axis current command value Iq ^* and d-axis current command value of the first embodiment). Id ^* ) is calculated.

Further, the voltage command value generator has a temporary voltage command value (temporary d-axis of the first embodiment) according to the motor model formula (formula (37.1) and formula (37.2) of the second embodiment) based on the current command value. The voltage command value Vd_m and the temporary q-axis voltage command value Vq_m) are calculated, and the voltage vector angle is calculated based on the temporary voltage command value.

By doing so, the voltage vector angle can be calculated by feedforward, and integration control becomes unnecessary, so that the occurrence of windup (saturation phenomenon) can be prevented.

Further, the voltage command value calculator (voltage command value calculator 14 of Example 2) has an error (Example 2) between the current command value and the motor current (d-axis current Id and q-axis current Iq of Example 2). The compensation voltage (d-axis compensation voltage Vd_p and q-axis compensation voltage Vq_p in Example 2) is obtained by multiplying the d-axis current error Id_p and the q-axis current error Iq_p) by the proportionality constant (proportional constant ka in Example 2). It has a proportional controller (current proportional controller 14b7 of the second embodiment) to be calculated, and calculates a voltage command value to which a compensation voltage is added.

By doing so, it is possible to suppress the harmonic current generated by the induced voltage distortion or the like only by proportional control and improve the followability to the current command value.

100 Motor control device 14 Voltage command value generator 14a Normal control area Voltage command value generator 14b Voltage saturation area Voltage command value generator 14b1 Output voltage limit command value generator 14b2 Induced voltage command value calculator 14b3 Current command value calculator 14b4 Temporary voltage command value calculator 14b5 Voltage vector angle calculator 14b6 Voltage command value calculator 14b7 Current proportional controller

Claims (4)

A voltage command value generator that generates a voltage command value for the motor from a torque command value based on the speed command value and the speed of the motor.
A control switching determination unit for determining whether or not the control region of the motor is in the voltage saturation region is provided.
When the control switching determination unit determines that the control region is in the voltage saturation region, the voltage command value generator determines that the control region is in the voltage saturation region.
The voltage command value is generated based on the torque command value, the limit value of the maximum voltage that can be output to the motor, and the voltage vector angle of the output voltage applied to the motor.
Motor control device. The voltage command value generator calculates the current command value of the motor based on the intersection of the constant torque curve of the motor of the torque command value and the constant induced voltage ellipse of the motor, and the voltage vector from the current command value. Find the angle and generate the voltage command value,
The motor control device according to claim 1. The voltage command value generator calculates a temporary voltage command value according to a motor model formula based on the current command value, and calculates the voltage vector angle based on the temporary voltage command value.
The motor control device according to claim 2. The voltage command value generator calculates the compensation voltage by multiplying the error between the current command value of the motor and the current of the motor by a proportional constant, and generates the voltage command value to which the compensation voltage is added. ,
The motor control device according to claim 1.

Images