Method and apparatus for controlling hybrid excitation type permanent magnet synchronous motor.

<b>BACKGROUND OF THE INVENTION</b>

1. Field of the Invention

The present invention relates to method and apparatus for controlling a hybrid excitation type permanent magnet synchronous motor.

2. Description of the Prior Art

Conventionally, a speed control of a permanent magnet synchronous motor has been implemented such that three-phase current or voltage for controlling the inverter is produced by detecting a rotation speed of the motor and compensating the rotation speed command by the feedback of the detected rotation speed. On the other hand, in a torque control the feedback of the detected rotation speed has been implemented to a current command instead of the rotation speed command. In case that a constant-output operation of such a motor is implemented by keeping the power source voltage constant, as a method for broadening operable speed area of the motor by raising the maximum rotation speed, a so-called demagnetization control has been proposed, in which an electric current is applied to an armature coil for canceling the magnetic flux of the permanent magnet, in order to equivalently reduce induced voltage. Fig. 36 shows a characteristic curve obtained from such a demagnetization control.

<b>SUMMARY OF THE INVENTION</b>

It is an object of the present invention to a control method and control system of a hybrid excitation type permanent magnet synchronous motor, which control method and system enable the operable speed range of the motor to be further broadened.

A motor control method according to the present invention is for a hybrid excitation type synchronous motor which has a field of a permanent magnet and a dc excitation coil. The motor is arranged such that magnetic flux of the field is controlled by adjusting dc electric current of the dc excitation coil. The method includes the steps of: controlling the magnetic flux of the field so as to keep constant when a rotation speed of the motor is smaller than a predetermined base speed and so as to change in reverse proportion with the rotation speed of the motor when the rotation speed is larger than the predetermined base speed; controlling electric current and voltage to the armature such that a current ratio of a direct-axis component and a quadrature-axis component which components are obtained by dividing an electric current to the armature coil into the direct-axis component and the quadrature-axis component relative to the voltage induced in the armature coil, is kept constant; and selecting at least one of the magnetic flux control and current-and-voltage control to control the motor.

Further according to the invention, there is provided a control system for a hybrid excitation type synchronous motor which has a field of a permanent magnet and a dc excitation coil. The motor is arranged such that magnetic flux of the field is controlled by adjusting dc electric current of the dc excitation coil. the control system is comprised of first control means which controls the magnetic flux of the field so as to keep constant when a rotation speed of the motor is smaller than a predetermined base speed and so as to change in reverse proportion with the rotation speed of the motor when the rotation speed is larger than the predetermined base speed. Second control means controls electric current and voltage to the armature such that a current ratio of a direct-axis component and a quadrature-axis component which components are obtained by dividing an electric current to the armature coil into the direct-axis component and the quadrature-axis component relative to the voltage induced in the armature coil, is kept constant. Selecting means selects at least one of the first and second control means to control the motor.

<b>BRIEF DESCRIPTION OF THE DRAWINGS</b>

In the drawings, like reference numeral designate like parts and elements throughout all figures, in which: <ul id="ul0001" list-style="none"><li>Fig. 1 is a structural view of a hybrid excitation type permanent magnet motor controlled by a control method according to the present invention;</li><li>Fig. 2A, 2B and 2C are structural views of a rotor of the motor of Fig. 1;</li><li>Fig. 3 is a perspective view of a rotor of the motor Fig. 1;</li><li>Fig. 4 is a structural view of an exiting coil of the motor of Fig. 1;</li><li>Fig. 5 is a block diagram of a first embodiment of a motor control system of a hybrid excitation type permanent magnet motor according to the present invention;</li><li>Fig. 6 is a characteristic curve between a torque and a value K of Fig. 1;</li><li>Fig. 7 is a characteristic curve between Λ and ω;</li><li>Fig. 8 is a block diagram of a second embodiment of a motor control system according to the present invention;</li><li>Fig. 9 is a block diagram of a speed control circuit;</li><li>Fig. 10 is a block diagram of a third embodiment of the motor control system according to the present invention;</li><li>Fig. 11 is a characteristic curve between Λ and ω of Fig. 10;</li><li>Fig. 12 is a block diagram of a fourth embodiment of the motor control system according to the present invention;</li><li>Fig. 13A a block diagram of a current control type inverter, and Fig. 13B is a block diagram of a voltage control type inverter;</li><li>Fig. 14 is characteristic curves which show differences among efficiencies of the embodiments;</li><li>Figs. 15A and 15B are equivalent circuits of d-axis and q-axis;</li><li>Fig. 16 is a block diagram of a HPM control system;</li><li>Fig. 17 is a another block diagram of a HPM control system;</li><li>Fig. 18 is a flow chart which indicates a method for calculating k-table and Λ-table;</li><li>Fig. 19 is an explanatory view for the calculation of k;</li><li>Fig. 20 is an efficiency map of a motor which is obtained by a simulation;</li><li>Fig. 21 is a block diagram of a practical system;</li><li>Fig. 22A and 22B are equivalent circuits including a transient condition;</li><li>Fig. 23 is a modification of a control section of Fig. 21;</li><li>Fig. 24 is a block diagram of a non-interference PI control;</li><li>Fig. 25 is a block diagram of the d-axis and q-axis of Fig. 24;</li><li>Fig. 26 is a block diagram of a non-interference IP control;</li><li>Fig. 27 is a block diagram of the d-axis and q-axis of Fig. 26;</li><li>Figs. 28A, 28B and 28C are explanatory views which show a difference between the PI control and the IP control;</li><li>Figs. 29A and 29B are equivalent circuits where iron loss is neglected and a transient condition is included;</li><li>Fig. 30 is a block diagram of a non-interference PI control;</li><li>Fig. 31 is a block diagram of a non-interference IP control;</li><li>Fig. 32 is a block diagram of a PI control on the basis of a voltage model taking into consideration with the iron loss;</li><li>Fig. 33 is a block diagram of a PI control on the basis of a voltage model neglecting iron loss;</li><li>Fig. 34 is a block diagram of a control system of PM;</li><li>Fig. 35 is an efficiency map obtained by a simulation of the PM; and</li><li>Fig. 36 is a characteristic curve in the torque-constant area and the output-constant area of a conventional motor control.</li></ul>

<b>DETAILED DESCRIPTION OF THE INVENTION</b>

Figs. 1 to 4 show a structure of a hybrid excitation type permanent magnet synchronous motor whose control is realized by a method and control system according to the present invention.

As shown in Fig. 1, an armature 1 of a stator is constituted by an armature core 2, an armature coil 3 and a cylindrical yoke 4. The armature core 2 is divided into two parts in the axial direction of the motor. One side of the armature core 2 is a N-pole core 2a, and the other side of the armature core 2 is a S-pole core 2b. A ring-shaped DC (direct current) excitation coil 5 as shown in Fig. 4 is disposed between the N-pole core 2a and the S-pole core 2b. The N-pole and S-pole cores 2a and 2b are magnetically connected with each other through the yoke 4 and mechanically supported by the yoke 4. The armature coil 3 is disposed so as to cross the N-pole and S-pole cores 2a and 2b. The excitation coil 5 is formed by winding a wire 5a and treating it by the insulation process. The number of turns of the excitation coil 5 is adjusted so as to generate a necessary magnetmotive force upon matching with a capacity of the electric source and the machinery dimension.

A rotor 11 is constituted by a rotor core 12 and a plurality of permanent magnets (PM) 13. The rotor core 12 is fixedly supported to a yoke 14 connected to a shaft 15. The rotor core 12 has salient pole portions 12a which project form the rotor core 12 and function as a salient pole, and the salient pole portions 12a are disposed at portions where the permanent magnets 13 are not located. As shown in Figs. 1, 2A and 2B, the permanent magnet 13 is embedded in the rotor 11 and covered with the rotor core 12. A slit 20 is formed between each permanent magnet 13 and each pole salient portion 12a in relation that the permanent magnets 13 and the salient pole portions 12a are fixed on the yoke. In order to reinforce the rotor 11, non-magnetic reinforcement plates 21 as shown in Fig. 2C are disposed in the several portions in the axial direction. The non-magnetic reinforcement plates 21, the salient pole portions 12a and the rotor cores 12 are integrally fitted with each other by inserting aluminum-alloy die-cast member or copper bar to slots 22. The inserted member and a pair of end rings installed at an axial end form an electric current passage. The salient pole portions 12a are divided into N-pole salient pole portions 12aN and S-pole salient pole portions 12aS which are dividedly disposed so as to be opposite to the N-pole and S-pole cores 2a and 2b, respectively. The permanent magnets 13 are fixedly attached on the rotor core 12, and the rotor core 12 is inserted into the yoke 12 to be supported thereby.

The N-pole and S-pole salient pole portions 12aN and 12aS are formed such that the length thereof corresponds to those of the N-pole core 2a and the S-pole core 2b, respectively, and that width thereof are constant along the circumferential direction. Further, the N-pole salient pole portions 12aN are arranged side by side with the N-pole side of the permanent magnets 13, as shown in Fig. 2A. The S-pole salient pole portions 12aS are arranged side by side with the S-pole sides of the permanent magnets 13, as shown in Fig. 2B. Furthermore, the N-pole salient pole portion 12aN and the S-pole side of the permanent magnets 13 are aligned along the axial direction of the motor at predetermined intervals. Similarly, the S-pole salient pole portion 12aS and the N-pole side of the permanent magnets 13 are aligned along the axial direction of the motor at predetermined intervals.

That is, the rotor 11 is formed as shown in Fig. 3, in which the N-pole salient pole portions 12aN and the N-pole side of the permanent magnets 13 are alternately arranged in the circumferential direction. The S-pole salient pole portions 12aS and the S-pole side of the permanent magnets 13 are alternately arranged in the circumferential direction. Furthermore, the arrangement units of the N-pole side and the S-pole side are separated with each other by a width of the excitation coil 5, and the salient pole portions 12a and the permanent magnets 13 are aligned in the axial direction. The number of the salient pole portions 12a is the same as that of the permanent magnets 13.

Although the motor in Figs. 1 to 3 has been shown and described such that the permanent magnet 13 is disposed at six poles, it will be understood that the number of the poles does not limited to six. Furthermore, while the motor of Figs. 1 to 3 has been shown and described such that the surface of the salient pole portion 12a and that of the permanent magnet 13 are located within a same circumferential surface, the salient pole portions 12a may be formed to further project toward the armature 1 so as to reduce the gap between the salient pole portions 12 and the armature 1. This arrangement will cause an effect that the magnetic flux passing through the salient pole portions 12 is increased. Because of the same reason, the width of the salient pole portions 12a may be broadened.

Referring to Figs. 6 to 35, there are shown embodiments of a control system of the above-mentioned motor, according to the present invention.

In general, a quadrature-axis reactance <maths id="math0001" num=""><math display="inline"><mrow><msub><mrow><mtext>X</mtext></mrow><mrow><mtext>q</mtext></mrow></msub><msub><mrow><mtext> = ωL</mtext></mrow><mrow><mtext>q</mtext></mrow></msub></mrow></math><img file="EP0631373A2_D0001.tif" /></maths> (ω: electric source angle frequency, Lq: a quadrature-axis self-inductance) and a direct-axis reactance <maths id="math0002" num=""><math display="inline"><mrow><msub><mrow><mtext>X</mtext></mrow><mrow><mtext>d</mtext></mrow></msub><msub><mrow><mtext> = ωL</mtext></mrow><mrow><mtext>d</mtext></mrow></msub></mrow></math><img file="EP0631373A2_D0002.tif" /></maths> (L_d: direct-axis self-inductance) exist in a motor. A voltage-equation and a generated torque thereof are represented by the following equations (1) and (2), respectively:<maths id="math0003" num=""><img file="EP0631373A2_D0003.tif" /></maths>    where V_d is a voltage on a direct-axis, V_q is a voltage on a quadrature-axis, i_d is a current on the direct-axis, i_q is a current on the quadrature-axis, R_l is a coil resistance, Λ is the number of an interlinkage magnetic-flux lines, Pole is the number of poles, and P is a differential operator <maths id="math0004" num=""><math display="inline"><mrow><mtext>P = d/dt</mtext></mrow></math><img file="EP0631373A2_D0004.tif" /></maths> (t is time).

In general, Λ can be generally expressed by a function of a direct-current excitation current I_f, that is, <maths id="math0005" num=""><math display="inline"><mrow><msub><mrow><mtext>Λ = f (I</mtext></mrow><mrow><mtext>f</mtext></mrow></msub><mtext>)</mtext></mrow></math><img file="EP0631373A2_D0005.tif" /></maths> . Accordingly, it is possible to control the field magnetic flux by the control of I_f.

Above-described equations are used as a basic equation in a wide-range variable speed control method of a motor having a DC exciting-current coil.

Control of a brushless DC motor

The manner of operation of the control of a brushless DC motor in a permanent magnet embedded type hybrid excitation type permanent magnet motor will be discussed hereinafter.

In a normal brushless DC motor, assuming that I_d is set at 0 (I_d = 0), only I_q is applied to the motor, and Λ is controlled such that Λ is kept constant at a speed lower than a base speed and Λ is varied in reverse-proportion with a rotation speed of the motor at a speed higher than a base speed to ensure a constant-output properties; a constant-torque output characteristic is obtained. In this case, the generated torque is represented by the following equation (3):<maths id="math0006" num=""><img file="EP0631373A2_D0006.tif" /></maths>    Fig. 5 shows a control system of a brushless DC motor in that I_d is kept at 0 (I_d = 0 control). The control system is of a speed control system whose main circuit comprises a torque calculator 30 for obtaining a torque T* from a command speed ω*, an electric current calculator 31 for obtaining a quadrature-axis electric current I_q form the torque T*, a three-phase electric current calculator 31 for obtaining a three-phase electric current from the quadrature-axis electric current by applying the equation (3) and the following equations (4), a current controller 33, and an inverter 35. <maths id="math0007" num=""><math display="inline"><mrow><msub><mrow><mtext>i</mtext></mrow><mrow><mtext>u</mtext></mrow></msub><msub><mrow><mtext> = √2I</mtext></mrow><mrow><mtext>l</mtext></mrow></msub><mtext> cosϑ</mtext></mrow></math><img file="EP0631373A2_D0007.tif" /></maths><maths id="math0008" num=""><math display="block"><mrow><msub><mrow><mtext>i</mtext></mrow><mrow><mtext>v</mtext></mrow></msub><msub><mrow><mtext> = √2I</mtext></mrow><mrow><mtext>l</mtext></mrow></msub><mtext> cos(ϑ - </mtext><mfrac><mrow><mtext>2</mtext></mrow><mrow><mtext>3</mtext></mrow></mfrac><mtext>π)   (4)</mtext></mrow></math><img file="EP0631373A2_D0008.tif" /></maths><maths id="math0009" num=""><math display="block"><mrow><msub><mrow><mtext>i</mtext></mrow><mrow><mtext>w</mtext></mrow></msub><msub><mrow><mtext> = √2I</mtext></mrow><mrow><mtext>l</mtext></mrow></msub><mtext> cos(ϑ + </mtext><mfrac><mrow><mtext>2</mtext></mrow><mrow><mtext>3</mtext></mrow></mfrac><mtext>π)</mtext></mrow></math><img file="EP0631373A2_D0009.tif" /></maths> where ϑ is a rotor position angle, <maths id="math0010" num=""><math display="inline"><mrow><msub><mrow><mtext>I</mtext></mrow><mrow><mtext>l</mtext></mrow></msub><mtext>= Iq/√3</mtext></mrow></math><img file="EP0631373A2_D0010.tif" /></maths> .

A rotation speed of the motor 35 is controlled by the inverter 34 and detected by a rotation speed detector 35. The detected signal by the rotation speed detector 35 is inputted as a rotational position to the three-phase current calculator 32 through a position detecting circuit 37. Further, the detected signal by the rotation speed detector 36 is modulated as a present speed in a speed detecting circuit 38 and is then added to the speed command ω*. The modulated signal at the speed detecting circuit 38 is transformed into a field-magnet magnetic flux Λ at an calculator 39. An excitation current calculator 40 obtains an electric current i_f from the field magnetic-flux, and a DC electric current control circuit 41 is controlled according to the obtained electric current i_f. Accordingly, the exciting current is varied according to the speed such that the field magnetic-flux Λ is kept constant when the speed is lower than a base speed and is gradually decreased in reverse proportion with the speed when the speed is higher than the base speed. This control extends an operable speed area of the motor 35.

If the control system of Fig. 5 is used as a torque control, the left hand side part divided by a broken line of Fig. 5 is not necessary in such a torque control system.

Control during low speed and/or low-load operation

In addition to the above-mentioned current control, this embodiment provides a further improved control by which the motor 30 is operable in a further high-speed a high-torque control by utilizing the characteristics that the quadrature-axis reactance of this motor 35 is large. Assuming that the motor 35 is put in a stationary condition in that R₁ is neglected since the R₁ is extremely small as compared with the reactance in the equations (1) and (2), the equations (1) and (2) are expressed by the following equations (5) and (6): <maths id="math0011" num=""><math display="inline"><mrow><msub><mrow><mtext>V</mtext></mrow><mrow><mtext>d</mtext></mrow></msub><msub><mrow><mtext> = -ωL</mtext></mrow><mrow><mtext>q</mtext></mrow></msub><msub><mrow><mtext>I</mtext></mrow><mrow><mtext>q</mtext></mrow></msub></mrow></math><img file="EP0631373A2_D0011.tif" /></maths><maths id="math0012" num=""><math display="inline"><mrow><msub><mrow><mtext>V</mtext></mrow><mrow><mtext>q</mtext></mrow></msub><msub><mrow><mtext> = ωL</mtext></mrow><mrow><mtext>d</mtext></mrow></msub><msub><mrow><mtext>I</mtext></mrow><mrow><mtext>d</mtext></mrow></msub><mtext> + ωΛ   (5)</mtext></mrow></math><img file="EP0631373A2_D0012.tif" /></maths><maths id="math0013" num=""><img file="EP0631373A2_D0013.tif" /></maths>    If the equation (5) is applied to a high-speed condition, ωΛ and I_d become considerably large. Therefore, the equation (5) is applied only for a low-speed and/or low-load operation. The voltage of the motor terminal is obtained from the equation (5) as follows:<maths id="math0014" num=""><img file="EP0631373A2_D0014.tif" /></maths>    When a possible output voltage of the inverter is V_max, the motor is operable within the range in that V ≦ V_max. That is, the operable range of the motor is under a condition V ≦ V_max.

It is assumed that the ratio between I_d and I_q is constant in order to uniformly determine the electric current relative to the torque command. That is, when the ratio between the I_d and I_q is K, the following equation (8) is obtained. <maths id="math0015" num=""><math display="inline"><mrow><msub><mrow><mtext>I</mtext></mrow><mrow><mtext>d</mtext></mrow></msub><msub><mrow><mtext> = KI</mtext></mrow><mrow><mtext>q</mtext></mrow></msub><mtext>   (8)</mtext></mrow></math><img file="EP0631373A2_D0015.tif" /></maths>    By applying the equation (8) to the equation (6), the following equation (9) is obtained.<maths id="math0016" num=""><img file="EP0631373A2_D0016.tif" /></maths>    As is clear from the equation (9), if the magnetic flux Λ is constant and the torque T is obtained, the electric current I_q is obtained.

In order to obtain an experimental data relative to the value K, as to a hybrid excitation type permanent magnet motor designed as a 45kW-8 pole motor, the relationship between the value K and the torque was searched in a condition that I_l is kept at 200A (I_l = 200A). Fig. 6 shows such examination data which indicates that the torque T is largely varied according to the value K. The result of Fig. 6 indicates that it is preferable to set the value K within a range from -0.3 to -1.2 as compared with the control at K= 0 (I_d = 0), in order to effectively obtain the output torque T. Although the torque T is determined from the valve K and the primary current, it is preferable to set K &lt; 0 in a K value constant control in order to obtain a high-torque relative by the small electric current.

If the value K is determined as mentioned above and the interlinkage magnetic flux Λ is varied within the torque-constant range and a speed higher than the base speed as shown in Fig. 7, the electric current value I_q is univocally determined relative to the torque T. In Fig. 7, Λ_max is a maximum value of Λ which is determined from the magnetic saturation, and Λ is a multiple of Λ_max and the speed ratio <maths id="math0017" num=""><math display="inline"><mrow><msub><mrow><mtext>ω0/ ω (Λ = Λ</mtext></mrow><mrow><mtext>max</mtext></mrow></msub><mtext> × ω0/ ω)</mtext></mrow></math><img file="EP0631373A2_D0017.tif" /></maths> . As a result, I_q is obtained from the determination of K and Λ in the equation (9), and I_d is obtained from the equation (8). Accordingly, each of phase currents is obtained as follows: <maths id="math0018" num=""><math display="inline"><mrow><msub><mrow><mtext>i</mtext></mrow><mrow><mtext>u</mtext></mrow></msub><msub><mrow><mtext> = √2I</mtext></mrow><mrow><mtext>l</mtext></mrow></msub><mtext> sinϑ</mtext></mrow></math><img file="EP0631373A2_D0018.tif" /></maths><maths id="math0019" num=""><math display="block"><mrow><msub><mrow><mtext>i</mtext></mrow><mrow><mtext>v</mtext></mrow></msub><msub><mrow><mtext> = √2I</mtext></mrow><mrow><mtext>l</mtext></mrow></msub><mtext> sin(ϑ + φ - </mtext><mfrac><mrow><mtext>2</mtext></mrow><mrow><mtext>3</mtext></mrow></mfrac><mtext>π)   (10)</mtext></mrow></math><img file="EP0631373A2_D0019.tif" /></maths><maths id="math0020" num=""><math display="block"><mrow><msub><mrow><mtext>i</mtext></mrow><mrow><mtext>w</mtext></mrow></msub><msub><mrow><mtext> = √2I</mtext></mrow><mrow><mtext>l</mtext></mrow></msub><mtext> sin(ϑ + φ + </mtext><mfrac><mrow><mtext>2</mtext></mrow><mrow><mtext>3</mtext></mrow></mfrac><mtext>π)</mtext></mrow></math><img file="EP0631373A2_D0020.tif" /></maths>    where Q is a rotor position angle, <maths id="math0021" num=""><math display="inline"><mrow><msub><mrow><mtext>φ = tan⁻¹(I</mtext></mrow><mrow><mtext>q</mtext></mrow></msub><msub><mrow><mtext>/I</mtext></mrow><mrow><mtext>d</mtext></mrow></msub><mtext>)</mtext></mrow></math><img file="EP0631373A2_D0021.tif" /></maths> , <maths id="math0022" num=""><math display="inline"><mrow><msub><mrow><mtext>I</mtext></mrow><mrow><mtext>l</mtext></mrow></msub><mtext> = √</mtext><mover accent="true"><mrow><msub><mrow><mtext>(I</mtext></mrow><mrow><mtext>d</mtext></mrow></msub><msub><mrow><mtext>² + I</mtext></mrow><mrow><mtext>q</mtext></mrow></msub><mtext>²)/3</mtext></mrow><mo>¯</mo></mover></mrow></math><img file="EP0631373A2_D0022.tif" /></maths> .

As a result, if the K value is set at minus, the direct-axis current I_d becomes minus. This lowers V_q and the operable range of the motor is extended.

Fig. 8 shows a block diagram of a I_d- I_q ratio constant control system of a torque control circuit. The quadrature-axis current I_q is obtained from the torque command T* and a feed-back magnetic flux command Λ in a calculator 50. a three-phase current calculator receives the I_q and the I_d obtained by the multiple of the I_q and K. The electric current values i_u, i_v, i_w of the equation (10) are outputted from the three-phase current calculator 51 and inputted to an inverter 54 through a current controller 53. Upon receipt of the electric current, the inverter 54 is operated.

The detected value by a rotation speed detector 56 for a motor 55 is inputted to the three-phase current calculator 50 after modulated as a rotational position in a position detecting section 57. A speed detecting section 58 receives the detected value at the rotation speed detector 56 and outputs a speed ω. The speed ω is inputted to a calculator 59, which has a property shown in Fig. 8, outputs the magnetic flux Λ. The output Λ from the calculator 59 is inputted to the I_q calculator 50 and an exiting current calculator 60. An electric current if is outputted from the exciting current calculator 60 and inputted to a DC current controller 61.

Fig. 9 shows a speed control circuit to which a speed command ω* is inputted. The speed control circuit is generally similar to the torque control circuit of Fig. 9 except that a torque calculator 62 for calculating a torque command T* from the speed command ω* is provided in front of the I_q calculator 50 and a feed back command from a speed detecting section 58 is added to the speed command ω*.

Control at high-speed range or heavy-loading

Although at low-speed range of light-loading the operating property of the motor is improved by changing the magnetic flux Λ while keeping K constant, in this high-speed or heavy-loading the voltage V in the equations (6) and (7) frequently exceeds the maximum voltage V_max of the inverter due to the increase of the speed ω. Accordingly, in this control the motor is operated by keeping the voltage V at V_max ( <maths id="math0023" num=""><math display="inline"><mrow><msub><mrow><mtext>V= V</mtext></mrow><mrow><mtext>max</mtext></mrow></msub></mrow></math><img file="EP0631373A2_D0023.tif" /></maths> ). From the equation (7) and <maths id="math0024" num=""><math display="inline"><mrow><msub><mrow><mtext>V= V</mtext></mrow><mrow><mtext>max</mtext></mrow></msub></mrow></math><img file="EP0631373A2_D0024.tif" /></maths> , the following equation (11) is derived.<maths id="math0025" num=""><img file="EP0631373A2_D0025.tif" /></maths>    Although the equation (11) shows I_q for keeping the terminal voltage at V_max, in order to facilitate the control, the root part of the equation (11) is set at 0<maths id="math0026" num=""><img file="EP0631373A2_D0026.tif" /></maths> That is, so as to satisfy the equation (12) the magnetic flux Λ is changed for the control of torque. <maths id="math0027" num=""><math display="inline"><mrow><msub><mrow><mtext>I</mtext></mrow><mrow><mtext>q</mtext></mrow></msub><msub><mrow><mtext> = V</mtext></mrow><mrow><mtext>max</mtext></mrow></msub><msub><mrow><mtext>/ωL</mtext></mrow><mrow><mtext>q</mtext></mrow></msub><mtext>   (12)</mtext></mrow></math><img file="EP0631373A2_D0027.tif" /></maths>    When I_q is controlled according to the equation (12), I_d is represented as follows: <maths id="math0028" num=""><math display="inline"><mrow><msub><mrow><mtext>I</mtext></mrow><mrow><mtext>d</mtext></mrow></msub><msub><mrow><mtext> = - Λ/L</mtext></mrow><mrow><mtext>d</mtext></mrow></msub><mtext>   (13)</mtext></mrow></math><img file="EP0631373A2_D0028.tif" /></maths>    By substituting the equations (12) and (13) into the equation (6), the following equation (14) is obtained.<maths id="math0029" num=""><img file="EP0631373A2_D0029.tif" /></maths>    As is clear from the equation (14), when the speed ω is kept constant, the torque T is in proportion with the number Λ of the interlinkage magnetic flux. Accordingly, the torque control is carried out by controlling Λ. In this control, From the equations (5) and (13), <maths id="math0030" num=""><math display="inline"><mrow><msub><mrow><mtext>V</mtext></mrow><mrow><mtext>q</mtext></mrow></msub><msub><mrow><mtext> = ωL</mtext></mrow><mrow><mtext>d</mtext></mrow></msub><msub><mrow><mtext>L</mtext></mrow><mrow><mtext>d</mtext></mrow></msub><mtext> + ωΛ = 0</mtext></mrow></math><img file="EP0631373A2_D0030.tif" /></maths> is obtained. Accordingly, when the torque command is applied, the following equation (15) is obtained.<maths id="math0031" num=""><img file="EP0631373A2_D0031.tif" /></maths>    As a result, V_max of the equation (15) is decided by the inverter, and by obtaining Λ after the detection of ω, I_q and I_d represented by the equations (12) and (13). Accordingly, it is possible to extend the operable range of the motor by controlling I_q and I_d so as to keep V_max constant, that is, by keeping V_q = 0.

Since Λ is saturated if too large, it should be kept to be smaller than Λ_max. That is, the maximum torque is limited by the following equation (16):<maths id="math0032" num=""><img file="EP0631373A2_D0032.tif" /></maths>    Fig. 10 shows a block diagram of the above-mentioned V_q=0 control system. As shown in Fig. 10, a Λ calculator 70 executes the calculation of Λ represented by the equation (15) upon receiving the torque command T* and the speed command ω*. Then, limiter 71 limits Λ within Λ_max. I_d-calculator 72 executes the calculation of I_d from the equation (13), and I_q-calculator 73 executes the calculation of I_q from the equation (12). Following this, a three-phase current calculator 74 calculates the current values i_u, i_v and i_w, and the obtained valves are used for driving the inverter 76 through the current control section 75.

The detected value by the rotation speed detector 78 for the motor 77 is modulated as a rotational position at a position detecting section 79 and inputted to a three-phase current calculator 74. Further, the detected value by the rotation speed detector 78 is inputted to a speed detecting section 80 where the speed ω is outputted. The outputted speed ω is feed-backed to the Λ calculator 70 and the I_q calculator 73.

The output Λ of the limiter 71 is inputted to the I_d calculator 72 and the exciting current calculator 81 for obtaining the current i_f. A DC electric current control circuit 82 is controlled according to I_d and i_f.

Although Fig. 10 shows a torque control circuit, by adding the speed feedback and the torque calculator 62 shown in a left side portion divided by the broken line in Fig. 9, the control circuit of Fig. 6 can be used as a speed control circuit.

Accordingly, by suppressing V_max so as to implement V_q= 0 control, the terminal voltage V is decided. This control becomes suitable for the heavy-loading and high-speed range control.

Wide-range speed control

In order to smoothly implement a wide-range speed control, it is necessary to interconnect the I_d-I_q ratio-constant control and the V_q= 0 control. Although it is possible to properly select the value K (I_d-I_q ratio), it is necessary to select the value K so as to avoid discontinuous of the electric current and the magnetic flux at a transient between the I_d-I_q ratio-constant control and the V_q= 0 control.

During the V_q = 0 control, <maths id="math0033" num=""><math display="inline"><mrow><msub><mrow><mtext>I</mtext></mrow><mrow><mtext>q</mtext></mrow></msub><msub><mrow><mtext> = V</mtext></mrow><mrow><mtext>max</mtext></mrow></msub><msub><mrow><mtext>/ωL</mtext></mrow><mrow><mtext>q</mtext></mrow></msub></mrow></math><img file="EP0631373A2_D0033.tif" /></maths> , <maths id="math0034" num=""><math display="inline"><mrow><msub><mrow><mtext>I</mtext></mrow><mrow><mtext>d</mtext></mrow></msub><msub><mrow><mtext> = 1L/L</mtext></mrow><mrow><mtext>d</mtext></mrow></msub></mrow></math><img file="EP0631373A2_D0034.tif" /></maths> is obtained from the equations (12) and (13). When the control is changed between the I_d-I_q ratio-constant control and the V_q= 0 control, <maths id="math0035" num=""><math display="inline"><mrow><msub><mrow><mtext>I</mtext></mrow><mrow><mtext>d</mtext></mrow></msub><msub><mrow><mtext> = KI</mtext></mrow><mrow><mtext>q</mtext></mrow></msub></mrow></math><img file="EP0631373A2_D0035.tif" /></maths> . Accordingly, the following equation (17) is obtained. <maths id="math0036" num=""><math display="inline"><mrow><msub><mrow><mtext>KI</mtext></mrow><mrow><mtext>q</mtext></mrow></msub><msub><mrow><mtext> = -Λ/L</mtext></mrow><mrow><mtext>d</mtext></mrow></msub><mtext>   (17)</mtext></mrow></math><img file="EP0631373A2_D0036.tif" /></maths>    By substituting the equation (17) to the above-mentioned equation <maths id="math0037" num=""><math display="inline"><mrow><msub><mrow><mtext>I</mtext></mrow><mrow><mtext>q</mtext></mrow></msub><msub><mrow><mtext> = V</mtext></mrow><mrow><mtext>max</mtext></mrow></msub><msub><mrow><mtext>/ωL</mtext></mrow><mrow><mtext>q</mtext></mrow></msub></mrow></math><img file="EP0631373A2_D0037.tif" /></maths> , the following equation (18) is obtained. <maths id="math0038" num=""><math display="inline"><mrow><msub><mrow><mtext>K·(V</mtext></mrow><mrow><mtext>max</mtext></mrow></msub><msub><mrow><mtext>/ωL</mtext></mrow><mrow><mtext>q</mtext></mrow></msub><msub><mrow><mtext>) = -Λ/L</mtext></mrow><mrow><mtext>d</mtext></mrow></msub><mtext>   (18)</mtext></mrow></math><img file="EP0631373A2_D0038.tif" /></maths>    Further the following equation (19) is obtained from the equation (18). <maths id="math0039" num=""><math display="inline"><mrow><msub><mrow><mtext>K = -(L</mtext></mrow><mrow><mtext>q</mtext></mrow></msub><msub><mrow><mtext>/L</mtext></mrow><mrow><mtext>d</mtext></mrow></msub><msub><mrow><mtext>)·(ωΛ/V</mtext></mrow><mrow><mtext>max</mtext></mrow></msub><mtext>)   (19)</mtext></mrow></math><img file="EP0631373A2_D0039.tif" /></maths>    Although the speed ω is variable in the equation (19), in the output-constant operation range during the I_d-I_q ratio constant control, Λω is kept at a constant value E₀ and therefore the induced voltage becomes constant, as shown in Fig. 11. As a result, if the motor is designed so that the switching of the control method is implemented within the constant-output operation range, the value K can be defined by the following equation (20): <maths id="math0040" num=""><math display="inline"><mrow><msub><mrow><mtext>K = -(L</mtext></mrow><mrow><mtext>q</mtext></mrow></msub><msub><mrow><mtext>/L</mtext></mrow><mrow><mtext>d</mtext></mrow></msub><msub><mrow><mtext>)·(E₀/V</mtext></mrow><mrow><mtext>max</mtext></mrow></msub><mtext>)   (20)</mtext></mrow></math><img file="EP0631373A2_D0040.tif" /></maths> Accordingly, it becomes possible to implement the K constant control within the constant-output operation range, and therefore the switching between the controls may be implemented when the speed of the motor is higher than the base-speed. That is, the switching from the I_d-I_q ratio constant control to the V_q= 0 control may be implemented when <maths id="math0041" num=""><math display="inline"><mrow><msub><mrow><mtext>I</mtext></mrow><mrow><mtext>q</mtext></mrow></msub><msub><mrow><mtext> ≧ V</mtext></mrow><mrow><mtext>max</mtext></mrow></msub><msub><mrow><mtext>/ωL</mtext></mrow><mrow><mtext>q</mtext></mrow></msub></mrow></math><img file="EP0631373A2_D0041.tif" /></maths> is satisfied.

When the motor is controlled such that I_d and I_q satisfy the equations (13) and (17), the following equations (21) and (22) are obtained:<maths id="math0042" num=""><img file="EP0631373A2_D0042.tif" /></maths>    By compounding the control systems of Figs. 8 and 10 in consideration with the equations (21) and (22), the control system of Fig. 12 is obtained. In Fig. 12, the switching between the I_d-I_q ratio constant control and the V_q = 0 control is implemented by the switching of switches SW₁, SW₂, and SW₃. When the I_q calculator 50 calculates the following equation (23) is implemented, the switches SW₁ and Sw₂ are connected to a-contact side, and the switch SW₃ is opened. On the other hand, when the I_q calculator 50 calculate the following equation (24), the switches SW₁ and SW₂ are connected to b-contact, and the switch SW₃ is closed.

When<maths id="math0043" num=""><img file="EP0631373A2_D0043.tif" /></maths>    When<maths id="math0044" num=""><img file="EP0631373A2_D0044.tif" /></maths>    If a speed feedback element is applied to the control system of Fig. 12, it becomes a speed control system.

Although the embodiments of the control system have been shown and described to comprise the inverter 34, 54, 76 of a current control type which is driven by the current control section 33, 53, 75 as shown in Fig. 13A, an inverter of a voltage control type as shown in Fig. 13B may be used in the control system since it is possible to calculate the voltage V_d and V_q by substituting I_d, I_q and Λ into the equation (1).

With the above-mentioned embodiments of the motor control system according to the present invention, it is possible to extend the operable range of the PM motor which includes a DC field magnet coil. Furthermore, it is possible to form the motor more smaller and to improve the efficiency due to the increase of the maximum torque and the raising of the maximum rotation speed. This further enables the current capacity of the electric source to be smaller.

In concrete, the experiment by the motor of 45kW-8P 3000rpm was implemented, and the data as shown in Fig. 14 was obtained. As is clear from the data of Fig. 14, the operable range ① of a field magnet control is extended as compared with the operable range ⑤ in case of no- field magnet control. By the Id-Iq ratio constant control, the operable range is further extended as shown by the line ②, and by the V_q = 0 control within a constant-output range the operable range thereof is extended as shown by the line ③. Further, maximum advantage is obtained by the combination of I_d-I_q ratio control and the V_q = 0 control as shown by the line ④.

Hereinbefore, the control method such as a high-speed range and heavy-loading control, and a wide-range speed control of the brushless DC motor, which control method is implemented by obtaining I or V on the basis of the voltage equation, has been mentioned.

Next, another control method will be discussed hereinafter. The control method is a method for deciding the current ratio and the magnetic flux by taking account the iron loss and using previously prepared a current ration (k) table and a magnet flux (Λ) table.

Fig. 15A and 15B show the equivalent circuits of this control system, in which the respective circuits of the q-axis and d-axis are described and where V_d is a d-axis component of an armature voltage, V_q is q-axis component of an armature voltage of the q-axis, I_d is a d-axis component of an armature current, I_q is a q-axis component of an armature current, R₁ is a resistance of the armature, R_c is an equivalent iron-loss resistance, L_d is a d-axis inductance, L_q is a q-axis inductance, ω is an angle frequency of an electric source, Λ is an interlinkage magnetic flux of the armature coil, I_cd and I_cq are current iron losses, and I_d' and I_q' are current loading.

As is clear from Figs. 15A and 15B, the current values I_d, I_q, I_cd and I_cq are obtained by the following equations (25) and (26): <maths id="math0045" num=""><math display="inline"><mrow><msub><mrow><mtext>I</mtext></mrow><mrow><mtext>d</mtext></mrow></msub><msub><mrow><mtext> = I</mtext></mrow><mrow><mtext>d</mtext></mrow></msub><msub><mrow><mtext>' + I</mtext></mrow><mrow><mtext>cd</mtext></mrow></msub></mrow></math><img file="EP0631373A2_D0045.tif" /></maths><maths id="math0046" num=""><math display="inline"><mrow><msub><mrow><mtext>I</mtext></mrow><mrow><mtext>q</mtext></mrow></msub><msub><mrow><mtext> = I</mtext></mrow><mrow><mtext>q</mtext></mrow></msub><msub><mrow><mtext>' + I</mtext></mrow><mrow><mtext>cq</mtext></mrow></msub><mtext>   (25)</mtext></mrow></math><img file="EP0631373A2_D0046.tif" /></maths><maths id="math0047" num=""><img file="EP0631373A2_D0047.tif" /></maths>    When Λ_m is an interlinkage magnetic flux only by the permanent magnet, M_f is a mutual inductance of a field magnet coil and an armature coil, and I_f is an exciting current, the following equation (27) is obtained. <maths id="math0048" num=""><math display="inline"><mrow><msub><mrow><mtext>Λ = Λ</mtext></mrow><mrow><mtext>m</mtext></mrow></msub><msub><mrow><mtext> + M</mtext></mrow><mrow><mtext>f</mtext></mrow></msub><msub><mrow><mtext>I</mtext></mrow><mrow><mtext>f</mtext></mrow></msub><mtext>   (27)</mtext></mrow></math><img file="EP0631373A2_D0048.tif" /></maths>    The following equations (28) and (29) are obtained by providing a voltage equation of Figs. 15A and 15B according to the equations (25) and (26).<maths id="math0049" num=""><img file="EP0631373A2_D0049.tif" /></maths>    where<maths id="math0050" num=""><img file="EP0631373A2_D0050.tif" /></maths>    Furthermore, the equation of the torque can be represented by the following equation (30):<maths id="math0051" num=""><img file="EP0631373A2_D0051.tif" /></maths>    where Pole is the number of poles.

As a result, when the hybrid excitation type permanent magnet synchronous motor is operated by the torque control, the operation is accomplished by controlling I_d', I_q' and Λ.

If it is defined that <maths id="math0052" num=""><math display="inline"><mrow><msub><mrow><mtext>I</mtext></mrow><mrow><mtext>d</mtext></mrow></msub><msub><mrow><mtext>' = kI</mtext></mrow><mrow><mtext>q</mtext></mrow></msub><mtext>'   (31)</mtext></mrow></math><img file="EP0631373A2_D0052.tif" /></maths> , the following equation (32) is derived from the equation (30).<maths id="math0053" num=""><img file="EP0631373A2_D0053.tif" /></maths>    As is clear from the equation (32), when the current ratio k of the d-axis component of current and the q-axis component of it and the magnetic flux Λ relative to the proper speed ωr and torque T can be obtained, it is possible to decide the current I_q'.

Hereinafter, the field magnet control method and the creating method of the control table will be discussed with reference to Figs. 16 to 35. Fig. 16 shows a block diagram of a control system of the hybrid excitation type permanent magnet synchronous motor of a current control type. In Fig. 16, the torque command T* is inputted to the k calculating circuit 101 and a Λ calculating circuit 102, and an appropriate values k* and Λ* corresponding to the speed ω are selected from the k-table and the Λ-table, respectively, in the circuit 101 and 102.

In a I_q' calculating circuit 103, the q-axis current component I_q' is obtained from k*, the torque command T* and the magnetic flux Λ*. In a multiplier 104, a d-axis current component is obtained by multiplying the current I_q' and k*. In an adder 106, an output I_q is obtained by adding I_q' and I_cq (I_q' + I_cq). In a I_cd calculator 107, an iron loss current I_cd is obtained by substituting the q-axis current component I_q' and the equivalent iron loss R_c into the equation (26). the output I_d is obtained by adding I_d' and I_cd (I_d' + I_cd) in an addor 109. In a three-phase current calculating circuit 110, three-phase current i_u*, i_v* and i_w* are obtained by substituting the outputs I_q and I_d, and the output q into the following equations (33): <maths id="math0054" num=""><math display="inline"><mrow><msub><mrow><mtext>i</mtext></mrow><mrow><mtext>u</mtext></mrow></msub><msub><mrow><mtext> = √2I</mtext></mrow><mrow><mtext>l</mtext></mrow></msub><mtext> sin(ϑ+ φ)</mtext></mrow></math><img file="EP0631373A2_D0054.tif" /></maths><maths id="math0055" num=""><math display="block"><mrow><msub><mrow><mtext>i</mtext></mrow><mrow><mtext>v</mtext></mrow></msub><msub><mrow><mtext> = √2I</mtext></mrow><mrow><mtext>l</mtext></mrow></msub><mtext> sin(ϑ + φ - </mtext><mfrac><mrow><mtext>2</mtext></mrow><mrow><mtext>3</mtext></mrow></mfrac><mtext>π)   (33)</mtext></mrow></math><img file="EP0631373A2_D0055.tif" /></maths><maths id="math0056" num=""><math display="block"><mrow><msub><mrow><mtext>i</mtext></mrow><mrow><mtext>w</mtext></mrow></msub><msub><mrow><mtext> = √2I</mtext></mrow><mrow><mtext>l</mtext></mrow></msub><mtext> sin(ϑ + φ + </mtext><mfrac><mrow><mtext>2</mtext></mrow><mrow><mtext>3</mtext></mrow></mfrac><mtext>π)</mtext></mrow></math><img file="EP0631373A2_D0056.tif" /></maths>    where ϑ is a rotor position angle and <maths id="math0057" num=""><math display="inline"><mrow><msub><mrow><mtext>tan⁻¹ (I</mtext></mrow><mrow><mtext>q</mtext></mrow></msub><msub><mrow><mtext>/I</mtext></mrow><mrow><mtext>d</mtext></mrow></msub><mtext>)</mtext></mrow></math><img file="EP0631373A2_D0057.tif" /></maths> , and <maths id="math0058" num=""><math display="inline"><mrow><msub><mrow><mtext>I</mtext></mrow><mrow><mtext>l</mtext></mrow></msub><mtext> is √</mtext><mover accent="true"><mrow><msub><mrow><mtext>(I</mtext></mrow><mrow><mtext>d</mtext></mrow></msub><msub><mrow><mtext>² + I</mtext></mrow><mrow><mtext>q</mtext></mrow></msub><mtext>²)/3</mtext></mrow><mo>¯</mo></mover></mrow></math><img file="EP0631373A2_D0058.tif" /></maths> .

In a current control circuit 112, the control output of the inverter 113 is obtained from the three-phase current values i_u, i_v and i_w. A rotation speed detector PS of the motor 114 is connected to a position detecting circuit 111 for obtaining the output Q and a speed detecting circuit 115. The speed detecting circuit 115 output ω which is send to the k calculating circuit 101, the Λ calculating circuit 102 and the R_c calculating circuit 108. Further, the exciting current calculating circuit 116 receives an output Λ from the Λ calculating circuit and outputs i_f* to the DC current control circuit 117.

Fig. 17 shows a block diagram of a motor control system of a voltage control type. The voltage control type control system is generally similar to the current control type of Fig. 16 expect that a three-phase voltage command calculating circuit 118 is used instead of the three-phase current calculating circuit 111 so as to output voltage output V_u*, V_v* and V_q* to the inverter 113.

Hereinafter, the calculation of the k-table of the k calculating circuit 101 and the Λ-table of the Λ calculating circuit 102 will be discussed with reference to Figs. 16 and 17.

When the torque command T is commanded and the speed ω is applied to the control system, k and Λ, which satisfy the inputted torque T and speed ω, infinitely exist. Accordingly, it is possible to implement the maximum torque control an the maximum efficiency control by the value k (the ratio between I_d' and I_q') and the magnetic flux Λ.

In case that the maximum efficiency control is implemented, the following equations (32), (26) and (25) are satisfied relative to arbitrary k and Λ.<maths id="math0059" num=""><img file="EP0631373A2_D0059.tif" /></maths>    where during operation I_q' &gt; 0, and k&lt; 0 (T&gt;0); and during regenerating I_q'&lt; 0, and k &gt; 0 (T&gt;0).<maths id="math0060" num=""><img file="EP0631373A2_D0060.tif" /></maths><maths id="math0061" num=""><math display="inline"><mrow><msub><mrow><mtext>I</mtext></mrow><mrow><mtext>d</mtext></mrow></msub><msub><mrow><mtext> = I</mtext></mrow><mrow><mtext>d</mtext></mrow></msub><msub><mrow><mtext>' + I</mtext></mrow><mrow><mtext>cd</mtext></mrow></msub></mrow></math><img file="EP0631373A2_D0061.tif" /></maths><maths id="math0062" num=""><math display="inline"><mrow><msub><mrow><mtext>I</mtext></mrow><mrow><mtext>q</mtext></mrow></msub><msub><mrow><mtext> = I</mtext></mrow><mrow><mtext>q</mtext></mrow></msub><msub><mrow><mtext>' + I</mtext></mrow><mrow><mtext>cq</mtext></mrow></msub><mtext>   (25)</mtext></mrow></math><img file="EP0631373A2_D0062.tif" /></maths>    The terminal voltage V₁ is represented by the following equation (34) from the equation (28), and the primary voltage I₁ is represented by the following equation (35). <maths id="math0063" num=""><math display="inline"><mrow><msub><mrow><mtext>V</mtext></mrow><mrow><mtext>l</mtext></mrow></msub><mtext> = √</mtext><mover accent="true"><mrow><msub><mrow><mtext>V</mtext></mrow><mrow><mtext>d</mtext></mrow></msub><mtext>²</mtext></mrow><mo>¯</mo></mover><mtext></mtext><mover accent="true"><mrow><mtext>-</mtext></mrow><mo>¯</mo></mover><mtext></mtext><mover accent="true"><mrow><msub><mrow><mtext>V</mtext></mrow><mrow><mtext>q</mtext></mrow></msub><mtext>²</mtext></mrow><mo>¯</mo></mover><mtext>   (34)</mtext></mrow></math><img file="EP0631373A2_D0063.tif" /></maths><maths id="math0064" num=""><math display="inline"><mrow><msub><mrow><mtext>   V</mtext></mrow><mrow><mtext>max</mtext></mrow></msub><mtext> = V₁</mtext></mrow></math><img file="EP0631373A2_D0064.tif" /></maths> : where V_max is an outputtable maximum voltage of the inverter. <maths id="math0065" num=""><math display="inline"><mrow><msub><mrow><mtext>I</mtext></mrow><mrow><mtext>l</mtext></mrow></msub><mtext> = √</mtext><mover accent="true"><mrow><msub><mrow><mtext>I</mtext></mrow><mrow><mtext>d</mtext></mrow></msub><mtext>²</mtext></mrow><mo>¯</mo></mover><mtext></mtext><mover accent="true"><mrow><mtext>-</mtext></mrow><mo>¯</mo></mover><mtext></mtext><mover accent="true"><mrow><msub><mrow><mtext>I</mtext></mrow><mrow><mtext>q</mtext></mrow></msub><mtext>²</mtext></mrow><mo>¯</mo></mover><mtext>/√3   (35)</mtext></mrow></math><img file="EP0631373A2_D0065.tif" /></maths>    Further, the various losses of the motor are represented as follows:    Inverter Loss; <maths id="math0066" num=""><math display="inline"><mrow><msub><mrow><mtext>W</mtext></mrow><mrow><mtext>INV</mtext></mrow></msub><mtext> = k₁I₁² + k₂ I₁</mtext></mrow></math><img file="EP0631373A2_D0066.tif" /></maths>    Motor Loss    Primary Loss; <maths id="math0067" num=""><math display="inline"><mrow><msub><mrow><mtext>W</mtext></mrow><mrow><mtext>cu</mtext></mrow></msub><msub><mrow><mtext> = R</mtext></mrow><mrow><mtext>l</mtext></mrow></msub><msub><mrow><mtext> (I</mtext></mrow><mrow><mtext>d</mtext></mrow></msub><msub><mrow><mtext>² + I</mtext></mrow><mrow><mtext>q</mtext></mrow></msub><mtext>²)</mtext></mrow></math><img file="EP0631373A2_D0067.tif" /></maths>    Iron Loss;<maths id="math0068" num=""><img file="EP0631373A2_D0068.tif" /></maths>    where<maths id="math0069" num=""><img file="EP0631373A2_D0069.tif" /></maths>    Machine Loss; <maths id="math0070" num=""><math display="inline"><mrow><msub><mrow><mtext>W</mtext></mrow><mrow><mtext>mc</mtext></mrow></msub><msup><mrow><mtext> = k₀ × n</mtext></mrow><mrow><mtext>1.6</mtext></mrow></msup></mrow></math><img file="EP0631373A2_D0070.tif" /></maths>    Stray Loss; <maths id="math0071" num=""><math display="inline"><mrow><msub><mrow><mtext>W</mtext></mrow><mrow><mtext>st</mtext></mrow></msub><msub><mrow><mtext> = 0.2W</mtext></mrow><mrow><mtext>cu</mtext></mrow></msub></mrow></math><img file="EP0631373A2_D0071.tif" /></maths>    As a result, a total loss can be represented as follows: <maths id="math0072" num=""><math display="inline"><mrow><msub><mrow><mtext>W</mtext></mrow><mrow><mtext>total</mtext></mrow></msub><msub><mrow><mtext> = W</mtext></mrow><mrow><mtext>INV</mtext></mrow></msub><msub><mrow><mtext> + 1.2W</mtext></mrow><mrow><mtext>cu</mtext></mrow></msub><msub><mrow><mtext> + W</mtext></mrow><mrow><mtext>fe</mtext></mrow></msub><msub><mrow><mtext> +W</mtext></mrow><mrow><mtext>mc</mtext></mrow></msub><msub><mrow><mtext> + W</mtext></mrow><mrow><mtext>st</mtext></mrow></msub></mrow></math><img file="EP0631373A2_D0072.tif" /></maths>    Accordingly, an efficiency η is represented as follows: <ul id="ul0002" list-style="none"><li>(1) during driving, <maths id="math0073" num=""><math display="inline"><mrow><msub><mrow><mtext>η = ω</mtext></mrow><mrow><mtext>r</mtext></mrow></msub><msub><mrow><mtext>T × 100/(ω</mtext></mrow><mrow><mtext>r</mtext></mrow></msub><msub><mrow><mtext>T + W</mtext></mrow><mrow><mtext>total</mtext></mrow></msub><mtext>)</mtext></mrow></math><img file="EP0631373A2_D0073.tif" /></maths></li><li>(2) during regenerating, <maths id="math0074" num=""><math display="inline"><mrow><msub><mrow><mtext>η = (ω</mtext></mrow><mrow><mtext>r</mtext></mrow></msub><msub><mrow><mtext>T + W</mtext></mrow><mrow><mtext>total</mtext></mrow></msub><msub><mrow><mtext>) × 100/ω</mtext></mrow><mrow><mtext>r</mtext></mrow></msub><mtext>T</mtext></mrow></math><img file="EP0631373A2_D0074.tif" /></maths></li></ul>    Therefore, by fixing the torque T and the speed ωr and varying k and Λ, it becomes possible to obtain the values of k and Λ by which the efficiency becomes maximum.

The derivation of the k-table and the Λ-table is implemented according to a flow chart shown in Fig. 18.

In a step S₁, the speed ωr is initialized (ωr = 0), in a step S2 the torque T is initialized (T = 0), in a step S3 the value k is set to k_min ( <maths id="math0075" num=""><math display="inline"><mrow><msub><mrow><mtext>k = k</mtext></mrow><mrow><mtext>min</mtext></mrow></msub></mrow></math><img file="EP0631373A2_D0075.tif" /></maths> ) and MAXEFF (maximum efficiency) is set at 0 (MAXEFF = 0), and in a step S4 the magnetic flux Λ is set at O.

In a step S5, on the basis of k and Λ, I_d and I_q are calculated. Following this, in a step S6 the total efficiency η is calculated. In a step S7, it is judged whether the total efficiency η is higher than the MAXEFF or not. When the judgement in the step S7 is "YES", the program jumps to a step S9. When the judgement in the step S7 is "NO", the program proceeds to a step S8 wherein MAXEFF = η and k and Λ are stored. In the step S9 it is judged whether Λ is smaller than Λ_max or not. When the judgement in the step S9 is "YES" (Λ&lt;Λ_max), the program proceeds to a step S13 wherein Λ is incremented by <maths id="math0076" num=""><math display="inline"><mrow><mtext>ΔΛ (Λ = Λ +ΔΛ)</mtext></mrow></math><img file="EP0631373A2_D0076.tif" /></maths> , and then returns to the step S5. When the judgement in the step S9 is "NO" (Λ≧ Λ_max), the program proceeds to a step S10 wherein it is judged whether k is smaller than k_max or not. When the judgement in the step S10 is "YES" (k &lt; k_max), the program proceeds to a step S14 wherein k is incremented by <maths id="math0077" num=""><math display="inline"><mrow><mtext>Δk (k= k+Δk)</mtext></mrow></math><img file="EP0631373A2_D0077.tif" /></maths> and returns to the step S4. When the judgement in the step S10 is "NO", the program proceeds to a step S11 wherein it is judged whether T is smaller than T_max or not. When the judgement in the step S11 is "YES" ( T &lt; T_max), the program proceeds to a step S15 wherein T is incremented by <maths id="math0078" num=""><math display="inline"><mrow><mtext>DT (T=T +ΔT)</mtext></mrow></math><img file="EP0631373A2_D0078.tif" /></maths> , and returns to the step S3. When the judgment in the step S11 is "NO", the program proceeds to a step S12 wherein it is judged whether ωr is smaller than ωr_max or not. When the judgement in the step S12 is "YES", the program proceeds to a step S16 wherein ωr is incremented by <maths id="math0079" num=""><math display="inline"><mrow><mtext>Δωr (ωr = ωr + Δωr)</mtext></mrow></math><img file="EP0631373A2_D0079.tif" /></maths> , and then returns to the step S2. When the judgement in the step S12 is "NO", the program proceeds to END.

In this flow chart, if ΔΛ and Δk is too small, the amount of the points to be calculated becomes too large. Accordingly, the ΔΛ and Δχ are decreased step by step.

For example, when the hybrid excitation type permanent magnet synchronous motor of 45kw is operated by varying the speed ωr and the torque T from 0 to maximum, respectively, ωr and the torque T is varied as follows: <maths id="math0080" num=""><math display="inline"><mrow><mtext>0 ≦ ωr ≦ 9000 (rpm) Δωr = 1000 (rpm)   10 steps</mtext></mrow></math><img file="EP0631373A2_D0080.tif" /></maths><maths id="math0081" num=""><math display="inline"><mrow><mtext>0 ≦ T</mtext></mrow></math><img file="EP0631373A2_D0081.tif" /></maths><maths id="math0082" num=""><math display="inline"><mrow><mtext>≦ 20 (kg m) ΔT = 2 (kg m)   11 steps</mtext></mrow></math><img file="EP0631373A2_D0082.tif" /></maths>    At the 110 (10 × 11) points of the operating points, the efficiency is calculated by varying k and Λ as follows: <maths id="math0083" num=""><math display="inline"><mrow><mtext>- 6≦ k≦   Δk =0.01</mtext></mrow></math><img file="EP0631373A2_D0083.tif" /></maths><maths id="math0084" num=""><math display="inline"><mrow><msub><mrow><mtext>0.3 Λ</mtext></mrow><mrow><mtext>max</mtext></mrow></msub><msub><mrow><mtext> ≦ Λ ≦ Λ</mtext></mrow><mrow><mtext>max   </mtext></mrow></msub><msub><mrow><mtext>ΔΛ =0.01Λ</mtext></mrow><mrow><mtext>max</mtext></mrow></msub><mtext>.</mtext></mrow></math><img file="EP0631373A2_D0084.tif" /></maths>

The varying of k is implemented as shown in Fig. 20, that is, first Δk = 1 (10 points), next Δk = 0.25 (7 points), third Δk =0.05 (9 points), and then Δk =0.01 (9 points). Accordingly, Λ may be varied from o.3 Λ_max to Λ_max relative to the 35 points of k.

The simulation data of the coefficient k and ΔΛ of the hybrid excitation type permanent magnet synchronous motor is shown as follows:    45kw - 8P-3000rpm    R₁ = 0.01582 (Ω)    L_d = 0.299 (mH) L_q = 0.192 (mH)    L₀ = 0.186    Inverter loss coefficient    k₁ = 0.0156/2 k₂= 9.02/2 The k-table and Λ-table were obtained as shown in the following Tables 1 and 2. <tables id="tabl0001" num="0001"><img file="EP0631373A2_D0085.tif" /></tables><tables id="tabl0002" num="0002"><img file="EP0631373A2_D0086.tif" /></tables>

As a result, the efficiency of the motor was obtained as shown in a map of Fig. 20, and the MAXEFF was 95% (MAXEFF=95%).

Although the maximum efficiency operation as mentioned above is implemented by making the k-table and the Λ-table, such a table making complicates the calculation thereof in the practical control system. Accordingly, in order to shorten the time for the calculation, the exciting current I_d, the torque current I_q and the magnetic flux Λ, by which the maximum efficiency is obtained relative to the arbitrary operating points, are made into a table. In this case, the algorithm for making tables of I_d and I_q is the same as that of the k-table and the Λ-table. That is, I_d and I_q obtained from k and Λ are made into a table, and the Λ-table which is the same as that of the previous case is used. Since the command current is directly made into table, it is necessary to make the table by using small steps.

For example, assuming that the parameter of the motor is the same as that of the previous case, the rotation speed ωr and the torque T are varied as follows: <maths id="math0085" num=""><math display="inline"><mrow><mtext>0 ≦ ωr ≦ 10400 (rpm) Δωr = 81.25 rpm   (128 steps)</mtext></mrow></math><img file="EP0631373A2_D0087.tif" /></maths><maths id="math0086" num=""><math display="inline"><mrow><mtext>0 ≦ T 165 (N, m) ΔT = 10.33 N, m   (16 steps)</mtext></mrow></math><img file="EP0631373A2_D0088.tif" /></maths>    Since such steps are calculated by a CPU, the number of steps is 2^X ( X is natural number). The k-table and Λ-table are obtained within the following range. <maths id="math0087" num=""><math display="inline"><mrow><mtext>-6 ≦ k ≦ 3</mtext></mrow></math><img file="EP0631373A2_D0089.tif" /></maths><maths id="math0088" num=""><math display="inline"><mrow><msub><mrow><mtext>0.3 Λ</mtext></mrow><mrow><mtext>max</mtext></mrow></msub><msub><mrow><mtext> ≦ Λ ≦ Λ</mtext></mrow><mrow><mtext>max</mtext></mrow></msub></mrow></math><img file="EP0631373A2_D0090.tif" /></maths>    Accordingly, the number of I_d, I_q, Λ data becomes 2193 ((128 +1) × (16 + 1) = 2193). The data is interpolated.

Fig. 21 shows a block diagram of the maximum efficiency control system which is practically used. The commands I_d*, I_q* and Λ* are outputted from the I_d-table, I_q-table and Λ-table. V_d* and V_q* are obtained therefrom and pass through a d-axis-q-axis coordinate transformation section and a 2φ - 3φ transformation section, and are used for the control of a PWM inverter, that is, the armature voltage. The command Λ* is transformed into V_f* and is used for the field magnet control at a PWM chopper.

A counterplan for disturbance is included in the control system shown in Fig. 21. That is, an enclosed area by a broken line in Fig. 21, which is disposed between the command I_d*, I_q* and the coordinate transformation (α-β transformation), includes an integral element (I element), K_Id/S, K_Iq/S, an adding proportion element (P element), K_pd, K_pq and adding elements ωL_qI_q, ω(L_dI_d +Λ). The integral element and the proportion element are provided for compensating a transient term and a normal term in the previous voltage equation and integral and adding in a feedback system. The adding element indicates an interference term and a compensation in a feed-forward system.

The disturbance in the equivalent circuit and the voltage equation will be discussed hereinafter. That is, I_q, I_d and Λ are previously made into tables as shown in Fig. 21 for obtaining maximum efficiency of k-table and Λ-table, and the current control is independently implemented on the respective d-axis and q-axis, in actual case. Accordingly, the disturbance is generated in the respective axes and degrades the performance of the control. In Fig. 22, taking into consideration with the iron loss of a typical disturbance, a voltage equation for an equivalent circuit of Fig. 22 including a transient condition is obtained as follows:<maths id="math0089" num=""><img file="EP0631373A2_D0091.tif" /></maths>    where R₁ is an armature resistance; R_c is a equivalent iron-loss resistance; L_d is a d-axis inductance; L_q is a q-axis inductance; Λ is an interlinkage magnetic flux of the armature coil; w is an electric source angle frequency; V_d is a d-axis component of the armature voltage; V_q is a q-axis component of the armature voltage; I_d is a d-axis component of the armature current; I_q is a q-axis component of the armature current; P is a deferential operator; and k₁ is <maths id="math0090" num=""><math display="inline"><mrow><msub><mrow><mtext>1+R₁/R</mtext></mrow><mrow><mtext>c</mtext></mrow></msub></mrow></math><img file="EP0631373A2_D0092.tif" /></maths> .

Furthermore, I_d, I_q, I_cd and I_cq in Figs. 22A and 22B are represented by the equations (25) and (26), and <maths id="math0091" num=""><math display="inline"><mrow><msub><mrow><mtext>Λ = Λm + M</mtext></mrow><mrow><mtext>f</mtext></mrow></msub><msub><mrow><mtext>I</mtext></mrow><mrow><mtext>f</mtext></mrow></msub></mrow></math><img file="EP0631373A2_D0093.tif" /></maths> . As a result, from the voltage equation (36) and the equivalent circuits of Figs. 22A and 22B, a feedback of I_d' and I_q' may be implemented for the control of V_d and V_q. Although in the control system of Fig. 21 a feedback of I_d and I_q is implemented from the output V_u and V_w of the PWM inverter, in order to obtain I_d' and I_q' from which iron-loss current has already eliminated, it is necessary to implement further transformation. That is, since only I_d and I_q are obtained from the coordinate transformation, it is necessary to implement the transformation represented by the following equation.<maths id="math0092" num=""><img file="EP0631373A2_D0094.tif" /></maths>    Fig. 23 shows a feedback system of a modulation of a feedback system shown by a hatching in Fig. 21 upon taking into consideration with the equation (37). That is, d'-q' transformation represented by the equation (37) is implemented after the 3φ-2φ transformation and the d-q transformation. An area enclosed by a broken line shows a non-interference IP control system which eliminates iron-loss.

In the voltage equation (36) and the equivalent circuits of Figs. 22A and 22B, a disturbance, which can not be controlled independently in d-axis and q-axis, that is, an interference term, is K_r(Λ_m + M_fI_f) at the d-axis which is of an induced voltage due to the q-axis inductance and the d-axis inductance. and ωkL_dI_q' and ωKr(Λ_m + M_fI_f) at the q-axis which is of an induced voltage generated by the change of the magnetic flux by the DC field magnet current of q-axis. Therefore, a feed-forward compensation is implemented relative to d-axis and q-axis as a control for eliminating the disturbance. That is, a feed-forward compensation of ω(L_dI_q + Λ) shown in Fig. 21 is implemented.

Hereinafter, a current control section for a feed-forward compensation and a feed-forward system will be discussed.

Figs. 24 and 25 show an example of a non-interference processing by a feed-forward compensation, which is of a proportion-integral (PI) control system. A current control section between the table command and α-β transformation includes a proportion integral element <maths id="math0093" num=""><math display="inline"><mrow><msub><mrow><mtext>K</mtext></mrow><mrow><mtext>pd</mtext></mrow></msub><msub><mrow><mtext>(1+1/ST</mtext></mrow><mrow><mtext>Id</mtext></mrow></msub><msub><mrow><mtext>) · K</mtext></mrow><mrow><mtext>pd</mtext></mrow></msub><msub><mrow><mtext>(1 + 1/ST</mtext></mrow><mrow><mtext>Iq</mtext></mrow></msub><mtext>)</mtext></mrow></math><img file="EP0631373A2_D0095.tif" /></maths> , and feed-forward compensation elements ωK_rL_d, ωKρK_q and ωK_r(Λ_m+M_fI_f). Command values I_d* and I_q*, and feedback values I_d' and I_q' are inputted thereto and the command values V_d* and V_q* are outputted therefrom upon eliminating the disturbance. Accordingly, with this non-interference processing, d-axis and q-axis are independently separated as shown in Fig. 25, and therefore the performance of the control is improved.

In this case, if the proportion gain of d-axis is <maths id="math0094" num=""><math display="inline"><mrow><msub><mrow><mtext>K</mtext></mrow><mrow><mtext>pd</mtext></mrow></msub><msub><mrow><mtext>=ωcL</mtext></mrow><mrow><mtext>pd</mtext></mrow></msub><mtext>/2π</mtext></mrow></math><img file="EP0631373A2_D0096.tif" /></maths> and an integral time constant is <maths id="math0095" num=""><math display="inline"><mrow><msub><mrow><mtext>T</mtext></mrow><mrow><mtext>Id</mtext></mrow></msub><msub><mrow><mtext>=L</mtext></mrow><mrow><mtext>d</mtext></mrow></msub><mtext>/R</mtext></mrow></math><img file="EP0631373A2_D0097.tif" /></maths> , a transfer function for I_d' and I_d*' are represented by the following equation (38).<maths id="math0096" num=""><img file="EP0631373A2_D0098.tif" /></maths>    Accordingly, a response characteristic is determined by the characteristic angle frequency ωc. Similarly, the same relationship as to the q-axis is obtained by assuming that <maths id="math0097" num=""><math display="inline"><mrow><msub><mrow><mtext>K</mtext></mrow><mrow><mtext>pd</mtext></mrow></msub><msub><mrow><mtext>=ωcL</mtext></mrow><mrow><mtext>pd</mtext></mrow></msub><mtext>/2π</mtext></mrow></math><img file="EP0631373A2_D0099.tif" /></maths> and <maths id="math0098" num=""><math display="inline"><mrow><msub><mrow><mtext>T</mtext></mrow><mrow><mtext>Id</mtext></mrow></msub><msub><mrow><mtext>=L</mtext></mrow><mrow><mtext>d</mtext></mrow></msub><mtext>/R₁.</mtext></mrow></math><img file="EP0631373A2_D0100.tif" /></maths>

Figs. 26 and 27 show another embodiment of the non-interference processing of a IP control system. The structure of the IP control system of Fig. 26 is the same as that of the current control section shown in Fig. 21. The IP control system incudes a proportion elements K_pd and K_pq and feed-forward compensations ωk_rL_d, ωK_rL_q and ωK_r(Λ_m +M_fI_f) in addition to integral elements K_Id/S and K_Iq/S. Accordingly, the control system of Fig. 26 also suppresses the disturbance. By the non-interference processing shown in Fig. 27, it becomes possible to separate d-axis and q-axis, and therefore the control performance is improved.

In this case, if the proportion gain is <maths id="math0099" num=""><math display="inline"><mrow><msub><mrow><mtext>K</mtext></mrow><mrow><mtext>pd</mtext></mrow></msub><msub><mrow><mtext> = 2ωcξL</mtext></mrow><mrow><mtext>d</mtext></mrow></msub><mtext> - R₁</mtext></mrow></math><img file="EP0631373A2_D0101.tif" /></maths> , and <maths id="math0100" num=""><math display="inline"><mrow><msub><mrow><mtext>K</mtext></mrow><mrow><mtext>Id</mtext></mrow></msub><mtext> = ωc²L₁</mtext></mrow></math><img file="EP0631373A2_D0102.tif" /></maths> at d-axis, a transfer function between I_d' and I_d*' is represented by the following equations (39) and (40).<maths id="math0101" num=""><img file="EP0631373A2_D0103.tif" /></maths>    Since this equation (39) is the same as a transfer function of a general secondary filter, it is represented as follows:<maths id="math0102" num=""><img file="EP0631373A2_D0104.tif" /></maths>    where ωc is a characteristic angle frequency; and ξ is a damping coefficient.

Accordingly, a current response is determined by the characteristic angle frequency ωc and the damping coefficient ξ. Similarly, the same relationship as to the q-axis is obtained by the similar assumption.

Figs. 24 and 25 show the PI control system which is implemented by using "primary delay" and Fig. 26 and 27 show IP control system which implemented by using "secondary delay". That is, in PI control the difference between the command value and the feedback value is inputted to both the integral term and the proportion term, and V_d* and V_q* are outputted so as to decrease the difference to 0. By the implementation of non-interference processing, it becomes possible to control the response of the motor current (I_d, I_q) as "primary delay", as shown in Fig. 28A and by changing the gains (K_pd, K_pq) as shown in Fig. 28B, it can be treated as an arbitrary primary delay.

On the other hand, In IP control the difference between the command value and the feed back value is inputted only to the integral term, and the feed back value is inputted to the proportion term. Further V_d* and V_q* are outputted so as to decrease the difference to O. By the non-interference processing, it becomes possible to control the motor current by the secondary delay as shown in Fig. 28C, and the secondary delay is arbitrarily determined by changing the gains (K_pd, K_Id, K_pq and K_Iq). Since the number of the gains to be changed is larger than that of the gain of the PI control, the degree of freedom of the response wave is increased. Further, it is possible to increase the response speed of the IP control rather than that of the PI control by the proper selection of the gain constituted by the same sample frequency.

When large quantity of noises are occurred in the feedback signals I_d and I_q in the control systems of Figs. 24 and 26, the current values I_d', I_q' and I_f of the interference term may be exchanged by the command I_d*, I_q* and I_f*. That is, ωK_rL_dI_d' to ωK_rL_dI_d'*, ωK_rL_qI_q' to ωK_rL_qI_q'*, and ωK_r(Λ_m + M_fI_f) to ωK_r(Λ_m+M_fI_f*). Where I_f* is a DC exiting current command.

In the above explanation, the equivalent circuit including a transient condition taking account the iron loss is applied as promise. Figs. 29A and 29B show equivalent circuits including the transient condition in case that the iron loss is neglected and corresponds to a case that the resistance R_c in the circuits of Figs. 22A and 22B is set infinity (∞). The voltage equation in this case is represented by the following equation (41).<maths id="math0103" num=""><img file="EP0631373A2_D0105.tif" /></maths>    When <maths id="math0104" num=""><math display="inline"><mrow><msub><mrow><mtext>R</mtext></mrow><mrow><mtext>c</mtext></mrow></msub><mtext>= ∞</mtext></mrow></math><img file="EP0631373A2_D0106.tif" /></maths> is substituted into the equation (26), it becomes that <maths id="math0105" num=""><math display="inline"><mrow><msub><mrow><mtext>I</mtext></mrow><mrow><mtext>d</mtext></mrow></msub><msub><mrow><mtext> = I</mtext></mrow><mrow><mtext>d</mtext></mrow></msub><mtext>'</mtext></mrow></math><img file="EP0631373A2_D0107.tif" /></maths> and the current control section shown in Fig. 21 controls V_d and V_q by the feedback of I_d and I_q. Accordingly, it becomes unnecessary to transform I_d' and I_q' in Fig. 23. Since the interference term in the equation (41) of Fig. 29 becomes -ωL_qI_q at the d-axis side and <maths id="math0106" num=""><math display="inline"><mrow><msub><mrow><mtext>ωL</mtext></mrow><mrow><mtext>d</mtext></mrow></msub><msub><mrow><mtext>I</mtext></mrow><mrow><mtext>d</mtext></mrow></msub><msub><mrow><mtext> + ω(Λ</mtext></mrow><mrow><mtext>m</mtext></mrow></msub><msub><mrow><mtext> - M</mtext></mrow><mrow><mtext>f</mtext></mrow></msub><msub><mrow><mtext>I</mtext></mrow><mrow><mtext>f</mtext></mrow></msub><mtext>)</mtext></mrow></math><img file="EP0631373A2_D0108.tif" /></maths> at the q-axis side, the value may be compensated by the feed-forward compensation. Furthermore, as similar to Figs. 25 and 27, the stationary term (R₁ + L_pd)I_d and the transient term (R₁ + L_qp)I_q are compensated by the PI control or IP control.

Fig. 30 shows a PI control system by which non-interference processing is implemented. Fig. 31 shows a IP control system. The difference between Fig. 24 and Fig. 26 is that a coefficient relating to the iron loss is omitted. The embodiments in Figs. 30 and 31 are equivalent with embodiments that I_d'*, I_q'*, I_d' and I_q' are replaced with I_d*, I_q*, I_d and I_q, respectively in the embodiments of Figs. 25 and 27. Further, the response characteristics of Figs. 30 and 31 are equivalent with the equations (38) and (39). Furthermore, it is a matter of fact that I_d, I_q and I_f of the interference term may be replaced with I_d*, I_q* and I_f*, respectively.

Fig. 32 shows a PI control system which uses a voltage model and to which non-interference is applied. Furthermore, stationary terms R₁I_d' and R₁I_q' are compensated by a feed-forward on the d-axis and the q-axis. That is, the PI control is an equation shown by the equation (36), and is that the error of the voltage command of either the transient term of each motor constant (L_d, L_q, R₁ and the like) is compensated thereby. Also, instead of the compensation of R₁ term and the interference term with a current feedback value (I_d', I_q' and I_f) as mentioned above, it may be compensated by the current command (I_d*, I_q* and I_f*).

Although a control method of Fig. 33 has been shown and described so as to be controlled on the basis of the voltage equation (36) taking into consideration with the iron loss, if a case of Fig. 30 is applied thereto, the control method of Fig. 33 may be replaced with that on the basis of the voltage equation (41), where R₁I_d' and R₁I_q' of R₁ term are changed into R₁I_d and R₁I_q, respectively. Thus, by treating the interference term by non-interference so as to independently control on the d-axis and the q-axis. That is, by applying the following interference term to the current control system by feed-forward, it becomes possible to independently control on the d-axis and the q-axis. and therefore the control response is improved.    d-axis side ωK_rL_qI_q'    In case that iron loss is taken account,    q-axis side ωK_eL_dI_d' + ωK_r(Λ_m +M_fI_f)    d-axis side ωL_qI_q    In case that iron loss is neglected,    q-axis side ωL_dI_d + ω(Λ_m +M_fI_f)    In case that noise is generated in the feedback signal (I_d, I_q and I_f), by calculating the interference term with the command value (I_d*, I_q* and I_f*, or I_d'*, I_q'* and I_f'), it becomes possible to accomplish non-interference processing. Furthermore, by constructing the voltage model on the basis of the voltage equation of the hybrid excitation type permanent magnet synchronous motor, the control performance is improved.

Although the explanation hereinbefore is relating to the hybrid excitation type permanent magnet synchronous motor, if the magnet flux Λ is kept constant, the above-mentioned control systems may be applied to a PM motor. As shown in Fig. 34, by eliminating a circuit relating to Λ from the control system of Fig. 16, a control system of the PM motor can be constituted.

As a result of the simulation of the maximum efficiency control of the PM motor, the constants of the embedded type PM motor are decided as follows and the k-table is set as shown in Table 3 and Λ-table is set at 1 in all region.    45kw-8P-3000rpm    R₁ = 0.0121 Ω, R_co = 28.8Ω    L_q = 0.450mH, L_d = 0. 161mH    k₀ = 2.3337 × 10⁻⁴, L₀ = 0.0893    k₁ = 0.0156/2, k₂ = 9.02/2 Therefore, the efficiency map of the embedded type PM motor was obtained as shown in Fig. 35, and MAXEFF = 95.28% was obtained. <tables id="tabl0003" num="0003"><img file="EP0631373A2_D0109.tif" /></tables>

In a maximum torque control, when the torque T and the speed ωr are applied, the primary current I₁ becomes minimum at the maximum torque. Accordingly, k and Λ or I_d, I_q and Λ may be obtained. That is, I_d, I_q of the equation (25) is obtained from the equation (32), and I₁ is obtained from the equation (35). Therefore, k and Λ, or I_d, i_q and Λ by which I₁ becomes minimum relative to arbitrary ωr and T may be obtained by repeat calculations.

With the thus arranged embodiments, by proper determination of k and Λ by k-table and Λ-table or the determination of I_d, I_q and Λ by I_d, I_q and Λ table, the operable area of the hybrid excitation type permanent magnet synchronous motor, d.c. machine and a PM motor is extended. Furthermore, it becomes possible to implement the maximum efficiency control and the maximum torque control. Additionally, it becomes possible to implement the non-interference control by removing the disturbance by the feed-forward control of the interference term.