WO2021153020A1 - モータ制御装置 - Google Patents
モータ制御装置 Download PDFInfo
- Publication number
- WO2021153020A1 WO2021153020A1 PCT/JP2020/045886 JP2020045886W WO2021153020A1 WO 2021153020 A1 WO2021153020 A1 WO 2021153020A1 JP 2020045886 W JP2020045886 W JP 2020045886W WO 2021153020 A1 WO2021153020 A1 WO 2021153020A1
- Authority
- WO
- WIPO (PCT)
- Prior art keywords
- command value
- voltage
- motor
- current
- axis
- Prior art date
Links
- 238000013016 damping Methods 0.000 abstract description 8
- 230000000694 effects Effects 0.000 abstract description 7
- BPTHCFBWOUNCBG-UHFFFAOYSA-N iodo-methyl-triphenoxy-$l^{5}-phosphane Chemical compound C=1C=CC=CC=1OP(I)(OC=1C=CC=CC=1)(C)OC1=CC=CC=C1 BPTHCFBWOUNCBG-UHFFFAOYSA-N 0.000 description 67
- 238000010586 diagram Methods 0.000 description 27
- 230000004044 response Effects 0.000 description 9
- 230000008859 change Effects 0.000 description 8
- 238000000034 method Methods 0.000 description 8
- 230000004907 flux Effects 0.000 description 6
- 230000000737 periodic effect Effects 0.000 description 6
- 230000003313 weakening effect Effects 0.000 description 3
- 230000010354 integration Effects 0.000 description 2
- 230000001360 synchronised effect Effects 0.000 description 2
- 238000013459 approach Methods 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000006835 compression Effects 0.000 description 1
- 238000007906 compression Methods 0.000 description 1
- 230000005347 demagnetization Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000006698 induction Effects 0.000 description 1
- 238000007562 laser obscuration time method Methods 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 239000003507 refrigerant Substances 0.000 description 1
- 229920006395 saturated elastomer Polymers 0.000 description 1
- 230000035939 shock Effects 0.000 description 1
- 230000001629 suppression Effects 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
Images
Classifications
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02P—CONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
- H02P21/00—Arrangements or methods for the control of electric machines by vector control, e.g. by control of field orientation
- H02P21/05—Arrangements or methods for the control of electric machines by vector control, e.g. by control of field orientation specially adapted for damping motor oscillations, e.g. for reducing hunting
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02P—CONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
- H02P21/00—Arrangements or methods for the control of electric machines by vector control, e.g. by control of field orientation
- H02P21/0003—Control strategies in general, e.g. linear type, e.g. P, PI, PID, using robust control
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02P—CONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
- H02P21/00—Arrangements or methods for the control of electric machines by vector control, e.g. by control of field orientation
- H02P21/14—Estimation or adaptation of machine parameters, e.g. flux, current or voltage
- H02P21/20—Estimation of torque
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02P—CONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
- H02P21/00—Arrangements or methods for the control of electric machines by vector control, e.g. by control of field orientation
- H02P21/22—Current control, e.g. using a current control loop
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02P—CONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
- H02P6/00—Arrangements for controlling synchronous motors or other dynamo-electric motors using electronic commutation dependent on the rotor position; Electronic commutators therefor
- H02P6/34—Modelling or simulation for control purposes
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02P—CONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
- H02P27/00—Arrangements or methods for the control of AC motors characterised by the kind of supply voltage
- H02P27/04—Arrangements or methods for the control of AC motors characterised by the kind of supply voltage using variable-frequency supply voltage, e.g. inverter or converter supply voltage
- H02P27/06—Arrangements or methods for the control of AC motors characterised by the kind of supply voltage using variable-frequency supply voltage, e.g. inverter or converter supply voltage using dc to ac converters or inverters
Definitions
- This disclosure relates to a motor control device.
- 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”).
- speed fluctuation causes motor vibration.
- torque control periodic disturbance suppression control
- 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).
- the low rotation region for example, the normal control region where the maximum torque / current control of the motor is performed.
- vibration occurs even in a high rotation region (for example, a voltage saturation region in which weakening magnetic flux control is performed), and the peak current of the motor increases due to the generation of vibration.
- the protection function of the inverter may operate to prevent demagnetization of the motor, and the motor may stop.
- 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.
- 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.
- a motor control device that controls 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
- PMSM Permanent Magnet Synchronous Motor
- a motor control device used as a storage device or the like will be described as an example.
- 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.
- 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.
- FIG. 1A shows the current vector locus when the estimated electrical angle angular velocity ⁇ e is constant.
- 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.
- 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.
- 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.
- FIG. 2 is a diagram showing a configuration example of the motor control device according to the first embodiment of the present disclosure.
- 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 dq / 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.
- the motor control device 100 has 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.
- 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.
- PLL Phase Locked Loop
- 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.
- 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.
- the control signal CONTROL_TYPE: B voltage saturation control
- 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 dq / u, v, w converter 23.
- 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.
- 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 dq / 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.
- the dq / 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.
- 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.
- 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.
- 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 / dq converter 29. Using the d-axis current Id and the q-axis current Iq, the axis error ⁇ (difference between the estimated rotation axis and the actual rotation axis) is calculated. 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 dq / u, v, w converter 23 and the u, v, w / dq 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 number of pole pairs 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 permissible value
- periodic rate is the variation of mechanical angle estimated angular velocity change
- * generates fluctuation torque command value ⁇ T for suppressing below, the resulting variations
- the torque command value ⁇ T is output to the adder 13.
- * 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 ⁇ .
- FIG. 3 is a diagram showing a configuration example of the control switching determination unit according to the first embodiment of the present disclosure.
- 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. 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.
- 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.
- 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.
- FIG. 4 is a diagram showing a configuration example of the correction torque generator according to the first embodiment of the present disclosure.
- 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.
- the fluctuation torque command value (correction torque) ⁇ T so that the speed fluctuation amplitude
- 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.
- ⁇ sin and ⁇ cos are values that are updated every mechanical angle period.
- the velocity fluctuation amplitude calculator 34b calculates the velocity fluctuation amplitude
- the subtractor 34c calculates the velocity fluctuation error
- * defines the permissible speed fluctuation amplitude
- the corrected torque amplitude calculator 34d adjusts the corrected torque amplitude
- 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. do.
- the ⁇ sin_i_old in the equation (5.1) is ⁇ sin_i in the previous mechanical angle period
- the ⁇ cos_i_old in the equation (5.2) is ⁇ cos_i in the previous mechanical angle period.
- 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.
- the orthogonal component separator 34f and the sin component ( ⁇ sin_i) of the velocity fluctuation correction phase ⁇ i are 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.
- 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).
- 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).
- 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.
- 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 of the first embodiment of the present disclosure.
- 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 an apparatus 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 estimated electrical angle 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.
- 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.
- 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 electric angle estimated 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.
- the p-term voltage is shown using the time derivative of the current change.
- 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.
- 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).
- equation (14.1) 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.
- 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. Adjust the average output voltage command value Va0 * 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.
- 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
- 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 is separated 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.
- 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).
- 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
- 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.
- 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.
- the MTPI voltage amplitude limiting processor 14b1-7 calculates the amplitude ratio scale of the output voltage fluctuation component according to the formulas (23.1) to (23.3), and based on the calculated amplitude ratio scale, the formula ( 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.
- 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 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).
- 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).
- 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).
- 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).
- the current command value calculator 14b3 calculates the torque T_M on the M point boundary according to the equation (35).
- 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.
- the temporary voltage command value calculator 14b4 has equations (37.1) and (37.2) based on the electric angle estimated 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.
- 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.
- 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.
- 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.
- the MTPI voltage amplitude limiting processor 14b1-7 sets the amplitude ratio scale of the output voltage fluctuation component to “1”.
- the output voltage limit command value Va * coincides with the MTPI assumed output voltage fluctuation component ⁇ Va_mtpi.
- 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.
- the output voltage limit command value Va * 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.
- 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.
- 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.
- the intersection of the constant torque curve and the MTPI curve is, for example, the motor torque equation shown in the equation (29) and the equation (17.1) showing the relationship between the d-axis current Id and the q-axis current Iq in the MTPI curve.
- the first term represents the magnet torque
- the second term represents the reluctance torque
- the magnet torque includes only the q-axis current Iq
- the reluctance torque represents the q-axis current Iq and the d-axis current. Includes both with 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.
- the equation (40) which is a quartic equation relating to the q-axis current Iq, can be obtained.
- 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). 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 *.
- 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 / dq 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 / dq 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.
- 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.
- 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 velocity, 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. 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 d-axis current error correction value ⁇ Id and the q-axis current error correction value ⁇ Iq are generated.
- 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).
- 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
- 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.
- 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.
- 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 has a sin component Iq_err_sin of the q-axis current fluctuation error Iq_err and a 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
- Iq_err_cos_i in the equation (46.2) is the integrated value of IQ_err_cos up to the current machine angle period.
- IQ_err_sin_i_old in the equation (46.1) is IQ_err_sin_i up to the previous machine angle period
- 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).
- 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.
- the d-axis current error integrator 37g is provided with 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 according to the equations (50.1) and (50.2), respectively.
- 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
- Id_err_cos_i in the equation (50.2) is an integrated value of Id_err_cos up to the current machine angle period.
- Id_err_sin_i_old in the equation (50.1) is Id_err_sin_i up to the previous machine angle period
- 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.
- Example 1 of the present disclosure is effective when the distortion of the induced voltage of the motor M is small.
- the actual induced voltage waveform has an induced voltage distortion due to the structure of the motor M, and when the induced voltage distortion is large, the harmonic of the current of the motor M becomes large.
- 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.
- 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.
- 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).
- the voltage saturation region voltage command value generator 14b further includes a current proportional controller 14b7, adders 14b8, 14b9, and subtractors 14b10, 14b11.
- 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.
- 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 *.
- 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.
- 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.
- 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.
- 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.
- 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.
- 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.
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Control Of Ac Motors In General (AREA)
Abstract
Description
<モータ制御装置の動作>
図1A及び図1Bは、本開示の実施例1のモータ制御装置の動作例の説明に供する図である。
図2は、本開示の実施例1のモータ制御装置の構成例を示す図である。図2において、モータ制御装置100は、減算器11,38と、速度制御器12と、加算器13と、電圧指令値生成器14と、制御切替判定部15と、d-q/u,v,w変換器23と、PWM(Pulse Width Modulation)変調器24と、IPM(Intelligent Power Module)25とを有する。IPM25は、モータMに接続される。モータMの一例としてPMSMが挙げられる。
図3は、本開示の実施例1の制御切替判定部の構成例を示す図である。図3において、制御切替判定部15は、電圧振幅算出器15aと、制御切替判定器15bとを有し、モータの現在の制御領域が通常制御領域か電圧飽和領域かを以下のようにして判定する。
図4は、本開示の実施例1の補正トルク生成器の構成例を示す図である。図4において、補正トルク生成器34は、速度変動成分分離器34aと、速度変動振幅算出器34bと、減算器34cと、補正トルク振幅算出器34dと、速度変動位相修正器34eと、直交成分分離器34fと、補正トルク復調器34gとを有する。
図5は、本開示の実施例1の電圧飽和領域電圧指令値生成器の構成例を示す図である。図5において、電圧飽和領域電圧指令値生成器14bは、出力電圧制限指令値生成器14b1と、誘起電圧指令値算出器14b2と、電流指令値算出器14b3と、仮電圧指令値算出器14b4と、電圧ベクトル角算出器14b5と、電圧指令値算出器14b6とを有する。
図9は、本開示の実施例1のMTPI電圧振幅制限処理器の動作例の説明に供する図である。
図11は、本開示の実施例1の通常制御領域電圧指令値生成器の構成例を示す図である。図11において、通常制御領域電圧指令値生成器14aは、電流指令値算出器14a1と、加算器16,17,21,22と、減算器18,19と、電圧指令値算出器20と、IIR(Infinite Impulse Response)フィルタ35a,35bと、非干渉化制御器36と、電流誤差補正値生成器37とを有する。
図12は、本開示の実施例1の電流誤差補正値生成器の構成例を示す図である。図12において、電流誤差補正値生成器37は、減算器37a,37eと、q軸電流誤差成分分離器37bと、q軸電流誤差積算器37cと、q軸電流誤差補正値復調器37dと、d軸電流誤差成分分離器37fと、d軸電流誤差積算器37gと、d軸電流誤差補正値復調器37hとを有する。
本開示の実施例1は、モータMの誘起電圧の歪みが小さい場合には有効である。しかし、実際の誘起電圧波形にはモータMの構造に起因した誘起電圧歪みがあり、この誘起電圧歪みが大きい場合、モータMの電流の高調波が大きくなる。その結果、モータMの電流制御の追従性が悪くなる場合がある。そこで、実施例2では、モータMの電流制御の追従性を高めて、制御の安定性をさらに向上させる。このため、実施例2では、電圧飽和領域電圧指令値生成器14bの構成が実施例1と一部異なる。
図13は、本開示の実施例2の電圧飽和領域電圧指令値生成器14bの構成例を示す図である。図13において、電圧飽和領域電圧指令値生成器14bが、出力電圧制限指令値生成器14b1と、誘起電圧指令値算出器14b2と、電流指令値算出器14b3と、仮電圧指令値算出器14b4と、電圧ベクトル角算出器14b5と、電圧指令値算出器14b6とを有する点は、実施例1(図5)と同一である。実施例2では、電圧飽和領域電圧指令値生成器14bは、さらに電流比例制御器14b7と、加算器14b8,14b9と、減算器14b10,14b11とを有する。
14 電圧指令値生成器
14a 通常制御領域電圧指令値生成器
14b 電圧飽和領域電圧指令値生成器
14b1 出力電圧制限指令値生成器
14b2 誘起電圧指令値算出器
14b3 電流指令値算出器
14b4 仮電圧指令値算出器
14b5 電圧ベクトル角算出器
14b6 電圧指令値算出器
14b7 電流比例制御器
Claims (4)
- 速度指令値とモータの速度とに基づいたトルク指令値から前記モータの電圧指令値を生成する電圧指令値生成器と、
前記モータの制御領域が電圧飽和領域にあるか否かを判定する制御切替判定部と、を具備し、
前記電圧指令値生成器は、前記制御切替判定部により前記制御領域が前記電圧飽和領域にあると判定された場合には、
前記トルク指令値と、前記モータへ出力可能な最大電圧の限界値と、前記モータに印加される出力電圧の電圧ベクトル角とに基づいて前記電圧指令値を生成する、
モータ制御装置。 - 前記電圧指令値生成器は、前記トルク指令値の前記モータの定トルク曲線と前記モータの定誘起電圧楕円との交点に基づいて前記モータの電流指令値を算出し、前記電流指令値から電圧ベクトル角を求めて電圧指令値を生成する、
請求項1に記載のモータ制御装置。 - 前記電圧指令値生成器は、前記電流指令値に基づいてモータモデル式に従って仮電圧指令値を算出し、前記仮電圧指令値に基づいて前記電圧ベクトル角を算出する、
請求項2に記載のモータ制御装置。 - 前記電圧指令値生成器は、前記モータの電流指令値と前記モータの電流との誤差に比例定数を乗算することにより補償電圧を算出し、前記補償電圧が加算された前記電圧指令値を生成する、
請求項1に記載のモータ制御装置。
Priority Applications (4)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
EP20916873.1A EP4099559A4 (en) | 2020-01-31 | 2020-12-09 | ENGINE CONTROL DEVICE |
US17/789,921 US12028002B2 (en) | 2020-01-31 | 2020-12-09 | Motor control device |
CN202080093047.0A CN114982123A (zh) | 2020-01-31 | 2020-12-09 | 电动机控制装置 |
AU2020426771A AU2020426771B2 (en) | 2020-01-31 | 2020-12-09 | Motor control device |
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2020-015699 | 2020-01-31 | ||
JP2020015699A JP6852823B1 (ja) | 2020-01-31 | 2020-01-31 | モータ制御装置 |
Publications (1)
Publication Number | Publication Date |
---|---|
WO2021153020A1 true WO2021153020A1 (ja) | 2021-08-05 |
Family
ID=75154785
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
PCT/JP2020/045886 WO2021153020A1 (ja) | 2020-01-31 | 2020-12-09 | モータ制御装置 |
Country Status (6)
Country | Link |
---|---|
US (1) | US12028002B2 (ja) |
EP (1) | EP4099559A4 (ja) |
JP (1) | JP6852823B1 (ja) |
CN (1) | CN114982123A (ja) |
AU (1) | AU2020426771B2 (ja) |
WO (1) | WO2021153020A1 (ja) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2023195273A1 (ja) * | 2022-04-06 | 2023-10-12 | 日立Astemo株式会社 | モータ制御装置並びに電気車 |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US11926221B2 (en) * | 2020-09-24 | 2024-03-12 | GM Global Technology Operations LLC | Open-loop control for transient operation of a rotary electric machine |
CN115913018A (zh) * | 2022-11-04 | 2023-04-04 | 珠海市伟高变频科技有限公司 | 一种基于转速控制的单转子压缩机振动抑制方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2017139892A (ja) * | 2016-02-04 | 2017-08-10 | 株式会社明電舎 | 埋込磁石同期モータの電流指令テーブル自動生成システムおよび電流指令テーブル自動生成方法 |
JP2017158415A (ja) | 2016-03-04 | 2017-09-07 | 株式会社富士通ゼネラル | モータ制御装置 |
JP2017158414A (ja) | 2016-03-04 | 2017-09-07 | 株式会社富士通ゼネラル | モータ制御装置 |
JP2019180173A (ja) * | 2018-03-30 | 2019-10-17 | 株式会社富士通ゼネラル | モータ制御装置 |
JP2019201471A (ja) * | 2018-05-15 | 2019-11-21 | 株式会社富士通ゼネラル | モータ制御装置 |
JP2019213309A (ja) * | 2018-06-01 | 2019-12-12 | 日産自動車株式会社 | 巻線界磁型同期モータの制御方法、及び、制御装置 |
Family Cites Families (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CA2708054C (en) * | 2007-12-04 | 2013-07-09 | Mitsubishi Electric Corporation | Control apparatus of alternating-current motor |
US9007004B2 (en) * | 2009-11-06 | 2015-04-14 | University Of Technology, Sydney | Sensorless AC motor controller |
KR101562419B1 (ko) * | 2011-07-05 | 2015-10-22 | 엘에스산전 주식회사 | 매입형 영구자석 동기 전동기의 구동장치 |
US9219439B2 (en) | 2011-12-09 | 2015-12-22 | Panasonic Intellectual Property Management Co., Ltd. | Electric motor control device |
JP6260502B2 (ja) | 2014-09-16 | 2018-01-17 | 株式会社デンソー | モータ制御装置 |
CN106961232B (zh) * | 2017-04-17 | 2019-12-06 | 上海大学 | 采用小容量薄膜电容的高功率因数二极管整流器永磁电机控制方法 |
JP6867267B2 (ja) * | 2017-10-13 | 2021-04-28 | ルネサスエレクトロニクス株式会社 | モータ制御装置およびモータシステム |
-
2020
- 2020-01-31 JP JP2020015699A patent/JP6852823B1/ja active Active
- 2020-12-09 EP EP20916873.1A patent/EP4099559A4/en active Pending
- 2020-12-09 CN CN202080093047.0A patent/CN114982123A/zh active Pending
- 2020-12-09 US US17/789,921 patent/US12028002B2/en active Active
- 2020-12-09 WO PCT/JP2020/045886 patent/WO2021153020A1/ja unknown
- 2020-12-09 AU AU2020426771A patent/AU2020426771B2/en active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2017139892A (ja) * | 2016-02-04 | 2017-08-10 | 株式会社明電舎 | 埋込磁石同期モータの電流指令テーブル自動生成システムおよび電流指令テーブル自動生成方法 |
JP2017158415A (ja) | 2016-03-04 | 2017-09-07 | 株式会社富士通ゼネラル | モータ制御装置 |
JP2017158414A (ja) | 2016-03-04 | 2017-09-07 | 株式会社富士通ゼネラル | モータ制御装置 |
JP2019180173A (ja) * | 2018-03-30 | 2019-10-17 | 株式会社富士通ゼネラル | モータ制御装置 |
JP2019201471A (ja) * | 2018-05-15 | 2019-11-21 | 株式会社富士通ゼネラル | モータ制御装置 |
JP2019213309A (ja) * | 2018-06-01 | 2019-12-12 | 日産自動車株式会社 | 巻線界磁型同期モータの制御方法、及び、制御装置 |
Non-Patent Citations (1)
Title |
---|
See also references of EP4099559A4 |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2023195273A1 (ja) * | 2022-04-06 | 2023-10-12 | 日立Astemo株式会社 | モータ制御装置並びに電気車 |
Also Published As
Publication number | Publication date |
---|---|
CN114982123A (zh) | 2022-08-30 |
US12028002B2 (en) | 2024-07-02 |
EP4099559A1 (en) | 2022-12-07 |
JP2021125898A (ja) | 2021-08-30 |
US20220399843A1 (en) | 2022-12-15 |
EP4099559A4 (en) | 2024-02-28 |
AU2020426771A1 (en) | 2022-07-21 |
AU2020426771B2 (en) | 2023-08-03 |
JP6852823B1 (ja) | 2021-03-31 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
WO2021153020A1 (ja) | モータ制御装置 | |
EP2034605B1 (en) | Electric motor driving device, and compressor driving device | |
JP3840905B2 (ja) | 同期電動機の駆動装置 | |
JP5413400B2 (ja) | 交流電動機の制御装置 | |
JP7225564B2 (ja) | モータ制御装置 | |
JP6672902B2 (ja) | モータ制御装置 | |
JPH1118496A (ja) | 電気車の制御装置および制御方法 | |
US20020171387A1 (en) | Apparatus and method for controlling permanent magnet electric machines | |
JP7081274B2 (ja) | モータ制御装置 | |
JP7020093B2 (ja) | モータ制御装置 | |
JP7363524B2 (ja) | センサレスモータ制御装置 | |
JP7318547B2 (ja) | モータ制御装置 | |
JP6769050B2 (ja) | モータ制御装置 | |
JP7020112B2 (ja) | モータ制御装置 | |
JP6984672B2 (ja) | モータ制御装置 | |
JP7009861B2 (ja) | モータ制御装置 | |
JP7127657B2 (ja) | モータ制御装置 | |
JP2023051559A (ja) | モータ制御装置 | |
JP2012175776A (ja) | モータ制御装置及びモータ駆動システム | |
WO2023013302A1 (ja) | 交流電動機の駆動制御装置および駆動制御方法 | |
JP5511531B2 (ja) | 同期電動機の制御装置 | |
WO2023195273A1 (ja) | モータ制御装置並びに電気車 | |
JP2023170542A (ja) | 回転電機の制御方法、及び、回転電機の制御装置 | |
JP2023047965A (ja) | モータ制御装置 | |
JP2023051557A (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: 20916873 Country of ref document: EP Kind code of ref document: A1 |
|
ENP | Entry into the national phase |
Ref document number: 2020426771 Country of ref document: AU Date of ref document: 20201209 Kind code of ref document: A |
|
NENP | Non-entry into the national phase |
Ref country code: DE |
|
ENP | Entry into the national phase |
Ref document number: 2020916873 Country of ref document: EP Effective date: 20220831 |