WO2013024780A1 - Ｐｍモータの位置センサレス制御装置 - Google Patents

Ｐｍモータの位置センサレス制御装置 Download PDF

Info

Publication number
WO2013024780A1
WO2013024780A1 PCT/JP2012/070312 JP2012070312W WO2013024780A1 WO 2013024780 A1 WO2013024780 A1 WO 2013024780A1 JP 2012070312 W JP2012070312 W JP 2012070312W WO 2013024780 A1 WO2013024780 A1 WO 2013024780A1
Authority
WO
WIPO (PCT)
Prior art keywords
current
speed
sensorless control
voltage
coordinate conversion
Prior art date
Application number
PCT/JP2012/070312
Other languages
English (en)
French (fr)
Inventor
昌司 滝口
夢樹 小野
Original Assignee
株式会社明電舎
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by 株式会社明電舎 filed Critical 株式会社明電舎
Priority to RU2014109115/07A priority Critical patent/RU2560090C1/ru
Priority to CN201280039669.0A priority patent/CN103733504B/zh
Priority to JP2013528991A priority patent/JP5733404B2/ja
Publication of WO2013024780A1 publication Critical patent/WO2013024780A1/ja

Links

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
    • H02P6/00—Arrangements for controlling synchronous motors or other dynamo-electric motors using electronic commutation dependent on the rotor position; Electronic commutators therefor
    • H02P6/14—Electronic commutators
    • H02P6/16—Circuit arrangements for detecting position
    • H02P6/18—Circuit arrangements for detecting position without separate position detecting elements
    • 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/24—Vector control not involving the use of rotor position or rotor speed sensors
    • 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
    • H02P2203/00—Indexing scheme relating to controlling arrangements characterised by the means for detecting the position of the rotor
    • H02P2203/09—Motor speed determination based on the current and/or voltage without using a tachogenerator or a physical encoder
    • 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
    • H02P2207/00—Indexing scheme relating to controlling arrangements characterised by the type of motor
    • H02P2207/05—Synchronous machines, e.g. with permanent magnets or DC excitation
    • 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
    • H02P27/08—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 with pulse width modulation

Definitions

  • the present invention relates to a position sensorless control device for a PM motor.
  • a PM motor When a PM motor is operated at a low speed, an induced voltage is estimated using differential information for a current, and a rotational speed and a magnetic pole position can be estimated using the estimated voltage.
  • This relates to a sensorless control device.
  • PM motors synchronous motors
  • sensors that detect this position information often have built-in electronic components such as semiconductor elements, and small ones have problems in environmental resistance and durability, such as low mechanical strength. Yes. Therefore, in applications where high position control accuracy and responsiveness are not required, reliability is improved over control performance by applying a position sensorless control method that estimates the magnetic pole position from inverter voltage and current information without using a position sensor. A lot of research has been done to date.
  • This position sensorless control method is roughly divided into the following two types.
  • One is a method for estimating the speed electromotive force component generated by the magnetic flux of the field magnet.
  • This method estimates the speed electromotive force due to the field magnetic flux of the motor from the fundamental wave component of the inverter output voltage and current.
  • PM motors can be classified in terms of non-salient pole machine characteristics in which the inductance components of the field axis (d axis) and its orthogonal axis (q axis) are equal, and salient pole machine characteristics in which the inductance of each axis is not equal.
  • this method can be applied to both, it can be used for PM motors in general.
  • the other is a method of estimating the magnetic salient pole axis while measuring the inductance by superimposing a high frequency component on the output voltage or current component of the inverter.
  • this method if the PM motor has a saliency with a difference in inductance between the d axis and the q axis, the phase of the field magnetic pole axis can be estimated by continuously measuring the inductance of each axis component. .
  • this method cannot be applied to a PM motor with non-salient pole machine characteristics, and the field magnetic pole axis (d-axis) has two types of polarities, N pole and S pole. It is also necessary to add an additional magnetic pole discrimination method using magnetic saturation to the control.
  • the former method of estimating the electromotive force cannot be accurately operated in the low speed region, and the latter method of superimposing the high frequency can operate in the low speed region. It cannot be applied to PM motors having mechanical characteristics or low magnetic saturation.
  • Non-Patent Document 1 As a low-speed position sensorless control method, and sensorless control is performed using a control block as shown in FIG.
  • the control method shown in Non-Patent Document 1 is basically a kind of electromotive force estimation method, and the principle of the estimation method is described in Non-Patent Document 2, but the estimation is a rotational speed estimation surrounded by a dotted line shown in FIG. Done in the department.
  • the induced voltage e is obtained by subtracting the voltage drop component of the resistance component of the winding and the armature reaction component, which is the product of the current command i ref and the motor constant Rs + j ⁇ 1 Ls, from the voltage command value Vs.
  • the estimated field magnetic pole phase is set as a control reference axis having the d axis, and this reference axis is regarded as a real axis.
  • the imaginary part component of the induced voltage e is obtained as the q-axis induced voltage eq, and the real part is obtained as the d-axis induced voltage ed.
  • the eq voltage component is corrected by multiplying the estimated angular velocity ⁇ 1 by the gain ⁇ s and ed, and this induced voltage is divided by the magnetic flux ⁇ m and used for velocity estimation.
  • the estimated angular velocity ⁇ 1 is estimated after applying a low-pass filter that limits the broadband so that it does not oscillate due to this loop. ing.
  • the feature of the method shown in FIG. 6 is that the estimated speed is automatically corrected to the normal rotation direction even if the forward rotation and the reverse rotation are erroneously estimated in the reverse rotation direction near zero speed.
  • voltage information for example, in the method of substituting the voltage command Vs input to the conventional PWM modulation, if measures against voltage error due to dead time error or switching delay of the switching element are not properly applied, rotation is performed. Even if there is a direction estimation error, it cannot be detected immediately, and it cannot return to the normal rotation direction unless a certain reverse rotation speed is reached.
  • the purpose of the present invention is to extend the lower limit of the speed at which it can return to the normal rotation direction to a lower speed even if there is a voltage error between the voltage command and the actual voltage.
  • it is to provide a position sensorless control device for a PM motor that enables accurate position estimation.
  • a current command is generated from the speed command and the estimated speed
  • a voltage command is generated from the current command and a current detection value of ⁇ - ⁇ coordinates detected through the rotation coordinate conversion unit
  • An induced voltage calculation unit for calculating the speed electromotive forces e ⁇ and e ⁇ by inputting the current detection values i ⁇ and i ⁇ and current differential information pi ⁇ and pi ⁇ in the ⁇ - ⁇ coordinates during the zero voltage vector period;
  • a rotational speed estimator that calculates the estimated speed ⁇ ⁇ using the speed electromotive forces e ⁇ and e ⁇ calculated by the induced voltage calculator, and an estimated phase ⁇ ⁇ by calculating the estimated speed over time. Is output to the rotation coordinate conversion unit and the reverse rotation coordinate conversion unit and used as a reference phase of the rotation coordinate.
  • the induced voltage calculation unit is characterized in that the speed electromotive forces e ⁇ and e ⁇ are calculated by the following equations.
  • the current differential detection Is provided on the detection current input side of the rotation coordinate conversion unit, the three-phase current detection generated during the zero voltage vector period of the PWM modulation pattern and the current differential amount is input to the rotation coordinate converter,
  • the induced voltage calculation unit calculates the speed electromotive force e ⁇ and e ⁇ by inputting the current detection values i ⁇ and i ⁇ and the current differential information pi ⁇ ′ and pi ⁇ ′ that have been subjected to the rotation coordinate conversion by the rotation coordinate conversion unit. It is characterized by doing.
  • the induced voltage calculation unit is characterized in that the speed electromotive forces e ⁇ and e ⁇ are calculated by the following equations.
  • the current differential information pi ⁇ ′, pi ⁇ ′ according to another aspect of the present invention is converted into current components i ⁇ , i ⁇ of orthogonal two axes ⁇ , ⁇ of fixed coordinates by a rotating coordinate conversion unit, and then a differential operation is performed. It is characterized by being obtained.
  • a voltage drop correction unit that corrects a voltage drop of a switching element is provided on an input side of the rotational coordinate conversion unit, and a three-phase generated during a zero voltage vector period of the PWM modulation pattern Input the voltage drop correction amount from the current detection and voltage drop correction unit to the rotating coordinate conversion unit, A speed electromotive force is calculated by inputting the detected current value, the current differential information, and the voltage drop correction amount to the induced voltage calculation unit.
  • the current differential information according to another aspect of the present invention is characterized in that it is based on either a rotating coordinate system or a fixed coordinate system.
  • the induced voltage calculation unit is characterized in that the speed electromotive forces e ⁇ and e ⁇ are calculated by the following equations.
  • vce ⁇ , vce ⁇ voltage correction amount
  • the induced voltage calculation unit according to another aspect of the present invention is characterized in that the speed electromotive forces e ⁇ and e ⁇ are calculated by the following equations.
  • the current differential information according to another aspect of the present invention is characterized in that the current differential information is input to the rotating coordinate conversion unit as a difference calculation value between two current sampling values before and after the carrier vertex in the zero voltage vector period.
  • the current differential information according to still another aspect of the present invention is characterized in that a value obtained by differentiating a moving average value of current sampling values between a carrier peak and a base in a zero voltage vector period is input to the rotating coordinate conversion unit. Is.
  • the voltage drop correction amount in the voltage drop correction unit according to another aspect of the present invention is obtained by using table data corresponding to the three-phase detection current value, and the table data is separately for each phase of the three phases and positive and negative currents. It is characterized by being set individually for each polarity.
  • the estimated phase ⁇ ⁇ is calculated using the current signals i ⁇ and i ⁇ in the ⁇ - ⁇ coordinates and the current differential information pi ⁇ and pi ⁇ during the zero voltage vector period. .
  • the block diagram of the control apparatus which shows embodiment of this invention.
  • the block diagram of the control apparatus which shows other embodiment of this invention.
  • the block diagram of the control apparatus which shows other embodiment of this invention.
  • Explanatory drawing of the current sampling of this invention The wave form diagram of the simulation result of this invention.
  • the current signal and current differential information detected by the induced voltage calculation unit are input to calculate the speed electromotive force, and the estimated magnetic pole phase is obtained based on this speed electromotive force.
  • FIG. 1 is a block diagram of a position sensorless control apparatus showing a first embodiment of the present invention.
  • 1 is an inverter controlled by PWM
  • 2 is a PM motor.
  • the N-pole field axis of the actual PM motor is defined as d-axis
  • the phase advanced by 90 ° in the forward rotation direction is defined as q-axis.
  • the imaginary N-pole axis based on the magnetic pole estimation is defined as the ⁇ -axis
  • the phase advanced by 90 ° in the rotation direction by the electrical angle is defined as ⁇ .
  • Reference numeral 3 denotes a rotating coordinate conversion unit which inputs the three-phase currents i u , i v and i w detected by the current sensor and converts the coordinates to i ⁇ and i ⁇ which are estimated axes.
  • the coordinate-converted current signals i ⁇ and i ⁇ are input to the current differential detection unit 4, the induced voltage calculation unit 5, and the current control unit 9, respectively.
  • the current differential detection unit 4 detects the changes pi ⁇ and pi ⁇ of the ⁇ - ⁇ axis current during the zero voltage vector period and inputs them to the induced voltage calculation unit 5.
  • the current signal input i gamma, i [delta] a current differential information pi gamma, speed electromotive force, as will be described later with reference to pi ⁇ e ⁇ . e Calculate ⁇ .
  • Reference numeral 6 denotes a rotational speed estimator
  • 7 denotes an integrator, which calculates an estimated magnetic pole position ⁇ ⁇ by integrating the estimated speed obtained by the rotational speed estimator 6 and outputs it to the rotational coordinate converter 3 and the reverse rotational coordinate converter 10.
  • 8 is a speed control unit
  • 9 is a current control unit.
  • the speed control unit 8 outputs the current command i ⁇ * of the ⁇ -axis component corresponding to the torque command from the input information of the speed command ⁇ * and the speed estimation ⁇ ⁇ .
  • This current command i ⁇ * and an arbitrary current command i ⁇ * on the ⁇ axis are input to the current control unit 9.
  • the current control unit 9 is supplied with current signals i ⁇ and i ⁇ converted into rotational coordinates via the rotational coordinate conversion unit 3, and the current control unit 9 receives these current commands i ⁇ * and i ⁇ *.
  • the current signals i ⁇ and i ⁇ are compared to perform feedback calculation, and voltage commands v ⁇ * and v ⁇ * on rotational coordinates ( ⁇ - ⁇ coordinates) based on the estimated magnetic pole axis are output.
  • the voltage commands v ⁇ * and v ⁇ * are obtained by performing reverse rotation conversion to fixed coordinates, two-phase three-phase conversion, and the like by the reverse rotation coordinate conversion unit 10 that performs the reverse operation of the rotation coordinate conversion unit 3.
  • the alternating current voltage commands v u , v v and v w are input to the inverter 1.
  • the three-phase voltage commands v u , v v , and v w are power-amplified by PWM modulation and a substantially equivalent voltage is output.
  • the above is a common control block whether or not there is a position sensor.
  • the current differential detection unit 4 includes a change amount pi ⁇ , ⁇ - ⁇ -axis current change during the zero voltage vector period. pi ⁇ is detected.
  • the ⁇ - ⁇ axis current can be obtained by performing rotational coordinate conversion from the three-phase current detection i u , i v , i w , but as a rotational coordinate conversion method, signal conversion is performed by an analog multiplier or an analog adder / subtractor.
  • the continuous signal method to be realized and the three-phase current detection i u , i v , i w are converted into a digital signal by an analog / digital converter (A / D converter) and then converted by a digital converter such as a CPU. There is a method using a discrete signal for executing the operation of the part.
  • any rotation coordinate conversion method may be used.
  • a signal after the rotation coordinate conversion is performed.
  • an analog differentiator By applying an analog differentiator to and sampling and holding the current differential component in the zero voltage vector period, the changes pi ⁇ and pi ⁇ are obtained.
  • an A / D converter when an A / D converter is used, current is sampled at a number of times during the zero voltage vector period and converted by a signal conversion unit, and a differential component is obtained by approximation from a plurality of discrete current values at a plurality of times. Obtainable.
  • Position sensorless control is configured using the current signals i ⁇ and i ⁇ obtained by either method and the current differential information (variations) pi ⁇ and pi ⁇ .
  • the voltage equation of the PM motor is expressed in an orthogonal coordinate system (dq coordinate) represented by an electrical angle with the N-pole axis of the actual machine as the d-axis.
  • R winding resistance
  • L winding inductance
  • ⁇ d magnetic flux linked to the stator winding by the magnetic flux generated by the magnet Component
  • ⁇ rotor magnetic pole position (electrical angle)
  • p differential operator (d / dt)
  • a PM motor having non-salient pole characteristics is a control target, so that the inductance components of the d-axis and the q-axis are regarded as being equal and expressed by a common coefficient.
  • equation (2) the axial error ⁇ e itself is assumed to have little time change, and the differential term of this axial error is ignored.
  • the current differential information pi ⁇ and pi ⁇ in the ⁇ - ⁇ coordinate in the zero voltage vector period is obtained by differentiating the detected currents i ⁇ and i ⁇ obtained by sequentially rotating the coordinates with the estimated phase ⁇ ⁇ . is there.
  • ⁇ 0 is a fixed value
  • ⁇ 1 and ⁇ 2 are positive gains
  • z ⁇ 1 is the previous sampling value
  • the rotational speed estimation unit 6 shown in FIG. 1 is estimated based on equations (4) and (5).
  • the speed ⁇ ⁇ is calculated, and the estimated phase ⁇ ⁇ is output via the integrator 7. That is, the induced voltage calculation unit 5 calculates the speed electromotive force in the zero voltage vector period from the current signals i ⁇ and i ⁇ in the ⁇ - ⁇ coordinate and the current differential information pi ⁇ and pi ⁇ in the zero voltage vector period according to the equation (3).
  • e ⁇ and e ⁇ are output.
  • the rotational speed estimator 6 performs calculations of equations (4) and (5) using the input speed electromotive forces e ⁇ and e ⁇ and outputs an estimated speed ⁇ ⁇ .
  • the integrator 7 calculates the estimated phase ⁇ ⁇ by time-integrating the estimated speed ⁇ ⁇ and outputs it to the rotating coordinate conversion unit 3 and the reverse rotating coordinate conversion unit 10 to be used for the reference phase of the rotating coordinate.
  • FIG. 2 shows the second embodiment.
  • the difference from the first embodiment shown in FIG. 1 is that the current differential detection unit 11 is provided in a fixed coordinate system and current differential information pi ′ ⁇ , pi ′. is that ⁇ is obtained.
  • the estimated phase ⁇ ⁇ used for rotational coordinate conversion in order to accurately obtain the current differential information, the estimated phase ⁇ ⁇ used for rotational coordinate conversion must also be a continuous value. Therefore, when the three-phase current signal is converted into a digital value by A / D conversion, the estimated phase ⁇ ⁇ used for the rotation coordinate conversion needs to be sequentially updated, which increases the amount of calculation.
  • the second embodiment takes this point into consideration, and it is not necessary to sequentially update the estimated phase ⁇ ⁇ by configuring to obtain the current differential information pi ′ ⁇ and pi ′ ⁇ in a fixed coordinate system. It is possible.
  • the current signals i ⁇ and i ⁇ and the current differential information pi ′ ⁇ and pi ′ ⁇ in the zero voltage vector period are sampled only at the time of calculating the speed electromotive forces e ⁇ and e ⁇ in the zero voltage vector period. Coordinate conversion need only be performed during this zero voltage vector period. As a result, the calculation time in the CPU or digital circuit for executing the calculation is shortened, and an advantage that the present invention can be applied even when the circuit operation is slow is obtained.
  • the three-phase current detections i u , i v , and i w are displayed so as to be directly differentiated.
  • an analog signal or a digital signal can be realized.
  • Any information corresponding to the current differential information pi ′ ⁇ , pi ′ ⁇ in the ⁇ - ⁇ coordinates necessary for obtaining the speed electromotive forces e ⁇ , e ⁇ during the zero voltage vector period may be used.
  • the three-phase two-phase conversion is first performed to convert the current signals i ⁇ and i ⁇ of the orthogonal two axes ( ⁇ coordinates) of the fixed coordinate system, and then the differential operation is performed, and then the rotation coordinate conversion is applied.
  • the operation order may be changed.
  • v ⁇ , v ⁇ ⁇ , ⁇ axis voltage, i ⁇ , i ⁇ : ⁇ , ⁇ axis current
  • R winding resistance
  • L Winding inductance
  • ⁇ d Interlinkage magnetic flux of magnet
  • ⁇ Rotor angular velocity (electrical angle)
  • ⁇ Rotor magnetic pole position
  • p Differential operator Since (1) was handled in the rotary coordinate system, The term ⁇ ⁇ L corresponding to the speed electromotive force exists in the impedance matrix of the first term on the right side due to the magnetic flux generated by the armature current, but this component is handled in the fixed coordinate system in equation (6). Does not exist, which makes the arithmetic expression easier.
  • t and t + ⁇ T indicate two current sampling times having a time difference ⁇ T in the same zero voltage vector period. Further, ⁇ i ⁇ and ⁇ i ⁇ indicate current differences in this time interval.
  • an estimated phase angle ⁇ ⁇ corresponding to an intermediate time between time t and time t + ⁇ T is calculated and used.
  • the rotational coordinate transformation does not have to consider the amount of change in the estimated phase angle ⁇ ⁇ . Good.
  • FIG. 2 shows a configuration in which the induced voltage is estimated using equation (9) and the rotational speed and the magnetic pole position are estimated in the same manner as in the first embodiment.
  • the fixed coordinate information was explained in the example of a discrete system when explaining the physical meaning of the differential component, but this clearly shows that there is no phase change in the discrete system than in the analog system.
  • the rotating coordinate conversion since the calculation corresponding to the current differentiation is performed in the rotating coordinate system, the rotating coordinate conversion may be performed only in the induced voltage calculating period. However, although it increases to double coordinate conversion of current detection and current differentiation, it can be simplified compared with sequential rotation coordinates.
  • the induced voltage is calculated using current differentiation in the zero voltage vector period. This is characterized in that since a change in current is detected during a period in which the switching element used in the inverter is not operating, it is not affected by dead time or delay in operation of the switching element.
  • the switching element is composed of semiconductor elements such as IGBTs and diodes, the inverter output voltage does not become zero even during the zero voltage vector period because there is a voltage drop of each element. .
  • FIG. 3 shows a case where it is applied to FIG.
  • reference numeral 12 denotes a voltage drop correction unit, which inputs three-phase current detections i u , i v , i w and detects the voltage correction amount v ceu (i u for the voltage drop of the semiconductor element from the detected value of the current component of each phase. ), vcev (iu), vcew (iu) are obtained from the table data. Further, it is converted into a rotational coordinate system of the estimated phase by the equation (10) and converted into voltage correction amounts vce ⁇ and vce ⁇ .
  • the induced voltage calculation unit 5 adds the corrections of the voltage correction amounts vce ⁇ and vce ⁇ to the equation (9), corrects them to the equation (11), and outputs the speed electromotive forces e ⁇ and e ⁇ .
  • the induced voltage calculation unit 5 executes the calculation of the equation (12) instead of the equation (3) to obtain the speed electromotive force e ⁇ , e Find ⁇ .
  • the electromotive force in consideration of the voltage drop component of the semiconductor switching element, it becomes possible to control normally to a lower speed. Further, the accuracy of the estimated speed and magnetic pole position is improved, and the estimated speed and 6f ripple component of the magnetic pole position caused by the voltage drop of the element can be reduced.
  • Example 3 the voltage drop component of the semiconductor switching element is corrected to enable stable operation up to a lower speed.
  • the current components for calculating the voltage drop components due to the voltage correction amounts vce ⁇ and vce ⁇ and the winding resistance R also cause a voltage error unless the temporal matching with the current differential information is accurate.
  • the voltage error factor is reduced in consideration of the error reduction method by the statistical processing and the current sampling time.
  • This embodiment is specified as a discrete system.
  • S0 to S7 are current sampling timings
  • TS is a current sampling period.
  • the differential calculation of the current is performed by calculating the difference after the rotation coordinate system is converted into the rotation coordinate system
  • the rotation coordinate conversion is performed after calculating the difference in the fixed coordinate system.
  • multipoint current sampling synchronized with the carrier is performed, and current difference calculation is performed using data between the two sampling points sandwiching the apex of the triangular wave carrier as necessary.
  • the current value for calculating the current control and the voltage drop component of R uses the PWM carrier synchronization current sampling value detected at the times of the vertices S0 and S4 of the triangular wave carrier. That is, a necessary portion is selected and utilized from multipoint current detection in which current detection and its differentiation are synchronized with the PWM carrier.
  • a normal three-phase inverter is configured with three arms by switching elements such as IGBTs having six reverse conducting diodes, and there are current paths that circulate through the switching elements inside the inverter during two types of zero voltage vector periods.
  • switching elements such as IGBTs having six reverse conducting diodes
  • a voltage error due to variations in switching elements occurs. Therefore, in order to average the two types, a moving average of two times of the current sampled at the apex portion and the bottom portion of the triangular wave carrier waveform is used in the calculation of the equations (11) and (12). Thereby, the voltage error component can be suppressed statistically.
  • the detected current value used for the voltage drop component and R voltage drop component obtained by the voltage drop correction unit 12 is the current sampling value at the time of the maximum vertices S0 and S3 when it is near the vertex of the triangular wave carrier waveform.
  • the current sampling value at the time of the minimum vertex S4 may be used.
  • the average value of the current sampling values of S3 and S5 can be used instead of the current sampling value of the time S4, and the average value of the current sampling values of S7 and S9 can be used instead of the current sampling value of S8. .
  • FIG. 5 shows a simulation result when the fourth embodiment is applied.
  • (A) is a speed
  • (b) is a biaxial current component
  • (c) is a torque
  • (d) is a difference between an actual magnetic pole position and an estimated magnetic pole position.
  • FIG. 5 shows the case where the load torque is changed at time t2 and the speed command is changed at time t6.
  • the axis error ⁇ e is very small as shown in (d), and when the speed is 5% or less or 0% speed passes from forward rotation to reverse rotation. It turns out that it can be estimated accurately even at times.
  • Temporal matching is performed in which differentiation is performed using the current sampled near the peak of the carrier and the peak as the detected current value.
  • the current value obtained by sampling the synchronous current at the average carrier vertex of the two current components before and after the vertex that approximates the current differentiation is substituted.
  • the current differential component and the voltage correction amount vce and the current used for the voltage drop of R can be temporally matched and the influence of noise can be reduced.
  • the voltage drop correction unit 12 in the third embodiment shown in FIG. 3 obtains the voltage drop corresponding to the detected value of the current component of each of the three phases using table data.
  • the voltage drop correction unit 12 also uses table data to calculate the voltage drop component.
  • the table data is individually obtained for each phase and for each positive and negative current polarity. It is configured so that it can be set. As a result, even if there is a variation in the characteristics of the semiconductor, it can be corrected more accurately.
  • the estimated phase ⁇ ⁇ is calculated using the current signals i ⁇ and i ⁇ in the ⁇ - ⁇ coordinates and the current differential information pi ⁇ and pi ⁇ during the zero voltage vector period. .

Landscapes

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

Abstract

【課題】ＰＭモータの位置センサレス制御装置における推定位相は、低速領域では正確に推定できていない。 【解決手段】零電圧ベクトル期間中のγーδ座標での前記電流検出値と電流微分情報を入力して速度起電力ｅγ，ｅδを演算する誘起電力演算部を設ける。算出された速度起電力ｅγ，ｅδを回転速度推定部に入力して推定速度ω＾を求め、この推定速度を時間積分することで推定位相θ＾を算出したものである。

Description

ＰＭモータの位置センサレス制御装置
　本発明は、ＰＭモータの位置センサレス制御装置に係わり、ＰＭモータの低速運転時に、電流に微分情報を利用して誘起電圧を推定し、この推定電圧を利用して回転速度・磁極位置推定を可能とするセンサレス制御装置に関するものである。
　磁石の性能向上に伴い、制動（ダンパ）巻線のない永久磁石を界磁源とする同期電動機（以下ＰＭモータという）が実用に供されている。サーボ用途などでは、このＰＭモータを可変速制御するために回転子磁極の位置情報が必要である。しかし、この位置情報を検出するセンサには半導体素子などの電子部品が内蔵されていることが多く、また小形なものは機械的強度が弱いなど、耐環境性や耐久性に課題を有している。
そのため、高い位置制御精度や応答性が要求されない用途では、位置センサを使用しないでインバータの電圧や電流情報から磁極位置を推定する位置センサレス制御法を適用することにより制御性能よりも信頼性を改善することが要望されており、現在まで多くの研究がなされてきた。
　この位置センサレス制御方式の原理を大別すると次の２種類がある。
一つは、界磁磁石の磁束によって発生する速度起電力成分を推定する方法である。この方法は、インバータ出力の電圧や電流の基本波成分からモータの界磁磁束による速度起電力を推定するものである。ＰＭモータは、界磁軸（ｄ軸）とその直交軸（ｑ軸）のインダクタンス成分が等しいという非突極機特性と、また、各軸のインダクタンスが等しくない突極機特性という観点で分類できるが、この方法はその両方に適用できるため、ＰＭモータ全般に利用できることが特長である。しかし、回転速度が低い領域では速度起電力が小さくなってしまう。通常のインバータは、ＰＷＭ変調方式を利用した電圧制御によりモータに電力を供給しているが、インバータ出力電圧を検出する場合には、高周波成分を含むＰＷＭ波形から正確に、かつ高速に基本波成分の電圧を検出することは難しいという問題を有している。
　ＰＷＭ変調を適用する前段階での電圧指令を、電圧検出情報の代用に適用する方法もある。この方法は、デットタイムやスイッチング素子の動作遅れ時間などによる電圧誤差の影響が存在する。そのためＰＭモータの速度が低下して速度起電力が小さくなってくるのに対し、この電圧誤差分は減少しないため相対的に無視できなくなり、正確な位置センサレス制御が困難になる。
　もう一つは、インバータの出力電圧や電流成分に高周波成分を重畳することにより、インダクタンスを計測しながら磁気的な突極軸を推定する方法である。
この方法は、ｄ軸とｑ軸のインダクタンスに差がある突極性を有するＰＭモータであれば、各軸成分のインダクタンスを逐次計測し続けることにより、界磁磁極軸の位相を推定することができる。しかし、この方法は非突極機特性のＰＭモータには適用できないほか，界磁磁極軸（ｄ軸）にはＮ極とＳ極の２種類の極性が存在するので、これらの判別を行うために磁気飽和などを利用した付加的な磁極判別法も制御に追加する必要がある。
　つまり、前者の起電力を推定する方法は低速域では正確に動作できなくなり、後者の高周波を重畳する方法では低速域での動作は可能であるが、適用できるモータに制限があり、非突極機特性を有する場合や磁気飽和が少ないＰＭモータには適用できない。
M.Janson,L.Harnefors,O.Wallmark,andM.Leksell：　　　　　　　　　 "Synchronization at Startup and Stable Rotation 　　　　　　　　　　 Reversal of Sensorless Nonsalient PMSM Drives"，　　　　　　　　　　IEEE Trans. IE, Vol.53, No.2,pp.379-387(2006) L.Harnefors, M.Janson,R.Ottersten, and K.Pietilainen,　　　　　　　　　 "Unified Sensorless Vector Control of Synchronous　　　　　　　　　　 and Induction Motors" IEEE Trans.Ind.Electron.,　　　　　　　　　　 Vol.50, No.1,pp.153-160,Feb.2003,
　低速の位置センサレス制御方式としては非特許文献１があり、図６で示すような制御ブロックを用いてセンサレス制御を行っている。非特許文献１で示す制御方式は基本的には起電力推定方式の一種で、推定方式の原理は非特許文献２に記載されているが、推定は図６で示す点線で囲んだ回転速度推定部で行われる。
　すなわち、電圧指令値Vsから電流指令ｉrefとモータ定数Rs＋jω1Lsの積である巻線の抵抗成分と電機子反作用成分の電圧降下成分を減算して誘起電圧ｅを求める。位置センサレス制御の内部では推定した界磁磁極位相をｄ軸とする制御基準軸に設定し、この基準軸を実軸とみなす。そして前述の誘起電圧ｅの虚部成分をｑ軸の誘起電圧ｅｑとして、また、実部をｄ軸誘起電圧ｅｄとして求める。この制御基準軸の各成分のうちeq電圧成分を、推定角速度ω1に対応したゲインλsとedを乗じたもので補正し、さらにこの誘起電圧を磁束Ψmで除算して速度の推定に利用する。実際には外乱や位置推定と回転座標変換を介してeq電圧成分にフィードバックする項が存在するので、このループにより発振しないように広帯域を制限するローパスフィルタをかけてから推定角速度ω1を推定している。
　図６で示す方式の特徴は、零速度付近で正転と逆転とを誤って逆回転方向に推定しても、自動的に推定速度が正常な回転方向に修正される点にある。しかし電圧情報には、例えば、従来のＰＷＭ変調に入力される電圧指令Vsを代用する方式では、デットタイム誤差やスイッチング素子のスイッチング遅れなどによる電圧誤差に関する対策が適格に適用されていないと、回転方向の推定誤りがあった場合でも直ぐには検出できず、或る程度の逆転速度まで達しないと正常回転方向に戻ることはできない。
　本発明が目的とするのは、このような電圧指令と実電圧に電圧誤差が存在していても、正常回転方向に復帰できる速度の下限を、より低い速度まで拡張することであり、ひいては低速でも正確な位置推定を可能としたＰＭモータの位置センサレス制御装置を提供することにある。
　本発明の１つの観点によれば、速度指令と推定速度から電流指令を生成し、この電流指令と回転座標変換部を通して検出されたγーδ座標の電流検出値から電圧指令を生成し、逆回転座標変換部およびＰＷＭ変調パターンにより制御されるインバータを介してＰＭモータを制御するものにおいて、
零電圧ベクトル期間中のγーδ座標での前記電流検出値ｉγ，ｉδと電流微分情報ｐｉγ，ｐｉδを入力して速度起電力ｅγ，ｅδを演算する誘起電圧演算部と、誘起電圧演算部により演算された速度起電力ｅγ，ｅδを用いて推定速度ω＾を求める回転速度推定部と、推定速度を時間積分して推定位相θ＾を算出し、この推定位相を前記回転座標変換部と逆回転座標変換部に出力して回転座標の基準位相に利用することを特徴としたものである。
　本発明の１つの観点による誘起電圧演算部は、次式により速度起電力ｅγ，ｅδを算出することを特徴としたものである。
Figure JPOXMLDOC01-appb-M000006
ただし、Ｒ：巻線抵抗、Ｌ：巻線インダクタンス、φｄ：鎖交磁束、ｐ：微分演算子（ｄ／ｄｔ）
　本発明の１つの観点による前記回転速度推定部は、次式により推定速度ω＾を算出することを特徴としたものである。
Figure JPOXMLDOC01-appb-M000007
ただし、α＝α0＋λ2｜ω＾｜、α0は固定、λ1，λ2は正のゲイン、ｚ-1は前回のサンプリング値
　本発明の別の観点によれば、前記電流微分検出部を回転座標変換部の検出電流入力側に設け、前記ＰＷＭ変調パターンの零電圧ベクトル期間中に発生する三相電流検出とその電流微分量を回転座標変換器に入力し、
前記誘起電圧演算部は、回転座標変換部により回転座標変換した前記電流検出値ｉγ，ｉδと電流微分情報ｐｉγ’，ｐｉδ’を入力して速度起電力ｅγ，ｅδを算出することを特徴としたものである。
　本発明の別の観点による前記誘起電圧演算部は、次式により速度起電力ｅγ，ｅδを算出することを特徴としたものである。
Figure JPOXMLDOC01-appb-M000008
　本発明の別の観点による前記電流微分情報ｐｉγ’，ｐｉδ’は、回転座標変換部によって固定座標の直交二軸α，βの電流成分ｉα，ｉβに変換後、微分演算を実行して得ることを特徴としたものである。
　本発明の更に別の観点によれば、前記回転座標変換部の入力側にスイッチング素子の電圧降下を補正する電圧降下補正部を設け、前記ＰＷＭ変調パターンの零電圧ベクトル期間中に発生する三相電流検出と電圧降下補正部からの電圧降下補正量を回転座標変換部に入力し、
前記誘起電圧演算部に、前記電流検出値と電流微分情報、及び電圧降下補正量を入力して速度起電力を演算することを特徴としたものである。
　本発明の別の観点による前記電流微分情報は、回転座標系又は固定座標系の何れか一方の座標系に基づいたものであることを特徴としたものである。
　本発明の別の観点による前記誘起電圧演算部は、次式により速度起電力ｅγ，ｅδを算出することを特徴としたものである。
Figure JPOXMLDOC01-appb-M000009
ただし、ｖceγ，ｖceδ：電圧補正量
　本発明の他の観点による前記誘起電圧演算部は、次式により速度起電力ｅγ，ｅδを算出することを特徴としたものである。
Figure JPOXMLDOC01-appb-M000010
　本発明の他の観点による前記電流微分情報は、零電圧ベクトル期間におけるキャリア頂点前後２点の電流サンプリング値の差分演算値として前記回転座標変換部に入力することを特徴としたものである。
　本発明の更に他の観点による前記電流微分情報は、零電圧ベクトル期間におけるキャリア頂点と底辺間の電流サンプリング値の移動平均値を微分した値を前記回転座標変換部に入力することを特徴としたものである。
　本発明の他の観点による前記電圧降下補正部における電圧降下補正量は、三相検出電流値に対応したテーブルデータを用いて求め、テーブルデータは三相の各相毎個別に、且つ正負の電流極性毎に個別に設定されたことを特徴としたものである。
　以上のとおり、本発明によれば、零電圧ベクトル期間中のγーδ座標の電流信号ｉγ，ｉδと電流微分情報ｐｉγ，ｐｉδを用いて推定位相θ＾を算出したものである。これにより、デットタイムの影響を排除した誘起電圧の推定が可能になり、より低速まで正常な位置センサレス制御が実現できるなどの効果を奏するものである。
本発明の実施形態を示す制御装置の構成図。 本発明の他の実施形態を示す制御装置の構成図。 本発明の他の実施形態を示す制御装置の構成図。 本発明の電流サンプリングの説明図。 本発明のシミュレーション結果の波形図。 従来のＰＭモータのセンサレス制御装置における構成図。
　本発明は、誘起電圧演算部に検出された電流信号と電流微分情報を入力して速度起電力を演算し、この速度起電力に基づいて推定磁極位相を求めるようにしたものである。以下、各実施例に基づいて詳述する。
　図１は、本発明の第１の実施例を示す位置センサレス制御装置のブロック図を示したものである。１はＰＷＭ制御されるインバータ、２はＰＭモータで、ここでは、実機のＰＭモータのＮ極の界磁軸をｄ軸と定義し、それに対して正転方向に９０゜進んだ位相をｑ軸とする。しかし、位置センサが無いことでこのｄ軸とｑ軸は直接検出することはできない。そこで、磁極推定により仮想したＮ極軸をγ軸，回転方向に電気角で９０゜進んだ位相をδと定義する。
　３は回転座標変換部で、電流センサにより検出された３相電流ｉu,ｉv,ｉwを入力して推定軸であるｉγ，ｉδに座標変換する。座標変換された電流信号ｉγ，ｉδは電流微分検出部４と誘起電圧演算部５、および電流制御部９にそれぞれ入力される。電流微分検出部４では、零電圧ベクトル期間中のγーδ軸電流の変化量ｐｉγ，ｐｉδを検出し、誘起電圧演算部５に入力する。誘起電圧演算部５では、入力された電流信号ｉγ，ｉδと電流微分情報ｐｉγ，ｐｉδを用いて後述のように速度起電力ｅγ．ｅδを演算する。６は回転速度推定部、７は積分器で、回転速度推定部６による推定速度を積分して推定磁極位置θ＾を演算し、回転座標変換部３と逆回転座標変換部１０に出力する。８は速度制御部、９は電流制御部である。
　次に基本的な動作について説明する。
速度制御部８では、速度指令ω*と速度推定ω＾の入力情報から、トルク指令に相当するδ軸成分の電流指令ｉδ *を出力する。この電流指令ｉδ *とγ軸の任意の電流指令ｉγ *が電流制御部９に入力される。電流制御部９には、回転座標変換部３を介して回転座標に変換された電流信号ｉγ，ｉδが入力されており、電流制御部９ではこれらの電流指令ｉδ *，ｉγ *および電流信号ｉγ，ｉδを比較してフィードバック演算を行い、推定磁極軸を基準とする回転座標（γ‐δ座標）上の電圧指令ｖγ *，ｖδ *を出力する。この電圧指令ｖγ *，ｖδ *は、回転座標変換部３とは逆の動作を行う逆回転座標変換部１０により、固定座標への逆回転変換や２相３相変換などを行って三相交流の電圧指令ｖu，ｖv，ｖwとしてインバータ１に入力される。インバータ１では、三相電圧指令ｖu，ｖv，ｖwをＰＷＭ変調により電力増幅して略等価な電圧を出力する。
以上は位置センサがある場合もない場合も共通の制御ブロックである。
　次に位置センサレス制御の部分について説明する。
電流微分検出部４は、零電圧ベクトル期間中のγーδ軸電流の変化量ｐｉγ，
ｐｉδを検出する。γーδ軸電流は三相電流検出ｉu,ｉv,ｉwから回転座標変換を行うことで得られるが、回転座標変換方法としては、アナログ乗算器やアナログの加減算器などにより信号変換を実現する連続信号による方法と、三相電流検出ｉu,ｉv,ｉwをアナログ／ディジタル変換器（Ａ／Ｄ変換器）によってディジタル信号に変換した後に、ＣＰＵなどのディジタル変換器により信号変換部の演算を実行する離散信号による方法がある。
　本実施例では変化量である電流微分情報ｐｉγ，ｐｉδが検出されればよく、回転座標変換方法は何れの方法でもよいが、例えば、アナログ変換の場合には、回転座標変換後の信号にアナログ微分器を適用しておき、零電圧ベクトル期間における電流微分成分をサンプリングホールドすることで変化量ｐｉγ，ｐｉδが得られる。また、Ａ／Ｄ変換器を利用した場合には、零電圧ベクトル期間中において多数の時刻で電流をサンプリングして信号変換部にて変換し、複数時刻の離散電流値から差分近似により微分成分を得ることができる。何れかの方法によって得られた電流信号ｉγ，ｉδと電流微分情報（変化量）ｐｉγ，ｐｉδを利用して位置センサレス制御が構成される。
　ここで、基本的なＰＭモータの特性を表現する方程式を定義すると、実機のＮ極軸をｄ軸とする電気角で表した直交座標系（ｄ－ｑ座標）において、ＰＭモータの電圧方程式は（１）式となる。
Figure JPOXMLDOC01-appb-M000011
ただし、ｖｄ，ｖｑ：ｄ，ｑ軸電圧、ｉｄ，ｉｑ：ｄ，ｑ軸電流、Ｒ：巻線抵抗、Ｌ：巻線インダクタンス、φｄ：磁石の発生磁束により固定子巻線に鎖交する磁束成分、ω：回転子角速度（電気角，ω＝ｄθ／ｄｔ）、θ：回転子磁極位置（電気角）、ｐ：微分演算子（ｄ／ｄｔ）
本発明では非突極特性を有するＰＭモータを制御対象としているため、ｄ軸とｑ軸のインダクタンス成分は等しいものとみなして共通の係数で表現している。
　次に、推定位相θ＾を基準とするγーδ座標では、実位相θとの誤差位相（軸誤差）θe＝θーθ＾を考慮して（２）式で取り扱う。ここで、モータのインピーダンスや永久磁石による鎖交磁束などは真値が得られていると仮定し、（１）式と同じ変数を使用している。また、速度推定ω＾と推定位相θ＾には、ω＾＝θ＾／ｄｔの関係がある。
Figure JPOXMLDOC01-appb-M000012
　（２）式では軸誤差θe自体の時間変化は少ないものと仮定して、この軸誤差の微分項を無視している。また、（２）式では零電圧ベクトル期間におけるγーδ座標の電流微分情報ｐｉγ，ｐｉδは検出電流ｉγ，ｉδを推定位相θ＾により逐次回転座標変換したものを微分したものである。
　零電圧ベクトル期間中におけるインバータの出力電圧はｖγ＝０，ｖδ＝０であるので、（２）式の左辺を零とおき、さらに永久磁石による鎖交磁束φｄによる速度起電力ｅγ，ｅδを求める式に変形すると（３）式となる。
Figure JPOXMLDOC01-appb-M000013
　この零電圧ベクトル期間中の速度起電力ｅγ，ｅδを示す（３）式から磁極位置を推定している。（４）式は速度推定式に離散系の一次遅れを適用したものである。また、この一次遅れの帯域は（５）式のように、速度推定ω＾＝θ＾／ｄｔにより可変している。
Figure JPOXMLDOC01-appb-M000014
Figure JPOXMLDOC01-appb-M000015
ここで、α0は固定値、λ1，λ2は正のゲイン、ｚ-1は前回のサンプリング値図１で示す回転速度推定部６は、（４），（５）式に基づいて推定速度ω＾を算出し、積分器７を介して推定位相θ＾を出力する。すなわち、誘起電圧演算部５では、（３）式により零電圧ベクトル期間におけるγーδ座標の電流信号ｉγ，ｉδと電流微分情報ｐｉγ，ｐｉδから零電圧ベクトル期間中の速度起電力ｅγ，ｅδを出力する。回転速度推定部６は、入力された速度起電力ｅγ，ｅδを用いて（４），（５）式の演算を実行し、推定速度ω＾を出力する。積分器７では、推定速度ω＾を時間積分して推定位相θ＾を算出し、回転座標変換部３と逆回転座標変換部１０に出力して回転座標の基準位相に利用する。
　この実施例によれば、ＰＷＭ変調を行うインバータで非突極特性を有するＰＭモータを駆動する場合に、零電圧ベクトル期間中のγーδ座標の電流信号ｉγ，ｉδと電流微分情報ｐｉγ，ｐｉδを用いて推定位相θ＾を算出したものである。これにより、デットタイムの影響を排除した誘起電圧の推定が可能になり、より低速まで正常な位置センサレス制御が実現できるものである。
　図２は第２の実施例を示したもので、図１で示す第１の実施例との相違点は、電流微分検出部１１を固定座標系に設けて電流微分情報ｐｉ'γ，ｐｉ'δを得ていることである。
実施例１では、電流微分情報を正確に得るためには回転座標変換に使用する推定位相θ＾も連続系の値でなくてはならない。そのため，三相電流信号をＡ／Ｄ変換してディジタル値に変換する場合，回転座標変換に使用する推定位相θ＾も逐次更新する必要があり演算量が多くなる。この実施例２はその点を考慮したもので、固定座標系で電流微分情報ｐｉ'γ，ｐｉ'δを得るよう構成することによって推定位相θ＾を逐次更新する必要がなくなり簡単な演算実行を可能としたものである。
　具体的には、零電圧ベクトル期間中の速度起電力ｅγ，ｅδを演算する時刻だけ、零電圧ベクトル期間における電流信号ｉγ，ｉδと電流微分情報ｐｉ'γ，ｐｉ'δをサンプリングすればよく、座標変換もこの零電圧ベクトル期間にだけ実行すればよい。これにより、演算を実行するＣＰＵやディジタル回路での演算時間が短くなり、回路動作が遅くても適用可能であるという利点が得られる。
　なお、図２では三相電流検出ｉu,ｉv,ｉwを直接微分するように表示しているが、この微分演算の構成方法としてはアナログ信号でもディジタル信号でも実現でき、要は最終的に零電圧ベクトル期間中の速度起電力ｅγ，ｅδを求めるために必要なγーδ座標の電流微分情報ｐｉ'γ，ｐｉ'δに相当する情報であればよい。例えば、先に３相２相変換して固定座標系の直交二軸（αβ座標）の電流信号ｉα，ｉβに変換してから微分演算を実行し、その後に回転座標変換を適用するように演算順序を入れ替えてもよい。
　次に図２における電流微分情報ｐｉ'γ，ｐｉ'δを適用した場合における速度推定について説明する。
固定子巻線を基準とする固定座標系におけるＰＭモータの電圧方程式は（６）式となる。
Figure JPOXMLDOC01-appb-M000016
ここで、ｖα，ｖβ：α，β軸電圧、ｉα，ｉβ：α，β軸電流、R：巻線抵抗、
L：巻線インダクタンス、φd：磁石の鎖交磁束、ω：回転子角速度（電気角）、θ：回転子磁極位置、p：微分演算子
　(１)式では回転座標系で取り扱っていたため，右辺第１項のインピーダンス行列内に電機子電流により発生した磁束によって、速度起電力に相当する項ω＾Lが存在していたが、(６)式では固定座標系で取り扱っているためこの成分が存在してなく、これにより、演算式がより簡単となる。
　厳密に固定座標系における微分成分を角速度ωで回転する回転座標に変換する場合、ベクトル解析上ではαβ座標系に直交な速度ベクトルを定義して下記のような変換式を適用しなければならない。つまり回転座標変換に利用する位相変化も考慮した厳密な微分系では電機子反作用磁束の速度起電力項に注意が必要である。
ｐi(d,q)＋ω×i(α,β)＝ｐi(α,β)
しかし、図２の回転座標上の電流微分情報ｐｉ'γ，ｐｉ'δは、固定座標上において微分したものをある時刻の推定位相角θ＾で回転座標変換したものであり、回転座標変換する時刻における位相の変化項までは含まれていない。
　上述の内容を離散系で計算する例を説明すると次のようになる。
離散系において、αβ軸電流の微分演算を差分で近似する場合には次式で計算できる。
Figure JPOXMLDOC01-appb-M000017
ここで、ｔ，ｔ＋△Tは同一の零電圧ベクトル期間において時間差△Tを有する２回の電流サンプリング時刻を示している。また、△ｉα，△ｉβは、この時間間隔の電流差分を示している。
　電流差分△ｉα，△ｉβを利用して近似した微分成分を回転座標変換するには、時刻ｔと時刻ｔ＋△Tの中間時刻に相当する推定位相角θ＾を計算し、これを用いて回転座標変換する。つまり、電流差分△ｉα，△ｉβにはω×i(α,β)の項が既に含まれているので，回転座標変換においては推定位相角θ＾の変化量まで考慮しなくてもよい。
　このように電流微分成分のｐｉ'γ，ｐｉ'δの物理的な意味は (２)式とは異なっているので、代わりに（８）式で取り扱う必要がある。推定位相角θ＾を基準とするγーδ座標では、実位相θとの誤差位相θe＝θ－θ＾を考慮して（８）式で取り扱う。
Figure JPOXMLDOC01-appb-M000018
　後は実施例１と同じであり、誤差位相θeの変化は小さいとしてこの微分項は無視する。
零電圧ベクトル期間中に（８）式の左辺を零とおき、誘起電圧式の計算式に変形すると（９）式となる。
Figure JPOXMLDOC01-appb-M000019
　(９)式を用いて誘起電圧を推定し、実施例１と同様に回転速度及び磁極位置を推定するように構成したものが図２である。
固定座標情報を微分成分の物理的な意味を説明するときに離散系の例で説明したが、これはアナログ系で説明するよりも離散系の方が、位相変化が無いことを明確に示すことができるために利用したものであり、実施例１と同じように電流微分や回転座標変換を実現するには、アナログ演算回路でもディジタル回路でもどちらの利用でもよいことは勿論である。したがって，実施例２では電流微分演算を離散系に限定するものではない。
　この実施例２によれば、回転座標系で電流微分に相当する演算を行なったことで、回転座標変換は誘起電圧演算周期のみ実施すればよい。ただし、電流検出と電流微分という２倍の座標変換に増えるが，逐次回転座標を行なうよりも簡素化できるものである。
　実施例１，２では、零電圧ベクトル期間における電流微分を利用して誘起電圧を算出している。これは、インバータに使用されるスイッチング素子が動作していない期間での電流変化を検出しているので、デッドタイムやスイッチング素子の動作遅れなどの影響を受けないという特徴がある。しかし、厳密にはスイッチング素子をＩＧＢＴやダイオードなどの半導体素子によりインバータを構成した場合には、各素子の電圧降下分が存在するため、零電圧ベクトル期間中においてもインバータ出力電圧は零にはならない。
　この実施例はこの点を考慮して図３のように構成したものである。
通常，スイッチング素子の電圧降下成分は電流と関係があるので、三相各相の電流成分の検出値から半導体素子の電圧降下分を推定して補正するものである。
この実施例は、図１，図２で示す何れの実施例にも適用できるが、図３では図２に適用した場合を示している。
　図３において、１２は電圧降下補正部で、三相電流検出ｉu,ｉv,ｉwを入力して各相の電流成分の検出値から半導体素子の電圧降下分の電圧補正量ｖceu(ｉu), ｖcev(ｉu), ｖcew(ｉu)をテーブルデータから求める。さらに、（１０）式にて推定位相の回転座標系に変換して電圧補正量ｖceγ,ｖceδに変換する。
Figure JPOXMLDOC01-appb-M000020
そして、誘起電圧演算部５では、（９）式にこの電圧補正量ｖceγ，ｖceδの補正を加えて（１１）式に修正し、速度起電力ｅγ，ｅδを出力する。
Figure JPOXMLDOC01-appb-M000021
　なお、図１で示す実施例１に電圧成分の補正を行う場合には、誘起電圧演算部５は（３）式の代わりに（１２）式の演算を実行して速度起電力ｅγ，ｅδを求める。
Figure JPOXMLDOC01-appb-M000022
　この実施例によれば、半導体スイッチング素子の電圧降下成分まで考慮して起電力推定を行うことにより、より低速まで正常に制御できるようになる。また、推定した速度・磁極位置の精度が向上し，素子の電圧降下に起因する推定速度、磁極位置の６ｆリプル成分も小さくすることができるものである。
　実施例３では、半導体スイッチング素子の電圧降下成分を補正することにより、より低速までの安定した動作を可能にした。しかし、インバータの主回路素子の特性にバラツキが存在すると、電圧降下補正部１２による正確なる電圧補正ができなくなる可能性がある。また、電圧補正量ｖceγ，ｖceδや巻線抵抗Ｒによる電圧降下成分を計算する電流成分についても、電流微分情報との時間的な整合が正確でないと電圧誤差の要因となる。
　そこで、実施例４では統計処理による誤差の低減の仕方と、電流サンプリング時刻などを考慮して電圧誤差要因を減少するようにしたものである。
なお、この実施例は離散系に特定されるものである。
図４において、Ｓ０～Ｓ７は電流のサンプリングタイミング、ＴSは電流のサンプリング周期である。
　キャリア１周期にＳ０～Ｓ７の８回の電流サンプリングが実行される場合を考える。低速である場合には電圧指令ｖu,ｖv,ｖwの振幅は小さく，三相インバータのスイッチング素子の動作は三角波キャリアの中間付近のＳ１からＳ３の期間、およびＳ５からＳ７の期間に集中している。そのため、Ｓ３とＳ５の間は三相ともＯＮ状態であり、また、Ｓ７からＳ９の間は三相ともＯＦＦ状態という２種類の零電圧ベクトル期間になっている。そこで、この２種類の零電圧ベクトル期間において、Ｓ３とＳ５，またはＳ７とＳ９のように２回の電流サンプリングおよびＡ／Ｄ変換を行う。
　例えば、Ｓ３とＳ５のタイミングで検出した電流をＩ３，Ｉ５とすると，この時の電流の微分に相当する差分演算値△Ｉは（１３）式となる。
Figure JPOXMLDOC01-appb-M000023
　この電流の微分演算は、実施例１では回転座標系に変換した後に差分演算し、実施例２では固定座標系において差分演算してから回転座標変換する。このように実施例４では，キャリアに同期した多点電流サンプリングを行っておき、必要に応じて三角波キャリアの頂点をはさむ２点のサンプリング間のデータを利用して電流の差分演算を行なう。また，電流制御やＲの電圧降下成分を演算する電流値は、三角波キャリアの頂点Ｓ０とＳ４の時刻で検出したＰＷＭキャリア同期電流サンプリング値を利用する。すなわち、電流検出およびその微分をＰＷＭキャリアに同期した多点電流検出から必要な部分を選択して利用する。
　さらに，通常の３相インバータは６個の逆導通ダイオードを内蔵したＩＧＢＴなどのスイッチング素子により３アームを構成しており、２種類の零電圧ベクトル期間ではインバータ内部のスイッチング素子を循環する電流路が異なり、スイッチング素子のバラツキによる電圧誤差が発生する。そこで、２種類を平均するために、(１１)式や(１２)式の計算に三角波キャリア波形の頂点部分と底部部分でサンプリングした電流の２回分の移動平均を用いる。これにより統計的に電圧誤差成分が抑制できる。
　また、電圧降下補正部１２で求まった電圧降下成分やＲの電圧降下成分に使用する電流の検出値を、三角波キャリア波形の頂点付近の場合には最大頂点Ｓ０，Ｓ３の時刻の電流サンプリング値を、底辺部分であれば最小頂点Ｓ４時刻の電流サンプリング値を利用すればよい。しかし，Ｓ４時刻の電流サンプリング値の代わりに、Ｓ３とＳ５の電流サンプリング値の平均値で、また、Ｓ８の電流サンプリング値の代わりにＳ７とＳ９の電流サンプリング値の平均値で代用することができる。
　図５は、この実施例４を適用した場合のシミュレーション結果を示したものである。同図（ａ）は速度、（ｂ）は二軸電流成分、（ｃ）はトルク、（ｄ）は実際の磁極位置と推定した磁極位置の差である。
図５は時刻ｔ２で負荷トルクを変動させ、時刻ｔ６で速度指令を変動させた場合を示したもので、通常は、起電力を用いて位置・速度推定を行う場合、定格の５～１０％程度の速度までしか精度良く推定することができないが、（ｄ）図で示すように軸誤差θeは非常に少なくなっており、５％以下の速度時や正転から逆転への０％速度通過時でも精度よく推定できることが分る。
　したがって、この実施例によれば、
（１）電流の検出値としてキャリアの頂点とその頂点付近でサンプリングした電流を用いて微分するという時間的な整合をとる。
（２）または、電流微分を近似差分する頂点前後の２点の電流成分の平均キャリア頂点で同期電流サンプリングした電流値を代用にする。
これにより、電流微分成分と電圧補正量ｖceやＲの電圧降下に使用する電流との時間的な整合がとれてノイズの影響を低減することが可能となるものである。
また、頂点部と底辺部の２種類の零電圧ベクトル期間の電流検出値の移動平均をとることにより半導体素子のバラツキを統計的に抑制することが可能となるものである。
　図３で示す第３実施例での電圧降下補正部１２は、三相各相の電流成分の検出値に対応する電圧降下分はテーブルデータを利用して得ている。
実施例５でも、電圧降下補正部１２においても電圧降下成分を計算するためにテーブルデータを使用しているが、このテーブルデータを各相ごとに個別に、且つ正負の電流極性毎にも個別に設定できるように構成したものである。これにより、半導体の特性バラツキがあっても、より正確に補正することができるようにしたものである。
　したがって、この実施例によれば、スイッチング素子の特性バラツキを個別に補正したことにより、より正確な電圧降下成分の補正が可能になるものである。
　以上のとおり、本発明によれば、零電圧ベクトル期間中のγーδ座標の電流信号ｉγ，ｉδと電流微分情報ｐｉγ，ｐｉδを用いて推定位相θ＾を算出したものである。これにより、デットタイムの影響を排除した誘起電圧の推定が可能になり、より低速まで正常な位置センサレス制御が実現できるなどの効果を奏するものである。

Claims (13)

  1. 速度指令と推定速度から電流指令を生成し、この電流指令と回転座標変換部を通して検出されたγーδ座標の電流検出値から電圧指令を生成し、逆回転座標変換部およびＰＷＭ変調パターンにより制御されるインバータを介してＰＭモータを制御するものにおいて、
    零電圧ベクトル期間中のγーδ座標での前記電流検出値ｉγ，ｉδと電流微分情報ｐｉγ，ｐｉδを入力して速度起電力ｅγ，ｅδを演算する誘起電圧演算部と、誘起電圧演算部により演算された速度起電力ｅγ，ｅδを用いて推定速度ω＾を求める回転速度推定部と、推定速度を時間積分して推定位相θ＾を算出し、この推定位相を前記回転座標変換部と逆回転座標変換部に出力して回転座標の基準位相に利用することを特徴としたＰＭモータの位置センサレス制御装置。
  2. 前記誘起電圧演算部は、次式により速度起電力ｅγ，ｅδを算出することを特徴とした請求項１記載のＰＭモータの位置センサレス制御装置。
    Figure JPOXMLDOC01-appb-M000001
     ただし、Ｒ：巻線抵抗、Ｌ：巻線インダクタンス、φｄ：鎖交磁束、ｐ：微分演算子（ｄ／ｄｔ）
  3. 前記回転速度推定部は、次式により推定速度ω＾を算出することを特徴とした請求項１又は２記載のＰＭモータの位置センサレス制御装置。
    Figure JPOXMLDOC01-appb-M000002
     ただし、α＝α0＋λ2｜ω＾｜、α0は固定、λ1，λ2は正のゲイン、ｚ-1は前回のサンプリング値
  4. 前記電流微分検出部を回転座標変換部の検出電流入力側に設け、前記ＰＷＭ変調パターンの零電圧ベクトル期間中に発生する三相電流検出とその電流微分量を回転座標変換器に入力し、
    前記誘起電圧演算部は、回転座標変換部により回転座標変換した前記電流検出値ｉγ，ｉδと電流微分情報ｐｉγ’，ｐｉδ’を入力して速度起電力ｅγ，ｅδを算出することを特徴とした請求項１記載のＰＭモータの位置センサレス制御装置。
  5. 前記誘起電圧演算部は、次式により速度起電力ｅγ，ｅδを算出することを特徴とした請求項４記載のＰＭモータの位置センサレス制御装置。
    Figure JPOXMLDOC01-appb-M000003
  6. 前記電流微分情報ｐｉγ’，ｐｉδ’は、回転座標変換部によって固定座標の直交二軸α，βの電流成分ｉα，ｉβに変換後、微分演算を実行して得ることを特徴とする請求項４又は５記載のＰＭモータの位置センサレス制御装置。
  7. 前記回転座標変換部の入力側にスイッチング素子の電圧降下を補正する電圧降下補正部を設け、前記ＰＷＭ変調パターンの零電圧ベクトル期間中に発生する三相電流検出と電圧降下補正部からの電圧降下補正量を回転座標変換部に入力し、
    前記誘起電圧演算部に、前記電流検出値と電流微分情報、及び電圧降下補正量を入力して速度起電力を演算することを特徴とした請求項１又は４記載のＰＭモータの位置センサレス制御装置。
  8. 前記電流微分情報は、回転座標系又は固定座標系の何れか一方の座標系に基づいたものであることを特徴とした請求項１又は４又は７記載の何れかであるＰＭモータの位置センサレス制御装置。
  9. 前記誘起電圧演算部は、次式により速度起電力ｅγ，ｅδを算出することを特徴とした請求項４又は７記載のＰＭモータの位置センサレス制御装置。
    Figure JPOXMLDOC01-appb-M000004
     ただし、ｖceγ，ｖceδ：電圧補正量
  10. 前記誘起電圧演算部は、次式により速度起電力ｅγ，ｅδを算出することを特徴とした請求項１又は７記載のＰＭモータの位置センサレス制御装置。
    Figure JPOXMLDOC01-appb-M000005
  11. 前記電流微分情報は、零電圧ベクトル期間におけるキャリア頂点前後２点の電流サンプリング値の差分演算値として前記回転座標変換部に入力することを特徴とした請求項７乃至１０記載の何れかであるＰＭモータの位置センサレス制御装置。
  12. 前記電流微分情報は、零電圧ベクトル期間におけるキャリア頂点と底辺間の電流サンプリング値の移動平均値を微分した値を前記回転座標変換部に入力することを特徴とした請求項７乃至１０記載の何れかであるＰＭモータの位置センサレス制御装置。
  13. 前記電圧降下補正部における電圧降下補正量は、三相検出電流値に対応したテーブルデータを用いて求め、テーブルデータは三相の各相毎個別に、且つ正負の電流極性毎に個別に設定されたことを特徴とした請求項７乃至１２記載の何れかであるＰＭモータの位置センサレス制御装置。
PCT/JP2012/070312 2011-08-15 2012-08-09 Ｐｍモータの位置センサレス制御装置 WO2013024780A1 (ja)

Priority Applications (3)

Application Number Priority Date Filing Date Title
RU2014109115/07A RU2560090C1 (ru) 2011-08-15 2012-08-09 Устройство управления электродвигателем с постоянными магнитами без использования датчиков положения
CN201280039669.0A CN103733504B (zh) 2011-08-15 2012-08-09 Pm马达的无位置传感器控制装置
JP2013528991A JP5733404B2 (ja) 2011-08-15 2012-08-09 Ｐｍモータの位置センサレス制御装置

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
JP2011177371 2011-08-15
JP2011-177371 2011-08-15

Publications (1)

Publication Number Publication Date
WO2013024780A1 true WO2013024780A1 (ja) 2013-02-21

Family

ID=47715099

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2012/070312 WO2013024780A1 (ja) 2011-08-15 2012-08-09 Ｐｍモータの位置センサレス制御装置

Country Status (4)

Country Link
JP (1) JP5733404B2 (ja)
CN (1) CN103733504B (ja)
RU (1) RU2560090C1 (ja)
WO (1) WO2013024780A1 (ja)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2015037376A (ja) * 2013-08-09 2015-02-23 有限会社シー・アンド・エス国際研究所 交流電動機のデジタル式回転子位相速度推定装置
JP2015039278A (ja) * 2013-08-18 2015-02-26 有限会社シー・アンド・エス国際研究所 交流電動機のデジタル式回転子位相速度推定装置

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE102015214839A1 (de) * 2015-08-04 2017-02-09 Schaeffler Technologies AG & Co. KG Verfahren zum zeitdiskreten Regeln eines elektronisch kommutierten Elektromotors
JP6939800B2 (ja) * 2016-09-07 2021-09-22 日本電産株式会社 モータ制御方法、モータ制御システムおよび電動パワーステアリングシステム
CN107547027B (zh) * 2017-10-25 2020-07-21 株洲易力达机电有限公司 一种单电阻电机电流采样方法
RU2683586C1 (ru) * 2018-03-20 2019-03-29 Дмитрий Валерьевич Хачатуров Способ управления синхронным электродвигателем на постоянных магнитах

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2001251889A (ja) * 2000-03-06 2001-09-14 Hitachi Ltd 同期モータの回転子位置推定方法、位置センサレス制御方法及び制御装置
JP2004343963A (ja) * 2003-05-19 2004-12-02 Honda Motor Co Ltd ブラシレスｄｃモータの制御装置
JP2006262590A (ja) * 2005-03-16 2006-09-28 Honda Motor Co Ltd ブラシレスｄｃモータの制御装置
JP2007097263A (ja) * 2005-09-27 2007-04-12 Denso Corp 同期モータの磁極位置推定方法
JP2008092781A (ja) * 2006-09-30 2008-04-17 C & S Kokusai Kenkyusho:Kk 永久磁石同期電動機の駆動制御方法

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE3026202A1 (de) * 1980-07-10 1982-02-04 Siemens AG, 1000 Berlin und 8000 München Drehfeldmaschinenantrieb mit einer umrichtergespeisten drehfeldmaschine und einer mit zwei wechselspannungsintegratoren und einer rechenmodellschaltung verbundenen umrichtersteuerung
RU2092967C1 (ru) * 1995-07-27 1997-10-10 Новосибирский государственный технический университет Электропривод с синхронным двигателем
RU2141719C1 (ru) * 1998-03-25 1999-11-20 Мищенко Владислав Алексеевич Способ векторного управления синхронным электродвигателем с постоянными магнитами на роторе и электропривод для осуществления этого способа
JP4988374B2 (ja) * 2007-02-15 2012-08-01 三洋電機株式会社 モータ制御装置
CN101630938A (zh) * 2009-07-28 2010-01-20 哈尔滨工业大学 无位置传感器永磁同步电机转子初始位置辨识方法
RU99911U1 (ru) * 2010-08-05 2010-11-27 Кочанов Юрий Иванович Система векторного бездатчикового управления многофазным электродвигателем

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2001251889A (ja) * 2000-03-06 2001-09-14 Hitachi Ltd 同期モータの回転子位置推定方法、位置センサレス制御方法及び制御装置
JP2004343963A (ja) * 2003-05-19 2004-12-02 Honda Motor Co Ltd ブラシレスｄｃモータの制御装置
JP2006262590A (ja) * 2005-03-16 2006-09-28 Honda Motor Co Ltd ブラシレスｄｃモータの制御装置
JP2007097263A (ja) * 2005-09-27 2007-04-12 Denso Corp 同期モータの磁極位置推定方法
JP2008092781A (ja) * 2006-09-30 2008-04-17 C & S Kokusai Kenkyusho:Kk 永久磁石同期電動機の駆動制御方法

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2015037376A (ja) * 2013-08-09 2015-02-23 有限会社シー・アンド・エス国際研究所 交流電動機のデジタル式回転子位相速度推定装置
JP2015039278A (ja) * 2013-08-18 2015-02-26 有限会社シー・アンド・エス国際研究所 交流電動機のデジタル式回転子位相速度推定装置

Also Published As

Publication number Publication date
RU2560090C1 (ru) 2015-08-20
CN103733504B (zh) 2016-08-17
CN103733504A (zh) 2014-04-16
JPWO2013024780A1 (ja) 2015-03-05
JP5733404B2 (ja) 2015-06-10

Similar Documents

Publication Publication Date Title
JP4674525B2 (ja) 磁極位置推定方法及びモータ制御装置
JP5733404B2 (ja) Ｐｍモータの位置センサレス制御装置
Ghaderi et al. Wide-speed-range sensorless vector control of synchronous reluctance motors based on extended programmable cascaded low-pass filters
JP2009136085A (ja) 交流モータの制御装置
KR101840509B1 (ko) 동기전동기 센서리스 벡터제어를 위한 회전각 추정장치
JP6015486B2 (ja) 同期電動機の可変速制御装置
JP5428202B2 (ja) 永久磁石形同期電動機の制御装置
JP2004048958A (ja) Ｄｃブラシレスモータの制御装置
JP4632157B2 (ja) 永久磁石モータの駆動システム
JP4670405B2 (ja) 同期電動機のベクトル制御方法
JP6687228B1 (ja) 交流回転電機の制御装置
JP4781933B2 (ja) 電動機の制御装置
JP2010035352A (ja) 同期電動機のロータ位置推定装置
JP7304891B2 (ja) 回転機の制御装置および電動車両の制御装置
JP2012186911A (ja) モータ制御装置
JP5106295B2 (ja) 同期電動機のロータ位置推定装置
JP6384199B2 (ja) 位置推定装置、モータ駆動制御装置、位置推定方法及びプログラム
CN111543003B (zh) 旋转机的控制装置
Zhou et al. High-order terminal sliding-mode observer for speed estimation of induction motors
JP5186352B2 (ja) 電動機の磁極位置推定装置
JP5228435B2 (ja) インバータ制御装置とその制御方法
JP2007082380A (ja) 同期モータ制御装置
Hanamoto et al. Sensorless control of synchronous reluctance motor using modified flux linkage observer with an estimation error correct function
Noguchi et al. Study of parameter variations compensation in sensorless control of PMSM
CN110061669A (zh) 永磁同步电机直接转矩控制方法

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: 12823427

Country of ref document: EP

Kind code of ref document: A1

ENP Entry into the national phase

Ref document number: 2013528991

Country of ref document: JP

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: 2014109115

Country of ref document: RU

Kind code of ref document: A

122 Ep: pct application non-entry in european phase

Ref document number: 12823427

Country of ref document: EP

Kind code of ref document: A1