WO2019142301A1 - 駆動装置、流体利用装置及び空気調和機 - Google Patents

駆動装置、流体利用装置及び空気調和機 Download PDF

Info

Publication number
WO2019142301A1
WO2019142301A1 PCT/JP2018/001458 JP2018001458W WO2019142301A1 WO 2019142301 A1 WO2019142301 A1 WO 2019142301A1 JP 2018001458 W JP2018001458 W JP 2018001458W WO 2019142301 A1 WO2019142301 A1 WO 2019142301A1
Authority
WO
WIPO (PCT)
Prior art keywords
current
synchronous motor
sub
pulsation
torque
Prior art date
Application number
PCT/JP2018/001458
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 KR1020207017409A priority Critical patent/KR102412170B1/ko
Priority to EP18901287.5A priority patent/EP3742598B1/en
Priority to JP2019565637A priority patent/JP6877592B2/ja
Priority to US16/960,113 priority patent/US11133762B2/en
Priority to PCT/JP2018/001458 priority patent/WO2019142301A1/ja
Priority to CN201880085523.7A priority patent/CN111630769B/zh
Publication of WO2019142301A1 publication Critical patent/WO2019142301A1/ja
Priority to JP2021072758A priority patent/JP7113939B2/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
    • H02P5/00Arrangements specially adapted for regulating or controlling the speed or torque of two or more electric motors
    • H02P5/46Arrangements specially adapted for regulating or controlling the speed or torque of two or more electric motors for speed regulation of two or more dynamo-electric motors in relation to one another
    • H02P5/50Arrangements specially adapted for regulating or controlling the speed or torque of two or more electric motors for speed regulation of two or more dynamo-electric motors in relation to one another by comparing electrical values representing the speeds
    • 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/14Estimation or adaptation of machine parameters, e.g. flux, current or voltage
    • H02P21/18Estimation of position or speed
    • 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/22Current control, e.g. using a current control loop
    • 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
    • H02P5/00Arrangements specially adapted for regulating or controlling the speed or torque of two or more electric motors
    • H02P5/46Arrangements specially adapted for regulating or controlling the speed or torque of two or more electric motors for speed regulation of two or more dynamo-electric motors in relation to one another

Definitions

  • the present invention relates to a drive device for driving two synchronous motors connected in parallel to one power converter, a fluid utilization device including the drive device, and an air conditioner including the fluid utilization device.
  • the synchronous motor can not generate torque properly and stop operation or stop operation without applying an appropriate voltage according to the rotational position using information on the rotor position, which is the rotational position of the rotor.
  • one of two synchronous motors connected to one power converter is a main synchronous motor, the other is a sub synchronous motor, and the driving device includes two synchronous motors;
  • a method of controlling using the rotational speed of a stand synchronous motor is disclosed.
  • the driving device vector-controls the main synchronous motor, and the d-axis current of the main synchronous motor using the difference in rotational speed and the difference in rotational position between the two synchronous motors.
  • Vector control is a control method of decomposing a current flowing through a synchronous motor into a current component generating a torque and a current component generating a magnetic flux, and controlling each current component independently.
  • Patent Document 1 describes that the method described in Patent Document 1 can also be applied to position sensorless control that does not use a sensor such as a speed detection unit or a position detection unit. Further, Non-Patent Document 1 and Non-Patent Document 2 disclose the results of a demonstration test of position sensorless control using the technology of Patent Document 1. Position sensorless control generally has a problem in operational stability in a low speed range. The reason is that the speed electromotive force of the motor decreases in the low rotation range, and the influence of the output voltage error of the power converter becomes relatively large. Even in the case of driving one synchronous motor with one power converter, the operation stability at low speed becomes a problem, so when driving two synchronous motors in parallel with one power converter, the low speed It becomes increasingly difficult to ensure operational stability. However, Patent Document 1, Non-Patent Document 1 and Non-Patent Document 2 do not specifically describe the driving method in the low speed range.
  • the present invention has been made in view of the above, and suppresses rotation instability in a low speed range even when two synchronous motors are driven with position sensorless using one power converter.
  • the object is to obtain a driving device that can
  • the drive device of the present invention comprises: a power converter for supplying power to a first synchronous motor and a second synchronous motor connected in parallel; A first current detector for detecting a first current flowing to the motor and a second current detector for detecting a second current flowing to the second synchronous motor.
  • the drive device estimates a first magnetic pole position of a rotor of the first synchronous motor using a voltage command for driving the first synchronous motor and a first current.
  • a second magnetic pole position estimating unit that estimates a second magnetic pole position of a rotor of the second synchronous motor using the voltage command and the second current.
  • the drive device includes a control unit that outputs a voltage command using the torque current command, the flux current command, the first current, and the first magnetic pole position, and torque current pulsation included in the torque current flowing to the second synchronous motor.
  • a pulsating component extraction unit for extracting at least one pulsating component of the component and the active power pulsating component included in the active power consumed by the second synchronous motor.
  • the driving device uses a first magnetic pole position and a second magnetic pole position to obtain an angular difference which is a difference between the magnetic pole positions of the rotors respectively possessed by the first synchronous motor and the second synchronous motor.
  • a magnetic flux current command determination unit that determines a magnetic flux current command using a subtractor, a pulsation component extracted by the pulsation component extraction unit, and an angle difference is characterized.
  • FIG. 7 is a diagram showing a first configuration example of a sub-side torque current pulsation component extraction unit shown in FIG. 1;
  • the figure which shows the 3rd structural example of the sub side torque electric current pulsating component extraction part shown in FIG. A first diagram for explaining the behavior of torque change of the sub-side synchronous motor shown in FIG.
  • the 3rd figure for demonstrating the behavior of the torque change of the sub side synchronous motor shown in FIG.
  • the 4th figure for demonstrating the behavior of the torque change of the sub side synchronous motor shown in FIG. The 5th figure for demonstrating the behavior of the torque change of the sub side synchronous motor shown in FIG.
  • the 6th figure for demonstrating the behavior of the torque change of the sub side synchronous motor shown in FIG. A diagram showing the main d-axis current shown in FIGS. 6, 7, 9 and 10, the sign of the angle difference, and the state of the torque of the sub-side synchronous motor in association with each other.
  • the drive device In the drive device according to the first embodiment, it is a true value of the rotational speed of the sub-side synchronous motor when the two synchronous motors are parallelly driven at low speed, virtually assuming that the output voltage error of the power converter is zero.
  • Diagram showing true velocity, d-axis current and q-axis current In the drive device according to the first embodiment, the output voltage error equivalent to that of the actual machine is entered, and the true speed value of the rotational speed of the sub-side synchronous motor when two synchronous motors are parallelly driven at low speed Diagram showing d-axis current and q-axis current In the drive device concerning Embodiment 1, the output voltage error of a power converter is virtually considered to be zero, and the figure showing the speed estimation waveform at the time of making two synchronous motors be parallelly driven at low speed.
  • FIG. 19 is a first diagram for illustrating code determination processing by the code determination unit shown in FIG. 19; The 2nd figure for demonstrating the code
  • FIG. 1 The 3rd figure for demonstrating the code
  • the block diagram of the fluid utilization apparatus which concerns on Embodiment 4 of this invention The block diagram of the air conditioner concerning Embodiment 5 of this invention
  • FIG. 1 is a diagram showing a configuration of a drive device according to Embodiment 1 of the present invention.
  • the synchronous motor utilizes a permanent magnet field type synchronous motor in which a permanent magnet is provided on the rotor, a winding field type synchronous motor in which field windings are wound on the rotor, and saliency of the rotor. It is roughly divided into a reluctance type synchronous motor to obtain a rotational torque. Of the types of these synchronous motors, two synchronous motors of the same type, for example, permanent magnet field type synchronous motors, are connected in parallel to the drive device 100 according to the first embodiment.
  • one of the two synchronous motors is referred to as the main synchronous motor 1a, and the other is referred to as the sub synchronous motor 1b.
  • the main synchronous motor 1a is a first synchronous motor
  • the sub synchronous motor 1b is a second synchronous motor.
  • a three-phase permanent magnet field type synchronous motor is used, but the motor constants of the two synchronous motors may be approximately the same, and a synchronous motor other than the permanent magnet field type Alternatively, a synchronous motor having a number of phases other than three, such as two or five, may be used.
  • Drive device 100 includes a power converter 2 for supplying power to main side synchronous motor 1a and sub side synchronous motor 1b connected in parallel, and a current detection unit 4a for detecting a first current flowing through main side synchronous motor 1a. And a current detection unit 4b that detects a second current flowing through the sub-side synchronous motor 1b, and a magnetic pole position estimation unit 5a that is a first magnetic pole position estimation unit. Further, drive device 100 includes a magnetic pole position estimation unit 5b which is a second magnetic pole position estimation unit, a current control unit 6 which is a control unit which outputs a voltage command, a pulsation component extraction unit 70, and a first subtractor. And a magnetic flux current command determination unit 9.
  • the pulsation component extraction unit 70 includes the sub torque current pulsation component extraction unit 7. Below, the sub-side torque current pulsation component extraction unit 7 may be simply referred to as "pulsating component extraction unit 7".
  • the power converter 2 converts DC power supplied from the DC voltage source 3 into AC power and outputs the AC power to the main synchronous motor 1a and the sub synchronous motor 1b.
  • a voltage source inverter is used for power converter 2.
  • the voltage type inverter is a device that switches the DC voltage supplied from the DC voltage source 3 and converts it into an AC voltage.
  • Power converter 2 is not limited to a voltage type inverter as long as it can output AC power for driving main side synchronous motor 1a and sub side synchronous motor 1b, and current type inverter, AC power amplitude is not limited. And, it may be a circuit such as a matrix converter that converts into alternating current power with different frequencies, and a multilevel converter in which the outputs of a plurality of converters are connected in series or in parallel.
  • a current detection unit 4a which is a first current detector, detects a phase current flowing from the power converter 2 to the main synchronous motor 1a, and outputs current information indicating the value of the detected phase current.
  • a current detection unit 4b which is a second current detector, detects a phase current flowing from the power converter 2 to the sub-side synchronous motor 1b, and outputs current information indicating the value of the detected phase current.
  • the current detection units 4a and 4b may be current sensors using current transformers for instruments called CT (Current Transformer), or may be current sensors using shunt resistors.
  • CT Current Transformer
  • the current detection units 4a and 4b may be a combination of these.
  • current is detected by current detection units 4a and 4b provided near the synchronous motor.
  • the phase current flowing through the synchronous motor is directly detected in the example shown in FIG. 1, the current detection method is not limited to the example of directly detecting the current flowing through the synchronous motor according to Kirchhoff's current law.
  • a single-shunt current detection method using a shunt resistor provided on the negative DC bus of power converter 2 a lower-arm shunt current detection method using a shunt resistor connected in series with the lower arm of power converter 2, etc.
  • the phase current flowing in the synchronous motor may be detected using
  • the lower arm shunt current detection method is also referred to as a three shunt current detection method because shunt resistances connected in series to each of the three lower arms are used.
  • the magnetic pole position estimation unit 5a uses the first current detected by the current detection unit 4a, that is, the phase current flowing through the main synchronous motor 1a, and the voltage command output from the current control unit 6 to perform main side synchronization. The first magnetic pole position of the motor 1a is estimated.
  • the magnetic pole position estimation unit 5b uses the second current detected by the current detection unit 4b, that is, the phase current flowing to the sub synchronous motor 1b, and the voltage command output from the current control unit 6, The second magnetic pole position of the motor 1b is estimated.
  • the magnetic pole position is determined using information indicating the speed electromotive force of the synchronous motor in the medium to high speed region of the rotational speed range of the rotor included in the synchronous motor. Is common.
  • the speed electromotive force is an induced power generated inside the synchronous motor as the rotor rotates, and is proportional to a field generated between the rotor and the stator included in the synchronous motor, and a rotational speed of the rotor. Details of the method of estimating the magnetic pole position will be described later.
  • the current control unit 6 sets the direction of the magnetic flux by the permanent magnet of the rotor included in the main synchronous motor 1a as d axis, and the axis orthogonal to the d axis as q axis
  • the vector controller coordinate-transforms the current detected by the current detector 4a into a current command value in the dq coordinate system.
  • current control is performed on dq coordinates based on the magnetic poles of the rotor.
  • the AC amount becomes a DC amount and control becomes easy.
  • the q-axis current is proportional to the magnet torque of the synchronous motor, the q-axis is called a torque axis, and the q-axis current is called a torque current.
  • the d-axis current produces a change in the magnetic flux generated at the stator with respect to the q-axis current, and changes the amplitude of the output voltage of the synchronous motor. It is called an excitation current or the like.
  • the types of synchronous motors include surface magnet type synchronous AC motors provided with permanent magnets on the outer peripheral surface of the rotor core, permanent magnet embedded motors with permanent magnets embedded inside the rotor core, etc. .
  • the q-axis current is often referred to as a torque current.
  • the estimated value of the magnetic pole position calculated by the magnetic pole position estimation unit 5a is used for coordinate conversion.
  • the current control unit 6 may use a polar coordinate system such as an ⁇ stator coordinate system or a ⁇ coordinate system in addition to the dq coordinate system in vector control.
  • direct torque control may be adopted for the current control unit 6 instead of vector control.
  • DTC direct torque control
  • the torque current and the magnetic flux current can be calculated more strictly.
  • This coordinate system is often referred to as an ft coordinate system, an nt coordinate system, etc., but since it is known it will not be described in detail.
  • the q-axis current may be referred to as a torque current
  • the d-axis current may be referred to as a magnetic flux current.
  • magnet torque is not generated in principle This is not the case when using a reluctance synchronous motor.
  • the current control unit 6 is controlled such that the torque current flowing through the main synchronous motor 1a matches the value of the torque current command, and the magnetic flux current flowing through the main synchronous motor 1a matches the value of the magnetic flux current command To be controlled.
  • the current control unit 6 is generally configured by a proportional integral controller and a non-interference controller.
  • the torque current command may be calculated as a result of the speed control in the magnetic flux current command determination unit 9, or may be input from a host controller. Details of the magnetic flux current command will be described later.
  • the synchronous motor when the synchronous motor is driven in an open loop, the synchronous motor may self-oscillate at the natural angular frequency ⁇ n and control may become unstable.
  • the natural angular frequency ⁇ n is expressed by the following approximate expression (1).
  • P m is the number of pole pairs
  • ⁇ a is the armature flux linkage number
  • L a is the armature inductance
  • J is the moment of inertia.
  • the electromechanical coupling vibration may be referred to as an electric spring resonance
  • the natural angular frequency ⁇ n represented by the above equation (1) is also referred to as an electric spring resonance angular frequency.
  • the stabilization compensator is added to the technique disclosed in the above-mentioned reference 1 in order to suppress the electrical spring resonance, the same stabilization compensation is required also in the drive device 100. For that purpose, it is necessary to check how much the torque current flowing through the sub-side synchronous motor 1b shown in FIG. 1 vibrates due to the electrical spring resonance.
  • the speed difference which is the difference between the rotational speeds of the rotors of the main synchronous motor and the sub synchronous motor is determined, and the speed difference is stabilized by using this speed difference. Compensation is being done. Since the main side synchronous motor is controlled stably by this, it can be said that the technique disclosed by patent document 1 calculated
  • the component by acceleration / deceleration torque and the component by load torque are superimposed on the torque current which flows into the sub side synchronous motor 1b.
  • the acceleration / deceleration torque is an inertial torque associated with acceleration / deceleration of the synchronous motor.
  • the load torque is a torque obtained by subtracting the acceleration / deceleration torque and the loss such as friction from the output torque.
  • the pulsation component extraction unit 7 shown in FIG. 1 extracts a pulsation component in the vicinity of the electrical machine spring resonance angular frequency included in the torque current of the sub-side synchronous motor 1b.
  • the pulsation component extraction unit 7 has two types, a method using a high pass filter and a method using a band pass filter, which will be described in order.
  • FIG. 2 is a diagram showing a first configuration example of the sub-side torque current pulsation component extraction unit shown in FIG.
  • FIG. 2 shows a configuration example of the sub-side torque current pulsation component extraction unit 7A using a first-order high-pass filter, and its transfer function is expressed by the following equation (2).
  • s is the Laplace transform operator
  • ⁇ c is the cutoff angular frequency.
  • n is an integer of 2 or more.
  • the cutoff angular frequency ⁇ c is preferably set to 1/3 or less of the electric motor resonance frequency, for example, 1/5 to 1/20 of the electric spring resonance frequency.
  • FIG. 3 is a diagram showing a second configuration example of the sub-side torque current pulsation component extraction unit shown in FIG.
  • FIG. 3 shows a configuration example of the sub-side torque current pulsation component extraction unit 7B using a second-order band pass filter, and its transfer function is expressed by the following equation (3).
  • s is the Laplace transform operator and ⁇ p is the peak angular frequency.
  • q is a quality factor, which is a factor that determines the passband width of the filter.
  • a band pass filter of order m may be used.
  • m is an integer of 3 or more.
  • the sub-side torque current pulsation component extraction unit 7B matches the peak angular frequency ⁇ p with the electrical machine resonance angle frequency.
  • the electrical machine spring resonance angular frequency has the property of fluctuating depending on the driving condition. Therefore, the pass band width of the band pass filter needs to be designed to be wider in anticipation of fluctuation of the electrical spring resonance angular frequency.
  • the sub-side torque current pulsation component extraction unit 7B is configured to track the peak angular frequency ⁇ p to the electrical spring resonant angular frequency by measuring the electrical spring resonant angular frequency, that is, to dynamically move the center frequency of the band pass filter. In this case, the pass band width can be narrowed.
  • the center frequency corresponds to the peak angular frequency ⁇ p .
  • FIG. 4 is a diagram showing a third configuration example of the sub-side torque current pulsation component extraction unit shown in FIG.
  • the sub-side torque current pulsation component extraction unit 7C shown in FIG. 4 includes a pulsation frequency measurement unit 71, a cosine wave generator 72, a sine wave generator 73, a Fourier cosine coefficient operation unit 74, a Fourier sine coefficient operation unit 75, and an AC restorer. 76 is provided.
  • the pulsation frequency measurement unit 71 measures the pulsation frequency included in the input signal that is the current detected by the current detection unit 4b, that is, the pulsation frequency included in the current detected by the current detection unit 4b.
  • the sub-side torque current pulsation component extraction unit 7 includes a coordinate converter (not shown) that performs coordinate conversion using the current detected by the current detection unit 4b and the magnetic pole position estimated by the magnetic pole position estimation unit 5b.
  • the coordinate converter converts the current of the three-phase coordinate system detected by the current detection unit 4b to the rotating orthogonal coordinate system using the magnetic pole position estimated by the magnetic pole position estimation unit 5b, and outputs it.
  • the above input signal corresponds to the current converted by the coordinate converter.
  • the cosine wave generator 72 generates a cosine wave signal oscillating at a pulsating frequency
  • the sine wave generator 73 generates a sinusoidal signal oscillating at a pulsating frequency.
  • the Fourier cosine coefficient operation unit 74 performs Fourier series expansion of the input signal which is the current detected by the current detection unit 4b using the cosine wave signal from the cosine wave generator 72, and the specific frequency included in the input signal Among component magnitudes, a Fourier cosine coefficient, which is the magnitude of the cosine component, is calculated.
  • the Fourier cosine coefficient is a coefficient obtained by expanding an even function having an arbitrary period into a series of cos.
  • the Fourier sine coefficient operation unit 75 performs Fourier series expansion of the input signal using the sine wave signal from the sine wave generator 73, and the magnitude of the sine component in the magnitude of the specific frequency component of the input signal. Calculate the Fourier sine coefficient which is The Fourier sine coefficient is a coefficient obtained by expanding an odd function having an arbitrary period into a series of sin.
  • the ac restorer 76 comprises the cosine wave signal from the cosine wave generator 72, the sine wave signal from the sine wave generator 73, the Fourier cosine coefficient obtained by the Fourier series expansion, and the Fourier obtained by the Fourier series expansion.
  • the alternating current is restored using the sine coefficient.
  • Fourier series expansion is to extract the magnitude and phase of a specific frequency component from an input signal.
  • the magnitude of the specific frequency component can be represented by the magnitude of the cosine component and the magnitude of the sine component.
  • the phase can be represented by the ratio of the magnitude of the cosine component to the magnitude of the sine component.
  • the characteristics of the band pass filter can be obtained by the Fourier series expansion and the inverse transformation. In the inverse transformation, an alternating current of a specific frequency is output based on the magnitude of the cosine component of the specific frequency component extracted in the Fourier series expansion and the magnitude and phase of the sine component.
  • the function of the driving device 100 When the function of the driving device 100 is mounted on a processing device such as a microcomputer, the function needs to be discretized and mounted.
  • the band pass filter of the formula (3) above the peak angular frequency ⁇ p
  • the calculation accuracy tends to decrease.
  • Fourier series expansion can be discretized to suppress a decrease in calculation accuracy even if the peak angular frequency ⁇ p is changed. Therefore, the method of extracting the pulsation component in the vicinity of the electrical machine spring resonance angular frequency included in the torque current of the sub-side synchronous motor 1b using Fourier series expansion is excellent in terms of mounting.
  • the method of extracting the pulsation component using Fourier series expansion is considered to be useful when changing the peak frequency of the band pass filter, but if the calculation accuracy can be ensured, the pulsation frequency shown in FIG.
  • a band pass filter obtained by combining the measurement unit 71 and the equation (3) may be used as the pulsation component extraction unit 7 shown in FIG.
  • the pulsation component extraction unit 7 can extract the pulsation component in the vicinity of the electrical machine spring resonance angular frequency included in the torque current of the sub-side synchronous motor 1b, it may be configured by any of the filters shown in FIGS. Good. Note that since the high pass filter is easier to design and to be mounted on a circuit than a band pass filter, it is preferable to select the high pass filter when importance is placed on the ease of design and mounting on a circuit. When it is desired to obtain sharp cutoff characteristics, it is preferable to select a band pass filter.
  • the subtractor 8 shown in FIG. 1 obtains an angular difference which is a difference between the magnetic pole positions of the rotors of the main synchronous motor 1a and the sub synchronous motor 1b.
  • the magnetic pole position is equal to the rotational position or rotational angle of the rotor of each of the main synchronous motor 1a and the sub synchronous motor 1b.
  • the voltage equation is expressed as the following equation (4).
  • a torque equation is represented like the following (5) Formula.
  • the first term of the right side of the following equation (5) represents the magnet torque, and the second term represents the reluctance torque.
  • the magnet torque is proportional to the q-axis current, and the reluctance torque is proportional to the product of the d-axis current and the q-axis current.
  • R a is an armature resistance
  • L d is a d-axis inductance
  • L q is a q-axis inductance
  • P m is a pole pair number
  • a a is an armature flux linkage number
  • ⁇ e represents angular velocity
  • i d is the d-axis current
  • i q is the q-axis current
  • v d is d-axis voltage
  • v q is q-axis voltage
  • t the generation torque.
  • the subscripts “ x ” of these coefficients are for identifying whether the synchronous motor is the main side or the sub side.
  • the subscript "x” is attached to, or the subscript "x" is omitted.
  • the subscript "m” is attached in place of the "x” subscript represents the main side
  • "s” is attached in place of the "x” subscript represents the sub side.
  • FIG. 5 is a first diagram for explaining the behavior of the torque change of the sub-side synchronous motor shown in FIG.
  • FIG. 6 is a second diagram for explaining the behavior of the torque change of the sub-side synchronous motor shown in FIG.
  • FIG. 7 is a third diagram for explaining the behavior of the torque change of the sub-side synchronous motor shown in FIG.
  • FIG. 8 is a fourth diagram for explaining the behavior of the torque change of the sub-side synchronous motor shown in FIG.
  • FIG. 9 is a fifth diagram for explaining the behavior of the torque change of the sub-side synchronous motor shown in FIG. FIG.
  • FIG. 10 is a sixth diagram for explaining the behavior of the torque change of the sub-side synchronous motor shown in FIG.
  • FIG. 11 is a diagram showing the d-axis current on the main side shown in FIGS. 6, 7, 9 and 10, the sign of the angle difference, and the state of the torque of the sub-side synchronous motor in association with each other.
  • FIGS. 5 to 10 The behavior of the torque change of the sub-side synchronous motor 1b by magnetic flux current compensation is shown in FIGS. 5 to 10, and the contents of FIGS. 5 to 10 are disclosed in Patent Document 1.
  • FIG. The technology disclosed in Patent Document 1 is used for controlling the driving of the sub-side synchronous motor 1b using the angle difference between the magnetic pole positions of the rotors of the main-side synchronous motor 1a and the sub-side synchronous motor 1b.
  • the drive device 100 according to the first embodiment is common to the first embodiment. However, the drive apparatus 100 according to the first embodiment differs from the technology disclosed in Patent Document 1 in that the speed difference between the main synchronous motor 1a and the sub synchronous motor 1b is not used.
  • the load of the main synchronous motor 1a is larger than the load of the sub synchronous motor 1b. That is, the main synchronous motor 1a has a heavy load.
  • the drive device 100 applies the same voltage to the two synchronous motors, but when the main synchronous motor 1a has a heavy load, the magnetic flux current of the sub synchronous motor 1b flows in the positive direction. This is apparent by solving the above equation (4).
  • the case where a positive magnetic flux current flows in the main synchronous motor 1a will be considered.
  • q-axis voltage of the main synchronous motor 1a is by increasing in the positive direction, the voltage command vector v ⁇ changes from dq * into v ⁇ dq **.
  • the q-axis voltage of the main synchronous motor 1a changes as described above, the d-axis voltage of the sub synchronous motor 1b decreases, and the q axis voltage of the sub synchronous motor 1b increases.
  • the torque of the sub-side synchronous motor 1b changes compared to the torque of the sub-side synchronous motor 1b shown in FIG.
  • the synchronous motor is a surface magnet type synchronous AC motor in which permanent magnets are provided on the outer peripheral surface of the rotor core, and there is no reluctance torque.
  • the torque of the sub synchronous motor 1b when the current of the sub synchronous motor 1b changes is reduced compared to the state of FIG.
  • FIG. 7 shows the torque state of the sub-side synchronous motor 1b when the negative magnetic flux current flows through the main-side synchronous motor 1a, contrary to the case of FIG.
  • the voltage command vector v ⁇ changes from dq * into v ⁇ dq **.
  • the d-axis voltage of the sub-side synchronous motor 1 b increases, and the q-axis voltage of the sub-side synchronous motor 1 b decreases.
  • the q-axis current of the sub-side synchronous motor 1b increases as the d-axis voltage of the sub-side synchronous motor 1b increases, and the d-axis current of the sub-side synchronous motor 1b decreases as the q-axis voltage of the sub-side synchronous motor 1b decreases. Decreases. In this case, the torque of the sub synchronous motor 1b is increased compared to the state of FIG.
  • the d-axis current of the main synchronous motor 1a is zero, and the load of the main synchronous motor 1a is larger than the load of the sub synchronous motor 1b, that is, the main synchronous motor 1a is heavier. It has become. Since the same voltage is applied to the main synchronous motor 1a and the sub synchronous motor 1b, the d-axis current of the sub synchronous motor 1b flows in the negative direction when the load on the sub synchronous motor 1b is heavy. .
  • the q-axis voltage of the sub-side synchronous motor 1b increases, the d-axis current of the sub-side synchronous motor 1b decreases. In this case, the torque of the sub synchronous motor 1b is increased compared to the state of FIG.
  • FIG. 10 shows the torque state of the sub-side synchronous motor 1b when the negative magnetic flux current flows through the main-side synchronous motor 1a, contrary to the case of FIG. In this case, the q-axis current of the sub synchronous motor 1b decreases. Therefore, the torque of the sub synchronous motor 1b is reduced compared to the state of FIG.
  • FIG. 11 shows the main d-axis current shown in FIGS. 6, 7, 9 and 10, the sign of the angle difference, and the state of the torque of the sub-side synchronous motor 1b in association with each other.
  • the angle difference ⁇ between the two synchronous motors is positive when the d-axis current of the main synchronous motor 1a is increased if the angle difference ⁇ between the two synchronous motors is determined as in the following equation (6): For example, the torque of the sub synchronous motor 1b decreases, and if the angular difference ⁇ is negative, the torque of the sub synchronous motor 1b increases.
  • ⁇ es in the following equation (6) represents the magnetic pole position of the sub-side synchronous motor 1b in electrical angle
  • ⁇ em represents the magnetic pole position of the main-side synchronous motor 1a in electrical angle.
  • the subtractor 8 is used to determine the angular difference ⁇ between the magnetic pole positions of the two synchronous motors.
  • a magnetic flux current command determination unit 9 shown in FIG. 1 determines a magnetic flux current command for stabilizing driving of the sub-side synchronous motor 1b.
  • changing the torque of the sub-side synchronous motor 1b by changing the magnetic flux current is as described above, it is important in the drive device 100 how to determine the magnetic flux current command.
  • the magnetic flux current is determined using the speed difference between the two synchronous motors.
  • the flux current is determined from the pulsating component of the torque current of the sub-side synchronous motor 1b. Decide. In order to explain the reason, it is necessary to describe the speed estimation method of the synchronous motor and the speed estimation error caused by the flux current change.
  • the arctangent method is the most primitive position estimation method and is widely known.
  • the following equation (7) is a voltage equation of the surface magnet type synchronous AC motor on the stator coordinates.
  • p is a differential operator
  • ⁇ e is a magnetic pole position (electrical angle)
  • R a is an armature resistance
  • L a is an armature inductance
  • v ⁇ and v ⁇ are voltages on stator coordinates
  • ⁇ a is an armature
  • the number of flux linkages, i ⁇ and i ⁇ are currents on stator coordinates.
  • the second term on the right side of the above equation (7) represents speed electromotive force.
  • the term of speed electromotive force can be expressed in the form of the following equation (8).
  • e ⁇ is an ⁇ -axis velocity electromotive force
  • e ⁇ is a ⁇ -axis velocity electromotive force
  • p is a differential operator
  • ⁇ ⁇ r is a rotor ⁇ -axis flux
  • ⁇ ⁇ r is a rotor ⁇ -axis flux
  • ⁇ a is an armature chain
  • the cross magnetic flux number, ⁇ e is the magnetic pole position (electrical angle), and ⁇ e is the angular velocity.
  • the symbol “ ⁇ ” represents an estimated value.
  • the rotor magnetic pole position can be estimated by calculating the above equation (10) to obtain the rotor magnetic flux and calculating the arc tangent as in the following equation (11). Since the angular velocity can be calculated by using the estimated rotor magnetic pole position, the estimated angular velocity ⁇ ⁇ e is calculated by the following equation (12). However, in order to avoid the influence of differential noise, it is usual to apply a low pass filter when using this estimated angular velocity ⁇ ⁇ e for control.
  • the estimated angular velocity ⁇ ⁇ e can also be calculated by estimating the velocity electromotive force and dividing the amplitude by the armature flux linkage number ⁇ a as shown in the following equation (13).
  • the error by the speed estimation method by the following equation (12) is smaller than the error by the method shown by the following equation (13).
  • FIG. 12 is a diagram showing an example in which the magnetic pole position estimation unit shown in FIG. 1 is configured by an adaptive flux observer.
  • the magnetic pole position estimation unit 5c shown in FIG. 12 is configured by the magnetic pole position estimation units 5a and 5b shown in FIG. 1 by an adaptive flux observer.
  • the synchronous motor 1c shown in FIG. 12 corresponds to each of the main synchronous motor 1a and the sub synchronous motor 1b shown in FIG.
  • the magnetic pole position estimation unit 5c is generated by the current detection units 4a and 4b, and a voltage vector that is a voltage command generated by the current control unit 6 and input to the magnetic pole position estimation units 5a and 5b by the adaptive magnetic flux observer method.
  • the rotation speed of the synchronous motor 1c is estimated using current information that is the dq axis current detection value on the dq coordinate axis, and the estimated angular velocity ⁇ ⁇ e and the estimated magnetic pole position ⁇ ⁇ e are output.
  • the magnetic pole position estimation unit 5c calculates a model deviation ⁇ based on the voltage vector of the synchronous motor 1c, the current vector, the primary angular frequency ⁇ 1 of the inverter, and the estimated angular velocity ⁇ ⁇ e , and the model And an angular velocity estimator 52 that calculates an estimated angular velocity ⁇ ⁇ e based on the deviation ⁇ .
  • the magnetic pole position estimation unit 5c calculates the estimated flux vector, the estimated current vector, by using the estimated angular velocity omega ⁇ e, the primary angular frequency arithmetic unit 53 for calculating the primary angular frequency omega 1, the primary angular frequency omega 1
  • an integrator 54 for integrating and outputting an estimated magnetic pole position ⁇ ⁇ e .
  • Model deviation calculation unit 51 includes a voltage vector of the synchronous motor 1c, a current vector, the primary angular frequency omega 1, based on the estimated angular velocity omega ⁇ e, the estimated flux vector and the current estimation calculating and outputting an estimated current vector And a subtractor 512 which calculates and outputs a current deviation vector by subtracting the current vector from the estimated current vector.
  • the model deviation calculator 51 further includes a deviation calculator 513 which receives the current deviation vector from the subtractor 512, extracts the orthogonal component of the estimated magnetic flux vector as a scalar quantity, and outputs this value as a model deviation ⁇ .
  • the current estimator 511 estimates the current and the magnetic flux from the state equation of the synchronous motor 1c.
  • the synchronous motor 1c is a general permanent magnet embedded synchronous AC motor, but the current estimator 511 is in the state even if it is a synchronous motor other than the permanent magnet embedded synchronous AC motor. If equations can be established, current estimation can be performed in the same manner.
  • the state equations are expressed as the following equations (14) and (15).
  • L d represents an inductance of d axis
  • L q represents an inductance of q axis
  • R a represents an armature resistance
  • ⁇ 1 represents a primary angular frequency
  • v d represents a d-axis voltage
  • v q represents a q-axis voltage
  • id represents a d-axis current
  • iq represents a q-axis current.
  • phi ds is d-axis stator flux
  • phi qs is the q-axis stator flux
  • phi dr is the d-axis rotor flux
  • omega e is the angular velocity
  • h 32 from h 11 represents a observer gain.
  • the symbol " ⁇ " represents an estimate.
  • h 41 and h 42 represent observer gains.
  • the above equations (14) and (15) are equations based on ordinary induced voltages, the same calculation can be made even if the above equations (14) and (15) are modified and expressed in the form of extended induced voltages Can. Since the estimated angular velocity ⁇ ⁇ e is included in the above equation (14), if the estimated angular velocity ⁇ ⁇ e and the actual angular velocity ⁇ e do not match, an error occurs in the current estimation.
  • the model deviation ⁇ is defined as the following equation (17), and the magnetic pole position estimating unit 5c adjusts the value of the estimated angular velocity ⁇ ⁇ e using the angular velocity estimator 52 so that the model deviation ⁇ becomes zero. .
  • the angular velocity estimator 52 one configured by serially connecting an integrator to a proportional integral controller is known.
  • the first-order angular frequency calculator 53 calculates the first-order angular frequency ⁇ 1 from the estimated magnetic flux vector, the estimated current vector, and the estimated angular velocity ⁇ ⁇ e based on the equation (16).
  • the integrator 54 estimates the magnetic pole position by integrating the primary angular frequency ⁇ 1 .
  • Adaptive magnetic flux observers are recognized as high-performance speed estimation methods because they are robust to fluctuations in the number of interlinkages and do not generate stationary speed estimation errors.
  • the output voltage error is an error between the value of the voltage command that the current control unit 6 gives to the power converter 2 and the actual voltage that the power converter 2 actually outputs.
  • the short circuit prevention time of the semiconductor element of the upper and lower arms in series which comprises the power converter 2, the ON voltage of a semiconductor element, etc. are known.
  • Many of the commercially available motor drive power converters are provided with an output voltage error compensation function, but it is difficult to compensate the output voltage error when the current flowing through the power converter is close to zero. .
  • relatively low cost power converters for motor drive usually have some output voltage error.
  • FIG. 13 is a drive apparatus according to the first embodiment, in which the output voltage error of the power converter is virtually zero, and the rotational speed of the sub-side synchronous motor when two synchronous motors are parallelly driven at low speed in parallel. It is a figure showing speed true value which is a value, d axis current, and q axis current.
  • FIG. 14 is a drive apparatus according to the first embodiment, in which an output voltage error equivalent to that of a real machine is applied, and two synchronous motors are parallelly driven at low speed, which is the true value of the rotational speed of the sub synchronous motor. It is a figure showing a true value, d axis current, and q axis current. In FIG. 13 and FIG. 14, the true value of the rotational speed, the d-axis current and the q-axis current are shown from the top to the bottom. Each horizontal axis is time.
  • this operating condition is a very severe condition for the control method disclosed in Patent Document 1, and in the control method disclosed in Patent Document 1, when this operating condition is given, sub-side synchronous motor 1b. Drive becomes unstable and parallel drive is difficult.
  • the above operating conditions are that the two synchronous motors are driven in parallel at low speed, and that the angular difference between the magnetic pole positions of the two synchronous motors is near zero.
  • FIGS. 13 and 14 show the analysis results under the condition that the difference in load torque between the two synchronous motors is extremely reduced and the angle difference between the magnetic pole positions of the two synchronous motors is close to zero. In the magnetic flux current compensation method, the torque of the sub synchronous motor 1b is changed using the angle difference between the two synchronous motors.
  • the true value of the rotational speed vibrates at the electrical machine spring resonance angular frequency.
  • the rotational speed shows the maximum value for each cycle of the electrical machine resonance angle frequency.
  • the q-axis current of the sub-side synchronous motor 1b also oscillates at the electrical spring resonance frequency with the true value of the rotational speed.
  • the d-axis current supplied to stabilize the drive of the sub-side synchronous motor is positive when the angle difference between the two synchronous motors is positive
  • the d-axis current is, as shown in FIG. 13 and FIG. It becomes a waveform that shows the maximum value.
  • the d-axis current includes the pulsation of the frequency component that is an integral multiple of the electric motor resonance angle frequency.
  • FIG. 15 is a diagram showing a speed estimation waveform when two synchronous motors are parallelly driven at low speed with the output voltage error of the power converter virtually taken as zero in the drive device according to the first embodiment.
  • the true value of the rotational speed is indicated by a solid line
  • the estimated value of the rotational speed by the arctangent method is indicated by an alternate long and short dashed line
  • the estimated value of the rotational speed of the adaptive flux observer is indicated by a broken line.
  • the horizontal axis represents time.
  • the vertical axis represents the number of revolutions of the rotor of the sub-side synchronous motor 1b.
  • FIG. 16 is a diagram showing a speed estimation waveform when two synchronous motors are driven in parallel at a low speed in parallel with an output voltage error equivalent to that of an actual machine in the drive device according to the first embodiment. Similar to FIG. 15, in FIG. 16, the true value of the rotational speed is shown by a solid line, the estimated value of the rotational speed by the arc tangent method is shown by an alternate long and short dashed line, and the estimated value of the rotational speed of the adaptive flux observer is a broken line It is indicated by.
  • the horizontal axis represents time.
  • the vertical axis represents the number of revolutions of the rotor of the sub-side synchronous motor 1b.
  • Some fluid utilization devices using a common drive device set the carrier frequency to 10 kHz or more in order to reduce electromagnetic noise. Such fluid utilization devices tend to have large output voltage errors.
  • FIG. 15 and FIG. 16 are compared, in FIG. 16 it is understood that the S / N ratio (Signal to Noise Radio) of the estimated velocity waveform is greatly deteriorated. This tendency is the same in the speed estimation errors of the arctangent method and the adaptive flux observer, respectively.
  • FIG. 17 is a diagram showing the result of FFT analysis of the estimated velocity waveform of FIG.
  • FIG. 17 shows, from top to bottom, the true value of the rotational speed, the estimated value by the arc tangent method, and the estimated value by the adaptive flux observer.
  • Each horizontal axis is frequency.
  • FIG. 17 it can be seen that, as shown in the dotted frame B, a speed estimation error occurs in the order of integral multiples of the electrical machine resonance angle frequency. The frequency of these speed estimation errors coincides with the frequency of the pulsating component of the flux current which is the d-axis current shown in FIG.
  • FIG. 18 is a diagram showing the result of FFT analysis of the velocity estimation waveform of FIG. Similar to FIG. 17, FIG. 18 shows, from top to bottom, the true value of the rotational speed, the estimated value by the arc tangent method, and the estimated value by the adaptive flux observer. Each horizontal axis is frequency. According to FIG. 18, it can be seen that there is an error of several times in the first-order component of the electrical machine spring resonance angular frequency as shown in the frame A of the dotted line. Further, as shown in the dotted frame B, it can be seen that a very large peak stands in the order of integral multiples of the electrical machine resonance angle frequency. Further, as shown in the dotted line frame C, it can be seen that a peak is raised on the high frequency side by the output voltage error component of the power converter 2.
  • the load of the fluid utilization device is often a square reduction torque load, and the load characteristic thereof is a light load on the low rotation side.
  • the reduced-square torque load is a load whose load torque increases in proportion to the square of the rotational speed of the motor.
  • the torque current decreases on the low rotation side, but the compensation accuracy of the output voltage error of the power converter 2 decreases in the region where the current is small. Furthermore, in the fluid utilization device, the speed electromotive force of the motor decreases in the low rotation range, so the influence of the output voltage error becomes relatively large. As a result, the above-described speed estimation error increases, and the speed difference between the two synchronous motors can not be accurately determined, leading to a control unstable state. Although the inventor of the present invention applied various filtering processes to the velocity difference signal containing a large amount of velocity estimation error to improve the stability, satisfactory performance was not obtained.
  • the estimated magnetic pole position signal is subjected to integration processing in its calculation process to remove high frequency components of the error included in the estimated speed.
  • the low frequency component of the estimated magnetic pole position signal an error signal due to the change of the magnetic flux current command remains, but the error is within several degrees.
  • the torque current of the sub-side synchronous motor 1b is estimated if the error of the magnetic pole position is about several degrees.
  • the error with the torque current in the dq axis is less than a few percent. This is self-evident considering that the cosine function can be approximated to 1 when it is close to zero.
  • the torque current of the sub synchronous motor 1b can be obtained with relatively high accuracy even when the magnetic flux current is changed.
  • a steady position estimation error occurs due to the influence of a change in motor constant, this becomes an estimation error of a direct current component, and this is not a problem when the pulsation component extraction unit 7 extracts a pulsation component.
  • the mechanical system may have a large moment of inertia.
  • the pulsation component appearing in the speed signal may be very small, even when the torque pulsation is large enough to cause the inverter to stop over current.
  • the torque current signal since the torque current signal has a better S / N ratio than the estimated speed signal, it can be said that it is better to use the torque current signal for stabilization compensation.
  • the magnetic flux current command determination unit 9 determines the magnetic flux current command based on the pulsation component of the torque current flowing through the sub-side synchronous motor 1 b.
  • FIG. 19 is a diagram showing a configuration example of the magnetic flux current command determination unit shown in FIG.
  • FIG. 20 is a first diagram for describing a code judging process by the code judging device shown in FIG.
  • FIG. 21 is a second diagram for illustrating the sign determination processing by the sign determination unit shown in FIG.
  • FIG. 22 is a third diagram for explaining the sign judging process by the sign judging unit shown in FIG.
  • the magnetic flux current command determination unit 9 shown in FIG. 19 includes a pulsation suppression control unit 91 and a compensation direction determination unit 92.
  • the magnetic flux current command determination unit 9 receives the pulsation component of the torque current of the sub-side synchronous motor 1b as an input, and determines the magnetic flux current command using the pulsation suppression control unit 91 and the compensation direction determination unit 92.
  • the pulsation suppression control unit 91 includes a gain multiplication unit 911 and a phase adjustment unit 912.
  • the gain multiplication unit 911 adjusts the gain of the torque current pulsation component which is the input signal.
  • the phase adjustment unit 912 adjusts the phase of the torque current pulsation component which is the input signal, and outputs the pulsation component whose amplitude is adjusted. If stability of the system can be ensured by only one of the gain multiplication unit 911 and the phase adjustment unit 912, the pulsation suppression control unit 91 needs to include both the gain multiplication unit 911 and the phase adjustment unit 912. There is no.
  • the gain multiplication unit 911 multiplies the torque current pulsation component which is the input signal by a specific gain and outputs it, and has a role of adjusting the stability and quick response of the system.
  • the gain may be changed according to the operating condition. For example, the gain may be increased in the low speed region and reduced in the high speed region.
  • the phase adjustment unit 912 includes, for example, a phase delay compensator, a low pass filter, an integration controller, and the like.
  • the phase delay compensator is intended to stabilize by lowering the gain by a fixed value in a high frequency region, and is generally used in the industry. Since the low pass filter and the integral controller also have the property of changing the signal phase in the high frequency region, the low pass filter or the integral controller can be used in the same manner as the phase delay compensator.
  • the cut-off angular frequency may be set to 1/3 or less of the electrical machine resonance angle frequency. If possible, the value is 1/10 to 1/20 of the electric motor resonance angle frequency. By setting in this manner, the phase can be delayed by approximately 90 degrees in the vicinity of the electrical machine spring resonance angular frequency, and control stability is enhanced.
  • a dead zone may be provided in any of the input and output of the pulsation suppression control unit 91.
  • the dead zone is useful for removing frequency components other than the electrical spring resonance that could not be removed by the above-described pulsation component extraction unit 7.
  • the compensation direction determination unit 92 is composed of a code determination unit 921 and a multiplier 922, and based on the operation principle described in FIGS. 5 to 11, from the angle difference between the magnetic pole positions of the rotors provided in each of the two synchronous motors. , Determine the compensation direction of the flux current command.
  • the code judging unit 921 performs the code judging process shown in FIG. 20 to FIG.
  • the horizontal axes in FIGS. 20 to 22 represent angle differences which are inputs of the code determination unit 921.
  • the angular difference indicates a positive or negative value as shown in FIG.
  • the vertical axes in FIG. 20 to FIG. 22 show the value of the output of the code determination unit 921.
  • the most basic code determination process is the method shown in FIG.
  • the sign determiner 921 outputs “1” when the angle difference indicates positive, and outputs “ ⁇ 1” when the angle difference indicates negative.
  • chattering may occur when the angle difference is close to zero. Therefore, as shown in FIG. 21, in the region where the angle difference is close to zero, the output of the code determination unit 921 is gradually switched from “1" to "-1", or the output of the code determination unit 921 is "-1" to It may be configured to gradually switch to “1”.
  • the output value of the code determination unit 921 may be gradually decreased as the angle difference increases.
  • the multiplier 922 multiplies the output of the sign determination unit 921 and the output of the pulsation suppression control unit 91 to generate a magnetic flux current command. That is, the magnetic flux current command determination unit 9 determines the magnetic flux current command based on the pulsation component suppressed by the pulsation suppression control unit 91 and the compensation direction determined by the sign determiner 921 of the compensation direction determination unit 92.
  • the effects of using the generated flux current command are as follows.
  • drive device 100 determines the magnetic flux current command using the pulsation component of the torque current of sub side synchronous motor 1b and the angle difference between the magnetic pole positions of the two synchronous motors. It is configured to As a result, the S / N ratio of the magnetic flux current command is improved, and the influence of the speed estimation error is reduced. As a result, noise and vibration increase, motor efficiency decreases, and problems such as step out are solved. Further, since the stability at the time of low speed driving is improved, the lower limit value of the rotational speed can be maintained at the same value as that of the existing synchronous motor drive device which drives one synchronous motor with one power converter. This makes it easy to replace the existing synchronous motor drive device with drive device 100 according to the first embodiment.
  • Second Embodiment In the second embodiment, a configuration example will be described in which the magnetic flux current is determined using the pulsating component of the active power consumed by the sub-side synchronous motor 1b.
  • the magnetic flux current is determined using the pulsating component of the active power consumed by the sub-side synchronous motor 1b.
  • One method for that is the method using the pulsation component of the torque current described in the first embodiment.
  • the pulsating component of the active power is used instead of the pulsating component of the torque current to generate the magnetic flux current.
  • FIG. 23 is a diagram showing a configuration of a drive device according to Embodiment 2 of the present invention.
  • Drive device 100A according to the second embodiment includes a pulsation component extraction unit 70A instead of pulsation component extraction unit 70 shown in FIG.
  • the pulsation component extraction unit 70A includes the sub-side active power pulsation component extraction unit 10.
  • the other configuration is the same as or equivalent to the configuration of the first embodiment, and the same or equivalent components are denoted by the same reference numerals and redundant description will be omitted.
  • Active power P x can be obtained by the following equation (18) using three-phase voltage commands v u * , v v * , v w * and phase currents i u * , i v * , i w * .
  • the subscript " x " is to distinguish between the main side and the sub side.
  • Ra is an armature resistance.
  • the three-phase voltage commands v u * , v v * , v w * are obtained from the current control unit 6.
  • the phase currents i u * , i v * , i w * are obtained from the current detection unit 4 b.
  • the second term on the right side of the equation (18) represents copper loss due to armature resistance. Strictly speaking, copper loss is also a part of effective power, but it is better to subtract copper loss since information to be obtained here corresponds to torque pulsation of the sub-side synchronous motor 1b. However, since the armature resistance may be small enough to be ignored, in this case, only the first term on the right side may be calculated.
  • sub side active power pulsation component extraction unit 10 performs the same arithmetic processing as sub side torque current pulsation component extraction unit 7 described in the first embodiment, and pulsation is generated from the active power of sub side synchronous motor 1b. Extract the ingredients. If the magnetic flux current command is determined using this information, the same effect as that of the first embodiment can be obtained.
  • the magnetic flux current command determination unit 9 of the second embodiment includes the pulsation suppression control unit 91 configured by the gain multiplication unit 911 and the phase adjustment unit 912 as in the first embodiment, the gain of the second embodiment is the same.
  • the multiplication unit 911 adjusts the gain of the active power pulsation component which is the input signal
  • the phase adjustment unit 912 of the second embodiment adjusts the phase of the active power pulsation component which is the input signal.
  • the pulsation suppression control unit 91 according to the second embodiment can always ensure the stability of the system if only one of the gain multiplication unit 911 and the phase adjustment unit 912 can ensure the system stability. It is not necessary to include both the multiplication unit 911 and the phase adjustment unit 912.
  • the second embodiment is useful when the inertia moment of the load connected to the synchronous motor is large, and the calculation amount is small because the calculation of the coordinate conversion is not used as compared with the first embodiment, so a simple configuration is provided. It is useful when using an arithmetic unit of Specifically, sub-side torque current pulsation component extraction unit 7 of the first embodiment rotates the current of the three-phase coordinate system detected by current detection unit 4b using a signal from magnetic pole position estimation unit 5b. Coordinate conversion is performed to an orthogonal coordinate system to obtain a torque current, and a pulsation component of the torque current is extracted.
  • the sub-side active power pulsation component extraction unit 10 is effective by using the current of the three-phase coordinate system detected by the current detection unit 4b as it is as in equation (18). The power is determined and the pulsating component of this active power is extracted. Then, the magnetic flux current command determination unit 9 of the second embodiment can determine the magnetic flux current command by using this pulsation component. Therefore, in the second embodiment, coordinate conversion is unnecessary, and the amount of calculation is reduced. In the case of an application with a large moment of inertia, by observing the pulsating component of the active power as described above, one coordinate conversion becomes unnecessary, so that it is possible to reduce the calculation load in the second embodiment.
  • FIG. 24 is a diagram showing a configuration of a drive device according to Embodiment 3 of the present invention.
  • the drive device 100B according to the third embodiment includes a pulsation component extraction unit 70B instead of the pulsation component extraction unit 70 illustrated in FIG.
  • the pulsation component extraction unit 70B includes a sub torque current pulsation component extraction unit 7, a main torque current pulsation component extraction unit 11, and a subtractor 8a that is a second subtractor.
  • the other configuration is the same as or equivalent to the configuration of the first embodiment, and the same or equivalent components are denoted by the same reference numerals and redundant description will be omitted.
  • the magnetic flux current command is determined from the difference between the pulsating component of the torque current of the sub-side synchronous motor 1b and the pulsating component of the torque current of the main side synchronous motor 1a.
  • Drive devices 100 and 100A of the first and second embodiments are configured on the premise that the pulsation component of the torque current of main side synchronous motor 1a is small in the steady state. Since the main synchronous motor 1a is vector controlled, if the torque current command value is a constant value, the torque current of the main synchronous motor 1a should follow the command value. However, in reality, the torque current of the main synchronous motor 1a pulsates due to various disturbance factors.
  • Disturbance factors include the short circuit prevention time of the semiconductor elements of the upper and lower arms in series constituting the power converter 2, offset of the current sensor, gain imbalance of the current sensor, distortion of the magnetic flux generated from the magnet provided in the rotor, etc. Conceivable.
  • the torque current pulsation due to these factors also occurs in the sub synchronous motor 1b in the same manner.
  • any AC component when any AC component is superimposed on the torque current command, this also causes disturbance to the magnetic flux current command determination unit 9.
  • the disturbance component superimposed on the torque current command generates torque current pulsation of that frequency also in the sub-side synchronous motor 1b, but the disturbance component is generated due to a different cause from self-oscillation due to electrical spring resonance. Therefore, it is not appropriate to feed back the disturbance component to the flux current command determination unit.
  • the sub-side torque current pulsation component extraction unit 7 of the first embodiment and the main side torque current pulsation component extraction unit 11 of the third embodiment are configured by a high pass filter, a band pass filter, and the like. In order to stabilize the sub synchronous motor 1b more accurately, it is desirable to eliminate the influence of the above-mentioned disturbance factor. However, in the high-pass filter, the rejection characteristics of disturbances are poor, and in order to improve the rejection characteristics of disturbances even in the band-pass filters, it is necessary to measure the electrical spring resonance frequency.
  • the third embodiment is characterized in that the pulsation component of the torque current generated in the main synchronous motor 1a is pulsated in the torque current of the sub synchronous motor 1b. It was configured to be subtracted from the ingredients.
  • the drive device 100B in addition to the sub-side torque current pulsation component extraction unit 7, includes the main-side torque current pulsation component extraction unit 11. Further, the drive device 100B includes a subtractor 8a that obtains a difference between the torque current pulsation component from the sub-side torque current pulsation component extraction unit 7 and the torque current pulsation component from the main torque current pulsation component extraction unit 11.
  • the main torque current pulsation component extraction unit 11 calculates a pulsation component of the torque current of the main synchronous motor 1 a.
  • the calculation method may be the same method as that of the sub-side torque current pulsation component extraction unit 7 described in the first embodiment.
  • the subtractor 8a calculates the difference between the pulsation components of the torque current generated in each of the two synchronous motors, and the magnetic flux current command determination unit 9 determines the magnetic flux current command using the difference.
  • the third embodiment the method of using the difference between the pulsating components of the torque current generated in each of the two synchronous motors has been described.
  • the difference between the pulsating components of the active power may be used.
  • one of sub-side torque current pulsation component extraction unit 7 and main-side torque current pulsation component extraction unit 11 calculates the pulsation of the torque current, and the other is the pulsation of the active power. After calculating ⁇ , these pulsations may be converted to the same scale to obtain a difference.
  • FIG. 25 is a configuration diagram of a fluid utilization device according to Embodiment 4 of the present invention.
  • a fluid utilization apparatus 300 will be described in which a propeller fan 300a is provided on the rotation shaft of the main synchronous motor 1a and a propeller fan 300b is provided on the rotation shaft of the sub synchronous motor 1b.
  • the fluid utilization device 300 shown in FIG. 25 includes the drive device 100 according to the first embodiment, and the drive device 100 includes the power converter drive device 200.
  • the power converter driver 200 includes a processor 201 and a memory 202.
  • the functions shown in FIG. 1, that is, the current control unit 6, the magnetic pole position estimation units 5a and 5b, the sub-side torque current pulsation component extraction unit 7, the subtractor 8 and the magnetic flux current command determination unit 9 use the processor 201 and the memory 202. Function is realized.
  • the processor 201 and the memory 202 are used as shown in FIG. 25, each of the functions described above is realized by software, firmware, or a combination thereof.
  • the software or firmware is written as a program and stored in the memory 202.
  • the processor 201 reads out and executes a program stored in the memory 202. It can also be said that these programs cause a computer to execute the procedures and methods performed by each of the above functions.
  • the memory 202 is a semiconductor memory such as a random access memory (RAM), a read only memory (ROM), a flash memory, an erasable programmable read only memory (EPROM), or an electrically erasable programmable read only memory (EEPROM) (registered trademark). .
  • the semiconductor memory may be non-volatile memory or volatile memory.
  • the memory 202 corresponds to a magnetic disk, a flexible disk, an optical disk, a compact disk, a mini disk, or a DVD (Digital Versatile Disc).
  • the processor 201 may output or store data such as the calculation result to the memory 202, or may store the data in an auxiliary storage device (not shown) via the memory 202.
  • the fluid utilization apparatus 300 may include the drive device 100A of the second embodiment or the drive device 100B of the third embodiment instead of the drive device 100.
  • current control unit 6 magnetic pole position estimation units 5a and 5b, sub-side active power pulsation component extraction unit 10, subtractor 8 and magnetic flux current command determination unit 9 shown in FIG.
  • the function is realized.
  • the determination unit 9 realizes its function using the processor 201 and the memory 202.
  • the power converter 2 may basically have any circuit configuration as long as it can supply any alternating current power to the main synchronous motor 1a and the sub synchronous motor 1b. .
  • Information on the current detected by the current detection units 4a and 4b is transmitted to the processor 201.
  • the two propeller fans 300a and 300b may have the same shape as each other, or may have different shapes. Further, the air flow paths of the two propeller fans 300a and 300b may not necessarily be the same.
  • the fluid utilization device 300 is an air conditioner
  • the two propeller fans 300a and 300b correspond to two blower fans provided in a blower chamber in the outdoor unit of the air conditioner, and the air flow path It corresponds to the air blowing chamber.
  • the air blowing chamber is a space formed by being surrounded by the side plate, the ceiling plate, the bottom plate, the heat exchanger, and the like of the outdoor unit. A flow of air is formed in the blowing chamber as the propeller fans 300a and 300b rotate.
  • two synchronous motors may be provided with fans of different shapes.
  • the cross-sectional area of the flow path in which the fan is provided may be smaller than the cross-sectional area of the flow path in which the other fan is provided.
  • the propeller fan may be driven by one of the synchronous motors, and the pump may be driven by the other synchronous motor, or the like, and fluid utilization devices of different specifications may be driven.
  • the fluid utilization device 300 includes a voltage detection unit that detects a voltage output from the power converter 2, and voltage information detected by the voltage detection unit is input to the processor 201. It may be configured as follows. Although not shown in FIG. 25, the fluid utilization device 300 may include a wind speed sensor that measures the wind speed of a fan, and the wind speed information detected by the wind speed sensor may be input to the processor 201. . Although not shown in FIG. 25, a temperature sensor may be provided to detect the temperature of an object cooled by a fan, and temperature information detected by the temperature sensor may be input to the processor 201. .
  • the fluid load of the fluid utilization device 300 has a damper characteristic, and in a high rotation range, the damper characteristic stabilizes the driving of the open-loop driven synchronous motor. However, in the low rotation range, the damper characteristics weaken and the drive of the synchronous motor becomes unstable, so the fluid utilization device 300 utilizes the parallel drive method described in the first, second, and third embodiments. Thus, in the fourth embodiment, parallel driving of synchronous motors can be realized in a wide speed range. In the fourth embodiment, since there is no need for advanced torque control, the cost increase can be suppressed by modifying the existing synchronous motor drive device that drives one synchronous motor with one power converter. However, it is possible to obtain a fluid utilization device 300 capable of driving the two propeller fans 300a and 300b.
  • FIG. 26 is a block diagram of an air conditioner according to Embodiment 5 of the present invention.
  • An air conditioner 400 according to Embodiment 5 includes a fluid utilization device 300, a refrigerant compressor 401, a condenser 403, a receiver 404, an expansion valve 405, and an evaporator 406.
  • the refrigerant compressor 401 and the condenser 403 are connected by piping.
  • the condenser 403 and the receiver 404 are connected by piping
  • the receiver 404 and the expansion valve 405 are connected by piping
  • the expansion valve 405 and the evaporator 406 are piping It connects, and between the evaporator 406 and the refrigerant compressor 401 is connected by piping.
  • the refrigerant circulates through the refrigerant compressor 401, the condenser 403, the liquid receiver 404, the expansion valve 405, and the evaporator 406.
  • the fluid utilization device 300 includes the current detection units 4a and 4b, the magnetic pole position estimation units 5a and 5c, etc. shown in FIG.
  • the steps of evaporation, compression, condensation, and expansion of the refrigerant are repeatedly performed, so that the refrigerant changes from liquid to gas and further from gas to liquid, so that it is between the refrigerant and the air outside the machine. Heat exchange takes place.
  • the evaporator 406 exerts a cooling function by evaporating the refrigerant liquid in a low pressure state and depriving the air around the evaporator 406 of heat.
  • the refrigerant compressor 401 is for compressing the refrigerant gas gasified by the evaporator 406 to condense the refrigerant into a high pressure gas.
  • the condenser 403 condenses the high-pressure refrigerant gas by releasing the heat of the refrigerant gas that has become high temperature in the refrigerant compressor 401, and converts the refrigerant gas into a refrigerant liquid.
  • the fluid utilization device 300 generates a wind by rotating the propeller fans 300 a and 300 b and cools the condenser 403 by passing the wind to the condenser 403.
  • the expansion valve 405 squeezes and expands the refrigerant liquid to convert the refrigerant liquid into a low pressure liquid in order to evaporate the refrigerant.
  • the receiver 404 is provided to adjust the amount of circulating refrigerant, and may be omitted in a small device.
  • the air conditioner 400 is highly required to reduce the cost, and on the other hand, since the energy saving regulations are being tightened year by year, high efficiency is also required. In recent energy saving regulations, not only the rated operating point but also the driving efficiency at the operating point of low output driving is regarded as important. Therefore, it is necessary to lower the lower limit of the operating speed of the cooling fan as much as possible.
  • the parallel drive device using the technology disclosed in Patent Document 1 is very excellent in cost, but when trying to configure it by position sensorless control, drive in a low rotation range There is a problem that becomes unstable.
  • the carrier frequency is often set high to 10 kHz or higher in order to reduce the carrier noise generated by the power converter 2 of the fluid utilization device 300, and the output voltage error increases, and Drive is prone to instability. Therefore, when the parallel drive device using the technique disclosed by patent document 1 is used for the air conditioner 400, the subject that the drive operation range of the air conditioner 400 became narrow occurred. Therefore, in the parallel drive device using the technology disclosed in Patent Document 1, it is difficult to achieve both cost reduction required for the cooling fan for the air conditioner 400 and high cooling performance.
  • the air conditioner 400 according to the fifth embodiment uses the parallel drive method described in the first to third embodiments, the drive in the low speed region does not become unstable, and the drivable range can be expanded. . Further, since the parallel driving method described in the first to third embodiments is premised on position sensorless control, the manufacturing cost of the air conditioner 400 can be reduced as compared with the case where a position sensor is used. Therefore, in the drive devices 100, 100A, and 100B according to the first to third embodiments, it is possible to achieve both cost reduction and high cooling performance required for the cooling fan for the air conditioner 400.
  • the configuration shown in the above embodiment shows an example of the contents of the present invention, and can be combined with another known technique, and one of the configurations is possible within the scope of the present invention. Parts can be omitted or changed.
  • Reference Signs List 1a main side synchronous motor, 1b sub side synchronous motor, 1c synchronous motor, 2 power converters, 3 DC voltage sources, 4a, 4b current detection units, 5a, 5b, 5c magnetic pole position estimation units, 6 current control units, 7, 7A, 7B, 7C Sub-side torque current pulsation component extraction unit, 8, 8a subtractor, 9 flux current command determination unit, 10 sub-side active power pulsation component extraction unit, 11 Main side torque current pulsation component extraction unit, 51 Model deviation Calculation unit, 52 angular velocity estimator, 53 primary angular frequency calculator, 54 integrator, 70, 70A, 70B pulsation component extraction unit, 71 pulsation frequency measurement unit, 72 cosine wave generator, 73 sine wave generator, 74 Fourier cosine Coefficient operation unit, 75 Fourier sine coefficient operation unit, 76 AC restorer, 91 pulsation suppression control unit, 92 compensation direction determination unit, 100, 1 0A, 100B drive unit, 200 power converter drive unit, 201

Landscapes

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

Abstract

駆動装置(100)は、メイン側同期電動機(1a)に流れる電流を検出する電流検出部(4a)と、サブ側同期電動機(1b)に流れる電流を検出する電流検出部(4b)とを備える。駆動装置(100)は、メイン側同期電動機(1a)が有する回転子の磁極位置を推定する磁極位置推定部(5a)と、サブ側同期電動機(1b)が有する回転子の磁極位置を推定する磁極位置推定部(5b)とを備える。駆動装置(100)は、電圧指令を出力する電流制御部(6)と、脈動成分を抽出するサブ側トルク電流脈動成分抽出部(7)とを備える。駆動装置(100)は、角度差を求める減算器(8)と、磁束電流指令を決定する磁束電流指令決定部(9)とを備えることを特徴とする。

Description

駆動装置、流体利用装置及び空気調和機
 本発明は、1台の電力変換器に並列接続される2台の同期電動機を駆動する駆動装置、駆動装置を備える流体利用装置及び流体利用装置を備える空気調和機に関する。
 同期電動機は、回転子の回転位置である回転子位置に関する情報を用いて回転位置に応じた適切な電圧を印加しなければ、トルクを適切に発生できずに脱調し又は動作を停止する可能性がある。そのため、複数台の同期電動機を駆動する場合、それぞれの同期電動機が備える回転子の回転位置に応じた電圧を印加するために、従来では、同期電動機の台数と同じ台数の電力変換器が用いられていた。しかしながら、このように構成した場合、同期電動機の台数が増加する程、電力変換器の台数が増加するため、コストの増加が課題となる。そこで近年では、同期電動機の制御技術の高度化により、2台の同期電動機を1台の電力変換器で駆動する試みがなされている。
 特許文献1には、1台の電力変換器に接続される2台の同期電動機の内、一方をメイン側同期電動機とし他方をサブ側同期電動機として、駆動装置が2台の同期電動機を、2台の同期電動機の回転速度を用いて制御する手法が開示されている。特許文献1に記載される手法では、駆動装置は、メイン側同期電動機をベクトル制御し、2台の同期電動機間の回転速度の差及び回転位置の差を用いてメイン側同期電動機のd軸電流指令を決定することによって、サブ側同期電動機の駆動を安定化させている。ベクトル制御は、同期電動機を流れる電流を、トルクを発生する電流成分と磁束を発生する電流成分とに分解し、それぞれの電流成分を独立に制御する制御方式である。
米国特許出願公開第2015/0229245号明細書
 特許文献1には、特許文献1に記載の手法を、速度検出手段、位置検出手段といったセンサを用いない位置センサレス制御においても適用可能である旨が記載されている。また、非特許文献1及び非特許文献2には、特許文献1の技術を用いた位置センサレス制御の実証試験結果が掲載されている。位置センサレス制御は一般に低速域での動作安定性に課題がある。その理由は、低回転域では電動機の速度起電力が低下し、電力変換器の出力電圧誤差の影響が相対的に大きくなるためである。1つの電力変換器で1台の同期電動機を駆動する場合でさえ低速での動作安定性が問題となるため、1つの電力変換器で2台の同期電動機を並列駆動する場合には低速域の動作安定性を確保することがますます難しくなる。しかしながら、特許文献1、非特許文献1及び非特許文献2には低速域での駆動方法について具体的な記載がなされていない。
 本発明は、上記に鑑みてなされたものであって、1台の電力変換器を用いて2台の同期電動機を位置センサレスで駆動する場合でも、低速域の回転が不安定になることを抑制できる駆動装置を得ることを目的とする。
 上述した課題を解決し、目的を達成するために、本発明の駆動装置は、並列接続される第1の同期電動機及び第2の同期電動機に電力を供給する電力変換器と、第1の同期電動機に流れる第1の電流を検出する第1の電流検出器と、第2の同期電動機に流れる第2の電流を検出する第2の電流検出器とを備える。駆動装置は、第1の同期電動機を駆動するための電圧指令と、第1の電流とを用いて、第1の同期電動機が有する回転子の第1の磁極位置を推定する第1の磁極位置推定部と、電圧指令と、第2の電流とを用いて、第2の同期電動機が有する回転子の第2の磁極位置を推定する第2の磁極位置推定部とを備える。駆動装置は、トルク電流指令と磁束電流指令と第1の電流と第1の磁極位置とを用いて電圧指令を出力する制御部と、第2の同期電動機に流れるトルク電流に含まれるトルク電流脈動成分と第2の同期電動機で消費される有効電力に含まれる有効電力脈動成分との少なくとも一方の脈動成分を抽出する脈動成分抽出部とを備える。駆動装置は、第1の磁極位置と第2の磁極位置とを用いて、第1の同期電動機及び第2の同期電動機がそれぞれ有する回転子の磁極位置の差である角度差を求める第1の減算器と、脈動成分抽出部で抽出された脈動成分と、角度差とを用いて、磁束電流指令を決定する磁束電流指令決定部とを備えることを特徴とする。
 本発明によれば、1台の電力変換器を用いて2台の同期電動機を位置センサレスで駆動する場合でも、低速域の回転が不安定になることを抑制できる、という効果を奏する。
本発明の実施の形態1に係る駆動装置の構成を示す図 図1に示すサブ側トルク電流脈動成分抽出部の第1の構成例を示す図 図1に示すサブ側トルク電流脈動成分抽出部の第2の構成例を示す図 図1に示すサブ側トルク電流脈動成分抽出部の第3の構成例を示す図 図1に示すサブ側同期電動機のトルク変化の挙動を説明するための第1の図 図1に示すサブ側同期電動機のトルク変化の挙動を説明するための第2の図 図1に示すサブ側同期電動機のトルク変化の挙動を説明するための第3の図 図1に示すサブ側同期電動機のトルク変化の挙動を説明するための第4の図 図1に示すサブ側同期電動機のトルク変化の挙動を説明するための第5の図 図1に示すサブ側同期電動機のトルク変化の挙動を説明するための第6の図 図6、図7、図9及び図10に示すメイン側のd軸電流と、角度差の符号と、サブ側同期電動機のトルクの状態とを対応付けて示す図 図1に示す磁極位置推定部を適応磁束オブザーバにより構成した例を示す図 実施の形態1に係る駆動装置において、仮想的に電力変換器の出力電圧誤差を零として、2台の同期電動機を低速で並列駆動させたときのサブ側同期電動機の回転速度の真値である速度真値とd軸電流とq軸電流とを表す図 実施の形態1に係る駆動装置において、実機同等の出力電圧誤差を入れて、2台の同期電動機を低速で並列駆動させたときのサブ側同期電動機の回転速度の真値である速度真値とd軸電流とq軸電流とを表す図 実施の形態1に係る駆動装置において、仮想的に電力変換器の出力電圧誤差を零として、2台の同期電動機を低速で並列駆動させたときの速度推定波形を表す図 実施の形態1に係る駆動装置において、実機同等の出力電圧誤差を入れて、2台の同期電動機を低速で並列駆動させたときの速度推定波形を表す図 図15の速度推定波形をFFT(Fast Fourier Transform)解析した結果を示す図 図16の速度推定波形をFFT解析した結果を示す図 図1に示す磁束電流指令決定部の構成例を示す図 図19に示す符号判定器による符号判定処理を説明するための第1の図 図19に示す符号判定器による符号判定処理を説明するための第2の図 図19に示す符号判定器による符号判定処理を説明するための第3の図 本発明の実施の形態2に係る駆動装置の構成を示す図 本発明の実施の形態3に係る駆動装置の構成を示す図 本発明の実施の形態4に係る流体利用装置の構成図 本発明の実施の形態5に係る空気調和機の構成図
 以下に、本発明の実施の形態に係る駆動装置、流体利用装置及び空気調和機を図面に基づいて詳細に説明する。なお、この実施の形態によりこの発明が限定されるものではない。
実施の形態1.
 図1は本発明の実施の形態1に係る駆動装置の構成を示す図である。同期電動機は、回転子に永久磁石が設けられる永久磁石界磁式同期電動機と、回転子に界磁巻線が巻かれている巻線界磁式同期電動機と、回転子の突極性を利用して回転トルクを得るリラクタンス式同期電動機とに大別される。実施の形態1に係る駆動装置100には、これらの同期電動機の種別の内、同種の同期電動機、例えば永久磁石界磁式同期電動機が2台並列に接続されている。実施の形態1では、2台の同期電動機の内、一方をメイン側同期電動機1aと称し、他方をサブ側同期電動機1bと称する。メイン側同期電動機1aは第1の同期電動機であり、サブ側同期電動機1bは第2の同期電動機である。
 なお、実施の形態1では、三相の永久磁石界磁式同期電動機が用いられるが、2台の同期電動機のそれぞれのモータ定数が同程度であればよく、永久磁石界磁式以外の同期電動機を用いてもよいし、二相、五相などの三相以外の相数の同期電動機を用いてもよい。
 駆動装置100は、並列接続されるメイン側同期電動機1a及びサブ側同期電動機1bに電力を供給する電力変換器2と、メイン側同期電動機1aに流れる第1の電流を検出する電流検出部4aと、サブ側同期電動機1bに流れる第2の電流を検出する電流検出部4bと、第1の磁極位置推定部である磁極位置推定部5aとを備える。また、駆動装置100は、第2の磁極位置推定部である磁極位置推定部5bと、電圧指令を出力する制御部である電流制御部6と、脈動成分抽出部70と、第1の減算器である減算器8と、磁束電流指令決定部9とを備える。脈動成分抽出部70は、サブ側トルク電流脈動成分抽出部7を備える。以下では、サブ側トルク電流脈動成分抽出部7を単に「脈動成分抽出部7」と称する場合がある。
 電力変換器2は、直流電圧源3から供給される直流電力を交流電力に変換してメイン側同期電動機1a及びサブ側同期電動機1bへ出力する。実施の形態1では電力変換器2に電圧形インバータが用いられる。電圧形インバータは、直流電圧源3から供給される直流電圧をスイッチングして交流電圧に変換する装置である。なお、電力変換器2は、メイン側同期電動機1a及びサブ側同期電動機1bを駆動するための交流電力を出力できるものであれば、電圧形インバータに限定されず、電流形インバータ、交流電力を振幅及び周波数が異なる交流電力に変換するマトリックスコンバータ、複数の変換器の出力を直列又は並列に接続したマルチレベル変換器などの回路でもよい。
 第1の電流検出器である電流検出部4aは、電力変換器2からメイン側同期電動機1aに流れる相電流を検出し、検出した相電流の値を示す電流情報を出力する。第2の電流検出器である電流検出部4bは、電力変換器2からサブ側同期電動機1bに流れる相電流を検出し、検出した相電流の値を示す電流情報を出力する。
 電流検出部4a,4bはCT(Current Transformer)と呼ばれる計器用変流器を用いた電流センサであっても良いし、シャント抵抗を用いた電流センサであっても良い。また電流検出部4a,4bは、これらを組み合わせたものでも良い。実施の形態1に係る駆動装置100では、同期電動機の近くに設けられた電流検出部4a,4bによって電流が検出される。図1に示した例では、同期電動機に流れる相電流を直接検出しているが、電流検出方式は、キルヒホッフの電流則によって同期電動機に流れる電流を演算できればよく、直接検出する例に限定されない。例えば、電力変換器2の負側直流母線に設けられるシャント抵抗を用いた1シャント電流検出方式、電力変換器2の下アームと直列に接続されるシャント抵抗を用いた下アームシャント電流検出方式などを用いて同期電動機に流れる相電流を検出してもよい。なお、三相の電力変換器2の場合、下アームシャント電流検出方式は、3つの下アームのそれぞれに直列に接続されるシャント抵抗を用いるため、3シャント電流検出方式とも呼ばれる。但し、1シャント電流検出方式又は3シャント電流検出方式では、メイン側同期電動機1a及びサブ側同期電動機1bのそれぞれに流れる電流の合計値のみ計測されるため、メイン側同期電動機1a及びサブ側同期電動機1bの内、何れか一方の同期電動機の近くに電流センサを設ける必要がある。また、言うまでもないが、三相の同期電動機の場合、同期電動機に接続される三相の配線の内の何れか二相の配線に電流センサを設ければ、残りの一相の電流はキルヒホッフの電流則で計算可能であるため、三相の配線の全てに電流センサを設ける必要がない。また電流検出部4a及び電流検出部4bの構成及び配置に関しては様々な方式が考えられるが、基本的にはどの方式を用いても構わない。
 磁極位置推定部5aは、電流検出部4aで検出される第1の電流、すなわちメイン側同期電動機1aに流れる相電流と、電流制御部6から出力される電圧指令とを用いて、メイン側同期電動機1aの第1の磁極位置を推定する。
 磁極位置推定部5bは、電流検出部4bで検出される第2の電流、すなわちサブ側同期電動機1bに流れる相電流と、電流制御部6から出力される電圧指令とを用いて、サブ側同期電動機1bの第2の磁極位置を推定する。
 磁極位置の推定方法には様々な方法が存在するが、同期電動機が備える回転子の回転速度全域の内、中高速域では、同期電動機の速度起電力を示す情報を利用して磁極位置を求めるのが一般的である。速度起電力は、回転子が回転することによって同期電動機内部に生じる誘起電力であり、同期電動機が備える回転子と固定子の間に生じる界磁と、回転子の回転速度とに比例する。磁極位置の推定方法の詳細は後述する。
 電流制御部6は、メイン側同期電動機1aに流れる電流を制御するために、メイン側同期電動機1aが備える回転子の永久磁石による磁束の方向をd軸とし、d軸に直交する軸をq軸として、ベクトル制御によって、電流検出部4aで検出された電流をdq座標系の電流指令値に座標変換するベクトル制御器である。一般的なベクトル制御器では、回転子の磁極を基準としたdq座標上での電流制御が行われる。相電流をdq座標上の値に変換すると、交流量が直流量となり制御が容易となるためである。同期電動機では、q軸電流と同期電動機のマグネットトルクとが比例するため、q軸はトルク軸と称され、q軸電流はトルク電流と称される。q軸電流に対してd軸電流は、固定子で発生する磁束の変化をもたらし、同期電動機の出力電圧の振幅を変化させるため、d軸は磁束軸と称され、d軸電流は磁束電流、励磁電流などと称される。なお、同期電動機の種類には、回転子鉄心の外周面に永久磁石が設けられた表面磁石型同期交流電動機、回転子鉄心の内部に永久磁石が埋め込まれた永久磁石埋込型電動機などがある。永久磁石埋込型電動機では、d軸電流によってリラクタンストルクが変化するため、q軸電流のみがトルクに作用するわけではないが、一般的にはq軸電流をトルク電流と呼ぶことが多い。
 座標変換には、磁極位置推定部5aで演算される磁極位置の推定値が用いられる。なお、電流制御部6には、ベクトル制御におけるdq座標系以外にも、αβ固定子座標系、γδ座標系などの極座標系を用いてもよい。また、電流制御部6には、ベクトル制御の代わりに直接トルク制御(Direct Torque Control:DTC)を採用してもよい。但し、DTCを採用する場合、電流指令を、磁束指令及びトルク指令に換算する必要がある。
 なお、dq座標系ではなく、固定子から発生する磁束を基準とした座標系で制御を行えば、トルク電流と磁束電流をより厳密に計算できる。この座標系は、f-t座標系、n-t座標系などと呼ばれることが多いが、公知であるため詳細については説明を割愛する。実施の形態1では、q軸電流をトルク電流と称し、d軸電流を磁束電流と称する場合があるが、dq座標系以外の座標系を使って制御する場合、マグネットトルクが原理的に発生しないリラクタンス同期電動機を用いる場合などは、この限りではない。
 なお、電流制御部6は、メイン側同期電動機1aに流れるトルク電流がトルク電流指令の値と一致するように制御され、またメイン側同期電動機1aに流れる磁束電流が磁束電流指令の値と一致するように制御される。電流制御部6の具体的な実現方法はどのような方法であってもよいが、電流制御部6は一般的には比例積分制御器及び非干渉化制御器により構成される。トルク電流指令は、磁束電流指令決定部9において速度制御の結果、算出されるものであっても良いし、上位のコントローラから入力されるものであってもよい。磁束電流指令の詳細は後述する。
 電流制御部6によってメイン側同期電動機1aがベクトル制御されたとき、サブ側同期電動機1bはメイン側同期電動機1aに連れ回り駆動されるため、サブ側同期電動機1bがオープンループ駆動している状態になる。同期電動機のオープンループ駆動に関する有名な論文として以下の参考文献1がある。
 (参考文献1)伊東淳一、豊崎次郎、大沢博著 「永久磁石同期電動機のV/f制御の高性能化」、電気学会論文誌D、122巻(2002年)3号 P253-259
 上記参考文献1によれば、同期電動機をオープンループ駆動すると、同期電動機が固有角周波数ωで自己発振して、制御が不安定になる場合があると述べられている。固有角周波数ωは、下記(1)式の近似式により表される。但し、Pは極対数、Φは電機子鎖交磁束数、Lは電機子インダクタンス、Jは慣性モーメントを表す。
Figure JPOXMLDOC01-appb-M000001
 機電連成振動は、電機ばね共振と呼ばれる場合があるため、上記(1)式によって表される固有角周波数ωは、電機ばね共振角周波数とも呼ばれる。上記参考文献1に開示される技術には、電機ばね共振を抑えるために安定化補償器が追加されているが、駆動装置100でも同様の安定化補償が必要となる。そのために、図1に示すサブ側同期電動機1bに流れるトルク電流が電機ばね共振によってどの程度振動しているかを調べる必要がある。
 なお、特許文献1に開示される技術では、メイン側同期電動機及びサブ側同期電動機がそれぞれ有する回転子の回転速度の差である速度差が求められ、この速度差を用いることによって速度差安定化補償が行われている。これによりメイン側同期電動機は安定に制御されているため、特許文献1に開示される技術は、サブ側同期電動機の速度脈動成分を求めて安定化補償を行っていたと言える。特許文献1に開示される技術と実施の形態1との相違点の詳細については後述する。
 サブ側同期電動機1bに流れるトルク電流には、加減速トルクによる成分と、負荷トルクによる成分とが重畳されている。加減速トルクは、同期電動機の加減速に伴う慣性トルクである。負荷トルクは、出力トルクから加減速トルクと、摩擦などの損失とを差し引いたトルクである。図1に示す脈動成分抽出部7では、サブ側同期電動機1bのトルク電流に含まれる電機ばね共振角周波数付近の脈動成分が抽出される。なお、脈動成分抽出部7は、ハイパスフィルタを用いる方法と、バンドパスフィルタを用いる方法の2種類があり、これらは順に説明する。
 図2は図1に示すサブ側トルク電流脈動成分抽出部の第1の構成例を示す図である。図2には、1次のハイパスフィルタを用いたサブ側トルク電流脈動成分抽出部7Aの構成例が示され、その伝達関数は下記(2)式で表される。但し、sはラプラス変換の演算子、ωはカットオフ角周波数である。
Figure JPOXMLDOC01-appb-M000002
 上記(2)式には、1次のハイパスフィルタを用いた場合の伝達関数が示されるが、より急峻なフィルタ特性を得たい場合、次数がnのハイパスフィルタを用いても良い。nは2以上の整数である。ハイパスフィルタを用いる場合、カットオフ角周波数ωは、電機ばね共振角周波数の1/3以下、例えば電機ばね共振角周波数の1/5から1/20の値に設定するのが好適である。
 図3は図1に示すサブ側トルク電流脈動成分抽出部の第2の構成例を示す図である。図3には、2次のバンドパスフィルタを用いたサブ側トルク電流脈動成分抽出部7Bの構成例が示され、その伝達関数は下記(3)式で表される。但し、sはラプラス変換の演算子、ωはピーク角周波数を表す。qは、クオリティファクタであり、フィルタの通過帯域幅を決定する係数である。
Figure JPOXMLDOC01-appb-M000003
 上記(3)式には、2次のバンドパスフィルタを用いた場合の伝達関数が示されるが、より急峻なフィルタ特性を得たい場合、次数がmのバンドパスフィルタを用いても良い。mは3以上の整数である。バンドパスフィルタを用いる場合、サブ側トルク電流脈動成分抽出部7Bは、ピーク角周波数ωと電機ばね共振角周波数とを一致させる。但し、上記の参考文献1では言及されていないが、電機ばね共振角周波数は、駆動条件により変動する性質がある。そのため、バンドパスフィルタの通過帯域幅は、電機ばね共振角周波数の変動分を見越して広めに設計する必要がある。なお、サブ側トルク電流脈動成分抽出部7Bは、電機ばね共振角周波数を実測してピーク角周波数ωを電機ばね共振角周波数にトラッキングするような構成、すなわちバンドパスフィルタの中心周波数を動的に変更する構成としても良く、その場合は通過帯域幅を狭くできる。中心周波数は、ピーク角周波数ωに相当する。
 なお、上記(3)式の計算を行う代わりに、図4に示すようにフーリエ級数展開を用いたバンドパスフィルタを用いても良い。図4は図1に示すサブ側トルク電流脈動成分抽出部の第3の構成例を示す図である。図4に示すサブ側トルク電流脈動成分抽出部7Cは、脈動周波数計測部71、余弦波発生器72、正弦波発生器73、フーリエ余弦係数演算部74、フーリエ正弦係数演算部75及び交流復元器76を備える。
 電流検出部4bで検出された電流である入力信号に含まれる脈動周波数、すなわち電流検出部4bで検出された電流に含まれる脈動周波数が、脈動周波数計測部71で計測される。なお、サブ側トルク電流脈動成分抽出部7は、電流検出部4bで検出された電流と磁極位置推定部5bで推定された磁極位置を用いて座標変換を行う不図示の座標変換器を備える。この座標変換器は、磁極位置推定部5bで推定された磁極位置を用いて、電流検出部4bで検出された三相座標系の電流を、回転する直交座標系に座標変換して出力する。上記の入力信号は、座標変換器で座標変換された電流に相当する。余弦波発生器72は、脈動周波数で振動する余弦波信号を発生し、正弦波発生器73は、脈動周波数で振動する正弦波信号を発生する。
 フーリエ余弦係数演算部74は、余弦波発生器72からの余弦波信号を用いて、電流検出部4bで検出された電流である入力信号のフーリエ級数展開を行い、当該入力信号に含まれる特定周波数成分の大きさのうち、余弦成分の大きさであるフーリエ余弦係数を演算する。フーリエ余弦係数は、任意の周期を持つ偶関数をcosの級数に展開したときの係数である。フーリエ正弦係数演算部75は、正弦波発生器73からの正弦波信号を用いて、当該入力信号のフーリエ級数展開を行い、当該入力信号の特定周波数成分の大きさのうち、正弦成分の大きさであるフーリエ正弦係数を演算する。フーリエ正弦係数は、任意の周期を持つ奇関数をsinの級数に展開したときの係数である。
 交流復元器76は、余弦波発生器72からの余弦波信号と、正弦波発生器73からの正弦波信号と、フーリエ級数展開によって得られたフーリエ余弦係数と、フーリエ級数展開によって得られたフーリエ正弦係数とを用いて、交流を復元する。フーリエ級数展開は、入力信号から特定の周波数成分の大きさと位相とを取り出すことである。特定の周波数成分の大きさは、余弦成分の大きさと正弦成分の大きさとで表すことができる。位相は、余弦成分の大きさと正弦成分の大きさの比で表すことができる。図4に示すサブ側トルク電流脈動成分抽出部7Cによれば、フーリエ級数展開と逆変換とによって、バンドパスフィルタの特性が得られる。逆変換では、フーリエ級数展開で取り出した特定周波数成分の余弦成分の大きさと正弦成分の大きさと位相とに基づいて、特定周波数の交流が出力される。
 マイコンなどの処理装置に駆動装置100の機能を実装する場合、当該機能を離散化して実装する必要があるが、上記(3)式のバンドパスフィルタを離散化して用いる場合、ピーク角周波数ωを変化させると計算精度が変動し、特にピーク角周波数ωを大きくすると計算精度が低下し易くなる。一方、フーリエ級数展開は、離散化して、ピーク角周波数ωを変化させても計算精度の低下を抑制することができる。従って、サブ側同期電動機1bのトルク電流に含まれる電機ばね共振角周波数付近の脈動成分を、フーリエ級数展開を用いて抽出する方法は、実装の面で優れている。このことから、フーリエ級数展開を用いて脈動成分を抽出する方法は、バンドパスフィルタのピーク周波数を変化させる場合に有用と考えられるが、計算精度を確保できる場合には、図4に示す脈動周波数計測部71と上記(3)式とを組み合わせたバンドパスフィルタを、図1に示す脈動成分抽出部7として利用してもよい。
 このように、脈動成分抽出部7は、サブ側同期電動機1bのトルク電流に含まれる電機ばね共振角周波数付近の脈動成分を抽出できれば、図2から図4に示す何れのフィルタで構成してもよい。なお、バンドパスフィルタに比べてハイパスフィルタは、設計及び回路への実装が簡便であるため、設計及び回路への実装の簡便さを重視する場合、ハイパスフィルタを選択すると良い。また、鋭い遮断特性を得たい場合には、バンドバスフィルタを選択すると良い。
 図1に示す減算器8は、メイン側同期電動機1a及びサブ側同期電動機1bがそれぞれ有する回転子の磁極位置の差である角度差を求める。磁極位置は、メイン側同期電動機1a及びサブ側同期電動機1bのそれぞれの回転子の回転位置又は回転角度に等しい。角度差を求める理由を説明するために、以下に永久磁石同期電動機の定常状態における電圧方程式とトルク方程式を示す。
 電圧方程式は下記(4)式のように表される。またトルク方程式は下記(5)式のように表される。下記(5)式の右辺の第一項はマグネットトルクを表し、第二項はリラクタンストルクを表す。マグネットトルクはq軸電流に比例し、リラクタンストルクはd軸電流とq軸電流の積に比例する。
Figure JPOXMLDOC01-appb-M000004
Figure JPOXMLDOC01-appb-M000005
 上記(4)式及び(5)式において、Rは電機子抵抗、Lはd軸インダクタンス、Lはq軸インダクタンス、Pは極対数、Φは電機子鎖交磁束数、ωは角速度、iはd軸電流、iはq軸電流、vはd軸電圧、vはq軸電圧、tは発生トルクを表す。これらの各係数の添字「」は、同期電動機がメイン側であるかサブ側であるかを識別するためのものである。例えば、メイン側とサブ側を識別する必要がない場合、添字に「」が付けられ、又は添字「」が省略される。また添字に「」の代わりに「」が付けられた場合にはメイン側を表し、添字に「」の代わりに「」が付けられた場合にはサブ側を表す。
 次に、図5から図11を用いて、磁束電流補償によるサブ側同期電動機1bのトルク変化の挙動について説明する。図5は図1に示すサブ側同期電動機のトルク変化の挙動を説明するための第1の図である。図6は図1に示すサブ側同期電動機のトルク変化の挙動を説明するための第2の図である。図7は図1に示すサブ側同期電動機のトルク変化の挙動を説明するための第3の図である。図8は図1に示すサブ側同期電動機のトルク変化の挙動を説明するための第4の図である。図9は図1に示すサブ側同期電動機のトルク変化の挙動を説明するための第5の図である。図10は図1に示すサブ側同期電動機のトルク変化の挙動を説明するための第6の図である。図11は図6、図7、図9及び図10に示すメイン側のd軸電流と、角度差の符号と、サブ側同期電動機のトルクの状態とを対応付けて示す図である。
 図5から図10には磁束電流補償によるサブ側同期電動機1bのトルク変化の挙動が示され、図5から図10の内容は特許文献1で開示されている。サブ側同期電動機1bの駆動を安定化する制御に、メイン側同期電動機1a及びサブ側同期電動機1bがそれぞれ有する回転子の磁極位置の角度差を用いる点に関しては、特許文献1で開示される技術と実施の形態1に係る駆動装置100は共通である。しかしながら、実施の形態1に係る駆動装置100では、メイン側同期電動機1aとサブ側同期電動機1bとの速度差を用いない点が、特許文献1で開示される技術と異なる。
 まず図5を基準として、メイン側同期電動機1aのd軸がサブ側同期電動機1bのd軸よりも遅れ位相となっているケースについて説明する。図5ではメイン側同期電動機1aの磁束電流idmが零であり、メイン側同期電動機1aのトルク電流iqmが正方向に流れている場合、電圧指令ベクトルv dq は第二象限の方向に発生する。電圧指令ベクトルv dq は回転する直交座標系(dq軸)の電圧指令を表わす。2台の同期電動機に異なる負荷トルクが発生したとき、2台の同期電動機のモータ定数が等しい場合には、重負荷の同期電動機の位相が遅れる。そのため、図5のケースでは、メイン側同期電動機1aの負荷がサブ側同期電動機1bの負荷よりも大きいと言える。すなわちメイン側同期電動機1aの方が重負荷である。駆動装置100は2台の同期電動機に同じ電圧を印加するが、メイン側同期電動機1aの方が重負荷である場合、サブ側同期電動機1bの磁束電流は正方向に流れる。これは上記(4)式を解くことにより明らかである。
 ここで、図6のようにメイン側同期電動機1aに正の磁束電流が流れた場合を考える。この場合、メイン側同期電動機1aのq軸電圧が正方向に増加することによって、電圧指令ベクトルがv dq からv dq **に変化する。このようにメイン側同期電動機1aのq軸電圧が変化することにより、サブ側同期電動機1bのd軸電圧が減少して、サブ側同期電動機1bのq軸電圧が増加する。サブ側同期電動機1bのd軸電圧が減少すると、サブ側同期電動機1bのq軸の電機子反作用であるωesqsqsが減少する。そのため、サブ側同期電動機1bのq軸電流が減少する。また、サブ側同期電動機1bのq軸電圧が増加することによって、サブ側同期電動機1bのd軸電流が増加する。このようにメイン側同期電動機1aの磁束電流idmであるd軸電流を変化させることによって、サブ側同期電動機1bの電流が変化する。このサブ側同期電動機1bの電流の変化によって、サブ側同期電動機1bのトルクは、図5に示すサブ側同期電動機1bのトルクと比べて変化する。ここでは説明を簡単にするために、同期電動機が、回転子鉄心の外周面に永久磁石が設けられた表面磁石型同期交流電動機であるとして、リラクタンストルクがないものとする。この場合、サブ側同期電動機1bの電流が変化したときのサブ側同期電動機1bのトルクは、図5の状態に比べて減少する。
 図7には、図6の場合とは逆に、メイン側同期電動機1aに負の磁束電流が流れた場合のサブ側同期電動機1bのトルク状態が示される。この場合、メイン側同期電動機1aのq軸電圧が減少することによって、電圧指令ベクトルがv dq からv dq **に変化する。これにより、サブ側同期電動機1bのd軸電圧が増加して、サブ側同期電動機1bのq軸電圧が減少する。サブ側同期電動機1bのd軸電圧が増加したことによりサブ側同期電動機1bのq軸電流は増加し、サブ側同期電動機1bのq軸電圧が減少することによりサブ側同期電動機1bのd軸電流は減少する。この場合、サブ側同期電動機1bのトルクは、図5の状態に比べて増加する。
 次に図8を基準として、メイン側同期電動機1aのd軸がサブ側同期電動機1bのd軸よりも進み位相となっているケースについて説明する。図8では、メイン側同期電動機1aのd軸電流が零であり、メイン側同期電動機1aの負荷がサブ側同期電動機1bの負荷よりも大きい状態、すなわちメイン側同期電動機1aの方が重負荷となっている。メイン側同期電動機1aとサブ側同期電動機1bには同じ電圧が印加されているため、サブ側同期電動機1bの方が重負荷となる場合、サブ側同期電動機1bのd軸電流は負方向に流れる。
 ここで、図9のようにメイン側同期電動機1aに正のd軸電流が流れた場合を考える。この場合、メイン側同期電動機1aのq軸電圧が増加することによって、電圧指令ベクトルがv dq からv dq **に変化する。このようにメイン側同期電動機1aのq軸電圧が変化することにより、サブ側同期電動機1bのd軸電圧が増加して、サブ側同期電動機1bのq軸電圧も増加する。サブ側同期電動機1bのd軸電圧が増加したことにより、サブ側同期電動機1bのq軸電流が増加する。また、サブ側同期電動機1bのq軸電圧が増加することによりサブ側同期電動機1bのd軸電流が減少する。この場合、サブ側同期電動機1bのトルクは、図8の状態に比べて増加する。
 図10には、図9の場合とは逆に、メイン側同期電動機1aに負の磁束電流が流れた場合のサブ側同期電動機1bのトルク状態が示される。この場合、サブ側同期電動機1bのq軸電流は減少する。従って、サブ側同期電動機1bのトルクは、図8の状態に比べて減少する。
 図6、図7、図9及び図10に示すメイン側のd軸電流と、角度差の符号と、サブ側同期電動機1bのトルクの状態とを対応付けて示すものが図11である。2台の同期電動機の角度差λを下記(6)式のように定めた場合、メイン側同期電動機1aのd軸電流を増加させた際、2台の同期電動機の角度差λが正であれば、サブ側同期電動機1bのトルクは減少し、角度差λが負であればサブ側同期電動機1bのトルクは増加する。但し、下記(6)式のθesは、サブ側同期電動機1bの磁極位置を電気角で表したものであり、θemは、メイン側同期電動機1aの磁極位置を電気角で表したものである。一方、メイン側同期電動機1aのd軸電流を減少させた際、2台の同期電動機の角度差λが正であればサブ側同期電動機1bのトルクは増加し、角度差λが負であればサブ側同期電動機1bのトルクは減少する。すなわち、メイン側同期電動機1aのd軸電流を変化させることによりサブ側同期電動機1bのトルクを変化させてサブ側同期電動機1bの駆動を安定化させようとする場合、2台の同期電動機の角度差λが正であるか負であるかによって、d軸電流の補償方向を決定する必要がある。このような理由から、減算器8を用いて、2台の同期電動機の磁極位置の角度差λを求めている。
Figure JPOXMLDOC01-appb-M000006
 図1に示す磁束電流指令決定部9は、サブ側同期電動機1bの駆動を安定化させるための磁束電流指令を決定する。磁束電流を変化させることでサブ側同期電動機1bのトルクを変えることはすでに述べたとおりであるが、駆動装置100では磁束電流指令をどのように決定するかが重要である。特許文献1では、2台の同期電動機の速度差を用いて磁束電流を決定していたが、実施の形態1に係る駆動装置100では、サブ側同期電動機1bのトルク電流の脈動成分から磁束電流を決定する。その理由を説明するには、同期電動機の速度推定法と磁束電流変化によって生じる速度推定誤差とについて述べる必要がある。
 先に述べたとおり、同期電動機の磁極位置の推定方法又は同期電動機が備える回転子の回転速度の推定方法には、様々な方法が検討されているが、同期電動機が備える回転子の回転速度全域の内、中高速域では、同期電動機の速度起電力情報を利用して磁極位置を求めるのが一般的である。ここではアークタンジェント法と適応磁束オブザーバとの2種類の方式について述べる。
 アークタンジェント法は最もプリミティブな位置推定法であり、広く知られている。下記(7)式は固定子座標上における表面磁石型同期交流電動機の電圧方程式である。但し、pは微分演算子、θは磁極位置(電気角)、Rは電機子抵抗、Lは電機子インダクタンス、vα,vβは固定子座標上の電圧、Φは電機子鎖交磁束数、iα,iβは固定子座標上の電流である。
Figure JPOXMLDOC01-appb-M000007
 上記(7)式の右辺第二項は速度起電力を表している。なお、速度起電力の項は下記(8)式のような形で表現可能である。但し、eαはα軸速度起電力、eβはβ軸速度起電力、pは微分演算子、Φαrは回転子α軸磁束、Φβrは回転子β軸磁束、Φは電機子鎖交磁束数、θは磁極位置(電気角)、ωは角速度である。
Figure JPOXMLDOC01-appb-M000008
 上記(7)式を見てわかるとおり、速度起電力に磁極位置情報であるθが含まれるため、上記(7)式を整理して磁極位置を演算する。まず、回転子磁束項を左辺に、それ以外の項を右辺にまとめると、下記(9)式が得られる。微分計算はノイズを増幅させるため、両辺積分して下記(10)式を得る。ここで、電圧センサなどに直流オフセットがある場合、純粋積分を使うと積分値が発散するため、下記(10)式の計算を行う際は近似積分を使い、直流分は積分しないようにするのが慣例である。
Figure JPOXMLDOC01-appb-M000009
Figure JPOXMLDOC01-appb-M000010
 上記(9)式及び(10)式において、記号「」は推定値を表す。上記(10)式を計算して、回転子磁束を求め、それを下記(11)式のようにアークタンジェント計算することで、回転子磁極位置を推定することができる。推定された回転子磁極位置を用いれば、角速度を計算可能なため、下記(12)式で推定角速度ω を計算する。但し、微分ノイズの影響を避けるため、この推定角速度ω を制御に使用する場合にはローパスフィルタを掛けるのが普通である。また、下記(13)式のように、速度起電力の推定を行い、その振幅を電機子鎖交磁束数Φで除算することにより、推定角速度ω を計算することもできる。しかしながら、永久磁石の磁束は温度変化により変動するため、下記(13)式の計算方法では温度変化によって定常的な速度推定誤差が生じる。そのため下記(12)式による速度推定方法による誤差は、下記(13)式に示す方法よる誤差よりも少ない。実施の形態1には下記(12)式による速度推定を行った場合について説明する。
Figure JPOXMLDOC01-appb-M000011
Figure JPOXMLDOC01-appb-M000012
Figure JPOXMLDOC01-appb-M000013
 アークタンジェント法以外にも様々な速度推定法が提案されている。以下では、図12を用いて、アークタンジェント法以外の速度推定法の代表例として、適応磁束オブザーバについて述べる。図12は図1に示す磁極位置推定部を適応磁束オブザーバにより構成した例を示す図である。図12に示す磁極位置推定部5cは、図1に示す磁極位置推定部5a,5bを適応磁束オブザーバにより構成したものである。図12に示す同期電動機1cは、図1に示すメイン側同期電動機1a及びサブ側同期電動機1bのそれぞれに対応する。磁極位置推定部5cは、適応磁束オブザーバの手法により、電流制御部6により生成されて磁極位置推定部5a,5bに入力される電圧指令である電圧ベクトルと、電流検出部4a,4bにより生成される電流情報、すなわちdq座標軸上におけるdq軸電流検出値である電流ベクトルとを用いて、同期電動機1cの回転速度を推定し、推定角速度ω と推定磁極位置θ とを出力する。
 磁極位置推定部5cは、同期電動機1cの電圧ベクトルと、電流ベクトルと、インバータの一次角周波数ωと、推定角速度ω とに基づきモデル偏差εを演算するモデル偏差演算部51と、モデル偏差εに基づき推定角速度ω を演算する角速度推定器52とを備える。また磁極位置推定部5cは、推定磁束ベクトルと、推定電流ベクトルと、推定角速度ω とを用いて、一次角周波数ωを演算する一次角周波数演算器53と、一次角周波数ωを積分して推定磁極位置θ を出力する積分器54とを備える。
 モデル偏差演算部51は、同期電動機1cの電圧ベクトルと、電流ベクトルと、一次角周波数ωと、推定角速度ω とに基づき、推定磁束ベクトル及び推定電流ベクトルを演算して出力する電流推定器511と、推定電流ベクトルから電流ベクトルを減算することによって電流偏差ベクトルを演算して出力する減算器512とを備える。またモデル偏差演算部51は、減算器512からの電流偏差ベクトルを入力とし、推定磁束ベクトルの直交成分をスカラ量として抽出し、この値をモデル偏差εとして出力する偏差演算器513を備える。推定磁束ベクトルの直交成分をスカラ量として抽出する手法としては、電流偏差ベクトルを回転二軸上に座標変換する手法と、電流偏差ベクトルと推定磁束ベクトルとの外積値の大きさを演算する手法とが公知である。
 電流推定器511は、同期電動機1cの状態方程式から電流と磁束を推定する。ここでは同期電動機1cは一般的な永久磁石埋込型の同期交流電動機であると仮定するが、永久磁石埋込型の同期交流電動機以外の同期電動機であっても、電流推定器511は、状態方程式が立式できれば、同様の方法で電流推定が可能である。
 同期電動機1cが永久磁石埋込型の同期交流電動機の場合、状態方程式は下記(14)式及び(15)式のように表される。但し、Lはd軸のインダクタンス、Lはq軸のインダクタンス、Rは電機子抵抗、ωは一次角周波数を表す。vはd軸電圧、vはq軸電圧、iはd軸電流、iはq軸電流を表す。また、φdsはd軸固定子磁束、φqsはq軸固定子磁束、φdrはd軸回転子磁束、ωは角速度、h11からh32はオブザーバゲインを表す。記号「」は推定値を表す。
Figure JPOXMLDOC01-appb-M000014
Figure JPOXMLDOC01-appb-M000015
 また、一次角周波数ωは下記(16)式のように与えられる。h41,42はオブザーバゲインを表す。
Figure JPOXMLDOC01-appb-M000016
 上記(14)式及び(15)式は通常の誘起電圧に基づく式であるが、上記(14)式及び(15)式に変形を加えて拡張誘起電圧の形式で表現しても同様の計算ができる。上記(14)式には推定角速度ω が含まれるため、推定角速度ω と実際の角速度ωとが一致していない場合、電流推定に誤差が生じる。ここではモデル偏差εを下記(17)式のように定義し、磁極位置推定部5cはモデル偏差εが零になるように、角速度推定器52を用いて推定角速度ω の値を調整する。角速度推定器52としては、比例積分制御器に積分器を直列接続して構成されたものが公知である。
Figure JPOXMLDOC01-appb-M000017
 一次角周波数演算器53は、上記(16)式を基に、推定磁束ベクトル、推定電流ベクトル及び推定角速度ω から一次角周波数ωを演算する。積分器54は、一次角周波数ωを積分することにより磁極位置を推定する。
 適応磁束オブザーバは、鎖交磁束数の変動にロバストであり、定常的な速度推定誤差が発生しない点で優れているため、高性能な速度推定法として世間に認知されている。
 磁極位置推定部5a,5bの構成例について詳しく説明したところで、次は磁束電流の変化が速度推定誤差に与える影響について述べる。ここでは課題を明確にするため、電力変換器2の出力電圧誤差が無い場合と有る場合との2種類の解析結果について述べる。当該課題は、低回転域で電力変換器2の出力電圧誤差の補償精度が低下すると、低回転域では電動機の速度起電力が低下して、電力変換器2の出力電圧誤差の影響が相対的に大きくなることを意味する。すなわち、当該課題は、出力電圧誤差がある場合、特許文献1の方式だけでは、低回転域で制御不安定に陥りやすいことである。なお、出力電圧誤差は、電流制御部6が電力変換器2へ与える電圧指令の値と電力変換器2が実際に出力する実電圧との誤差のことである。出力電圧誤差の要因としては、電力変換器2を構成する直列の上下アームの半導体素子の短絡防止時間、半導体素子のオン電圧などが知られている。市販される電動機駆動用の電力変換器の多くには、出力電圧誤差の補償機能が設けられているが、当該電力変換器に流れる電流が零に近い場合、出力電圧誤差を補償することが難しい。従って、比較的低価格な電動機駆動用の電力変換器では多少の出力電圧誤差が生じるのが普通である。
 図13は実施の形態1に係る駆動装置において、仮想的に電力変換器の出力電圧誤差を零として、2台の同期電動機を低速で並列駆動させたときのサブ側同期電動機の回転速度の真値である速度真値とd軸電流とq軸電流とを表す図である。図14は実施の形態1に係る駆動装置において、実機同等の出力電圧誤差を入れて、2台の同期電動機を低速で並列駆動させたときのサブ側同期電動機の回転速度の真値である速度真値とd軸電流とq軸電流とを表す図である。図13及び図14には、上から順位に回転速度の真値、d軸電流及びq軸電流が示される。それぞれの横軸は時間である。
 なお、この動作条件は、特許文献1に開示される制御方式にとっては非常に厳しい条件であり、特許文献1に開示される制御方式では、この動作条件が与えられた場合、サブ側同期電動機1bの駆動が不安定になり、並列駆動が困難である。上記の動作条件は、2台の同期電動機を低速で並列駆動させることと、2台の同期電動機の磁極位置の角度差が零近くであることである。図13及び図14には、2台の同期電動機の負荷トルクの差を極めて小さくし、2台の同期電動機の磁極位置の角度差が零近くとなる条件での解析結果が示される。磁束電流補償法では、2台の同期電動機の角度差を利用してサブ側同期電動機1bのトルクを変化させるため、角度差が零に近い場合、サブ側同期電動機1bの速度リプルを完全に零にすることは難しい。そのため、図13及び図14の上から1段目に示すように、回転速度の真値が電機ばね共振角周波数で振動している。電機ばね共振角周波数の1周期毎に回転速度が最大値を示す。また、サブ側同期電動機1bのq軸電流も、図13及び図14の上から3段目に示すように、回転速度の真値が電機ばね共振角周波数で振動している。磁束電流指令の決定法の詳細については後述するが、サブ側同期電動機の駆動を安定化させるために流すd軸電流は、前述したように、2台の同期電動機の角度差が正であるか負であるかによって、d軸電流の補償方向を変える必要があるため、d軸電流は、図13及び図14の上から2段目に示すように、電機ばね共振角周波数の1周期毎に最大値を示すような波形となる。このように、角度差の符号変化の影響によって、d軸電流には、電機ばね共振角周波数の整数倍の周波数成分の脈動が含まれる。なお、電力変換器2の出力電圧に誤差がある場合、dq軸電流には電気角周波数の6倍の脈動成分が重畳される。そのことを除けば、図13に示される回転速度の真値、d軸電流及びq軸電流のそれぞれの波形と、図14に示される回転速度の真値、d軸電流及びq軸電流のそれぞれの波形との間には大きな相違はない。
 図15は実施の形態1に係る駆動装置において、仮想的に電力変換器の出力電圧誤差を零として、2台の同期電動機を低速で並列駆動させたときの速度推定波形を表す図である。図15には、回転数の真値が実線で示され、アークタンジェント法による回転数の推定値が一点鎖線で示され、適応磁束オブザーバにおる回転数の推定値が破線で示される。横軸は時間を表す。縦軸はサブ側同期電動機1bの回転子の回転数を表す。
 実機同等の出力電圧誤差を入れて2台の同期電動機を低速で並列駆動させた場合、磁束電流の変化に伴って過渡的な速度推定誤差が発生する。図15によれば、適応磁束オブザーバによる速度推定誤差は、アークタンジェント法に比べて小さいが、どちらの方式でも速度推定誤差が生じていることが分かる。
 図16は実施の形態1に係る駆動装置において、実機同等の出力電圧誤差を入れて、2台の同期電動機を低速で並列駆動させたときの速度推定波形を表す図である。図16には、図15と同様に、回転数の真値が実線で示され、アークタンジェント法による回転数の推定値が一点鎖線で示され、適応磁束オブザーバにおる回転数の推定値が破線で示される。横軸は時間を表す。縦軸はサブ側同期電動機1bの回転子の回転数を表す。
 一般的な駆動装置を用いた、流体利用装置の中には、電磁騒音を減らすために、キャリア周波数を10kHz以上に設定するものがある。このような流体利用装置では、出力電圧誤差が大きくなる傾向がある。図15と図16を比べると、図16では、推定速度波形のS/N比(Signal to Noise Radio)が大きく悪化していることがわかる。この傾向は、アークタンジェント法及び適応磁束オブザーバのそれぞれの速度推定誤差において同じである。
 図17は図15の速度推定波形をFFT解析した結果を示す図である。図17には、上から順位に回転速度の真値と、アークタンジェント法による速度推定値と、適応磁束オブザーバによる速度推定値とが示される。それぞれの横軸は周波数である。図17によれば、点線の枠B内に示すように、電機ばね共振角周波数の整数倍の次数の速度推定誤差が生じていることがわかる。これらの速度推定誤差の周波数は、図13に示したd軸電流である磁束電流の脈動成分の周波数と一致する。また、点線の枠A内に示すように、電機ばね共振角周波数の1次成分に関しても、真値と推定値との間に振幅の誤差がある。
 図18は図16の速度推定波形をFFT解析した結果を示す図である。図18には、図17と同様に、上から順位に回転速度の真値と、アークタンジェント法による速度推定値と、適応磁束オブザーバによる速度推定値とが示される。それぞれの横軸は周波数である。図18によれば、点線の枠A内に示すように、電機ばね共振角周波数の1次成分に数倍の誤差があることが分かる。また、点線の枠B内に示すように、電機ばね共振角周波数の整数倍の次数に極めて大きなピークが立っていることが分かる。また、点線の枠C内に示すように、電力変換器2の出力電圧誤差成分によって高周波側にピークが立っていることが分かる。
 図17及び図18に示す解析結果は、磁束電流の変化が既存の速度推定器にとって外乱となることを示している。また図17に示すように電力変換器2の出力電圧誤差がない場合、点線の枠B内に示すように、磁束電流の脈動周波数と同じ周波数で速度推定誤差が生じているため、適切に速度推定を行うことができない。また、図17の点線の枠B内に示す周波数成分の速度推定誤差は、また図18のように電力変換器2の出力電圧誤差がある場合、数倍から数十倍に増加することがあり、駆動装置にとって大きな問題となる。特許文献1では、2台の同期電動機の速度差を用いて磁束電流指令を演算し、サブ側同期電動機1bの駆動を安定化させているが、図17の点線の枠B内に示すような速度推定誤差がある場合、サブ側同期電動機1bの安定化を図ることができない。その理由は以下のとおりである。
 まず、サブ側同期電動機1bの駆動を安定化させるために磁束電流指令を脈動させると、それによって速度推定に予期せぬ誤差成分が現れる。その誤差成分を抑え込むために磁束電流指令を周波数で脈動させてしまうと、サブ側同期電動機1bが加振され、サブ側同期電動機1bの振動が大きくなる。サブ側同期電動機1bの振動が大きくなると、サブ側同期電動機1bの駆動の安定化を図るためにメイン側同期電動機1aの磁束電流指令をより大きく変化させなければならない。これにより速度推定誤差がさらに増大するという悪循環が生じる。その結果として、同期電動機には、騒音及び振動が増加し、モータ効率が低下するなどの様々な現象が発生する。また、同期電動機は、トルクを適切に発生できずに脱調し、又は動作を停止する可能性がある。
 特に流体利用装置に、特許文献1に開示される技術を用いた場合、低速域の回転が不安定になるという現象が顕著になる。流体利用装置の負荷は、二乗低減トルク負荷であることが多く、その負荷特性は、低回転側では軽負荷である。二乗低減トルク負荷は、電動機の回転速度の二乗に比例して負荷トルクが増加する負荷である。
 従って、流体利用装置では、低回転側でトルク電流が小さくなるが、電流が小さい領域では電力変換器2の出力電圧誤差の補償精度が低下する。さらに、流体利用装置では、低回転域で電動機の速度起電力が低下するため、出力電圧誤差の影響が相対的に大きくなる。これにより、上述した速度推定誤差が増大して、2台の同期電動機の速度差を正確に求めることができなくなり、制御不安定な状態に至る。本願発明者は、この速度推定誤差を多く含んだ速度差信号に対して様々なフィルタリング処理を施し、安定性の改善を図ったが、満足する性能が得られなかった。
 出力電圧誤差により低周波側で速度推定誤差が生じることは一般に知られているが、磁束電流を変化させることにより速度推定誤差が生じることは、本願発明者の検討によって発見されたものであり、公知ではない。通常、磁束電流指令の変化は緩やかであるため、このような問題が生じないためである。しかしながら、駆動装置では、角度差が零に近いときに磁束電流を急激に変動させる必要がある。本願発明者は、そのような事例を詳しく検証する中で、磁束電流を変化させることにより速度推定誤差が生じることを発見した。そして、本願発明者は、このような速度推定誤差がある状態でも、2台の同期電動機を安定に並列駆動するために、過渡的な速度推定誤差の影響を排除する手法を確立する必要があると考え、サブ側同期電動機1bのトルク電流の脈動成分から磁束電流指令を演算する手法を考案するに至った。本願発明者の調査の結果、サブ側トルク電流の脈動成分から磁束電流指令を決定した場合、特許文献1に開示される方式に比べて、磁束電流指令のS/N比が大きく改善されることが判明した。その理由は以下のとおりである。
 前述のとおり、磁束電流変化により推定速度信号には多くの誤差成分が発生する。この影響を回避する手段として推定磁極位置信号に着目する。推定磁極位置信号は、その計算プロセスにおける積分処理によって、推定速度に含まれる誤差のうちの高周波成分が除去されている。推定磁極位置信号の低周波成分に関しては磁束電流指令の変化による誤差信号が残存するが、その誤差は数度程度に収まる。
 ここで、サブ側同期電動機1bのトルク電流について考えてみる。相電流を磁極位置の真値で座標変換した場合と、磁極位置の推定値で座標変換した場合とでは、磁極位置の誤差が数度程度であれば、真のdq軸におけるトルク電流と推定されるdq軸におけるトルク電流との誤差は、数%未満である。これは、cosine関数が零に近いときには1に近似できることを考えれば、自明である。
 このように、サブ側同期電動機1bのトルク電流は、磁束電流を変化させた場合でも比較的高精度に求めることができる。モータ定数の変動の影響によって定常的な位置推定誤差が生じる場合はあるが、これは直流分の推定誤差となるため、脈動成分抽出部7により脈動成分を抽出する場合には問題とならない。
 速度差ではなく、サブ側トルク電流を用いて安定化補償を行った方が良い理由は他にもある。ファン、ブロワーなどの流体利用装置では、機械系の慣性モーメントが大きい場合がある。このような場合、インバータが過電流停止するほどにトルク脈動が大きいときでも、速度信号に現れる脈動成分はごく小さいことがある。この場合、速度脈動が観測できるほど大きくなってから安定化を図るのではなく、トルク電流の脈動がある程度大きくなった段階で安定化補償を行った方が良い。このような事例では、推定速度信号よりもトルク電流信号の方がS/N比が良いので、トルク電流信号を安定化補償に用いた方が良いと言える。
 以上の理由により、実施の形態1では、磁束電流指令決定部9により、サブ側同期電動機1bに流れるトルク電流の脈動成分を基に磁束電流指令が決定される。
 図19は図1に示す磁束電流指令決定部の構成例を示す図である。図20は図19に示す符号判定器による符号判定処理を説明するための第1の図である。図21は図19に示す符号判定器による符号判定処理を説明するための第2の図である。図22は図19に示す符号判定器による符号判定処理を説明するための第3の図である。
 図19に示す磁束電流指令決定部9は、脈動抑制制御部91及び補償方向決定部92を備える。磁束電流指令決定部9は、サブ側同期電動機1bのトルク電流の脈動成分を入力とし、脈動抑制制御部91と補償方向決定部92を用いて、磁束電流指令を決定する。脈動抑制制御部91は、ゲイン乗算部911及び位相調整部912により構成される。
 ゲイン乗算部911は、入力信号であるトルク電流脈動成分のゲインを調整する。位相調整部912は、入力信号であるトルク電流脈動成分の位相を調整し、振幅が調整された脈動成分を出力する。なお、ゲイン乗算部911と位相調整部912の何れか一方だけで系の安定性を確保できるのであれば、脈動抑制制御部91は、必ずしもゲイン乗算部911及び位相調整部912の双方を備える必要は無い。
 ゲイン乗算部911は、入力信号であるトルク電流脈動成分に特定のゲインを乗算して出力するものであり、系の安定性と即応性を調節する役割を持つ。ゲインは動作条件に応じて変更しても良い。例えば、低速域ではゲインを高くし、高速域ではゲインを低くしても良い。位相調整部912は例えば、位相遅れ補償器、ローパスフィルタ、積分制御器などで構成される。位相遅れ補償器は高周波域でゲインを一定値下げて安定化を図るものであり、産業界で一般的に用いられている。ローパスフィルタ及び積分制御器にも高周波域の信号位相を変化させる性質があるため、位相遅れ補償器と同じようにローパスフィルタ又は積分制御器を用いることができる。
 1次ローパスフィルタによる近似積分器を位相調整部912として使用する場合、そのカットオフ角周波数は、電機ばね共振角周波数の1/3以下に設定すると良い。可能であれば電機ばね共振角周波数の1/10から1/20の値とする。このように設定すると、電機ばね共振角周波数付近で位相を90度前後遅らせることができ、制御安定性が高まる。
 図19には示されていないが、脈動抑制制御部91の入出力の何れかに不感帯を設けても良い。この不感帯は前述の脈動成分抽出部7で除去しきれなかった電機ばね共振以外の周波数成分を除去するのに役立つ。
 補償方向決定部92は、符号判定器921及び乗算器922により構成され、図5から図11で説明した動作原理に則り、2台の同期電動機のそれぞれが備える回転子の磁極位置の角度差から、磁束電流指令の補償方向を決定する。符号判定器921は図20から図22に示される符号判定処理を行う。図20から図22の横軸は、符号判定器921の入力である角度差を表す。角度差は図11に示すように正又は負の値を示す。図20から図22の縦軸は符号判定器921の出力の値を示す。
 最も基本的な符号判定処理は図20に示す方法である。符号判定器921は、角度差が正を示す場合には「1」を出力し、角度差が負を示す場合には「-1」を出力する。但し、図20の方法では角度差が零に近いときにチャタリングが発生するおそれがある。そのため、図21のように角度差が零に近い領域では、符号判定器921の出力を「1」から「-1」へ徐々に切り替え、又は符号判定器921の出力を「-1」から「1」へ徐々に切り替えるように構成してもよい。
 また、角度差が大きい条件では、角度差が小さい条件に比べて、メイン側同期電動機1aの磁束電流の変化量がサブ側同期電動機1bのトルクの変化量に及ぼす影響が大きくなる。そのため、図22のように角度差が大きい領域では、角度差が大きくなるに従って、符号判定器921の出力値を徐々に下げるようにしても良い。
 乗算器922は符号判定器921の出力と脈動抑制制御部91の出力とを掛け合わせ、磁束電流指令を生成する。すなわち、磁束電流指令決定部9では、脈動抑制制御部91により抑制された脈動成分と、補償方向決定部92の符号判定器921により決定された補償方向とにより、磁束電流指令を決定する。このように生成された磁束電流指令を用いることによる効果は以下の通りである。
 ここまで説明してきたとおり、位置センサレス制御による駆動装置では、低速駆動時に磁束電流指令の変化によって過渡的な速度推定誤差が生じる。特許文献1に開示されるように、速度差を用いてサブ側同期電動機1bの駆動を安定化させる方法では、速度推定誤差の影響を直接的に受けてしまうため、制御不安定な状態となり、騒音及び振動が増加し、モータ効率が低下するなどの様々な問題が生じる。そのため、特許文献1では、1つの電力変換器で1台の同期電動機を駆動する既存の同期電動機駆動装置に比べて、回転数の下限を高くしなければならないという課題があった。そのため、既存の同期電動機駆動装置から、特許文献1で開示される技術を用いた並列駆動装置への置き換えは困難である。
 これに対して、実施の形態1に係る駆動装置100は、サブ側同期電動機1bのトルク電流の脈動成分と、2台の同期電動機の磁極位置の角度差とを用いて、磁束電流指令を決定するように構成されている。これにより、磁束電流指令のS/N比が改善され、速度推定誤差の影響を受け難くなる。その結果、騒音及び振動が増加し、モータ効率が低下し、脱調するなどの問題が解決される。また、低速駆動時の安定性が改善されるため、回転数の下限値は、1つの電力変換器で1台の同期電動機を駆動する既存の同期電動機駆動装置と同等の値を維持できる。これにより、既存の同期電動機駆動装置を実施の形態1に係る駆動装置100に置き換えることが容易となる。
実施の形態2.
 実施の形態2では、サブ側同期電動機1bで消費される有効電力の脈動成分を利用して磁束電流を決定する構成例について説明する。騒音及び振動が大きくなり、モータ効率が低下するなどの課題を解決するには、電機ばね共振によるサブ側同期電動機1bの自己発振現象を、磁束電流が大きく変化する条件でも正確に検出する必要がある。そのための1つの手法が実施の形態1で述べたトルク電流の脈動成分を用いる手法である。但し、メイン側同期電動機1a及びサブ側同期電動機1bに接続される負荷である機械系の慣性モーメントが比較的大きい場合、トルク電流の脈動成分の代わりに、有効電力の脈動成分を用いて磁束電流を決定しても良い。前述したとおり、ファン、ブロワーなどの流体利用装置では機械系の慣性モーメントが大きい場合があるため、このような場合、推定速度信号を観測するよりも、有効電力の脈動成分を観測した方が良い。
 図23は本発明の実施の形態2に係る駆動装置の構成を示す図である。実施の形態2に係る駆動装置100Aは、図1に示す脈動成分抽出部70の代わりに、脈動成分抽出部70Aを備える。脈動成分抽出部70Aは、サブ側有効電力脈動成分抽出部10を備える。その他の構成については、実施の形態1の構成と同一又は同等であり、同一又は同等の構成部には同一の符号を付して、重複する説明は省略する。
 有効電力Pは三相電圧指令v ,v ,v と相電流i ,i ,i を用いて下記(18)式により求めることができる。添字「」はメイン側とサブ側を区別するためのものである。Rは電機子抵抗である。三相電圧指令v ,v ,v は、電流制御部6から得られる。相電流i ,i ,i は電流検出部4bから得られる。
Figure JPOXMLDOC01-appb-M000018
 上記(18)式の右辺の第二項は電機子抵抗による銅損を表す。厳密には銅損も有効電力の一部であるが、ここで知りたい情報はサブ側同期電動機1bのトルク脈動に相当するものであるから、銅損を差し引いておいた方が良い。但し、電機子抵抗が無視できる程度に小さい場合もあるため、その場合は右辺第一項のみを計算するようにしても構わない。
 機械系の慣性モーメントが大きく、かつ、速度脈動が微小である場合、有効電力の脈動はトルクの脈動によって生じたものであると考えられる。従って、この場合、サブ側有効電力脈動成分抽出部10は、実施の形態1で説明したサブ側トルク電流脈動成分抽出部7と同様の演算処理を行い、サブ側同期電動機1bの有効電力から脈動成分を抽出する。この情報を用いて磁束電流指令を決定すれば、実施の形態1と同様の効果が得られる。
 なお、実施の形態2の磁束電流指令決定部9は、実施の形態1と同様にゲイン乗算部911及び位相調整部912により構成される脈動抑制制御部91を備えるが、実施の形態2のゲイン乗算部911は、入力信号である有効電力脈動成分のゲインを調整し、実施の形態2の位相調整部912は、入力信号である有効電力脈動成分の位相を調整する。また、実施の形態2の脈動抑制制御部91は、実施の形態1と同様に、ゲイン乗算部911と位相調整部912の何れか一方だけで系の安定性を確保できるのであれば、必ずしもゲイン乗算部911及び位相調整部912の双方を備える必要は無い。
 実施の形態2は、同期電動機に接続される負荷の慣性モーメントが大きい場合に有用であり、実施の形態1に比べて、座標変換の演算を用いない分、計算量が少ないため、簡易な構成の演算装置を用いる場合には有用である。具体的には、実施の形態1のサブ側トルク電流脈動成分抽出部7は、電流検出部4bで検出された三相座標系の電流を、磁極位置推定部5bからの信号を用いて、回転する直交座標系に座標変換して、トルク電流を求め、そのトルク電流の脈動成分を抽出している。これに対して、実施の形態2では、サブ側有効電力脈動成分抽出部10が、上記(18)式のように、電流検出部4bで検出された三相座標系の電流をそのまま用いて有効電力を求め、この有効電力の脈動成分を抽出している。そして実施の形態2の磁束電流指令決定部9は、この脈動成分を利用して磁束電流指令を決定することができる。そのため実施の形態2では、座標変換が不要となり、計算量が少なくなる。慣性モーメントが大きいアプリケーションの場合、上記のように有効電力の脈動成分を観測することで、座標変換が1回不要になるため、実施の形態2では演算負荷を減らすことが可能となる。
実施の形態3.
 図24は本発明の実施の形態3に係る駆動装置の構成を示す図である。実施の形態3に係る駆動装置100Bは、図1に示す脈動成分抽出部70の代わりに、脈動成分抽出部70Bを備える。脈動成分抽出部70Bは、サブ側トルク電流脈動成分抽出部7と、メイン側トルク電流脈動成分抽出部11と、第2の減算器である減算器8aとを備える。その他の構成については、実施の形態1の構成と同一又は同等であり、同一又は同等の構成部には同一の符号を付して、重複する説明は省略する。
 実施の形態3では、サブ側同期電動機1bのトルク電流の脈動成分とメイン側同期電動機1aのトルク電流の脈動成分の差から磁束電流指令を決定する構成例について説明する。実施の形態1,2の駆動装置100,100Aは、メイン側同期電動機1aのトルク電流の脈動成分が定常状態では微小であることを前提とした構成となっている。メイン側同期電動機1aはベクトル制御されているので、トルク電流指令値が一定値であればメイン側同期電動機1aのトルク電流は指令値に追従するはずである。しかしながら、現実には様々な外乱要因によってメイン側同期電動機1aのトルク電流は脈動する。外乱要因としては、電力変換器2を構成する直列の上下アームの半導体素子の短絡防止時間、電流センサのオフセット、電流センサのゲインアンバランス、回転子に設けられる磁石から発生する磁束の歪みなどが考えられる。これらの要因によるトルク電流脈動は、サブ側同期電動機1bでも同様に発生する。また、実施の形態1,2では、トルク電流指令に何らかの交流成分が重畳されている場合、これも磁束電流指令決定部9にとっては外乱となる。トルク電流指令に重畳された外乱成分によって、サブ側同期電動機1bにもその周波数のトルク電流脈動が発生するが、当該外乱成分は、電機ばね共振による自己発振とは異なる原因で発生したものであるため、当該外乱成分を磁束電流指令決定部にフィードバックするのは適切でない。
 実施の形態1のサブ側トルク電流脈動成分抽出部7と実施の形態3のメイン側トルク電流脈動成分抽出部11はハイパスフィルタ、バンドパスフィルタなどにより構成される。サブ側同期電動機1bをより的確に安定化させるためには、上記の外乱要因の影響を排除することが望ましい。しかしながら、ハイパスフィルタでは外乱の除去特性が悪く、バンドパスフィルタでも外乱の除去特性を良くするためには電機ばね共振角周波数の計測が必要となる。こういった事情から、より簡易な方法で外乱の影響を除去するため、実施の形態3は、メイン側同期電動機1aで発生するトルク電流の脈動成分を、サブ側同期電動機1bのトルク電流の脈動成分から差し引くよう構成したのである。
 そのため実施の形態3に係る駆動装置100Bは、サブ側トルク電流脈動成分抽出部7に加えて、メイン側トルク電流脈動成分抽出部11を備える。また駆動装置100Bは、サブ側トルク電流脈動成分抽出部7からのトルク電流脈動成分とメイン側トルク電流脈動成分抽出部11からのトルク電流脈動成分との差分を求める減算器8aを備える。
 メイン側トルク電流脈動成分抽出部11は、メイン側同期電動機1aのトルク電流の脈動成分を計算する。計算法は実施の形態1で述べたサブ側トルク電流脈動成分抽出部7と同様の方法で良い。減算器8aでは、2台の同期電動機のそれぞれに発生するトルク電流の脈動成分の差分が計算され、磁束電流指令決定部9は、当該差分を用いて磁束電流指令を決定する。
 このように構成することで、より安定に2台の同期電動機を並列駆動することが可能となる。なお、実施の形態3では、2台の同期電動機のそれぞれに発生するトルク電流の脈動成分の差分を用いる方法を説明したが、その代わりに有効電力の脈動成分の差分を用いても良いことは言うまでもない。また、実施の形態3の駆動装置100Bは、サブ側トルク電流脈動成分抽出部7及びメイン側トルク電流脈動成分抽出部11の内、一方がトルク電流の脈動を計算し、他方が有効電力の脈動を計算した後、これらの脈動を同一のスケールに換算してから差分を求めるように構成してもよい。
実施の形態4.
 実施の形態4では、実施の形態1,2,3に係る駆動装置100,100A,100Bを用いた流体利用装置の構成例について説明する。図25は本発明の実施の形態4に係る流体利用装置の構成図である。実施の形態4では、メイン側同期電動機1aの回転軸にプロペラファン300aが設けられ、サブ側同期電動機1bの回転軸にプロペラファン300bが設けられている流体利用装置300について説明する。
 図25に示す流体利用装置300は、実施の形態1の駆動装置100を備え、駆動装置100は電力変換器駆動装置200を備える。電力変換器駆動装置200は、プロセッサ201及びメモリ202を備える。図1に示す各機能、すなわち電流制御部6、磁極位置推定部5a,5b、サブ側トルク電流脈動成分抽出部7、減算器8及び磁束電流指令決定部9は、プロセッサ201及びメモリ202を用いてその機能が実現される。
 図25に示すようにプロセッサ201及びメモリ202を利用する場合、上記の各機能のそれぞれは、ソフトウェア、ファームウェア又はこれらの組合せにより実現される。ソフトウェア又はファームウェアはプログラムとして記述され、メモリ202に記憶される。プロセッサ201はメモリ202に記憶されたプログラムを読み出して実行する。またこれらのプログラムは、上記の各機能のそれぞれが実行する手順及び方法をコンピュータに実行させるものであるとも言える。メモリ202は、RAM(Random Access Memory)、ROM(Read Only Memory)、フラッシュメモリー、EPROM(Erasable Programmable Read Only Memory)、又はEEPROM(Electrically Erasable Programmable Read Only Memory)(登録商標)といった半導体メモリが該当する。半導体メモリは不揮発性メモリでもよいし揮発性メモリでもよい。またメモリ202は、半導体メモリ以外にも、磁気ディスク、フレキシブルディスク、光ディスク、コンパクトディスク、ミニディスク又はDVD(Digital Versatile Disc)が該当する。なお、プロセッサ201は、演算結果等のデータをメモリ202に出力しても記憶させても良いし、メモリ202を介して不図示の補助記憶装置に当該データを記憶させてもよい。
 なお、流体利用装置300は、駆動装置100の代わりに実施の形態2の駆動装置100A又は実施の形態3の駆動装置100Bを備えてもよい。この場合、図23に示す電流制御部6、磁極位置推定部5a,5b、サブ側有効電力脈動成分抽出部10、減算器8及び磁束電流指令決定部9は、プロセッサ201及びメモリ202を用いてその機能が実現される。また、図24に示す電流制御部6、磁極位置推定部5a,5b、サブ側トルク電流脈動成分抽出部7、メイン側トルク電流脈動成分抽出部11、減算器8、減算器8a及び磁束電流指令決定部9は、プロセッサ201及びメモリ202を用いてその機能が実現される。
 実施の形態1でも述べたように、電力変換器2はメイン側同期電動機1a及びサブ側同期電動機1bに任意の交流電力を供給できるものであれば、基本的にどのような回路構成でも構わない。電流検出部4a,4bで検出された電流の情報はプロセッサ201へ送信される。
 2つのプロペラファン300a,300bは、互いに同一形状のものであっても良いし、異なる形状のものであっても良い。また、2つのプロペラファン300a,300bの空気の流路は必ずしも同じでなくとも良い。例えば流体利用装置300が空気調和機の場合、2つのプロペラファン300a,300bは、当該空気調和機の室外機内の送風室に設けられる2つの送風ファンに相当し、上記の空気の流路は、当該送風室に相当する。送風室は、室外機の側面板、天井板、底板、熱交換器などに囲まれることで形成される空間である。送風室には、プロペラファン300a,300bが回転することによって空気の流れが形成される。
 2つのプロペラファン300a,300bの回転数と負荷トルクとの特性は、異なっていた方が安定に並列駆動しやすいため、2台の同期電動機に異なる形状のファンを設けてもよいし、一方のファンが設けられる流路の断面積を、他方のファンが設けられる流路の断面積よりも小さくしても良い。また、一方の同期電動機でプロペラファンを駆動し、他方の同期電動機でポンプを駆動するなど、それぞれ異なる仕様の流体利用装置を駆動する構成としても良い。
 なお図25には示されていないが、流体利用装置300は、電力変換器2が出力する電圧を検出する電圧検出部を備え、電圧検出部で検出された電圧情報がプロセッサ201へ入力されるように構成してもよい。また図25には示されていないが、流体利用装置300は、ファンの風速を計測する風速センサを備え、風速センサで検出された風速情報がプロセッサ201へ入力されるように構成してもよい。また図25には示されていないが、ファンによって冷却される対象物の温度を検出する温度センサを備え、温度センサで検出された温度情報がプロセッサ201へ入力されるように構成してもよい。
 流体利用装置300の流体負荷は、ダンパ特性を持っており、高回転域では、そのダンパ特性がオープンループ駆動された同期電動機の駆動を安定化させる。しかしながら、低回転域では、そのダンパ特性が弱まり、同期電動機の駆動が不安定になるため、流体利用装置300は実施の形態1,2,3で述べた並列駆動法を利用する。これにより、実施の形態4では、幅広い速度範囲で同期電動機の並列駆動を実現することができる。また、実施の形態4では、高度なトルク制御の必要がないため、1つの電力変換器で1台の同期電動機を駆動する既存の同期電動機駆動装置を改修することによって、コストの増加を抑制しながら2つのプロペラファン300a,300bを駆動可能な流体利用装置300を得ることができる。
実施の形態5.
 実施の形態5では、実施の形態4に係る流体利用装置300を用いた空気調和機の構成例について説明する。図26は本発明の実施の形態5に係る空気調和機の構成図である。実施の形態5に係る空気調和機400は、流体利用装置300、冷媒圧縮機401、凝縮器403、受液器404、膨張弁405及び蒸発器406を備える。冷媒圧縮機401と凝縮器403との間は配管で接続される。同様に、凝縮器403と受液器404との間は配管で接続され、受液器404と膨張弁405との間は配管で接続され、膨張弁405と蒸発器406との間は配管で接続され、蒸発器406と冷媒圧縮機401との間は配管で接続される。これにより、冷媒圧縮機401、凝縮器403、受液器404、膨張弁405及び蒸発器406には冷媒が循環する。なお、図26では図示が省略されているが、流体利用装置300は、図1などに示す電流検出部4a,4b、磁極位置推定部5a,5cなどを備える。
 空気調和機400では冷媒の蒸発、圧縮、凝縮、膨張という工程が繰り返し行われるため、冷媒は液体から気体へ変化し、さらに気体から液体へ変化することにより、冷媒と機外空気との間で熱交換が行われる。
 蒸発器406は、低圧の状態で冷媒液を蒸発させ、蒸発器406の周囲の空気から熱を奪うことによって、冷却作用を発揮するものである。冷媒圧縮機401は、冷媒を凝縮するために蒸発器406でガス化された冷媒ガスを圧縮して、高圧のガスにするものである。凝縮器403は、冷媒圧縮機401で高温になった冷媒ガスの熱を放出することで、高圧の冷媒ガスを凝縮し、冷媒液に変換するものである。流体利用装置300は、プロペラファン300a,300bを回転することによって、風を発生させて、この風を凝縮器403に通過させることにより、凝縮器403を冷却する。膨張弁405は、冷媒を蒸発させるために、冷媒液を絞り膨張して、冷媒液を低圧の液に変換するものである。受液器404は循環する冷媒量の調節のために設けられるもので、小型の装置では省略しても良い。
 空気調和機400の大出力化に伴って凝縮器403が大型化すると、凝縮器403を冷却するための冷却装置として機能する流体利用装置300の冷却性能を増加させる必要が生じる。但し、凝縮器403の寸法を大きくするのに合わせて、冷却装置として機能する流体利用装置300の仕様変更を行うのは煩雑である。また、流体利用装置300の冷却性能を増加させるために、流体利用装置300を大出力化するためには、流体利用装置300を量産するための製造ラインの変更が必要となる場合もあり、製造ラインを構築するための初期投資がかさむ。そのため、大型の空気調和機400では、複数の冷却ファンを備えた流体利用装置300を使用することで冷却性能を向上させている。
 また、空気調和機400には、低コスト化の要求が高く、その一方で省エネルギー規制が年々強化されているため高効率化も要求されている。近年の省エネルギー規制では、定格動作点だけでなく、低出力駆動の動作点での駆動効率も重要視される。そのために冷却ファンの動作回転数の下限値を極力引き下げる必要がある。
 ここまで述べてきたとおり、特許文献1で開示される技術を用いた並列駆動装置は、コスト面では非常に優れているものの、これを位置センサレス制御で構成しようとすると、低回転域での駆動が不安定となる課題がある。特に空気調和機400では、流体利用装置300の電力変換器2で発生するキャリア騒音を減らすため、キャリア周波数を10kHz以上に高く設定する場合が多く、出力電圧誤差が増加して、低回転域での駆動が不安定になりやすい。従って、空気調和機400に、特許文献1で開示される技術を用いた並列駆動装置を用いた場合、空気調和機400の駆動動作範囲が狭まるという課題があった。そのため、特許文献1で開示される技術を用いた並列駆動装置では、空気調和機400用の冷却ファンに要求される低コスト化と高冷却性能とを両立させることが困難である。
 実施の形態5に係る空気調和機400は、実施の形態1から3で述べた並列駆動法を利用しているため、低速域の駆動が不安定になることがなく、駆動可能範囲を拡大できる。また実施の形態1から3で述べた並列駆動法は位置センサレス制御を前提としているため、位置センサを用いる場合に比べて、空気調和機400の製造コストを低減できる。従って、実施の形態1から3に係る駆動装置100,100A,100Bでは、空気調和機400用の冷却ファンに要求される低コスト化と高冷却性能とを両立させることが可能である。
 以上の実施の形態に示した構成は、本発明の内容の一例を示すものであり、別の公知の技術と組み合わせることも可能であるし、本発明の要旨を逸脱しない範囲で、構成の一部を省略、変更することも可能である。
 1a メイン側同期電動機、1b サブ側同期電動機、1c 同期電動機、2 電力変換器、3 直流電圧源、4a,4b 電流検出部、5a,5b,5c 磁極位置推定部、6 電流制御部、7,7A,7B,7C サブ側トルク電流脈動成分抽出部、8,8a 減算器、9 磁束電流指令決定部、10 サブ側有効電力脈動成分抽出部、11 メイン側トルク電流脈動成分抽出部、51 モデル偏差演算部、52 角速度推定器、53 一次角周波数演算器、54 積分器、70,70A,70B 脈動成分抽出部、71 脈動周波数計測部、72 余弦波発生器、73 正弦波発生器、74 フーリエ余弦係数演算部、75 フーリエ正弦係数演算部、76 交流復元器、91 脈動抑制制御部、92 補償方向決定部、100,100A,100B 駆動装置、200 電力変換器駆動装置、201 プロセッサ、202 メモリ、300 流体利用装置、300a,300b プロペラファン、400 空気調和機、401 冷媒圧縮機、403 凝縮器、404 受液器、405 膨張弁、406 蒸発器、511 電流推定器、512 減算器、513 偏差演算器、911 ゲイン乗算部、912 位相調整部、921 符号判定器、922 乗算器。

Claims (15)

  1.  並列接続される第1の同期電動機及び第2の同期電動機に電力を供給する電力変換器と、
     前記第1の同期電動機に流れる第1の電流を検出する第1の電流検出器と、
     前記第2の同期電動機に流れる第2の電流を検出する第2の電流検出器と、
     前記第1の同期電動機を駆動するための電圧指令と、前記第1の電流とを用いて、前記第1の同期電動機が有する回転子の第1の磁極位置を推定する第1の磁極位置推定部と、
     前記電圧指令と、前記第2の電流とを用いて、前記第2の同期電動機が有する回転子の第2の磁極位置を推定する第2の磁極位置推定部と、
     トルク電流指令と磁束電流指令と前記第1の電流と前記第1の磁極位置とを用いて前記電圧指令を出力する制御部と、
     前記第2の同期電動機に流れるトルク電流に含まれるトルク電流脈動成分と前記第2の同期電動機で消費される有効電力に含まれる有効電力脈動成分との少なくとも一方の脈動成分を抽出する脈動成分抽出部と、
     前記第1の磁極位置と前記第2の磁極位置とを用いて、前記第1の同期電動機及び前記第2の同期電動機がそれぞれ有する回転子の磁極位置の差である角度差を求める第1の減算器と、
     前記脈動成分抽出部で抽出された脈動成分と、前記角度差とを用いて、前記磁束電流指令を決定する磁束電流指令決定部と、
     を備えることを特徴とする駆動装置。
  2.  前記脈動成分抽出部は、
     前記第1の同期電動機に流れるトルク電流に含まれるトルク電流脈動成分と、前記第2の同期電動機に流れるトルク電流に含まれるトルク電流脈動成分との差分を求めて、前記脈動成分として出力する第2の減算器を備えることを特徴とする請求項1に記載の駆動装置。
  3.  前記脈動成分抽出部は、
     前記第1の同期電動機で消費される有効電力に含まれる有効電力脈動成分と、前記第2の同期電動機で消費される有効電力に含まれる有効電力脈動成分との差分を求めて、前記脈動成分として出力する第2の減算器を備えることを特徴とする請求項1に記載の駆動装置。
  4.  前記磁束電流指令決定部は、
     前記角度差を用いて前記磁束電流指令の補償方向を決定する補償方向決定部と、
     前記脈動成分を用いて前記脈動成分を抑制する脈動抑制制御部と、
     を備え、
     前記脈動抑制制御部により抑制された前記脈動成分と前記補償方向決定部により決定された前記補償方向とにより、前記磁束電流指令を決定することを特徴とする請求項1から3の何れか一項に記載の駆動装置。
  5.  前記脈動抑制制御部は、
     前記脈動成分の位相を調整し、位相が調整された脈動成分を出力する位相調整部と、前記脈動成分にゲインを乗じて、振幅が調整された脈動成分を出力するゲイン乗算部との少なくとも一方を備えることを特徴とする請求項4に記載の駆動装置。
  6.  前記位相調整部は、前記脈動成分の位相を遅らせる位相遅れ補償器であることを特徴とする請求項5に記載の駆動装置。
  7.  前記位相調整部は、1次ローパスフィルタによる近似積分器であり、
     前記近似積分器のカットオフ角周波数は、電機ばね共振角周波数の1/3以下に設定されることを特徴とする請求項5に記載の駆動装置。
  8.  前記脈動成分抽出部は、
     前記第1の電流又は前記第2の電流をdq座標系に座標変換する座標変換器と、
     座標変換後のq軸電流に含まれる直流成分を除去することにより、前記トルク電流脈動成分を演算するハイパスフィルタと、
     を備えることを特徴とする請求項1から7の何れか一項に記載の駆動装置。
  9.  前記脈動成分抽出部は、
     前記第1の電流又は前記第2の電流をdq座標系に座標変換する座標変換器と、
     座標変換後のq軸電流に含まれる電機ばね共振角周波数成分を抽出することにより、電機ばね共振角周波数を中心周波数とする前記トルク電流脈動成分を演算するバンドパスフィルタと、
     を備えることを特徴とする請求項1から7の何れか一項に記載の駆動装置。
  10.  前記脈動成分抽出部は、
     前記電機ばね共振角周波数を計測する共振周波数計測手段を備え、
     前記バンドパスフィルタの中心周波数を動的に変更することを特徴とする請求項9に記載の駆動装置。
  11.  前記脈動成分抽出部は、
     前記有効電力を計算する有効電力計算手段と、
     前記有効電力に含まれる直流分を除去することにより、前記有効電力の脈動成分を演算するハイパスフィルタと
     を備えることを特徴とする請求項1から7の何れか一項に記載の駆動装置。
  12.  前記脈動成分抽出部は、
     前記有効電力を計算する有効電力計算手段と、
     前記有効電力に含まれる電機ばね共振角周波数成分を抽出することにより、電機ばね共振角周波数を中心周波数とする前記有効電力脈動成分を演算するバンドパスフィルタと、
     を備えることを特徴とする請求項1から7の何れか一項に記載の駆動装置。
  13.  前記脈動成分抽出部は、
     前記電機ばね共振角周波数を計測する共振周波数計測手段を備え、
     前記バンドパスフィルタの中心周波数を動的に変更することを特徴とする請求項12に記載の駆動装置。
  14.  請求項1から13の何れか一項に記載の駆動装置を備えることを特徴とする流体利用装置。
  15.  請求項14に記載の流体利用装置を備えることを特徴とする空気調和機。
PCT/JP2018/001458 2018-01-18 2018-01-18 駆動装置、流体利用装置及び空気調和機 WO2019142301A1 (ja)

Priority Applications (7)

Application Number Priority Date Filing Date Title
KR1020207017409A KR102412170B1 (ko) 2018-01-18 2018-01-18 구동 장치 및 공기 조화기
EP18901287.5A EP3742598B1 (en) 2018-01-18 2018-01-18 Drive device, fluid utilization device, and air conditioner
JP2019565637A JP6877592B2 (ja) 2018-01-18 2018-01-18 駆動装置、流体利用装置及び空気調和機
US16/960,113 US11133762B2 (en) 2018-01-18 2018-01-18 Drive system, fluid utilization apparatus, and air conditioner
PCT/JP2018/001458 WO2019142301A1 (ja) 2018-01-18 2018-01-18 駆動装置、流体利用装置及び空気調和機
CN201880085523.7A CN111630769B (zh) 2018-01-18 2018-01-18 驱动装置、流体利用装置以及空调机
JP2021072758A JP7113939B2 (ja) 2018-01-18 2021-04-22 駆動装置及び空気調和機

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
PCT/JP2018/001458 WO2019142301A1 (ja) 2018-01-18 2018-01-18 駆動装置、流体利用装置及び空気調和機

Publications (1)

Publication Number Publication Date
WO2019142301A1 true WO2019142301A1 (ja) 2019-07-25

Family

ID=67302110

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2018/001458 WO2019142301A1 (ja) 2018-01-18 2018-01-18 駆動装置、流体利用装置及び空気調和機

Country Status (6)

Country Link
US (1) US11133762B2 (ja)
EP (1) EP3742598B1 (ja)
JP (2) JP6877592B2 (ja)
KR (1) KR102412170B1 (ja)
CN (1) CN111630769B (ja)
WO (1) WO2019142301A1 (ja)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP3998702B1 (en) * 2019-07-09 2024-01-03 Mitsubishi Electric Corporation Driving device, fluid utilization device, and air conditioner
US11349416B2 (en) * 2020-04-22 2022-05-31 Steering Solutions Ip Holding Corporation Feedback compensation of parameter imbalance induced current harmonics in synchronous motor drives
CN113312745B (zh) * 2021-04-13 2022-08-30 国网西藏电力有限公司经济技术研究院 受控电压源的静态负荷电磁暂态建模电路和建模仿真方法
CN114257109B (zh) * 2021-12-21 2023-07-14 郑州轻工业大学 一种npc变换器电流传感器故障模型预测控制方法
US20230327583A1 (en) * 2022-04-08 2023-10-12 Abb Schweiz Ag Parallel synchronous machines with single motor drive

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2003116293A (ja) * 2001-10-05 2003-04-18 Fuji Electric Co Ltd Dcブラシレスモータの並列駆動回路
US20070273310A1 (en) * 2006-03-24 2007-11-29 Airbus France Power system comprising several synchronous machines synchronously self-controlled by a converter and control method for such a system
JP2012196091A (ja) * 2011-03-17 2012-10-11 Fujitsu Telecom Networks Ltd モータ駆動電流アンプとモータ駆動装置とその駆動方法とモータ駆動システム
JP2013141345A (ja) * 2011-12-28 2013-07-18 Mitsubishi Heavy Ind Ltd モータ制御装置及び空気調和機
JP2014217226A (ja) * 2013-04-26 2014-11-17 東洋電機製造株式会社 電動機制御装置
US20150229245A1 (en) 2014-02-11 2015-08-13 Samsung Electronics Co., Ltd. Power apparatus, controlling method thereof and motor driving apparatus therein

Family Cites Families (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP4032370B2 (ja) * 1998-12-02 2008-01-16 株式会社名機製作所 同期電動機及び同期電動機制御装置
JP4958431B2 (ja) 2005-12-09 2012-06-20 日立アプライアンス株式会社 電動機制御装置
JP5326429B2 (ja) * 2008-08-26 2013-10-30 株式会社明電舎 電動機の脈動抑制装置
KR102273139B1 (ko) * 2014-02-11 2021-07-05 삼성전자주식회사 동력 장치 및 동력 장치에 포함되는 전동기 구동 장치
JP6364971B2 (ja) 2014-06-03 2018-08-01 株式会社リコー 画像形成装置
JP6308894B2 (ja) * 2014-07-08 2018-04-11 三菱電機株式会社 三相交流回転機の制御装置
KR101623652B1 (ko) * 2014-10-29 2016-05-23 엘지전자 주식회사 모터제어장치
US9634591B2 (en) * 2015-08-17 2017-04-25 Infineon Technologies Ag System and method for motor control using position sensors
US10389286B2 (en) * 2016-03-07 2019-08-20 Schlumberger Technology Corporation Methods and systems for control of parallel connected motor drives
CN106533298A (zh) * 2016-12-24 2017-03-22 天津工业大学 双永磁同步电机驱动系统转速同步控制方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2003116293A (ja) * 2001-10-05 2003-04-18 Fuji Electric Co Ltd Dcブラシレスモータの並列駆動回路
US20070273310A1 (en) * 2006-03-24 2007-11-29 Airbus France Power system comprising several synchronous machines synchronously self-controlled by a converter and control method for such a system
JP2012196091A (ja) * 2011-03-17 2012-10-11 Fujitsu Telecom Networks Ltd モータ駆動電流アンプとモータ駆動装置とその駆動方法とモータ駆動システム
JP2013141345A (ja) * 2011-12-28 2013-07-18 Mitsubishi Heavy Ind Ltd モータ制御装置及び空気調和機
JP2014217226A (ja) * 2013-04-26 2014-11-17 東洋電機製造株式会社 電動機制御装置
US20150229245A1 (en) 2014-02-11 2015-08-13 Samsung Electronics Co., Ltd. Power apparatus, controlling method thereof and motor driving apparatus therein

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
JUNICHI ITOJIRO TOYOSAKIHIROSHI OSAWA: "High performance V/f control method for PM Motor", IEEJ TRANSACTIONS ON INDUSTRY APPLICATIONS, vol. 122, no. 3, 2002, pages 253 - 259
See also references of EP3742598A4
YONGJAE LEEJUNG-IK HA: "Analysis and Control of Mono Inverter Dual Parallel SPMSM Drive System", 2014, IEEE
YONGJAE LEEJUNG-IK HA: "Control Method for Mono Inverter Dual Parallel Surface-Mounted Permanent-Magnet Synchronous Machine Drive System", 2015, IEEE

Also Published As

Publication number Publication date
KR20200090197A (ko) 2020-07-28
US11133762B2 (en) 2021-09-28
CN111630769B (zh) 2023-04-25
EP3742598A1 (en) 2020-11-25
JPWO2019142301A1 (ja) 2020-09-24
KR102412170B1 (ko) 2022-06-23
JP7113939B2 (ja) 2022-08-05
JP2021106499A (ja) 2021-07-26
EP3742598B1 (en) 2021-08-18
EP3742598A4 (en) 2020-12-16
JP6877592B2 (ja) 2021-05-26
CN111630769A (zh) 2020-09-04
US20210067066A1 (en) 2021-03-04

Similar Documents

Publication Publication Date Title
WO2019142301A1 (ja) 駆動装置、流体利用装置及び空気調和機
JP6343037B1 (ja) モータ駆動装置および冷凍機器
WO2011105355A1 (ja) 回転電気機械のトルクリプル抑制制御装置および制御方法
JP6537725B2 (ja) 交流電動機の速度推定装置、交流電動機の駆動装置、冷媒圧縮機および冷凍サイクル装置
TW200524265A (en) Controller for synchronous motor, electric appliance and module
CN111212980B (zh) 用于压缩机应用的闭环扭矩补偿
JP4644010B2 (ja) 同期リラクタンス電動機のベクトル制御方法及び同装置
JP5488043B2 (ja) モータのトルク制御装置
JP6833071B2 (ja) 駆動装置、流体利用装置及び空気調和機
WO2020202655A1 (ja) 永久磁石同期機の駆動装置、永久磁石同期機のトルク補償方法、および電気車
JP5644203B2 (ja) 回転電気機械のトルクリプル抑制制御装置および制御方法
JP4906819B2 (ja) 圧縮機およびトルク制御装置並びに空気調和機
JP2015042010A (ja) モータ駆動装置、モータ駆動モジュール、圧縮機及びこれを備えた冷凍装置
CN107425761B (zh) 永磁同步电机低频控制方法、装置、压缩机以及空调器
JP3968688B2 (ja) 交流電動機のベクトル制御方法及び同装置
JP7118274B2 (ja) 駆動装置、流体利用装置及び空気調和機
JP3968689B2 (ja) 同期リラクタンス電動機のベクトル制御方法及び同装置
WO2021144856A1 (ja) 交流電動機の駆動装置、圧縮機駆動装置及び冷凍サイクル装置
TWI662782B (zh) 馬達驅動裝置、及具備其之冷凍循環裝置、以及馬達驅動方法
JP2023089904A (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: 18901287

Country of ref document: EP

Kind code of ref document: A1

ENP Entry into the national phase

Ref document number: 2019565637

Country of ref document: JP

Kind code of ref document: A

ENP Entry into the national phase

Ref document number: 20207017409

Country of ref document: KR

Kind code of ref document: A

NENP Non-entry into the national phase

Ref country code: DE

ENP Entry into the national phase

Ref document number: 2018901287

Country of ref document: EP

Effective date: 20200818