The present invention relates to a control system for a motor, and more particularly to a control system for a motor that determines a ripple compensating current.
The output torque of a permanent magnet synchronous motor (PMSM) (either a surface permanent magnet (SPM) or an interior permanent magnet (IPM) motor) may be determined by a voltage command and a phase advance angle. A specific output torque of the PMSM is determined by first selecting a specific quadrature axis (also referred to as the q-axis) reference current and a direct axis (also referred to as the d-axis) reference current, and then determining the voltage command and the phase advance angle based on the selected quadrature axis reference current and the direct axis reference current.
Interior permanent magnet synchronous motors (IPMSMs) are relatively inexpensive, but typically produce a relatively high torque ripple during operation. Electric motors utilized in electric power steering (EPS) applications are generally required to produce relatively low torque ripple. Thus, the torque ripple produced by an IPMSM or an SPM may need to be reduced before being used in an EPS application. However, the effects of torque ripple compensation at relatively high motor velocities may be diminished.
In one embodiment, a motor control system for determining a ripple compensating current is provided. The motor control system includes a motor having a plurality of motor harmonics and a motor frequency, and a bandwidth compensation controller that is in communication with the motor. The bandwidth compensation controller is configured to determine a magnitude response compensation value and a phase compensation value for each of the plurality of motor harmonics. The magnitude response compensation value and the phase compensation value are both based on the motor frequency. The bandwidth compensation controller is configured to determine the ripple compensating current based on the magnitude response compensation value and the phase compensation value.
In another embodiment, a method of determining a ripple compensating current for a motor having a plurality of motor harmonics is provided. The method includes determining a magnitude response compensation value and a phase compensation value for each of the plurality of motor harmonics by a bandwidth compensation controller. The magnitude response compensation value and the phase compensation value are both based on a motor frequency. The method also includes determining the ripple compensating current based on the magnitude response compensation value and the phase compensation value.
These and other advantages and features will become more apparent from the following description taken in conjunction with the drawings.
The subject matter which is regarded as the invention is particularly pointed out and distinctly claimed in the claims at the conclusion of the specification. The foregoing and other features, and advantages of the invention are apparent from the following detailed description taken in conjunction with the accompanying drawings in which:
Referring now to the Figures, where the invention will be described with reference to specific embodiments, without limiting same,
In the embodiment as shown in
For feedback control purposes, the phase currents ia and ib transmitted to the motor 32 by lines 50 and 52 may be detected to determine the instantaneous current flow to the motor 32. Specifically, a transducer 56 may be used to monitor the phase current ia on the line 50, and a transducer 58 may be used to monitor the phase current ib on the line 52. It should be noted that although transducer 56 and transducer 58 are illustrated, only one of the lines 50 or 52 may be monitored to measure either phase current ia or phase current ib. A control signal 60 representing the measured phase current ia may be sent to the a-axis current amplifier 40 from the transducer 56, and a control signal 62 representing the measured phase current ib may be sent to the b-axis current amplifier 42 from the transducer 58. An augmented or amplified value of the phase current ia is then sent to the a-axis ADC 44 from the a-axis current amplifier 40, and an amplified value of the phase current ib 62 is sent to the b-axis ADC 46 from the b-axis current amplifier 42. The a-axis ADC 44 converts the amplified value of the phase current ia into a digital value 64. The digital value 64 represent the magnitude of the phase current ia. The b-axis ADC 46 converts the amplified value of the phase current ib into a digital value 66. The digital value 66 represents the magnitude of the phase current ib.
The transform controller 38 receives as input the digital value 64 from the ADC 44 and the digital value 66 from the ADC 46. In one embodiment, the transform controller 38 is a three-phase to two-phase transformation controller where measured values for the AC current (e.g., the digital value 64 representing the phase current ia and the digital value 66 representing the phase current ib) are converted into equivalent measured currents, which are a measured d-axis current IdMEASURED and a measured q-axis current IqMEASURED. The measured d-axis current IdMEASURED is sent to a subtractor 70 and the measured q-axis current IqMEASURED is sent to a subtractor 72.
The command current controller 20 receives as input a torque reference command Te, an angular speed ωm, and the control signal 53 representing the bridge voltage Vecu from the transducer 51. The torque reference command Te represents a commanded torque value, and may be derived from another controller (not shown), or may correspond to a torque value generated by an operator. The angular speed ωm, is measured by the speed sensor 36. The speed sensor 36 may include, for example, an encoder and a speed calculation circuit for calculating the angular speed of a rotor (not shown) of the motor 32 based on a signal received by the encoder. The command current controller 20 calculates a reference d-axis current Id
The torque ripple compensation controller 48 receives as input the reference d-axis current Id
The adder 69 receives the reference d-axis current Id
The adder 71 receives the reference q-axis current Iq
The d-axis PI controller 22 receives as input the d-axis error signal 74 from the subtractor 70. The d-axis PI controller 22 calculates a d-axis voltage signal VD. The d-axis voltage signal VD is based on a d-axis proportional gain KP, and a d-axis integral gain Ki. Likewise, the q-axis PI controller 23 receives as input the q-axis error signal 76 from the subtractor 72. The q-axis PI controller 23 calculates a q-axis voltage signal VQ. The q-axis voltage signal VQ is based on a q-axis proportional gain KP, and a q-axis integral gain Ki.
The polar conversion controller 24 receives as input the d-axis voltage signal VD from the d-axis PI controller 22 and the q-axis voltage signal VQ from the q-axis PI controller 23. Based on the inputs, the polar conversion controller 24 determines a voltage command Vcmd and a phase advance angle δ. The PWM inverter controller 26 receives as inputs the voltage command Vcmd and the phase advance angle δ from the polar conversion controller 24. The PWM inverter controller 26 also receives the rotor angle value θr measured by the motor position sensor 34. In one exemplary embodiment, the PWM inverter controller 26 may include an over-modulation space vector PWM unit to generate three respective duty cycle values Da, Db, and Dc. The duty cycle values Da, Db, and Dc are used to drive gate drive circuits (not shown) of the inverter 28 that energize phases of the of the motor 32.
Determining the q-axis ripple compensating current Iq
The torque ripple compensation module 80 receives as input the reference d-axis current Id
The torque ripple compensation module 80 may first determine the torque ripple magnitude 82 and the torque ripple phase 83 for each motor harmonic using the compensation look-up tables 84. The compensation look-up tables 84 each contain specific values for either the torque ripple magnitude 82 or the torque ripple phase 83 for a specific motor harmonic. The values for the torque ripple magnitude 82 or the torque ripple phase 83 in each of the compensation look-up tables 84 are based on the reference d-axis current Id
Once the values for the torque ripple magnitude 82 and the torque ripple phase 83 for a specific motor harmonic are found in the compensation look-up tables 84, the torque ripple compensation module 80 may then calculate a current harmonic q-axis magnitude 86 and a current q-axis phase output 88 for the specific motor harmonic. Alternatively, or in addition to the current harmonic q-axis magnitude 86 and the current q-axis phase output 88, the torque ripple compensation module 80 may calculate a current harmonic d-axis magnitude 90 and a current d-axis phase output 92 for the specific motor harmonic. In one embodiment, the current harmonic q-axis magnitude 86 for a specific motor harmonic may be determined by either Equations 1 or 2, the current q-axis phase output 88 for a specific motor harmonic may be determined by Equation 3, the current harmonic d-axis magnitude 90 for a specific motor harmonic may be determined by Equation 4, and the current d-axis phase output 92 for a specific motor harmonic may be determined by Equation 5:
where MtrCurrQax_iMag is the current harmonic q-axis magnitude 86, t_MtrTrq_ripple_iMag(n) is the torque ripple magnitude 82 found in one of the compensation look-up tables 84 for an ith harmonic (where i represents a specific harmonic such as, for example, the eighth harmonic), MtrCurrQax_iPhase is the current q-axis phase output 88, t_MtrTrq_ripple_iPh(n) is the torque ripple phase 83 found in one of the compensation look-up tables 84 for the ith harmonic, MtrCurrDax_iMag is the current harmonic d-axis magnitude 90, and MtrCurrDax_iPhase is the current d-axis phase output 92.
The torque ripple compensation module 80 calculates the current harmonic q-axis magnitude 86 and the current q-axis phase output 88 for each motor harmonic. For example, as seen in
The bandwidth compensation module 49 includes a magnitude compensation module 91, a phase compensation module 93, and a ripple compensation block 94. The magnitude compensation module 91 receives as input from the torque ripple compensation module 80 the m number of values of the current harmonic q-axis magnitude 86. Alternatively, or in addition to the m number of values of the current harmonic q-axis magnitude 86, the magnitude compensation module 91 receives as input from the torque ripple compensation module 80 the m number of values of the current harmonic d-axis magnitude 90. The magnitude compensation module 91 also receives as input the angular speed ωm, (from the speed sensor 36 shown in
The magnitude compensation module 91 may then determine a PI compensating magnitude based on the frequency compensation value n. The PI compensating magnitude is a unitless gain value for a specific motor harmonic. In one embodiment, the magnitude compensation module 91 includes a magnitude look-up table 95, where a specific value of the PI compensating magnitude is selected based on the frequency compensation value n. The magnitude compensation module 91 may then determine the magnitude response compensation value 96 using Equations 6 and 7, and the d-axis magnitude response compensation value 98 using Equations 6 and 8:
MtrCurrPI—iMag_Qax=t_MtrCurrPIComp—iMag(n)*MtrCurrQax_IMag Equation 7
MtrCurrPI—iMag_Dax=t_MtrCurrPIComp—iMag(n)*MtrCurrDax—iMag Equation 8
where MtrCurrPI_iMag_Qax is the q-axis magnitude response compensation value 96, MtrCurrQax_iMag is the current harmonic q-axis magnitude 86 for an ith harmonic, t_MtrCurrPIComp_iMag(n) is the PI compensating magnitude an ith harmonic, MtrCurrPI_iMag_Dax is the d-axis magnitude response compensation value 98, MtrCurrDax_iMag is the current harmonic d-axis magnitude 90 for an ith harmonic.
The magnitude compensation module 91 calculates the q-axis magnitude response compensation value 96, the d-axis magnitude response compensation value 98, or both the q-axis magnitude response compensation value 96 and the d-axis magnitude response compensation value 98 for each motor harmonic (e.g., the m number of motor harmonics).
The phase compensation module 93 receives as input from the torque ripple compensation module 80 the m number of values for the current q-axis phase output 88 and the m number of values for the torque ripple phase 83. Alternatively, or in addition to the m number of values of the current q-axis phase output 88, the phase compensation module 93 receives as input from the torque ripple compensation module 80 the m number of values for the current d-axis phase output 92. The phase compensation module 93 also receives as input the angular speed ωm. The phase compensation module 93 determines a q-axis phase compensation value 100 for an ith harmonic. Alternatively, or in addition to the q-axis phase compensation value 100, the phase compensation module 93 determines a d-axis phase compensation value 102 for an ith harmonic. Specifically, the phase compensation module 93 first determines a frequency compensation value n for an ith harmonic using Equation 6, which is described above.
The phase compensation module 93 may then determine a PI compensating phase based on the frequency compensation value n. The PI compensating phase is a phase in degrees for a specific motor harmonic. In one embodiment, the phase compensation module 93 includes a phase look-up table 97, where a specific value of the PI compensating phase is selected based on the frequency compensation value n. The phase compensation module 93 may then determine the q-axis phase compensation value 100 using Equations 8, and the d-axis phase compensation value 102 using Equation 9:
MtrCurrPI—iPh_Qax=t_MtrCurrPIComp—iPh(n)+MtrCurrQax_IPh Equation 8
MtrCurrPI—iPh_Dax=t_MtrCurrPIComp—iPh(n)+MtrCurrDax_IPh Equation 9
where MtrCurrPI_iPh_Qax is the q-axis phase compensation value 100, MtrCurrQax_IPh is the q-axis current harmonic q-axis phase 88 for an ith harmonic, t_MtrCurrPIComp_iPh(n) is the PI compensating phase an ith harmonic, MtrCurrPI_iPh_Dax is the d-axis phase compensation value 102, MtrCurrDax_IPh is the current d-axis phase output 92 for an ith harmonic.
The phase compensation module 93 calculates the q-axis phase response compensation value 100, the d-axis phase response compensation value 102, or both the q-axis phase response compensation value 100 and the d-axis phase response compensation value 102 for each motor harmonic (e.g., the m number of motor harmonics).
The ripple compensation block 94 receives as input the m number of values of the magnitude response compensation value 96, the d-axis magnitude response compensation value 98, or both, from the magnitude compensation module 91. The ripple compensation block 94 also receives as input the m number of values of the q-axis phase response compensation value 100, the d-axis phase response compensation value 102, or both the q-axis phase response compensation value 100 and the d-axis phase response compensation value 102 from the from the phase compensation module 93. The ripple compensation module 94 also receives the rotor angle value θr measured by the motor position sensor 34. The ripple compensation block 94 determines the q-axis ripple compensating current Iq
The ripple compensation block 94 determines the d-axis ripple compensating current Id
The q-axis ripple compensating current Iq
In block 304, the torque ripple compensation controller 48 determines the torque ripple magnitude and the torque ripple phase for each motor harmonic using the compensation look-up tables 84 based on the values of the reference d-axis current Id
In block 306, the torque ripple compensation controller 48 may then calculate the current harmonic q-axis magnitude 86 and the current q-axis phase output 88 for each motor harmonic using Equations 1-3. Alternatively, or in addition to the current harmonic q-axis magnitude 86 and the current q-axis phase output 88, the torque ripple compensation controller 48 determines the current harmonic d-axis magnitude 90 and the current d-axis phase output 92 for each motor harmonic using Equations 4-5. Method 300 may then proceed to block 308.
In block 308, the bandwidth compensation module 49 receives as input the m number of values of the current harmonic q-axis magnitude 86, the q-axis phase output 88 the rotor angle value θr measured by the motor position sensor 34, and the angular speed ωm. Alternatively, or in addition to the m number of values of the current harmonic q-axis magnitude 86 and the q-axis phase output 88, the bandwidth compensation module 49 receives as input the m number of values of the current harmonic d-axis magnitude 90 and the d-axis phase output 92. Method 300 may then proceed to block 310.
In block 310, the magnitude compensation module 91 determines the q-axis magnitude response compensation value 96 for an ith harmonic. Alternatively, or in addition to the q-axis magnitude response compensation value 96, the magnitude compensation module 91 determines the d-axis magnitude response compensation value 98 for an ith harmonic. In one embodiment, the magnitude response compensation value 96 is determined using Equations 6 and 7 (described above), and the d-axis magnitude response compensation value 98 is determined using Equations 6 and 8 (described above). Method 300 may then proceed to block 312.
In block 312, the phase compensation module 93 determines the q-axis phase compensation value 100 for an ith harmonic. Alternatively, or in addition to the q-axis phase compensation value 100, the phase compensation module 93 determines the d-axis phase compensation value 102 for an ith harmonic. In one embodiment, the q-axis phase compensation value 100 is determined using Equations 6 and 8 (described above), and the d-axis phase compensation value 102 is determined using Equation 6 and Equation 9 (described above). Method 300 may then proceed to block 314.
In block 314, the ripple compensation block 94 determines the q-axis ripple compensating current Iq
The bandwidth compensation controller 49 as described above determines the q-axis ripple compensating current Iq
While the invention has been described in detail in connection with only a limited number of embodiments, it should be readily understood that the invention is not limited to such disclosed embodiments. Rather, the invention can be modified to incorporate any number of variations, alterations, substitutions or equivalent arrangements not heretofore described, but which are commensurate with the spirit and scope of the invention. Additionally, while various embodiments of the invention have been described, it is to be understood that aspects of the invention may include only some of the described embodiments. Accordingly, the invention is not to be seen as limited by the foregoing description.
Number | Name | Date | Kind |
---|---|---|---|
4713596 | Bose | Dec 1987 | A |
4773149 | Kip et al. | Sep 1988 | A |
5196778 | Hayashida | Mar 1993 | A |
5223775 | Mongeau | Jun 1993 | A |
5410234 | Shibata et al. | Apr 1995 | A |
5652495 | Narazaki et al. | Jul 1997 | A |
5962999 | Nakamura et al. | Oct 1999 | A |
6152254 | Phillips | Nov 2000 | A |
6161068 | Kurishige et al. | Dec 2000 | A |
6222334 | Tamagawa et al. | Apr 2001 | B1 |
6288515 | Hiti et al. | Sep 2001 | B1 |
6370459 | Phillips | Apr 2002 | B1 |
6465975 | Naidu | Oct 2002 | B1 |
6499559 | Mc Cann et al. | Dec 2002 | B2 |
6605912 | Bharadwaj et al. | Aug 2003 | B1 |
6700342 | Hampo et al. | Mar 2004 | B2 |
6900607 | Kleinau et al. | May 2005 | B2 |
7034493 | Yoshimoto et al. | Apr 2006 | B2 |
7071649 | Shafer et al. | Jul 2006 | B2 |
7145310 | Ihm et al. | Dec 2006 | B2 |
7199549 | Kleinau et al. | Apr 2007 | B2 |
7207412 | Uryu | Apr 2007 | B2 |
7394214 | Endo et al. | Jul 2008 | B2 |
7576506 | Kleinau et al. | Aug 2009 | B2 |
7952308 | Schulz et al. | May 2011 | B2 |
8633766 | Khlat et al. | Jan 2014 | B2 |
8896244 | Kleinau | Nov 2014 | B2 |
20020175649 | Reutlinger | Nov 2002 | A1 |
20030146041 | Kanda | Aug 2003 | A1 |
20060100766 | Schwarz et al. | May 2006 | A1 |
20070043490 | Yokota et al. | Feb 2007 | A1 |
20070046126 | Sagoo et al. | Mar 2007 | A1 |
20070103105 | Endo et al. | May 2007 | A1 |
20080167779 | Suzuki | Jul 2008 | A1 |
20090026994 | Namuduri et al. | Jan 2009 | A1 |
20090027000 | Gallegos-Lopez et al. | Jan 2009 | A1 |
20090069979 | Yamashita et al. | Mar 2009 | A1 |
20090114470 | Shimizu et al. | May 2009 | A1 |
20090234538 | Ta et al. | Sep 2009 | A1 |
20090267555 | Schulz et al. | Oct 2009 | A1 |
20100153162 | Tam et al. | Jun 2010 | A1 |
20110175556 | Tobari et al. | Jul 2011 | A1 |
20120112549 | Perisic et al. | May 2012 | A1 |
20120221208 | Kojo et al. | Aug 2012 | A1 |
20130261896 | Gebregergis et al. | Oct 2013 | A1 |
20140191699 | Dixon | Jul 2014 | A1 |
20140239860 | Kleinau | Aug 2014 | A1 |
20140265961 | Gebregergis et al. | Sep 2014 | A1 |
20140375239 | Kim et al. | Dec 2014 | A1 |
Number | Date | Country |
---|---|---|
1675099 | Sep 2005 | CN |
1741368 | Mar 2006 | CN |
101218146 | Jul 2008 | CN |
101399516 | Apr 2009 | CN |
101456429 | Jun 2009 | CN |
101981804 | Feb 2011 | CN |
2003010 | Dec 2008 | EP |
2000108916 | Apr 2000 | JP |
2001247049 | Sep 2001 | JP |
Entry |
---|
F. Briz, et al., “Analysis and Design of Current Regulators Using Complex Vectors”, IEEE Industry Applications Society Annual Meeting, New Orleans, Louisiana; Oct. 5-9, 1997, pp. 1504-1511. |
European Search Report for related European Application No. 12196930.7, dated: Mar. 22, 2013, 7 pages. |
L Harnefors, et al., “Model-Based Current Control of AC Machines Using the Internal Model Control Method”, IEEE Transactions on Industry Applications, vol. 34, No. 1, Jan./Feb. 1998, pp. 133-141. |
J. Kirtley, “6.061 Introduction to Electric Power Systems, Class Notes Chapter 12 Permanent Magnet ‘Brushless DC’” Motors, Massachusetts Institute of Technology, Department of Electrical Engineering and Computer Science, Spring 2011. |
Chinese Office Action for Chinese Patent Application No. 201310104183.7 issued on Jan. 6, 2015. |
Number | Date | Country | |
---|---|---|---|
20140265962 A1 | Sep 2014 | US |