WO2020152785A1 - モータ制御方法、及び、モータ制御装置 - Google Patents

モータ制御方法、及び、モータ制御装置 Download PDF

Info

Publication number
WO2020152785A1
WO2020152785A1 PCT/JP2019/001902 JP2019001902W WO2020152785A1 WO 2020152785 A1 WO2020152785 A1 WO 2020152785A1 JP 2019001902 W JP2019001902 W JP 2019001902W WO 2020152785 A1 WO2020152785 A1 WO 2020152785A1
Authority
WO
WIPO (PCT)
Prior art keywords
axis
command value
axis current
voltage
current
Prior art date
Application number
PCT/JP2019/001902
Other languages
English (en)
French (fr)
Inventor
弘征 小松
純 本杉
Original Assignee
日産自動車株式会社
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by 日産自動車株式会社 filed Critical 日産自動車株式会社
Priority to PCT/JP2019/001902 priority Critical patent/WO2020152785A1/ja
Publication of WO2020152785A1 publication Critical patent/WO2020152785A1/ja

Links

Images

Classifications

    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02PCONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
    • H02P21/00Arrangements or methods for the control of electric machines by vector control, e.g. by control of field orientation
    • H02P21/05Arrangements 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
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02PCONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
    • H02P21/00Arrangements or methods for the control of electric machines by vector control, e.g. by control of field orientation
    • H02P21/24Vector control not involving the use of rotor position or rotor speed sensors
    • H02P21/26Rotor flux based control
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02PCONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
    • H02P25/00Arrangements or methods for the control of AC motors characterised by the kind of AC motor or by structural details
    • H02P25/02Arrangements or methods for the control of AC motors characterised by the kind of AC motor or by structural details characterised by the kind of motor
    • H02P25/022Synchronous motors

Definitions

  • the present invention relates to a motor control method and a motor control device.
  • the stator current is formed after the increase of the rotor current (field current) with the increase of the torque command value.
  • a technique is disclosed that controls the axial current to improve the responsiveness of the field current at the timing when the drive torque is increased such as when the motor starts to be driven.
  • the d-axis current is controlled on condition that the increase of the field current accompanying the increase of the torque command value is started. Therefore, the field current when decreasing the drive torque is controlled. Can not improve the responsiveness of.
  • the present invention aims to provide a technique capable of improving the responsiveness of a field current in a field winding type synchronous motor regardless of the driving torque increasing or decreasing direction.
  • a motor control method relates to a field winding type synchronous motor including a rotor having a rotor winding and a stator having a stator winding, and a stator current flowing through the stator winding. And a motor current controlling a rotor current flowing through the rotor winding.
  • a torque command value for generating a desired torque in the winding field type synchronous motor is set, and the stator is set based on the torque command value and the operating condition of the winding field type synchronous motor.
  • the d-axis current command value and the q-axis current command value for the current and the f-axis current command value for the rotor current are calculated, and the voltage is applied to the rotor winding, and the f-axis current command value is calculated.
  • the d-axis current command value is corrected according to the insufficient amount, and the d-axis current command value and the q-axis are calculated.
  • the winding field type synchronous motor is controlled based on the current command value and the f-axis current command value.
  • FIG. 1 is a block diagram of a motor control system to which a motor control method according to an embodiment is applied.
  • FIG. 2 is a block diagram showing the configuration of the current control unit according to the embodiment.
  • FIG. 3 is a block diagram showing the configuration of the f-axis pseudo non-interference voltage calculation unit according to the embodiment.
  • FIG. 4 is a block diagram showing the configuration of the d-axis model of one embodiment.
  • FIG. 5 is a block diagram showing the configuration of the d-axis current F/B model of one embodiment.
  • FIG. 6 is a block diagram showing the configuration of the d-axis pseudo non-interference voltage model of one embodiment.
  • FIG. 7 is a block diagram showing the configuration of the f-axis undervoltage calculating unit according to the embodiment.
  • FIG. 1 is a block diagram of a motor control system to which a motor control method according to an embodiment is applied.
  • FIG. 2 is a block diagram showing the configuration of the current control unit according to the embodiment.
  • FIG. 8 is a block diagram showing the configuration of the f-axis model of one embodiment.
  • FIG. 9 is a block diagram showing the configuration of the f-axis current F/B model of one embodiment.
  • FIG. 10 is a block diagram showing the configuration of the f-axis undervoltage calculation model according to the embodiment.
  • FIG. 11 is a block diagram showing the configuration of the d-axis current command value correction unit according to the embodiment.
  • FIG. 12 is a flowchart showing one control cycle of the control calculation according to the embodiment.
  • FIG. 13 is a time chart illustrating the effect of the motor control method according to the embodiment when the drive torque is decreasing.
  • FIG. 14 is a time chart illustrating the effect of the motor control method according to the embodiment when the driving torque is in the increasing direction.
  • FIG. 1 is a block diagram showing a configuration example of a motor control system 100 when a motor control method according to an embodiment of the present invention is applied to a winding field type synchronous motor.
  • the motor control system 100 is applied to, for example, an electric vehicle.
  • the motor control system 100 can be applied to a hybrid vehicle and a system other than the automobile, such as a railway.
  • the motor control system 100 of this embodiment includes a winding field type synchronous motor 101 to be controlled, a PWM converter 102, an inverter 103, a power supply voltage 104, an f-axis current output unit 105, and a current sensor 106. , A/D converter 107, coordinate converter 108, magnetic pole position detector 109, pulse counter 110, angular velocity calculator 111, prefetch compensator 112, current command value calculator 113, and current controller. 114, a non-interference control unit 115, a voltage command value calculation unit 116, and a coordinate converter 117.
  • a winding field type synchronous motor 101 (hereinafter, simply referred to as “motor 101”) includes a rotor having a rotor winding (field winding, rotor coil) and a stator winding (armature winding, stator). And a stator having a coil).
  • the motor 101 serves as a drive source for the vehicle.
  • the motor 101 is controlled by controlling a rotor current (field current) flowing through the rotor winding and a stator current flowing through the stator winding.
  • the PWM converter 102 drives PWM_Duty to a switching element (for example, an IGBT) included in the inverter 103 based on three-phase voltage command values v * u , v * v , and v * w output from a coordinate converter 117 described later.
  • Signals strong electric element drive signals
  • D * uu , D * ul , D * vu , D * vl , D * wu , and D * wl are generated and output to the inverter 103.
  • the inverter 103 is a three-phase voltage type inverter which is composed of three-phase six-arms and has two switching elements in total for each phase.
  • the inverter 103 converts the DC voltage of the power supply voltage 104 into AC voltages v u , v v , and v w based on the high-power element drive signal generated by the PWM converter 102, and supplies the AC voltages v u , v v , and v w to the motor 101.
  • the power supply voltage (battery) 104 is a power source for the inverter 3 and the f-axis current output unit 105, and is, for example, a laminated lithium ion battery.
  • the power supply voltage 104 supplies DC power to the inverter 3 and the f-axis current output unit 105.
  • the f-axis current output unit 105 uses the electric power supplied from the power supply voltage 104 to output the f-axis voltage v f for controlling the f-axis current if that flows through the rotor winding of the motor 101.
  • f-axis voltage v f is calculated in accordance with the second f-axis voltage command value v * f outputted from the voltage command value computing section 116 to be described later.
  • the f-axis current output unit 105 adjusts the f-axis voltage v f to be applied to the rotor winding to match the second f-axis voltage command value v * f.
  • Current sensor 106 of the three-phase alternating current supplied from the inverter 103 to the motor 101 at least two phases of current, for example, detecting the u-phase current i u, v-phase current i v.
  • the detected two-phase currents i u and iv are converted into digital signals (currents i us and i vs ) by the A/D (analog/digital) converter 107 and input to the coordinate converter 108.
  • the current sensor 106 detects only the two-phase current, the remaining one-phase current i ws can be obtained by the following equation (1).
  • the current sensor 106 detects a rotor current (f-axis current if ) supplied from the f-axis current output unit 105 to the motor 101.
  • the detected f-axis current if is converted into a digital signal by the A/D converter 107 and output to the f-axis current control unit 116.
  • the magnetic pole position detector 109 outputs, to the pulse counter 110, A-phase, B-phase, and Z-phase pulses (ABZ pulse) according to the position (angle) of the stator of the motor 101 in order to acquire the electrical angle of the stator. To do.
  • the pulse counter 110 calculates the electrical angle ⁇ re of the motor 101 based on the ABZ pulse, and outputs it to the angular velocity calculator 111.
  • the angular velocity calculator 111 calculates the electrical angular velocity ⁇ re and the mechanical angular velocity ⁇ rm from the time change rate of the input electrical angle ⁇ re .
  • the mechanical angular velocity ⁇ rm is obtained by dividing the electrical angular velocity ⁇ re by the motor pole pair number p.
  • the mechanical angular velocity ⁇ rm is output to the current command value calculation unit 113.
  • the electrical angular velocity ⁇ re is output to the non-interference control unit 115 and the look-ahead compensation unit 112.
  • the pre-reading compensation unit 112 inputs the electrical angle ⁇ re and the electrical angular velocity ⁇ re, and adds the product of the electrical angular velocity ⁇ re and the dead time of the control system to the electrical angle ⁇ re so that the pre-reading compensation is performed. Calculate the electrical angle ⁇ re '.
  • the pre-reading-compensated electrical angle ⁇ re ′ is output to the coordinate converter 117.
  • the coordinate converter 108 performs conversion from a three-phase AC coordinate system (uvw axis) to an orthogonal two-axis DC coordinate system (dq axes). Specifically, the coordinate converter 108 calculates the following equation from the input u-phase current i us , v-phase current i vs , the electrical angle ⁇ re, and the w-phase current i ws obtained by the above equation (1). The coordinate conversion process is performed using (2) to calculate the d-axis current i d and the q-axis current i q .
  • the current command value calculation unit 113 calculates the d-axis current command value i * d and the q-axis current command value i * q for the stator current and the rotor current (field) based on the torque command value and the operating condition of the motor 101.
  • the f-axis current command value i * f for the magnetic current is calculated.
  • the current command value calculation unit 113 receives the torque command value T * , the motor rotation speed (mechanical angular velocity ⁇ rm ), and the power supply voltage (DC voltage) V dc as input, and receives the d-axis current command value i * d , q-axis current command value i * q , and f-axis current command value i * f are calculated.
  • the d-axis current command value i * d , the q-axis current command value i * q , and the f-axis current command value i * f are the previously stored torque command value T * , motor rotation speed (mechanical angular velocity ⁇ rm ), By referring to the map data that defines the relationship between the power supply voltage V dc and each of the d-axis current command value i * d , the q-axis current command value i * q , and the f-axis current command value i * f. Desired.
  • the torque command value T * is a command value for generating a desired torque in the motor 101, and is set by a known method according to the accelerator operation of the driver.
  • the current control unit 114 the actual current (actual current) d-axis voltage command value of the d-axis current i d is a measured value in order to follow in steady-state error without desired response to d-axis current command value i * d of The (first d-axis voltage command value) is calculated and output to the voltage command value calculation unit 116.
  • the current control unit 114 causes the q-axis current i q , which is a measured value of the actual current (actual current), to follow the q-axis current command value i * q with a desired responsiveness without a steady deviation.
  • a command value (first q-axis voltage command value) is calculated and output to the voltage command value calculator 116.
  • the current control unit 114 the actual current (actual current) measurements in which f axis current i f the f-axis voltage to follow a steady deviation without the desired response in the f-axis current command value i * f
  • a command value (first f-axis current command value) is calculated and output to the voltage command value calculation unit 116.
  • Decoupling control unit 115 and the electrical angular velocity omega re inputted, the d-axis current i d and the q-axis current i q, and a f-axis current i f, d-axis, q-axis, and interference between the f-axis non-interference voltage v D_dcpl necessary to offset voltage, v q_dcpl, v calculates the F_dcpl.
  • the non-interference control section 115 is a measure of the actual current (actual current) d-axis current i d, instead of the q-axis current i q, and f-axis current i f, d-axis current command value i *
  • a d-axis current reference response i d_ref , a q-axis current reference response i q_ref , and an f-axis current reference response i which will be described later, are calculated based on d , the q-axis current reference value i * q , and the f-axis current reference value i * f.
  • the non-interference voltages v d_dcpl , v q_dcpl , and v f_dcpl may be calculated using f_ref (actual current).
  • the voltage command value calculation unit 116 outputs the first d-axis voltage command value v d_dsh , the first q-axis voltage command value v q_dsh , and the first f-axis voltage command value v f_dsh, which are the outputs of the current control unit 114. Is corrected (added in this embodiment) by using the non-interference voltages v d_dcpl , v q_dcpl , and v f_dcpl which are the outputs of the non-interference control unit 115.
  • the voltage command value calculation unit 116 outputs the second d-axis voltage command value v * d and the second q-axis voltage command value v * q obtained by the correction to the coordinate converter 117 and , And outputs the second f-axis voltage command value v * f to the f-axis current output unit 105.
  • the coordinate converter 117 performs conversion from the orthogonal two-axis DC coordinate system (dq axes) rotating at the electrical angular velocity ⁇ re to the three-phase AC coordinate system (uvw phase). Specifically, the coordinate converter 117 calculates the second d-axis voltage command value v * d , the second q-axis voltage command value v * q , and the pre-reading-compensated electrical angle ⁇ re ′ from the input second d-axis voltage command value v * d .
  • the voltage command values v * u , v * v , and v * w for each phase of uvw are calculated by performing coordinate conversion processing using the following equation (3).
  • the current controller 114, the non-interference controller 115, the voltage command value calculator 116, and the coordinate converter 117 are configured as one functional unit included in the controller 10 including at least one controller.
  • the controller 10 includes, for example, a central processing unit (CPU), a read only memory (ROM), a random access memory (RAM), and an input/output interface (I/O interface).
  • each parameter of the above equation (4) is as follows.
  • s in the formula is a Laplace operator.
  • i d d-axis current
  • i q q-axis current if : f-axis current v d : d-axis voltage
  • v q q-axis voltage
  • v f f-axis voltage
  • L d d-axis inductance
  • L q q-axis inductance
  • L f f-axis inductance
  • M mutual inductance L d between the stator / rotor ': d-axis dynamic inductance
  • L q' q-axis dynamic inductance
  • L f ' f-axis dynamic inductance
  • M' a stator / rotor Dynamic mutual inductance
  • R a Stator winding resistance
  • the voltage equation of the above equation (4) can be diagonalized as shown in the following equation (5).
  • the characteristics from the voltage to the current on the d-axis, the q-axis, and the f-axis are first-order lags as shown in the following equations (6), (7), and (8), respectively. ..
  • FIG. 2 is a block diagram showing the current control unit 114 of this embodiment.
  • the current control unit 114 includes an f-axis pseudo non-interference voltage calculation unit 201, an f-axis undervoltage calculation unit 202, a d-axis current command value correction unit 203, a d-axis current control unit 204, and a q-axis current control unit 205. And an f-axis current control unit 206.
  • the f-axis pseudo non-interference voltage calculation unit 201 receives the d-axis current command value i * d output from the current command value calculation unit 113 as an input and calculates the f-axis pseudo non-interference voltage v f_dcpl_dsh . Details will be described with reference to FIG.
  • FIG. 3 is a block diagram showing the f-axis pseudo non-interference voltage calculation unit 201.
  • the f-axis pseudo non-interference voltage calculation unit 201 includes a d-axis model 301, a d-axis current feedback (F/B) model 302, and an f-axis pseudo non-interference voltage model 303.
  • F/B d-axis current feedback
  • the d-axis model 301 receives the output value of a d-axis current F/B model 302 described later as an input, and performs a filtering process using a d-axis model that models the characteristics from the d-axis voltage v d to the d-axis current i d. As a result, the d-axis current reference response i d_ref and the differential value s ⁇ i d_ref of the d-axis current reference response are output. Details of the d-axis model 301 will be described with reference to FIG.
  • FIG. 4 is a block diagram showing the d-axis model 301.
  • the d-axis model 301 includes a multiplier 401, a subtractor 402, a divider 403, and an integrator 404.
  • the multiplier 401 multiplies the d-axis current reference response i d_ref , which is the output value of the integrator 404, by the stator winding resistance R a , which is the resistance value of the stator winding included in the motor 101, and calculates the calculated value. Output to the subtractor 402.
  • the subtractor 402 subtracts the output value of the multiplier 401 from the output value of the d-axis current F/B model 302, and outputs the calculated value to the divider 403.
  • the divider 403 divides the output value of the subtractor 402 by the d-axis dynamic inductance L d ′ to calculate the differential value s ⁇ id_ref of the d-axis current reference response.
  • the differential value s ⁇ id_ref of the d-axis current reference response is output to the integrator 404 and the f-axis pseudo non-interference voltage model 303 described later.
  • the integrator 404 calculates the d-axis current nominal response i d_ref by integration of the differential value s ⁇ i d_ref the d-axis current nominal response.
  • the d-axis current reference response i d_ref is output to the multiplier 401 and the d-axis current F/B model 302. In this way, the output of the integrator 404 is used as the current reference response, and the input of the integrator 404 is configured to be obtainable as the differential value of the current reference response, so that the differential value of the current reference response can be easily obtained. it can.
  • the d-axis current F/B model 302 receives the d-axis current command value i * d and the d-axis current normative response i d_ref , which is the output of the d-axis model 301 described above, as input.
  • the d-axis current F/B model output value (d-axis voltage estimated value) for causing the d-axis current normative response i d_ref to follow the value i * d with a desired responsiveness without a steady deviation is calculated, and the d-axis model 301 Output to. Details will be described with reference to FIG.
  • FIG. 5 is a block diagram showing details of the d-axis current F/B model 302.
  • the d-axis current F/B model 302 includes multipliers 501 and 502 and a subtractor 503.
  • the multiplier 501 multiplies the d-axis current command value i * d output from the current command value calculation unit 113 by the gain G ad , and outputs the calculated value to the subtractor 503.
  • the gain G ad is expressed by the following equation (9).
  • the multiplier 502 multiplies the d-axis current reference response id_ref output from the d-axis model 301 by the gain G bd , and outputs the calculated value to the subtractor 503.
  • the gain G bd is expressed by the following equation (9).
  • ⁇ m in the equation (9) is a d-axis and q-axis current control reference response time constant, that is, a time constant corresponding to a desired response characteristic to be realized.
  • the subtractor 503 calculates the d-axis current F/B model output value by subtracting the output value of the multiplier 502 from the output value of the multiplier 501, and outputs it to the above-mentioned d-axis model 301.
  • the f-axis pseudo non-interference voltage model 303 (see FIG. 3) inputs the differential value s ⁇ id_ref of the d-axis current reference response, which is the output of the d-axis model 301, and outputs the f-axis pseudo non-interference voltage v f_dcpl_dsh . To do. Details of the f-axis pseudo non-interference voltage model 303 will be described with reference to FIG.
  • FIG. 6 is a block diagram showing the f-axis pseudo non-interference voltage model 303.
  • the f-axis pseudo non-interference voltage model 303 is composed of a multiplier 601.
  • the multiplier 601 multiplies the differential value s ⁇ id_ref of the d-axis current reference response, which is the output of the d-axis model 301, by the dynamic mutual inductance M′ between the stator and the rotor, thereby f-axis pseudo non-interference.
  • the voltage v f_dcpl_dsh is calculated.
  • the f-axis pseudo non-interference voltage v f_dcpl_dsh is output to the f-axis undervoltage calculator 202 (see FIG. 2 ).
  • the f-axis undervoltage calculator 202 inputs the f-axis pseudo non-interference voltage v f_dcpl_dsh , the power supply voltage V dc, and the f-axis current command value i * f, and calculates the f-axis undervoltage ⁇ v f . Details of the f-axis undervoltage calculation unit 202 will be described with reference to FIG. 7.
  • FIG. 7 is a block diagram showing the f-axis undervoltage calculation unit 202.
  • the f-axis undervoltage calculation unit 202 includes an f-axis model 701, an f-axis current F/B model 702, and an f-axis undervoltage calculation model 703.
  • the f-axis model 701 is a model in which the characteristics from the f-axis voltage to the f-axis current are modeled by inputting the f-axis current F/B model output value output from the f-axis current F/B model 702, which will be described later. perform f-axis model) by filtering processing of the motor 101, and outputs the f-axis current nominal response i F_REF, and a differential value s ⁇ i F_REF of f-axis current nominal response. Details will be described with reference to FIG.
  • FIG. 8 is a block diagram showing details of the f-axis model 701.
  • the f-axis model 701 includes a multiplier 801, a subtractor 802, and an integrator 804.
  • the multiplier 801 multiplies the f-axis current reference response if_ref output from the integrator 804 by the stator winding resistance Rf , and outputs the obtained value to the subtractor 802.
  • the subtractor 802 subtracts the output value of the multiplier 801 from the output value of the f-axis current F/B model output from the f-axis current F/B model 702, and outputs the calculated value to the divider 803.
  • the divider 803 divides the output value of the subtractor 802 by the f-axis dynamic inductance L f ′ to calculate the differential value s ⁇ if_ref of the f-axis current reference response.
  • the differential value s ⁇ if_ref of the f-axis current reference response is output to the integrator 804.
  • the integrator 804 calculates an f axis current nominal response i F_REF by integration of the differential value s ⁇ i F_REF of f-axis current nominal response.
  • the f-axis current reference response if_ref is output to the multiplier 801 and the f-axis current F/B model 702.
  • the f-axis current F/B model 702 receives the f-axis current command value i * f and the f-axis current normative response if_ref , which is the output of the f-axis model 701 described above, as input.
  • An f-axis current F/B model output value (f-axis voltage estimated value) for causing the f-axis current reference response if_ref to follow the value i * f with desired responsiveness without steady deviation is calculated. Details of the f-axis current F/B model 702 will be described with reference to FIG. 9.
  • FIG. 9 is a block diagram showing the f-axis current F/B model 702.
  • the f-axis current F/B model 702 is configured to include multipliers 901 and 902 and a subtractor 903.
  • the multiplier 901 multiplies the f-axis current command value i * f output from the current command value calculation unit 113 by the gain G af , and outputs the calculated value to the subtractor 903.
  • the gain G af is expressed by the following equation (10).
  • the multiplier 902 multiplies the f-axis current reference response i * f_ref output from the f-axis model 701 by the gain G bf , and outputs the calculated value to the subtractor 903.
  • the gain G bf is expressed by the following equation (10).
  • ⁇ m in the expression (10) is a current control reference response time constant of the f-axis, that is, a time constant corresponding to a desired response characteristic to be realized.
  • the subtractor 903 calculates the f-axis current F/B model output value by subtracting the output value of the multiplier 902 from the output value of the multiplier 901, and calculates the f-axis model 701 and the f-axis undervoltage calculation model 703. And output to.
  • the f-axis current F/B model output value is an estimated value of the f-axis voltage v f (f-axis voltage estimated value).
  • the f-axis undervoltage calculation model 703 receives the f-axis current F/B model output value, the power supply voltage V dc, and the f-axis pseudo non-interference voltage v f_dcpl_dsh , and inputs the f-axis current F/B.
  • the f-axis undervoltage ⁇ v f obtained by limiting the model output value with a predetermined value (limit value) is output. More specifically, the f-axis undervoltage calculation model 703 corresponds to the f-axis voltage vf required to control the f-axis current if at a desired response speed with respect to the f-axis current command value i * f .
  • the shortage of the outputtable voltage with respect to the shaft current F/B model output value is calculated as the f-axis undervoltage ⁇ v f .
  • the outputtable voltage here is a voltage that is determined according to the power supply voltage Vdc as described later, and is a voltage that the f-axis current output unit 105 can output. Details of the f-axis undervoltage calculation model 703 will be described with reference to FIG. 10.
  • FIG. 10 is a block diagram showing the f-axis undervoltage calculation model 703.
  • the f-axis undervoltage calculation model 703 is configured to include a comparator 1001, an inverter 1002, a comparator 1003, and subtractors 1004, 1005, 1006.
  • the subtractor 1005 provided in the preceding stage of the comparator 1001 obtains a subtraction value obtained by subtracting the f-axis pseudo non-interference voltage v f_dcpl_dsh from the power supply voltage V dc of the power supply voltage 104. Then, the comparator 1001 compares the output value from the f-axis current pseudo F/B model 702 with the output value of the subtractor 1005, and outputs a smaller value to the comparator 1003.
  • the inverter 1002 inverts the sign of the power supply voltage V dc .
  • a subtractor 1006 is provided in the preceding stage of the comparator 1003, and the subtractor 1006 obtains a subtraction value obtained by subtracting the f-axis pseudo non-interference voltage v f_dcpl_dsh from the output of the inverter 1002. Then, the comparator 1003 compares the output value of the comparator 1001 with the output value of the subtractor 1006, and outputs a larger value to the subtractor 1004.
  • the f-axis undervoltage calculation model 703 can execute limit processing on the output value of the f-axis current F/B model 702 (subtractor 903) based on a predetermined value (limit value).
  • the voltage between the upper limit value and the lower limit value corresponds to the above-mentioned outputtable voltage.
  • the calculated value is the power supply voltage V It is possible to avoid saturation with respect to dc . Further, by executing the limit processing as described above, there is no need for a means for determining the presence/absence of saturation of the f-axis voltage, a filter switching means for switching the filter depending on the presence/absence of the saturation of the f-axis voltage, and the like. An appropriate f-axis voltage command value can be calculated.
  • the f-axis pseudo non-interference voltage v f_dcpl_dsh is calculated based on the d-axis current command value i * d before being corrected by the d-axis current command value correction unit 203 in the f-axis pseudo non-interference voltage calculation unit 201. It is calculated. That is, the predetermined upper and lower limit values (limit values) in the present embodiment are the f-axis pseudo non-interference voltage v f_dcpl_dsh calculated based on the power supply voltage V dc and the uncorrected d-axis current command value i * d. It is calculated based on
  • the subtractor 1004 calculates the f-axis undervoltage ⁇ v f by subtracting the output value of the comparator 1003 from the f-axis current F/B model output value. That is, the f-axis undervoltage ⁇ v f depends on the difference between the f-axis current F/B model output value and the value obtained by limiting the f-axis current F/B model output value with predetermined upper and lower limits. It is determined. Thereby, the shortage of the outputtable voltage with respect to the f-axis voltage vf for controlling the f-axis current with a desired response can be calculated.
  • the f-axis undervoltage ⁇ v f is output to the d-axis current command value correction unit 203.
  • the gain (901, 902) of the f-axis current F/B model 702 is stabilized by calculating the f-axis undervoltage ⁇ v f with such a configuration and controlling the f-axis current with a desired responsiveness. Since it is not necessary to set the gains (901, 902) to be small in consideration of the above, it is possible to set the gains (901, 902) to the gains that satisfy the desired responsiveness (see Formula (10)).
  • the d-axis current command value correction unit 203 receives the d-axis current command value i * d and the f-axis undervoltage ⁇ v f as input, and based on the f-axis undervoltage ⁇ v f , the d-axis current command value i * Correct d .
  • the corrected d-axis current command value i * d is output to the d-axis current control unit 204 as a corrected d-axis current command value. Details of the d-axis current command value correction unit 203 will be described with reference to FIG. 11.
  • FIG. 11 is a block diagram showing the d-axis current command value correction unit 203.
  • the d-axis current command value correction unit 203 includes a filter 1101 and a subtractor 1102.
  • the filter 1101 is composed of a function of H d — comp (s).
  • the filter 1101 performs filtering processing using H d_comp (s) on the input f-axis undervoltage ⁇ v f , and outputs the d-axis current correction value i d_comp to the subtractor 1102. Details of H d_comp (s) will be described later.
  • the subtractor 1102 subtracts the d-axis current correction value i d_comp from the d-axis current command value i * d to calculate the corrected d-axis current command value.
  • the d-axis current command value i * d is corrected using H d — comp (s).
  • Equation (11) in the current (d-axis current i d, f-axis current i f) substitution with a nominal response f-axis voltage v f is expressed by the following equation (12).
  • the f-axis voltage saturation can be suppressed by canceling out the f-axis undervoltage ⁇ v f in the second term on the right side of the equation (13).
  • the d-axis current correction value id_comp at this time is given by the following equation (14).
  • the d-axis current command value i is determined according to the f-axis undervoltage ⁇ v f , that is, according to the shortage of the outputtable voltage with respect to the f-axis voltage v f . * d can be corrected.
  • the response characteristic of the f-axis current if can be improved, the f-axis current if can be operated with a desired response (response speed), and the torque response can be improved. Can be made.
  • the d-axis current command value When the f-axis undervoltage ⁇ vf is positive, the d-axis current command value is corrected to the negative side (decreasing direction), and when the f-axis undervoltage ⁇ vf is negative, the d-axis current command value is increased (increase). Since it is configured to correct the torque command value, the f-axis current if can be operated with a desired responsiveness (response speed) regardless of the increasing or decreasing direction of the torque command value.
  • the d-axis current command value i * d is corrected only when the f-axis voltage is saturated.
  • the correction term since the correction term has the integral characteristic (1/s), the d-axis current command value i * d is corrected even when the f-axis voltage is not saturated. Therefore, as shown in the following expression (15), the integration characteristic is canceled by adding a high-pass filter to the above expression (14). Accordingly, H d — comp (s) can be configured to correct the d-axis current command value i * d only when the f-axis voltage is saturated.
  • the time constant ⁇ hr of the high-pass filter in Expression (15) is a tuning element and is appropriately adjusted to obtain a desired response.
  • the above is the details of H d — comp (s).
  • the corrected d-axis current command value calculated using such a function is output to the d-axis current control unit 204 (see FIG. 2).
  • the d-axis current control unit 204 is a first unit for causing the d-axis current i d , which is a measured value of the actual current (actual current), to follow the d-axis current command value i * d with a desired responsiveness without a steady deviation.
  • the d-axis voltage command value v d_dsh is calculated and output to the voltage command value calculation unit 116.
  • the q-axis current control unit 205 is a first unit for causing the q-axis current i q , which is a measured value of an actual current (actual current), to follow the q-axis current command value i * q with a desired responsiveness without a steady deviation.
  • the q-axis voltage command value v q_dsh is calculated and output to the voltage command value calculation unit 116.
  • the actual current (actual current) a measure of the f-axis current i f the f-axis current command value i * f the steady-state error without first for follow a desired response of
  • the f-axis voltage command value v f_dsh is calculated and output to the voltage command value calculation unit 116.
  • the above-mentioned d-axis current control unit 204, q-axis current control unit 205, and f-axis current control unit 206 ideally perform control (non-interference control) for canceling the interference voltage between the d-axis, the q-axis, and the f-axis. Functioning normally, it has a simple characteristic of 1 input and 1 output (controlled object characteristic) (see the above equations (6), (7), and (8)). Therefore, the d-axis current control unit 204, the q-axis current control unit 205, and the f-axis current control unit 206 of the present embodiment are known as simple PI feedback compensators or so-called robust model matching compensators. It can be realized by a compensator.
  • FIG. 12 is a flowchart showing a motor control method according to an embodiment.
  • the one control cycle from the start to the end shown in FIG. 12 is programmed in the controller 10 so as to be constantly executed at regular intervals while the motor control system 100 is activated.
  • step S11 the controller 10 (pulse counter 110, A/D converter 107) calculates the electrical angle ⁇ re based on the ABZ pulse, and obtains the acquired u and v phase currents i u , iv , and from f-axis current i f that flows through the rotor winding, for detecting a current i us, i vs, i fs as a digital signal.
  • step S12 the controller 10 (angular velocity calculator 111) calculates the electrical angular velocity ⁇ re and the mechanical angular velocity ⁇ rm from the temporal change rate of the input electrical angle ⁇ re .
  • step S13 the controller 10 (look-ahead compensation unit 112) calculates a look-ahead-compensated electrical angle ⁇ re ′ in which the dead time of the control system is taken into consideration, from the electrical angle ⁇ re and the electrical angular velocity ⁇ re .
  • step S14 the controller 10 (the coordinate converter 108), u, v, w of each phase current i us, i vs, i ws , and from the electric angle ⁇ re, d, q-axis current i d, a i q calculate.
  • step S15 the controller 10 (current command value calculation unit 113) determines the d-axis current command values i * d , q from the torque command value T * , the motor rotation speed (mechanical angular velocity ⁇ rm ), and the power supply voltage V dc.
  • the axis current command value i * q and the f-axis current command value i * f are calculated.
  • step S16 the controller 10 (current control unit 114) sets the d, q, f-axis currents id , iq , if to the d, q, f-axis current command values i * d , i * q , i * f.
  • the first voltage command values v d_dsh , v q_dsh , v f_dsh for each phase of the d, q, and f axes are calculated so as to follow the desired response with no steady deviation.
  • step S17 the controller 10 (non-interference control unit 115) determines between the d-axis, the q-axis, and the f-axis from the d, q, f-axis currents i d , i q , if, and the electrical angular velocity ⁇ re .
  • the non-interference voltages v d_dcpl , v q_dcpl , and v f_dcpl necessary to cancel the interference voltage are calculated.
  • step S18 the controller 10 (voltage command value calculation unit 116) sets the first voltage command values v d_dsh , v q_dsh , v f_dsh for the d, q, and f axes to the non-interference voltages v d_dcpl , v q_dcpl , v.
  • the second voltage command values v * d , v * q , and v * f for the d, q, and f axes are calculated by performing correction using f_dcpl .
  • step S19 the controller 10 (coordinate converter 117) causes the second d-axis voltage command value v * d , the second q-axis voltage command value v * q , and the pre-reading compensated electrical angle ⁇ re ′. From this, the voltage command values v * u , v * v , v * w for each phase of uvw are calculated.
  • FIG. 13 is a time chart showing a control result by the motor control method of this embodiment.
  • the horizontal axis represents time, and the vertical axis represents d-axis current [A], q-axis current [A], and f-axis current [A] in order from the left side, and f-axis in order from the right side.
  • the voltage [V] and the torque [N] are shown.
  • the broken lines in the figure relating to the d-axis current [A], the q-axis current [A], and the f-axis current [A] indicate current command values.
  • the broken line in the figure relating to the torque [N] indicates the torque command value.
  • the response of the actual current or the actual torque to the current command value or the torque command value is shown by a solid line in this embodiment, and the control by the prior art (conventional example) is shown by a chain line.
  • FIG. 13 shows a comparison result between the present embodiment and the conventional example when the torque command value is in the decreasing direction.
  • the d-axis current command value i * d the q-axis current command value i * q , and the f-axis current command value i * f are changed as the torque command value T* is applied in a decreasing direction in steps. Applied in steps.
  • the d-axis current i d , the q-axis current i q , and the f-axis current if are the applied d-axis current command value i * d , q-axis current command value i * q , and f-axis current command value i. * f converges to the goal.
  • the f-axis voltage v f is saturated by the output enable lower limit voltage (-Vdc)
  • the response of the f-axis current i f (response speed) is lowered.
  • the torque response also deteriorates.
  • a voltage (f-axis F/B model output value) for controlling the f-axis current if with a desired response and an outputtable voltage based on the power supply voltage V dc. depending on a difference f axis undervoltage Delta] v f, since the d-axis current i d transiently corrected, saturated f-axis voltage v f is suppressed, can be quickly responses f axis current i f it can. As a result, as shown in the figure, the torque response can be made faster than in the conventional example.
  • FIG. 14 shows that the motor control of the present embodiment when the torque command value is increasing and the function of correcting the d-axis current id according to the f-axis undervoltage ⁇ v f is turned off. It is a comparison result with the motor control of the embodiment (comparative example) in the case.
  • the horizontal axis represents time
  • the vertical axis represents the d-axis current [A], the q-axis current [A], and the f-axis current [A] in order from the left side.
  • the f-axis voltage [V] and the torque [N] are shown in order from the top on the right side.
  • the broken lines in the figure relating to the d-axis current [A], the q-axis current [A], and the f-axis current [A] indicate current command values.
  • the broken line in the figure relating to the torque [N] indicates the torque command value.
  • the response of the actual current or the actual torque to the current command value or the torque command value is shown by a solid line in this embodiment, and the control according to the comparative example is shown by a chain line.
  • the d-axis current command value i * d the q-axis current command value i * q , and the f-axis current command value i * f are increased as the torque command value T * is applied in the increasing direction in steps. Applied in steps.
  • the d-axis current i d , the q-axis current i q , and the f-axis current if are the applied d-axis current command value i * d , q-axis current command value i * q , and f-axis current command value i. * f converges to the goal.
  • the f-axis voltage v f is saturated by the output enable upper limit voltage (V dc)
  • the response of the f-axis current i f (response speed) is lowered.
  • the torque response also deteriorates.
  • a voltage (f-axis F/B model output value) for controlling the f-axis current if with a desired response and an outputtable voltage based on the power supply voltage V dc. depending on a difference f axis undervoltage Delta] v f, since the d-axis current i d transiently corrected, saturated f-axis voltage v f is suppressed, can be quickly responses f axis current i f it can. As a result, as shown in the figure, the torque response can be made faster than in the conventional example.
  • the motor control method that corrects the d-axis current command value according to the shortage of the f-axis voltage, the f-axis voltage saturation can be suppressed and the f-axis current can be reduced to a desired value. It can be operated responsively. As a result, the torque responsiveness can be improved.
  • the d-axis current command value is continuously corrected according to the shortage of the f-axis voltage, a non-continuous parameter by a switch or the like depending on the increasing/decreasing direction of the torque command value or It is possible to improve the torque responsiveness regardless of the increase/decrease direction of the torque command value without switching the specifications.
  • the motor control method is a stator that flows through a stator winding in a winding field type synchronous motor 101 that includes a rotor having a rotor winding and a stator having a stator winding.
  • the torque command value T * for generating a desired torque in the winding field type synchronous motor 101 is set, and the torque command value T * and the operating condition of the winding field type synchronous motor 101 are set.
  • the d-axis current command value i * d and the q-axis current command value i * q for the stator current and the f-axis current command value i * f for the rotor current are calculated.
  • the voltage applied to the rotor winding which is the f-axis voltage vf (f-axis current F/B model) required to control the rotor current at a predetermined response speed with respect to the f-axis current command value.
  • the d-axis current command value i * d is corrected, and the winding field type synchronous motor 101 is based on the d-axis current command value i * d , the q-axis current command value i * q , and the f-axis current command value i * f. To control.
  • the d-axis current command value i * d is corrected according to the shortage of the outputtable voltage with respect to the f-axis voltage v f, so that the f-axis voltage saturation is suppressed and the f-axis current if is desired responsiveness. Can be operated with. As a result, the torque responsiveness can be improved regardless of the increasing or decreasing direction of the torque command value.
  • the d-axis current command value i * d is corrected to the negative side to obtain the f-axis voltage.
  • the d-axis current command value i * d is corrected to the positive side.
  • the d-axis current command value i * d is appropriately corrected according to the increasing/decreasing direction of the torque command value, and the f-axis current if is calculated regardless of the increasing/decreasing direction of the torque command value. It can be operated with desired responsiveness.
  • the f-axis model 701 that models the characteristics of up to f axis current i f that constitute the rotor current from the f-axis voltage v f, f-axis current command value i *
  • the f- axis current feedback (F/B) model 702 to which f and the output of the f-axis model 701 are input forms a pseudo F/B system, and the output value of the f-axis current F/B model (f When the value obtained by limiting the output value of the shaft current F/B model (f-axis voltage estimated value) with the predetermined upper and lower limits is insufficient for the output value of the f-axis current F/B model.
  • the f-axis current F/B model 702 is considered in consideration of stability. Since it is not necessary to set the gains (901, 902) of (1) to be small, these gains (901, 902) can be set to gains that satisfy the desired responsiveness (see formula (10)).
  • the f-axis pseudo non-interference voltage v f_dcpl_dsh for decoupling the interference voltage between the d-axis of the stator current and the f-axis of the rotor current is corrected. It is calculated based on the previous d-axis current command value i * d , and the above predetermined upper and lower limit values are calculated based on the power supply voltage V dc and the f-axis pseudo non-interference voltage v f_dcpl_dsh .
  • the shortage of the power supply voltage V dc with respect to the f-axis voltage v f can be appropriately calculated.
  • the calculated value is saturated with respect to the power supply voltage V dc .
  • the embodiment of the present invention has been described above, the above embodiment merely shows a part of the application example of the present invention, and the technical scope of the present invention is limited to the specific configuration of the above embodiment. Absent.
  • the present invention is not limited to the above-described embodiments, and various modifications and applications are possible.
  • the pulse counter 110 does not need to be configured as one functional unit of the controller 10, and may be provided separately from the controller 10 as a device having a pulse counting function.

Landscapes

  • Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Control Of Ac Motors In General (AREA)

Abstract

モータ制御方法は、回転子巻線を有する回転子と固定子巻線を有する固定子とを備える巻線界磁型同期モータにおいて、固定子巻線に流れる固定子電流と回転子巻線に流れる回転子電流とを制御するモータ制御方法である。このモータ制御方法によれば、巻線界磁型同期モータに所望のトルクを発生させるトルク指令値を設定し、トルク指令値と巻線界磁型同期モータの動作状況とに基づいて、固定子電流に対するd軸電流指令値及びq軸電流指令値と、回転子電流に対するf軸電流指令値とを算出する。そして、回転子巻線に印加される電圧であって、f軸電流指令値に対して回転子電流を所定の応答速度で制御するのに必要なf軸電圧を推定し、推定されたf軸電圧に対して電源電圧に基づいて定まる出力可能電圧が不足した場合には、当該不足分に応じてd軸電流指令値を補正し、d軸電流指令値、q軸電流指令値、及びf軸電流指令値に基づいて巻線界磁型同期モータを制御する。

Description

モータ制御方法、及び、モータ制御装置
 本発明は、モータ制御方法、及び、モータ制御装置に関する。
 JP2016-158429Aでは、界磁巻線型の同期モータをベクトル制御で駆動するモータ制御装置において、トルク指令値の増加に伴う回転子電流(界磁電流)の増加の開始後に固定子電流を構成するd軸電流を制御して、モータの駆動開始時等の駆動トルクを増加させるタイミングにおける界磁電流の応答性を改善する技術が開示されている。
 ここで、界磁巻線型の同期モータでは、駆動トルクを増加させる時だけではなく、駆動トルクを減少させる時にも界磁電流の応答性が悪化するという課題がある。しかしながら、JP2016-158429Aに開示された技術では、トルク指令値の増加に伴う界磁電流の増加の開始後を条件としてd軸電流を制御しているので、駆動トルクを減少させる際における界磁電流の応答性を改善することはできない。
 本発明は、界磁巻線型の同期モータにおいて、駆動トルクの増加方向、減少方向に関わらず、界磁電流の応答性を改善することができる技術を提供することを目的とする。
 本発明の一態様におけるモータ制御方法は、回転子巻線を有する回転子と、固定子巻線を有する固定子とを備える巻線界磁型同期モータにおいて、固定子巻線に流れる固定子電流と回転子巻線に流れる回転子電流とを制御するモータ制御方法である。このモータ制御方法によれば、巻線界磁型同期モータに所望のトルクを発生させるトルク指令値を設定し、トルク指令値と巻線界磁型同期モータの動作状況とに基づいて、固定子電流に対するd軸電流指令値及びq軸電流指令値と、回転子電流に対するf軸電流指令値とを算出し、回転子巻線に印加される電圧であって、f軸電流指令値に対して回転子電流を所定の応答速度で制御するのに必要なf軸電圧を推定する。そして、推定されたf軸電圧に対して電源電圧に基づいて定まる出力可能電圧が不足した場合には、当該不足分に応じてd軸電流指令値を補正し、d軸電流指令値、q軸電流指令値、及びf軸電流指令値に基づいて巻線界磁型同期モータを制御する。
 本発明の実施形態については、添付された図面とともに以下に詳細に説明する。
図1は、一実施形態のモータ制御方法が適用されるモータ制御システムのブロック図である。 図2は、一実施形態の電流制御部の構成を示すブロック図である。 図3は、一実施形態のf軸疑似非干渉電圧算出部の構成を示すブロック図である。 図4は、一実施形態のd軸モデルの構成を示すブロック図である。 図5は、一実施形態のd軸電流F/Bモデルの構成を示すブロック図である。 図6は、一実施形態のd軸疑似非干渉電圧モデルの構成を示すブロック図である。 図7は、一実施形態のf軸不足電圧算出部の構成を示すブロック図である。 図8は、一実施形態のf軸モデルの構成を示すブロック図である。 図9は、一実施形態のf軸電流F/Bモデルの構成を示すブロック図である。 図10は、一実施形態のf軸不足電圧算出モデルの構成を示すブロック図である。 図11は、一実施形態のd軸電流指令値補正部の構成を示すブロック図である。 図12は、一実施形態の制御演算の一制御周期を示すフローチャートである。 図13は、駆動トルクが減少方向にある場合における一実施形態のモータ制御方法の効果を説明するタイムチャートである。 図14は、駆動トルクが増加方向にある場合における一実施形態のモータ制御方法の効果を説明するタイムチャートである。
 [一実施形態]
 図1は、本発明の一実施形態に係るモータ制御方法を巻線界磁型同期モータに適用した場合のモータ制御システム100の構成例を示すブロック図である。モータ制御システム100は、例えば電気自動車に適用される。なお、モータ制御システム100は、電気自動車以外に、ハイブリッド車両や、自動車以外の例えば鉄道等のシステムに適用することも可能である。
 本実施形態のモータ制御システム100は、制御対象の巻線界磁型同期モータ101と、PWM変換器102と、インバータ103と、電源電圧104と、f軸電流出力部105と、電流センサ106と、A/D変換器107と、座標変換器108と、磁極位置検出器109と、パルスカウンタ110と、角速度演算器111と、先読み補償部112と、電流指令値演算部113と、電流制御部114と、非干渉制御部115と、電圧指令値演算部116と、座標変換器117と、を備える。
 巻線界磁型同期モータ101(以下、単に「モータ101」という)は、回転子巻線(界磁巻線、ロータコイル)を有する回転子と、固定子巻線(電機子巻線、ステータコイル)を有する固定子とを備える巻線界磁型同期モータである。本実施形態のモータ制御システム100が車両に搭載される場合、モータ101は車両の駆動源となる。詳細は後述するが、モータ101は、回転子巻線を流れる回転子電流(界磁電流)と、固定子巻線を流れる固定子電流とが制御されることによって制御される。
 PWM変換器102は、後述の座標変換器117から出力される三相電圧指令値v* u、v* v、v* wに基づいて、インバータ103が備えるスイッチング素子(例えばIGBT)へのPWM_Duty駆動信号(強電素子駆動信号)D* uu、D* ul、D* vu、D* vl、D* wu、D* wlを生成し、インバータ103に出力する。
 インバータ103は、3相6アームで構成され、相ごとに2つずつ計6つのスイッチング素子を備えた三相電圧型インバータである。インバータ103は、PWM変換器102が生成する強電素子駆動信号に基づいて、電源電圧104の直流電圧を交流電圧vu、vv、vwに変換し、モータ101に供給する。
 電源電圧(バッテリ)104は、インバータ3とf軸電流出力部105の電力源であって、例えば積層型リチウムイオンバッテリである。電源電圧104は、インバータ3とf軸電流出力部105とに直流電力を供給する。
 f軸電流出力部105は、電源電圧104から供給される電力を用いて、モータ101の回転子巻線に流れるf軸電流ifを制御するためのf軸電圧vfを出力する。f軸電圧vfは、後述する電圧指令値演算部116から出力される第2のf軸電圧指令値v* fに応じて算出される。換言すると、f軸電流出力部105は、回転子巻線に印加するf軸電圧vfを第2のf軸電圧指令値v* fと一致するように調整する。
 電流センサ106は、インバータ103からモータ101に供給される三相交流電流のうち、少なくとも2相の電流、例えば、u相電流iu、v相電流ivを検出する。検出された2相の電流iu、ivは、A/D(アナログ/デジタル)変換器107でデジタル信号(電流ius、ivs)に変換され、座標変換器108に入力される。なお、電流センサ106が2相の電流のみを検出する場合、残りの1相の電流iwsは、次式(1)により求めることができる。
Figure JPOXMLDOC01-appb-M000001
 また、電流センサ106は、f軸電流出力部105からモータ101に供給される回転子電流(f軸電流if)を検出する。検出されたf軸電流ifは、A/D変換器107でデジタル信号に変換され、f軸電流制御部116に出力される。
 磁極位置検出器109は、モータ101の固定子の電気角度を取得するために、当該固定子の位置(角度)に応じたA相B相Z相のパルス(ABZパルス)をパルスカウンタ110に出力する。
 パルスカウンタ110は、ABZパルスに基づいてモータ101の電気角度θreを算出して、角速度演算器111に出力する。
 角速度演算器111は、入力される電気角度θreの時間変化率から、電気角速度ωreと、機械角速度ωrmとを算出する。機械角速度ωrmは、電気角速度ωreをモータ極対数pで除算することにより求められる。機械角速度ωrmは、電流指令値演算部113に出力される。電気角速度ωreは、非干渉制御部115と先読み補償部112とに出力される。
 先読み補償部112は、電気角度θreと電気角速度ωreとを入力して、電気角速度ωreと制御系が持つ無駄時間との乗算値を電気角度θreに加算することにより、先読み補償後電気角θre'を算出する。先読み補償後電気角θre'は、座標変換器117に出力される。
 座標変換器108は、3相交流座標系(uvw軸)から直交2軸直流座標系(d-q軸)への変換を行う。具体的には、座標変換器108は、入力されるu相電流ius、v相電流ivs、及び電気角度θreと、上記式(1)で求まるw相電流iwsとから、以下式(2)を用いて座標変換処理を行うことによって、d軸電流idとq軸電流iqを算出する。
Figure JPOXMLDOC01-appb-M000002
 電流指令値演算部113は、トルク指令値とモータ101の動作状況とに基づいて、固定子電流に対するd軸電流指令値i* d及びq軸電流指令値i* qと、回転子電流(界磁電流)に対するf軸電流指令値i* fとを算出する。より具体的には、電流指令値演算部113は、トルク指令値T*、モータ回転数(機械角速度ωrm)、及び、電源電圧(直流電圧)Vdcを入力とし、d軸電流指令値i* d、q軸電流指令値i* q、及びf軸電流指令値i* fを算出する。d軸電流指令値i* d、q軸電流指令値i* q、及びf軸電流指令値i* fは、予め記憶された、トルク指令値T*、モータ回転数(機械角速度ωrm)、及び、電源電圧Vdcと、d軸電流指令値i* d、q軸電流指令値i* q、及びf軸電流指令値i* fの各々との関係を定めたマップデータを参照することにより求められる。なお、トルク指令値T*は、モータ101に所望のトルクを発生させるための指令値であって、ドライバのアクセル操作等に応じて公知の方法により設定される。
 電流制御部114は、実際の電流(実電流)の計測値であるd軸電流idをd軸電流指令値i* dに定常偏差なく所望の応答性で追従させるためのd軸電圧指令値(第1のd軸電圧指令値)を算出して、電圧指令値演算部116に出力する。
 また、電流制御部114は、実際の電流(実電流)の計測値であるq軸電流iqをq軸電流指令値i* qに定常偏差なく所望の応答性で追従させるためのq軸電圧指令値(第1のq軸電圧指令値)を算出して、電圧指令値演算部116に出力する。
 また、電流制御部114は、実際の電流(実電流)の計測値であるf軸電流ifをf軸電流指令値i* fに定常偏差なく所望の応答性で追従させるためのf軸電圧指令値(第1のf軸電流指令値)を算出して、電圧指令値演算部116に出力する。
 電流制御部114の詳細は図2を参照して後述する。
 非干渉制御部115は、入力される電気角速度ωreと、d軸電流idと、q軸電流iqと、f軸電流ifとから、d軸、q軸、及びf軸間の干渉電圧を相殺するために必要な非干渉電圧vd_dcpl、vq_dcpl、vf_dcplを算出する。ただし、非干渉制御部115は、実際の電流(実電流)の計測値であるd軸電流id、q軸電流iq、及びf軸電流ifの代わりに、d軸電流指令値i* d、q軸電流指令値i* q、及びf軸電流指令値i* fに基づいて算出される後述のd軸電流規範応答id_ref、q軸電流規範応答iq_ref、f軸電流規範応答if_ref(実電流)を用いて非干渉電圧vd_dcpl、vq_dcpl、vf_dcplを算出してもよい。
 電圧指令値演算部116は、電流制御部114の出力である第1のd軸電圧指令値vd_dsh、第1のq軸電圧指令値vq_dsh、及び、第1のf軸電圧指令値vf_dshを、非干渉制御部115の出力である非干渉電圧vd_dcpl、vq_dcpl、vf_dcplを用いて補正(本実施形態では加算)する。そして、電圧指令値演算部116は、当該補正により得た、第2のd軸電圧指令値v* d、及び、第2のq軸電圧指令値v* qを座標変換器117に出力するとともに、第2のf軸電圧指令値v* fをf軸電流出力部105に出力する。
 座標変換器117は、電気角速度ωreで回転する直交2軸直流座標系(d‐q軸)から3相交流座標系(uvw相)への変換を行う。具体的には、座標変換器117は、入力される第2のd軸電圧指令値v* d、第2のq軸電圧指令値v* q、及び、先読み補償後電気角θre'から、以下式(3)を用いて座標変換処理を行うことによって、uvw各相の電圧指令値v* u、v* v、v* wを算出する。
Figure JPOXMLDOC01-appb-M000003
 以上が本実施形態のモータ制御システム100が備える構成の概要である。なお、本実施形態では、上述した構成のうち、PWM変換器102、A/D変換器107、座標変換器108、パルスカウンタ110、角速度演算器111、先読み補償部112、電流指令値演算部113、電流制御部114、非干渉制御部115、電圧指令値演算部116、座標変換器117は、少なくとも一つ以上のコントローラで構成されるコントローラ10が備える一機能部として構成される。コントローラ10は、例えば、中央演算装置(CPU)、読み出し専用メモリ(ROM)、ランダムアクセスメモリ(RAM)、及び、入出力インタフェース(I/Oインタフェース)から構成される。
 [電流制御部]
 以下では、電流制御部114の詳細について説明する。まず、本実施形態の制御対象である巻線界磁型同期モータ101の電圧方程式について説明する。巻線界磁型同期モータ101の電圧方程式は、以下式(4)で表される。
Figure JPOXMLDOC01-appb-M000004
 ただし、上記式(4)の各パラメータは、以下のとおりである。なお、式中のsはラプラス演算子である。
d    :     d軸電流
q    :     q軸電流
f    :     f軸電流
d    :     d軸電圧
q    :     q軸電圧
f    :     f軸電圧
d    :     d軸インダクタンス
q    :     q軸インダクタンス
f    :     f軸インダクタンス
M     :     固定子/回転子間の相互インダクタンス
d'   :     d軸動的インダクタンス
q'   :     q軸動的インダクタンス
f'   :     f軸動的インダクタンス
M'    :     固定子/回転子間の動的相互インダクタンス
a    :     固定子巻線抵抗
f    :     回転子巻線抵抗
ωre   :     電気角速度
 ここで、非干渉制御部115による非干渉制御が理想的に機能すれば、上記式(4)の電圧方程式は、以下式(5)に示すように対角化することができる。
Figure JPOXMLDOC01-appb-M000005
 上記式(5)によれば、d軸、q軸、及びf軸の電圧から電流までの特性は、それぞれ以下式(6)、(7)、及び(8)に示すとおりの一次遅れとなる。
Figure JPOXMLDOC01-appb-M000006
Figure JPOXMLDOC01-appb-M000007
Figure JPOXMLDOC01-appb-M000008
 続いて、図2を参照して、電流制御部114の構成の詳細を説明する。
 図2は、本実施形態の電流制御部114を示すブロック図である。電流制御部114は、f軸疑似非干渉電圧算出部201と、f軸不足電圧算出部202と、d軸電流指令値補正部203と、d軸電流制御部204と、q軸電流制御部205と、f軸電流制御部206と、を含んで構成される。以下、各構成の詳細について図2から図11を参照して説明する。
 f軸疑似非干渉電圧算出部201は、電流指令値演算部113から出力されるd軸電流指令値i* dを入力として、f軸疑似非干渉電圧vf_dcpl_dshを算出する。詳細は図3を用いて説明する。
 図3は、f軸疑似非干渉電圧算出部201を示すブロック図である。f軸疑似非干渉電圧算出部201は、d軸モデル301と、d軸電流フィードバック(F/B)モデル302と、f軸疑似非干渉電圧モデル303とを含む。
 d軸モデル301は、後述するd軸電流F/Bモデル302の出力値を入力として、d軸電圧vdからd軸電流idまでの特性をモデル化したd軸モデルを用いてフィルタリング処理することにより、d軸電流規範応答id_ref及びd軸電流規範応答の微分値s・id_refを出力する。d軸モデル301の詳細は図4を用いて説明する。
 図4は、d軸モデル301を示すブロック図である。d軸モデル301は、乗算器401と、減算器402と、除算器403と、積分器404とを含む。
 乗算器401は、積分器404の出力値であるd軸電流規範応答id_refに、モータ101が備える固定子巻線の抵抗値である固定子巻線抵抗Raを乗算して、算出値を減算器402に出力する。
 減算器402は、d軸電流F/Bモデル302の出力値から乗算器401の出力値を減算して、算出値を除算器403に出力する。
 除算器403は、減算器402の出力値に、d軸動的インダクタンスLd'を除算してd軸電流規範応答の微分値s・id_refを算出する。d軸電流規範応答の微分値s・id_refは、積分器404と、後述のf軸疑似非干渉電圧モデル303とに出力される。
 積分器404は、d軸電流規範応答の微分値s・id_refを積分処理することによりd軸電流規範応答id_refを算出する。d軸電流規範応答id_refは、乗算器401と、d軸電流F/Bモデル302とに出力される。このように、積分器404の出力を電流規範応答とするとともに、積分器404の入力を電流規範応答の微分値として取得可能に構成することにより、電流規範応答の微分値を容易に求めることができる。
 d軸電流F/Bモデル302(図3参照)は、d軸電流指令値i* dと、上述のd軸モデル301の出力であるd軸電流規範応答id_refを入力として、d軸電流指令値i* dにd軸電流規範応答id_refを定常偏差なく所望の応答性で追従させるためのd軸電流F/Bモデル出力値(d軸電圧推定値)を算出して、d軸モデル301に出力する。詳細は、図5を参照して説明する。
 図5は、d軸電流F/Bモデル302の詳細を示すブロック図である。d軸電流F/Bモデル302は、乗算器501、502と、減算器503と、を含んで構成される。
 乗算器501は、電流指令値演算部113から出力されるd軸電流指令値i* dにゲインGadを乗算し、算出した値を減算器503に出力する。ゲインGadは、下記式(9)で表される。
 乗算器502は、d軸モデル301から出力されるd軸電流規範応答id_refにゲインGbdを乗算し、算出した値を減算器503に出力する。ゲインGbdは、下記式(9)で表される。
Figure JPOXMLDOC01-appb-M000009
 ただし、式(9)中のτmは、d軸及びq軸の電流制御規範応答時定数、すなわち、実現したい所望の応答特性に相当する時定数である。
 そして、減算器503は、乗算器501の出力値から乗算器502の出力値を減算することによりd軸電流F/Bモデル出力値を算出して、上述のd軸モデル301に出力する。
 f軸疑似非干渉電圧モデル303(図3参照)は、d軸モデル301の出力であるd軸電流規範応答の微分値s・id_refを入力して、f軸疑似非干渉電圧vf_dcpl_dshを出力する。f軸疑似非干渉電圧モデル303の詳細は図6を用いて説明する。
 図6は、f軸疑似非干渉電圧モデル303を示すブロック図である。f軸疑似非干渉電圧モデル303は、乗算器601で構成される。
 乗算器601では、d軸モデル301の出力であるd軸電流規範応答の微分値s・id_refに、固定子/回転子間の動的相互インダクタンスM'を乗算することによりf軸疑似非干渉電圧vf_dcpl_dshを算出する。f軸疑似非干渉電圧vf_dcpl_dshは、f軸不足電圧算出部202(図2参照)に出力される。
 f軸不足電圧算出部202は、f軸疑似非干渉電圧vf_dcpl_dshと、電源電圧Vdcと、f軸電流指令値i* fとを入力とし、f軸不足電圧Δvfを算出する。f軸不足電圧算出部202の詳細は図7を用いて説明する。
 図7は、f軸不足電圧算出部202を示すブロック図である。f軸不足電圧算出部202は、f軸モデル701と、f軸電流F/Bモデル702と、f軸不足電圧算出モデル703とを含む。
 f軸モデル701は、後述する、f軸電流F/Bモデル702から出力されるf軸電流F/Bモデル出力値を入力として、f軸電圧からf軸電流までの特性をモデル化したモデル(モータ101のf軸モデル)によるフィルタリング処理を行い、f軸電流規範応答if_refと、f軸電流規範応答の微分値s・if_refとを出力する。詳細は図8を参照して説明する。
 図8は、f軸モデル701の詳細を示すブロック図である。f軸モデル701は、乗算器801と、減算器802と、積分器804とを含んで構成される。
 乗算器801は、積分器804の出力であるf軸電流規範応答if_refに固定子巻線抵抗Rfを乗算して、得た値を減算器802に出力する。
 減算器802は、f軸電流F/Bモデル702から出力されるf軸電流F/Bモデル出力値から乗算器801の出力値を減算して、算出した値を除算器803に出力する。
 除算器803は、減算器802の出力値をf軸動的インダクタンスLf'で除算してf軸電流規範応答の微分値s・if_refを算出する。f軸電流規範応答の微分値s・if_refは、積分器804に出力される。
 積分器804は、f軸電流規範応答の微分値s・if_refを積分処理することによりf軸電流規範応答if_refを算出する。f軸電流規範応答if_refは、乗算器801と、f軸電流F/Bモデル702とに出力される。
 f軸電流F/Bモデル702(図7参照)は、f軸電流指令値i* fと、上述のf軸モデル701の出力であるf軸電流規範応答if_refを入力として、f軸電流指令値i* fにf軸電流規範応答if_refを定常偏差なく所望の応答性で追従させるためのf軸電流F/Bモデル出力値(f軸電圧推定値)を算出する。f軸電流F/Bモデル702の詳細は、図9を参照して説明する。
 図9は、f軸電流F/Bモデル702を示すブロック図である。f軸電流F/Bモデル702は、乗算器901、902と、減算器903と、を含んで構成される。
 乗算器901は、電流指令値演算部113から出力されるf軸電流指令値i* fにゲインGafを乗算し、算出した値を減算器903に出力する。ゲインGafは、下記式(10)で表される。
 乗算器902は、f軸モデル701から出力されるf軸電流規範応答i* f_refにゲインGbfを乗算し、算出した値を減算器903に出力する。ゲインGbfは、下記式(10)で表される。
Figure JPOXMLDOC01-appb-M000010
 ただし、式(10)中のτmは、f軸の電流制御規範応答時定数、すなわち、実現したい所望の応答特性に相当する時定数である。
 減算器903は、乗算器901の出力値から乗算器902の出力値を減算することによりf軸電流F/Bモデル出力値を算出して、f軸モデル701と、f軸不足電圧算出モデル703とに出力する。なお、このf軸電流F/Bモデル出力値は、f軸電圧vfの推定値(f軸電圧推定値)である。
 f軸不足電圧算出モデル703(図7参照)は、f軸電流F/Bモデル出力値と、電源電圧Vdcと、f軸疑似非干渉電圧vf_dcpl_dshとを入力として、f軸電流F/Bモデル出力値を所定の値(リミット値)でリミット処理することにより得たf軸不足電圧Δvfを出力する。より具体的には、f軸不足電圧算出モデル703は、f軸電流指令値i* fに対してf軸電流ifを所望の応答速度で制御するために必要なf軸電圧vfに相当するf軸電流F/Bモデル出力値に対する出力可能電圧の不足分をf軸不足電圧Δvfとして算出する。なお、ここでの出力可能電圧は、後述するように電源電圧Vdcに応じて定まる電圧であって、f軸電流出力部105が出力可能な電圧である。f軸不足電圧算出モデル703の詳細は、図10を参照して説明する。
 図10は、f軸不足電圧算出モデル703を示すブロック図である。f軸不足電圧算出モデル703は、比較器1001と、反転器1002と、比較器1003と、減算器1004、1005、1006とを含んで構成される。
 比較器1001の前段に設けられる減算器1005においては、電源電圧104の電源電圧Vdcからf軸疑似非干渉電圧vf_dcpl_dshを差し引いた減算値が求められる。そして、比較器1001は、f軸電流疑似F/Bモデル702からの出力値と、減算器1005の出力値とを比較し、より小さな値を比較器1003へと出力する。
 反転器1002は、電源電圧Vdcの符号を反転させる。
 比較器1003の前段には減算器1006が設けられており、減算器1006においては、反転器1002の出力から、f軸疑似非干渉電圧vf_dcpl_dshを差し引いた減算値が求められる。そして、比較器1003は、比較器1001の出力値と、減算器1006の出力値とを比較し、より大きな値を減算器1004に出力する。
 このような構成により、f軸不足電圧算出モデル703は、f軸電流F/Bモデル702(減算器903)の出力値に対して所定の値(リミット値)に基づくリミット処理を実行することができる。本実施形態における所定の値(上下限値)は以下の通りである。すなわち、上限値には、電源電圧Vdcからf軸疑似非干渉電圧vf_dcpl_dshを減算した値が設定される(上限値=Vdc-vf_dcpl_dsh)。下限値には、-電源電圧Vdcからf軸非干渉電圧vf_dcpl_dshを減算した値が設定される(下限値=-Vdc-vf_dcpl_dsh)。なお、これら上限値と下限値との間の電圧が上述の出力可能電圧に相当する。
 上下限値がこのように設定されることにより、後段の電圧指令値演算部116において第1のf軸電圧指令値に非干渉電圧vf_dcplが加算されても、算出された値が電源電圧Vdcに対して飽和してしまうことを回避することができる。また、上述のようなリミット処理を実行することにより、f軸電圧の飽和の有無を判定する手段や、f軸電圧の飽和の有無に応じてフィルタを切替えるフィルタ切替え手段等を要さずに、適切なf軸電圧指令値を算出することができる。
 なお、f軸疑似非干渉電圧vf_dcpl_dshは、f軸疑似非干渉電圧算出部201において、後述するd軸電流指令値補正部203で補正される前のd軸電流指令値i* dに基づいて算出される。すなわち、本実施形態における所定の上下限値(リミット値)は、電源電圧Vdcと、補正前のd軸電流指令値i* dとに基づいて算出されたf軸疑似非干渉電圧vf_dcpl_dshとに基づいて算出される。
 そして、減算器1004は、f軸電流F/Bモデル出力値から比較器1003の出力値を差し引くことによりf軸不足電圧Δvfを算出する。すなわち、f軸不足電圧Δvfは、f軸電流F/Bモデル出力値と、f軸電流F/Bモデル出力値を所定の上下限値でリミット処理して得た値との差分に応じて決定される。これにより、f軸電流を所望の応答で制御するためのf軸電圧vfに対する出力可能電圧の不足分を算出することができる。f軸不足電圧Δvfは、d軸電流指令値補正部203に出力される。なお、このような構成によりf軸不足電圧Δvfを算出し、f軸電流を所望の応答性で制御することで、f軸電流F/Bモデル702のゲイン(901、902)を、安定性を考慮して小さく設定する必要がないので、これらゲイン(901、902)を所望の応答性を満足するゲインに設定することができる(数式(10)参照)。
 d軸電流指令値補正部203(図2参照)は、d軸電流指令値i* dと、f軸不足電圧Δvfを入力とし、f軸不足電圧Δvfに基づいてd軸電流指令値i* dを補正する。補正されたd軸電流指令値i* dは、補正後d軸電流指令値としてd軸電流制御部204に出力される。d軸電流指令値補正部203の詳細について図11を用いて説明する。
 図11は、d軸電流指令値補正部203を示すブロック図である。d軸電流指令値補正部203は、フィルタ1101と、減算器1102とを含んで構成される。
 フィルタ1101は、Hd_comp(s)なる関数で構成される。フィルタ1101は、入力されるf軸不足電圧Δvfに対して、Hd_comp(s)を用いたフィルタリング処理を行い、d軸電流補正値id_compを減算器1102に出力する。Hd_comp(s)の詳細については後述する。
 そして、減算器1102は、d軸電流指令値i* dからd軸電流補正値id_compを減算して、補正後d軸電流指令値を算出する。これにより、d軸電流指令値i* dは、Hd_comp(s)を用いて補正される。
 以下、Hd_comp(s)の導出過程について説明する。まず、上記式(4)で示した電圧方程式より、f軸電圧vfは次式(11)で表される。
Figure JPOXMLDOC01-appb-M000011
 式(11)中の電流(d軸電流id、f軸電流if)を規範応答で置換するとf軸電圧vfは次式(12)で表される。
Figure JPOXMLDOC01-appb-M000012
 ここで、図11を用いて説明したように、式(12)にd軸電流補正値i* dを補正するd軸電流補正値id_compを導入すると、f軸電圧vfは次式(13)で表される。
Figure JPOXMLDOC01-appb-M000013
 そして、(13)式の右辺第2項においてf軸不足電圧Δvfを相殺することで、f軸電圧飽和を抑制することができる。この時のd軸電流補正値id_compは次式(14)となる。
Figure JPOXMLDOC01-appb-M000014
 d軸電流補正値id_compがこのように構成されることにより、f軸不足電圧Δvfに応じて、すなわち、f軸電圧vfに対する出力可能電圧の不足分に応じてd軸電流指令値i* dを補正することができる。その結果、f軸電圧飽和が抑制されるので、f軸電流ifの応答特性を改善し、f軸電流ifを所望の応答性(応答速度)で動作させることができ、トルクの応答性を向上させることができる。また、f軸不足電圧Δvfが正の場合は、d軸電流指令値を負側(減少方向)に補正し、f軸不足電圧Δvfが負の場合は、d軸電流指令値を正側(増加方向)に補正するように構成されるので、トルク指令値の増加方向、減少方向に関わらず、f軸電流ifを所望の応答性(応答速度)で動作させることができる。
 なお、実際には、f軸電圧が飽和したときのみd軸電流指令値i* dが補正されることが好ましい。しかしながら、上記式(14)では補正項が積分特性(1/s)を持つため、f軸電圧飽和時以外においてもd軸電流指令値i* dを補正してしまう。そこで、次式(15)で示すように、上記式(14)にハイパスフィルタを追加することで積分特性を相殺する。これにより、Hd_comp(s)をf軸電圧が飽和したときのみd軸電流指令値i* dを補正するように構成することができる。
Figure JPOXMLDOC01-appb-M000015
 ただし、式(15)中のハイパスフィルタの時定数τhrは、チューニング要素であり、所望の応答性を得るために適宜調整される。
 以上がHd_comp(s)の詳細である。このような関数を用いて算出された補正後d軸電流指令値は、d軸電流制御部204(図2参照)に出力される。
 d軸電流制御部204は、実際の電流(実電流)の計測値であるd軸電流idをd軸電流指令値i* dに定常偏差なく所望の応答性で追従させるための第1のd軸電圧指令値vd_dshを算出して、電圧指令値演算部116に出力する。
 q軸電流制御部205は、実際の電流(実電流)の計測値であるq軸電流iqをq軸電流指令値i* qに定常偏差なく所望の応答性で追従させるための第1のq軸電圧指令値vq_dshを算出して、電圧指令値演算部116に出力する。
 f軸電流制御部206は、実際の電流(実電流)の計測値であるf軸電流ifをf軸電流指令値i* fに定常偏差なく所望の応答性で追従させるための第1のf軸電圧指令値vf_dshを算出して、電圧指令値演算部116に出力する。
 上記のd軸電流制御部204、q軸電流制御部205、及び、f軸電流制御部206は、d軸、q軸、及びf軸間の干渉電圧を相殺する制御(非干渉制御)が理想的に機能すれば、通常、1入力1出力の単純な特性(制御対象特性)となる(上記式(6)、(7)、及び(8)を参照)。従って、本実施形態のd軸電流制御部204、q軸電流制御部205、及び、f軸電流制御部206は、簡単なPIフィードバック補償器、又は、いわゆるロバストモデルマッチング補償器のような公知の補償器により実現することができる。
 以下では、これまで説明した一実施形態のモータ制御方法を適用した一制御周期の流れを、図12を参照して説明する。
 図12は、一実施形態のモータ制御方法を示すフローチャートである。図12で示す開始から終了までにかかる一制御周期は、モータ制御システム100が起動している間、一定の間隔で常時実行するように上記のコントローラ10にプログラムされている。
 ステップS11では、コントローラ10(パルスカウンタ110、A/D変換器107)は、ABZパルスに基づいて電気角度θreを算出するとともに、取得したu、v相の電流iu、iv、及び、回転子巻線を流れるf軸電流ifから、デジタル信号としての電流ius、ivs、ifsを検出する。
 ステップS12では、コントローラ10(角速度演算器111)は、入力される電気角度θreの時間変化率から、電気角速度ωreと機械角速度ωrmとを算出する。
 ステップS13では、コントローラ10(先読み補償部112)は、電気角度θreと電気角速度ωreとから、制御系が持つ無駄時間が考慮された先読み補償後電気角θre'を算出する。
 ステップS14では、コントローラ10(座標変換器108)は、u、v、w各相の電流ius、ivs、iws、及び電気角度θreから、d、q軸電流id、iqを算出する。
 ステップS15では、コントローラ10(電流指令値演算部113)は、トルク指令値T*、モータ回転数(機械角速度ωrm)、及び、電源電圧Vdcから、d軸電流指令値i* d、q軸電流指令値i* q、及びf軸電流指令値i* fを算出する。
 ステップS16では、コントローラ10(電流制御部114)は、d、q、f軸電流id、iq、ifをd、q、f軸電流指令値i* d、i* q、i* fに定常偏差なく所望の応答性で追従させるためのd、q、f軸各相の第1の電圧指令値vd_dsh、vq_dsh、vf_dshを算出する。
 ステップS17では、コントローラ10(非干渉制御部115)は、d、q、f軸電流id、iq、ifと、電気角速度ωreとから、d軸、q軸、及びf軸間の干渉電圧を相殺するために必要な非干渉電圧vd_dcpl、vq_dcpl、vf_dcplを算出する。
 ステップS18では、コントローラ10(電圧指令値演算部116)は、d、q、f各軸の第1の電圧指令値vd_dsh、vq_dsh、vf_dshを、非干渉電圧vd_dcpl、vq_dcpl、vf_dcplを用いて補正することにより、d、q、f各軸の第2の電圧指令値v* d、v* q、v* fを算出する。
 そして、ステップS19では、コントローラ10(座標変換器117)は、第2のd軸電圧指令値v* d、第2のq軸電圧指令値v* q、及び、先読み補償後電気角θre'から、uvw各相の電圧指令値v* u、v* v、v* wを算出する。
 以上が本実施形態のモータ制御方法による制御フローの概要である。以下では、図13及び図14を参照して、本実施形態のモータ制御方法による作用効果について説明する。
 図13は、本実施形態のモータ制御方法による制御結果を示すタイムチャートである。横軸は時間を表し、縦軸は、左側の上から順に、d軸電流[A]、q軸電流[A]、及び、f軸電流[A]を表し、右側の上から順に、f軸電圧[V]、トルク[N]を表している。d軸電流[A]、q軸電流[A]、及び、f軸電流[A]に係る図中の破線は電流指令値を示している。トルク[N]に係る図中の破線はトルク指令値を示している。そして、この電流指令値又はトルク指令値に対する実電流又は実トルクの応答を、本実施形態は実線で、従来技術による制御(従来例)は一点鎖線で示している。
 図13で示されるのは、トルク指令値が減少方向にある場合における本実施形態と従来例との比較結果である。
 時刻t0では、トルク指令値T*がステップで減少方向に印加されることに伴い、d軸電流指令値i* d、q軸電流指令値i* q、及びf軸電流指令値i* fがステップで印加される。
 時刻t1では、d軸電流id、q軸電流iq、f軸電流ifが、印加されたd軸電流指令値i* d、q軸電流指令値i* q、f軸電流指令値i* fを目標に収束していく。ただし、従来例では、f軸電圧vfが出力可能下限電圧(-Vdc)で飽和するため、f軸電流ifの応答(応答速度)が低下してしまう。その結果、トルクの応答も低下してしまう。
 一方、本実施形態のモータ制御方法によれば、f軸電流ifを所望の応答で制御するための電圧(f軸F/Bモデル出力値)と、電源電圧Vdcに基づく出力可能電圧との差分であるf軸不足電圧Δvfに応じて、d軸電流idを過渡的に補正するので、f軸電圧vfの飽和が抑制され、f軸電流ifの応答を早くすることができる。その結果、図で示すとおり、トルクの応答を従来例に比べて早くすることができている。
 続いて、図14で示されるのは、トルク指令値が増加傾向にある場合における本実施形態のモータ制御と、f軸不足電圧Δvfに応じてd軸電流idを補正する機能をOFFにした場合の実施形態(比較例)のモータ制御との比較結果である。
 図14も図13と同様に横軸は時間を表し、縦軸は、左側の上から順に、d軸電流[A]、q軸電流[A]、及び、f軸電流[A]を表し、右側の上から順に、f軸電圧[V]、トルク[N]を表している。d軸電流[A]、q軸電流[A]、及び、f軸電流[A]に係る図中の破線は電流指令値を示している。トルク[N]に係る図中の破線はトルク指令値を示している。そして、この電流指令値又はトルク指令値に対する実電流又は実トルクの応答を、本実施形態は実線で、比較例による制御は一点鎖線で示している。
 時刻t0では、トルク指令値T*がステップで増加方向に印加されることに伴い、d軸電流指令値i* d、q軸電流指令値i* q、及びf軸電流指令値i* fがステップで印加される。
 時刻t1では、d軸電流id、q軸電流iq、f軸電流ifが、印加されたd軸電流指令値i* d、q軸電流指令値i* q、f軸電流指令値i* fを目標に収束していく。ただし、比較例では、f軸電圧vfが出力可能上限電圧(Vdc)で飽和するため、f軸電流ifの応答(応答速度)が低下してしまう。その結果、トルクの応答も低下してしまう。
 一方、本実施形態のモータ制御方法によれば、f軸電流ifを所望の応答で制御するための電圧(f軸F/Bモデル出力値)と、電源電圧Vdcに基づく出力可能電圧との差分であるf軸不足電圧Δvfに応じて、d軸電流idを過渡的に補正するので、f軸電圧vfの飽和が抑制され、f軸電流ifの応答を早くすることができる。その結果、図で示すとおり、トルクの応答を従来例に比べて早くすることができている。
 このように、本実施形態では、f軸電圧の不足分に応じてd軸電流指令値を補正するモータ制御方法を適用することにより、f軸電圧飽和が抑制されえ、f軸電流を所望の応答性で動作させることができる。これにより、トルクの応答性を向上させることができる。
 また、本実施形態によれば、f軸電圧の不足分に応じてd軸電流指令値を連続的に補正するため、トルク指令値の増加/減少方向に応じたスイッチ等による非連続なパラメータ或いは仕様の切り替え等を要さずに、トルク指令値の増加/減少方向に関わらずにトルクの応答性を向上させることができる。
 以上、一実施形態のモータ制御方法は、回転子巻線を有する回転子と、固定子巻線を有する固定子とを備える巻線界磁型同期モータ101において、固定子巻線に流れる固定子電流と回転子巻線に流れる回転子電流とを制御するモータ制御方法である。このモータ制御方法によれば、巻線界磁型同期モータ101に所望のトルクを発生させるトルク指令値T*を設定し、トルク指令値T*と巻線界磁型同期モータ101の動作状況とに基づいて、固定子電流に対するd軸電流指令値i* d及びq軸電流指令値i* qと、回転子電流に対するf軸電流指令値i* fとを算出する。そして、回転子巻線に印加される電圧であって、f軸電流指令値に対して回転子電流を所定の応答速度で制御するのに必要なf軸電圧vf(f軸電流F/Bモデル出力値(f軸電圧推定値))を推定(算出)し、推定されたf軸電圧に対して電源電圧Vdcに基づいて定まる出力可能電圧が不足した場合には、当該不足分に応じてd軸電流指令値i* dを補正し、d軸電流指令値i* d、q軸電流指令値i* q、及びf軸電流指令値i* fに基づいて巻線界磁型同期モータ101を制御する。
 これにより、f軸電圧vfに対する出力可能電圧の不足分に応じてd軸電流指令値i* dが補正されるので、f軸電圧飽和が抑制され、f軸電流ifを所望の応答性で動作させることができる。その結果、トルク指令値の増加傾向、減少方向に関わらず、トルクの応答性を向上させることができる。
 また、一実施形態のモータ制御方法によれば、f軸電圧vfが出力可能電圧の上限値よりも大きい場合には、d軸電流指令値i* dを負側に補正し、f軸電圧vfが出力可能電圧の下限値よりも小さい場合には、d軸電流指令値i* dを正側に補正する。これにより、トルク指令値の増加傾向、減少方向に応じて、d軸電流指令値i* dを適切に補正して、トルク指令値の増加傾向、減少方向に関わらず、f軸電流ifを所望の応答性で動作させることができる。
 また、一実施形態のモータ制御方法によれば、f軸電圧vfから回転子電流を構成するf軸電流ifまでの特性をモデル化したf軸モデル701と、f軸電流指令値i* fとf軸モデル701の出力とが入力されるf軸電流フィードバック(F/B)モデル702と、により疑似的なF/B系を構成し、f軸電流F/Bモデルの出力値(f軸電流F/Bモデル出力値(f軸電圧推定値))を所定の上下限値でリミット処理することにより得た値が、f軸電流F/Bモデルの出力値に対して不足する場合には、当該不足分(Δvf)に応じてd軸電流指令値i* dを補正する。f軸電圧vfに対する出力可能電圧の不足分に応じたd軸電流指令値i* dの補正をこのような構成により実現することにより、安定性を考慮してf軸電流F/Bモデル702のゲイン(901、902)を小さく設定する必要がないので、これらゲイン(901、902)を所望の応答性を満足するゲインに設定することができる(数式(10)参照)。
 また、一実施形態のモータ制御方法によれば、固定子電流のd軸及び回転子電流のf軸間の干渉電圧を非干渉化するためのf軸疑似非干渉電圧vf_dcpl_dshを補正がなされる前のd軸電流指令値i* dに基づいて算出し、上記の所定の上下限値は、電源電圧Vdcとf軸疑似非干渉電圧vf_dcpl_dshとに基づいて算出される。このように補正前のd軸電流指令値i* dに基づいて上下限値を設定することにより、f軸電圧vfに対する電源電圧Vdcの不足分を適切に算出することができる。その結果、後段の電圧指令値演算部116において第1のf軸電圧指令値vf_dshに非干渉電圧vf_dcplが加算されても、算出された値が電源電圧Vdcに対して飽和してしまうことを回避して、f軸電流ifを所望の応答性で制御することができる。
 以上、本発明の実施形態について説明したが、上記実施形態は本発明の適用例の一部を示したに過ぎず、本発明の技術的範囲を上記実施形態の具体的構成に限定する趣旨ではない。本発明は、上述した実施形態に限定されることはなく、様々な変形や応用が可能である。例えば、図1で示す構成の一部はコントローラ10の一機能部として構成される旨上述したが、同様の機能を発揮可能な限り、コントローラ10の一機能部として構成される必要は必ずしもない。例えば、パルスカウンタ110は、コントローラ10の一機能部として構成される必要はなく、パルスカウント機能を備えるデバイスとしてコントローラ10とは別個に設けられても良い。

Claims (5)

  1.  回転子巻線を有する回転子と、固定子巻線を有する固定子とを備える巻線界磁型同期モータにおいて、前記固定子巻線に流れる固定子電流と前記回転子巻線に流れる回転子電流とを制御するモータ制御方法であって、
     前記巻線界磁型同期モータに所望のトルクを発生させるトルク指令値を設定し、
     前記トルク指令値と前記巻線界磁型同期モータの動作状況とに基づいて、前記固定子電流に対するd軸電流指令値及びq軸電流指令値と、前記回転子電流に対するf軸電流指令値とを算出し、
     前記回転子巻線に印加される電圧であって、前記f軸電流指令値に対して前記回転子電流を所定の応答速度で制御するのに必要なf軸電圧を推定し、
     推定された前記f軸電圧に対して電源電圧に基づいて定まる出力可能電圧が不足した場合には、当該不足分に応じて前記d軸電流指令値を補正し、
     前記d軸電流指令値、前記q軸電流指令値、及び前記f軸電流指令値に基づいて前記巻線界磁型同期モータを制御する、
    モータ制御方法。
  2.  請求項1に記載のモータ制御方法において、
     前記f軸電圧が前記出力可能電圧の上限値よりも大きい場合には、前記d軸電流指令値を負側に補正し、
     前記f軸電圧が前記出力可能電圧の下限値よりも小さい場合には、前記d軸電流指令値を正側に補正する、
    モータ制御方法。
  3.  請求項1又は2に記載のモータ制御方法において、
     前記f軸電圧から前記回転子電流を構成するf軸電流までの特性をモデル化したf軸モデルと、
     前記f軸電流指令値と前記f軸モデルの出力とが入力されるf軸電流フィードバック(F/B)モデルと、により疑似的なF/B系を構成し、
     前記f軸電流F/Bモデルの出力値を所定の上下限値でリミット処理することにより得た値が前記f軸電流F/Bモデルの出力値に対して不足する場合には、当該不足分に応じて前記d軸電流指令値を補正する、
    モータ制御方法。
  4.  請求項3に記載のモータ制御方法において、
     前記固定子電流のd軸及び前記回転子電流のf軸間の干渉電圧を非干渉化するためのf軸疑似非干渉電圧を前記補正がなされる前のd軸電流指令値に基づいて算出し、
     前記所定の上下限値は、前記電源電圧と前記f軸疑似非干渉電圧とに基づいて算出される、
    モータ制御方法。
  5.  回転子巻線を有する回転子と、固定子巻線を有する固定子とを備える巻線界磁型同期モータと、前記固定子巻線に流れる固定子電流と前記回転子巻線に流れる回転子電流とを制御するコントローラとを備えるモータ制御装置であって、
     前記コントローラは、
     前記巻線界磁型同期モータに所望のトルクを発生させるトルク指令値を設定し、
     前記トルク指令値と前記巻線界磁型同期モータの動作状況とに基づいて、前記固定子電流に対するd軸電流指令値及びq軸電流指令値と、前記回転子電流に対するf軸電流指令値とを算出し、
     前記回転子巻線に印加される電圧であって、前記f軸電流指令値に対して前記回転子電流を所定の応答速度で制御するのに必要なf軸電圧を推定し、
     推定された前記f軸電圧に対して電源電圧に基づいて定まる出力可能電圧が不足した場合には、当該不足分に応じて前記d軸電流指令値を補正し、
     前記d軸電流指令値、前記q軸電流指令値、及び前記f軸電流指令値に基づいて前記巻線界磁型同期モータを制御する、
    モータ制御装置。
PCT/JP2019/001902 2019-01-22 2019-01-22 モータ制御方法、及び、モータ制御装置 WO2020152785A1 (ja)

Priority Applications (1)

Application Number Priority Date Filing Date Title
PCT/JP2019/001902 WO2020152785A1 (ja) 2019-01-22 2019-01-22 モータ制御方法、及び、モータ制御装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
PCT/JP2019/001902 WO2020152785A1 (ja) 2019-01-22 2019-01-22 モータ制御方法、及び、モータ制御装置

Publications (1)

Publication Number Publication Date
WO2020152785A1 true WO2020152785A1 (ja) 2020-07-30

Family

ID=71736619

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2019/001902 WO2020152785A1 (ja) 2019-01-22 2019-01-22 モータ制御方法、及び、モータ制御装置

Country Status (1)

Country Link
WO (1) WO2020152785A1 (ja)

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS6115392B2 (ja) * 1977-09-30 1986-04-23 Toyo Jihoki Mfg
JPH11313498A (ja) * 1998-04-28 1999-11-09 Denso Corp 界磁巻線式同期機の駆動制御装置
WO2017014249A1 (ja) * 2015-07-21 2017-01-26 株式会社デンソー 回転電機の制御装置

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS6115392B2 (ja) * 1977-09-30 1986-04-23 Toyo Jihoki Mfg
JPH11313498A (ja) * 1998-04-28 1999-11-09 Denso Corp 界磁巻線式同期機の駆動制御装置
WO2017014249A1 (ja) * 2015-07-21 2017-01-26 株式会社デンソー 回転電機の制御装置

Similar Documents

Publication Publication Date Title
JP5130716B2 (ja) モータ制御装置および電気式動力舵取装置
JP5130031B2 (ja) 永久磁石モータの位置センサレス制御装置
JP5198332B2 (ja) 永久磁石同期電動機のトルク制御装置
JP5877733B2 (ja) 電動モータの制御装置
US20150222210A1 (en) Generating motor control reference signal with control voltage budget
JP2010200430A (ja) 電動機の駆動制御装置
WO2014084009A1 (ja) 電動機の制御装置
CN111418146B (zh) 电动机的控制方法以及电动机的控制装置
JP5412820B2 (ja) 交流電動機の制御装置及び制御方法
CN113574792B (zh) 永磁同步机的控制装置
JP2006340530A (ja) 電動機の制御装置
JP3637209B2 (ja) 速度センサレスベクトル制御を用いた電力変換装置
JP5397664B2 (ja) モータ制御装置
JP4680754B2 (ja) Dcブラシレスモータのロータ角度推定方法及びdcブラシレスモータの制御装置
KR20080019131A (ko) 전압제어기를 이용한 유도전동기 및 그 제어방법
WO2020100268A1 (ja) 電動車両の制御方法、及び、制御装置
JP5996485B2 (ja) モータの駆動制御装置
WO2020152785A1 (ja) モータ制御方法、及び、モータ制御装置
JP7035818B2 (ja) 巻線界磁型同期モータの制御方法、及び、制御装置
JP7251424B2 (ja) インバータ装置及びインバータ装置の制御方法
JP6680104B2 (ja) モータの制御装置、及び、制御方法
JP2009268183A (ja) 三相交流モータの駆動装置
JP2012175776A (ja) モータ制御装置及びモータ駆動システム
JP7225561B2 (ja) モータ制御方法、及び、モータ制御装置
WO2015019905A1 (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: 19912083

Country of ref document: EP

Kind code of ref document: A1

NENP Non-entry into the national phase

Ref country code: DE

122 Ep: pct application non-entry in european phase

Ref document number: 19912083

Country of ref document: EP

Kind code of ref document: A1

NENP Non-entry into the national phase

Ref country code: JP