The present invention relates to a power conversion device.
As a technology for controlling load torque pulsation of a compressor system or the like by position sensorless control, Patent Document 1 is known. Patent Document 1 discloses correction of a torque current command value that is an output of speed control by using a pulsating component extracted by a resonance filter that extracts a periodic pulsating component of a difference between a motor speed and a speed command.
The technology described in Patent Document 1 can suppress pulsation in the motor speed by adding an output of a resonance filter and speed control, but does not mention suppression of pulsation in torque.
An object of the invention is to provide a power conversion device that realizes stable, highly accurate, and highly efficient control characteristics by suppressing pulsation in a speed of a driven motor, and pulsation in torque.
According to an aspect of the invention, there is provided a power conversion device that drives a magnet motor. The power conversion device includes: a first filter unit that includes a primary system transfer function filter that varies in correspondence with a pulsating component of the magnet motor, and computes a second speed estimation value from a speed command value and a first speed estimation value;
Power to be supplied to the magnet motor is controlled on the basis of the second torque current command value.
According to the invention, stable, highly accurate, and highly efficient control characteristics are realized by suppressing pulsation in a speed of a driven motor and pulsation in torque.
Hereinafter, examples will be described in detail with reference to the accompanying drawings. In addition, in the drawings, the same reference numeral is given to a common configuration. In addition, the following examples are not limited to configurations illustrated in the drawings.
A magnet motor 1 outputs motor torque that is a combination of a torque component due to a magnetic flux of a permanent magnet and a torque component due to inductance of armature winding.
A power converter 2 outputs voltage values proportional to three-phase AC voltage command values vu*, vv*, vw*, and varies an output voltage value and an output frequency value of the magnet motor 1.
A DC power supply 3 supplies a DC voltage to the power converter 2.
A current detector 4 outputs iuc, ivc, iwc which are detection values of three-phase AC currents iu, iv, iw of the magnet motor 1. The current detector 4 may also detect AC currents of two phases among the three phases of the magnet motor 1, for example, the u-phase and the w-phase, and a v-phase AC current may be obtained as iv=−(iu+iw) from an AC condition (iu+iv+iw=0). In this embodiment, the current detector 4 is provided inside the power conversion device as an example, but may also be provided outside the power conversion device.
A coordinate conversion unit 5 outputs d-axis and q-axis current detection values idc and iqc from the detection values iuc, ivC, and iwc of the three-phase AC currents iu, iv, and iw, and a phase estimation value θdc.
A speed control computation unit 6 computes a difference between a speed command value ωr* and a new speed estimation value ωdc′, and computes and outputs a first torque current command value iq0* based on the difference (deviation).
A secondary system transfer function filter computation unit 7 outputs the first torque current command value iq0*, the speed command value ωr*, and a second torque current command value iq* calculated based on a secondary system transfer function filter of a Laplace operator s.
A primary system transfer function filter computation unit 8 outputs the new speed estimation value ωdc′ based on a speed estimation value ωdc and a primary system transfer function filter of the Laplace operator s.
A vector control computation unit 9 outputs d-axis and q-axis voltage command values vdc** and vqc** computed based on the d-axis current command value id*, the q-axis current command value (second torque current command value) iq*, the current detection values idc and iqc, the new speed estimation value ωdc′, and electrical circuit parameters of the magnet motor 1.
A phase error estimation computation unit 10 uses voltage command values vdc** and vqc** of a dc-axis and a qc-axis which are the control axes, the speed estimation value ωdc, the current detection values idc and iqc, and the electrical circuit parameters of the magnet motor 1 to output an estimation value Δθc of a phase error AO, which is a deviation between the phase θdc of the control axis and a phase θd of a magnet of the magnet motor 1.
A frequency and phase estimation computation unit 11 outputs a speed estimation value ωdc and a phase estimation value θdc based on the phase error estimation value Δθc.
A coordinate conversion unit 12 outputs three-phase AC voltage command values vu*, vv*, and vw* based on the dc-axis and qc-axis voltage command values vdc** and vqc**, and the phase estimation value ac.
First, description will be given of a basic operation of a sensorless vector control method using the secondary system transfer function filter computation unit 7 and the primary system transfer function filter computation unit 8 which are characteristics of this example.
The speed control computation unit 6 computes a first q-axis current command value iq0* in accordance with (Mathematical Formula 1) by proportional control and integral control so that the new speed estimation value ωdc′ to be described later conforms to the speed command value ωr*.
Here,
Next, description will be given of the secondary system transfer function filter computation unit 7 that is one characteristic of the invention.
In the secondary system transfer function filter computation unit 7, the second torque current command value iq* is computed in accordance with (Mathematical Formula 2) by using the first q-axis current command value iq0*, a pulsating component ωn that is a mechanical angular velocity of the magnetic motor 1, and attenuation ratios ζa and ζb in a secondary system transfer function filter 71 of the Laplace operator s.
A mechanical angular velocity conversion unit 72 computes a pulsating component ωn that is a value of the number of pole pairs Pm of the magnet motor 1 and is set to the secondary system transfer function filter 71 in accordance with (Mathematical Formula 3).
A control mode selection unit 73 determines a low-speed region and a high-speed region in accordance with the magnitude of the speed command value ωr*. When it is determined as the low-speed region, “control mode for suppressing pulsation in a speed” of the magnetic motor 1 is selected, and when it is determined as the high-speed region, “control mode for suppressing pulsation in torque” is selected.
An attenuation ratio selection unit 74 sets the attenuation ratios ζa and ζb relating to the control mode for suppressing pulsation in a speed when the “control mode for suppressing pulsation in a speed” is selected, and sets the attenuation ratios ζa and ζb relating to the control mode for suppressing pulsation in a speed when the “control mode for suppressing pulsation in torque” is selected. The secondary system transfer function filter 71 includes a pulsating component ωn as a parameter, and varies in correspondence with the pulsating component.
In the primary system transfer function filter computation unit 8, the new speed estimation value ωdc′ is computed in accordance with (Mathematical Formula 4) by using the speed command value ωr*, the speed estimation value ωdc, a progress time constant Ta, and a delay time constant Tb in the primary system transfer function filter 81 of the Laplace operator s.
In the primary system transfer function filter computation unit 8, the progress time constant Ta is computed in accordance with (Mathematical Formula 5).
progress time constant Ta may be a value that is inversely proportional to the pulsating component ωn that is a mechanical angular velocity. The delay time constant Tb may be sufficiently shorter than Ta. For example, the delay time constant Tb may be a time constant (hereinafter, referred to as a response time constant of current control) corresponding to a response frequency of current control. From (Mathematical Formula 5), the progress time constant Ta in a primary system transfer function filter 81 varies in correspondence with a pulsating component of the magnet motor. The primary progress time constant Ta may be rewritten based on the speed command value.
A “control mode for suppressing pulsation in a speed” is set so that a gain becomes maximum at the pulsating component ωn, and a “control mode for suppressing pulsation in torque” is set so that a gain becomes minimum at the pulsating component ωn. The characteristics can be obtained in accordance with setting of the attenuation ratios ζa and ζb.
In the vector control computation unit 9, first, dc-axis and qc-axis voltage reference values vdc* and vqc* are output in accordance with (Mathematical Formula 6) by using a setting value R* of coil resistance that is an electrical circuit parameter of the permanent magnet motor 1, a setting value Ld* of d-axis inductance, a setting value Lq* of q-axis inductance, a setting value Ke* of an induced voltage coefficient, dc-axis and qc-axis current command values id* and iq*, and the speed estimation value ωdc.
Here, Tacr: response time constant of current control. Second, dc-axis and qc-axis voltage correction values Δvdc and Δvqc are computed in accordance with (Mathematical Formula 7) by proportional control and integral control so that the current detection values idc and iqc of respective components conform to the dc-axis and qc-axis current command values id* and iq*.
Here,
In the phase error estimation computation unit 10, the phase error estimation value Δθc is computed in accordance with an extended induced voltage formula (Mathematical Formula 9) on the basis of the dc-axis and qc-axis voltage command values vdc** and vqc**, the current detection values idc and iqc, the speed estimation value ωdc and the electrical circuit parameters (R* and Lq*) of the magnet motor 1.
The frequency and phase estimation computation unit 11 computes the speed estimation value ωdc in accordance with (Mathematical Formula 10) by P (proportional)+I (integral) control so that the above-described phase error estimation value Δθc conforms to a command value Δθc* (=0) thereof. In addition, the frequency and phase estimation computation unit 11 computes the phase estimation value θdc in accordance with (Mathematical Formula 11) by I (integral) control.
Here, Kppll: proportional gain of PLL control, Kipll: integral gain of PLL control.
Next, description will be given of the principle by which the invention provides highly stable, highly accurate, and highly efficient control characteristics.
In addition,
In the gain characteristic at the upper stage, the gain is maximized at a point E for the pulsating component ωn, and at a crossover frequency that is a frequency point F where the gain is 0 dB (gain=1), the phase characteristic at the lower stage is −180 degrees near a point G. According to this, it is clear that a phase margin of the control system is insufficient and the speed control system is unstable.
In the gain characteristic at the upper stage, the gain is maximized at a point H for the pulsating component ωn, and at a crossover frequency that is a frequency point I where the gain is 0 dB (gain=1), and in the phase characteristic at the lower stage, the phase margin near a point J is improved. That is, at 50% of the base speed, it is clear that the speed control system is stable due to an effect of the primary system transfer function filter computation unit 8.
That is, when the motor speed increases, it is necessary to improve the frequency characteristics of the speed control system, and thus the primary system transfer function filter computation unit 8 is essential. The reason for this is that the frequency characteristics of the speed control system are affected by a control band that varies depending on the response frequency FPLL set in the frequency and phase estimation computation unit 11.
When the magnitude |ωdc| of the speed estimation value ωdc and the response frequency FPLL have a relationship of (Mathematical Formula 12), the primary system transfer function filter computation unit 8 may be used.
Due to an influence of the pulsating load torque, the maximum value of the motor torque at time L in the same drawing is approximately 72%. At this time, an output P of the motor is generated by a motor speed ωr and torque τm in accordance with a relationship of (Mathematical Formula 13).
Therefore, in the high-speed region, the pulsation of the motor output P is large, and the motor efficiency is not good.
It can be seen that the effect of the invention is clear by selecting the “control mode for suppressing pulsation in a speed” when the motor is a low-speed region and the “control mode for controlling pulsation in torque” when the motor is a high-speed region.
In this example, as an example, it is assumed that a region where the magnitude of the speed command value ωr* is 10% or 50% of the base speed is the low-speed region, and a region where the magnitude is 50% or more of the base speed is the high-speed region, but a region where the magnet motor 1 steps out due to load pulsating torque may be set as the low-speed region, and a region equal to or faster than the region may be set as the high-speed region. In addition, the low-speed region and the high-speed region may be uniquely determined by the magnitude of the speed command value. In addition, the pulsating width of the motor speed ωr and the motor torque τm can be intentionally controlled by changing the attenuation ratios ζa and ζb in correspondence with the magnitude of the speed command value ωr* or the speed estimation value ωdc.
According to this example, for example, in a case where a compressor system is driven by the magnet motor, it is possible to realize a power conversion device with highly stable, highly accurate, and highly efficient control characteristics from the low-speed region to the high-speed region by selecting the “control mode for suppressing pulsation in a speed” in the low-speed region of the magnet motor and the “control mode for suppressing pulsation in torque” in the high-speed region.
In addition, according to this example, when the compressor system is driven by a motor, pulsation in torque and speed of the motor, which occur in association with load torque pulsation that varies per one mechanical angular rotation, can be freely suppressed.
Description will be given of a verification method when this example is employed with reference to
Three-phase AC voltage detection values (vuc, vvc, and vwc) which are outputs of the voltage detector 21, three-phase AC current detection values (iuc, ivc, and iwc), and a position θ that is an output of the encoder are input to a vector component voltage/current calculation unit 24 to compute vector voltage components vdc and vqc, vector current components idc and iqc, and a speed detection value ωrc obtained by differentiating the position θ.
The speed command value ωr* applied to a controller of the power converter 2 is set to, for example, approximately 50% of the base speed, and the compressor embedded with the magnet motor 1 is driven.
When the magnitude |ωrc| of the speed detection value orc and the response frequency FPLL set in the power conversion device 20 have a relationship of (Mathematical Formula 14), there is a high possibility that the primary system transfer function filter computation unit 8 is being used.
A in the
B in the
When conducting the same test by setting the magnitude of the speed command value to for every 10% of the base speed, a state of a control mode including the “control mode for suppressing pulsation in a speed” and the “control mode for suppressing pulsation in torque” becomes clear. Although the waveform of the vector current component iqc has been observed, but a waveform of the speed detection value orc may be observed.
In this example, two secondary system transfer function filters 7a1 and 7a5 are provided in the secondary system transfer function filter computation unit 7a, and the two filters are switched between the low-speed region and the high-speed region, respectively. A magnet motor 1, a speed control computation unit 6, and a primary system transfer function filter computation unit 8 to a coordinate conversion unit 12 in
In the secondary system transfer function filter computation unit 7a, a second torque current command value iq1* is computed in accordance with (Mathematical Formula 15) by using a first q-axis current command value iq0*, a pulsating component ωn, and attenuation ratios ζa1 and ζb1 of the “control mode for suppressing pulsation in a speed” in the secondary system transfer function filter 7a1 of the Laplace operator s.
In addition, a second torque current command value iq2* is computed in accordance with (Mathematical Formula 16) by using the first q-axis current command value iq0*, the pulsating component ωn, and attenuation ratios ζa2 and ζb2 of the “control mode for suppressing pulsation in torque” in the secondary system transfer function filter 7a5.
A mechanical speed conversion unit 7a2 computes the pulsating component ωn that is a value of the number of pole pairs Pm of the magnet motor 1 and is set to the secondary system transfer function filter 7a1 in accordance with (Mathematical Formula 3) described above.
A control mode selection unit 7a3 determines the low-speed region and the high-speed region in accordance with the magnitude of the speed command value ωr*. The control mode selection unit 7a3 selects the “control mode for suppressing pulsation in a speed” In a case where it is determined as the low-speed region, and the “control mode for suppressing pulsation in torque” in a case where it is determined as the high-speed region.
In a case where the control mode selection unit 7a3 selects the “control mode for suppressing pulsation in a speed”, a filter switching unit 7a4 outputs iq1* that is an output of the secondary system transfer function filter 7a1 as a q-axis current command values iq*. In addition, in a case where the “control mode for suppressing pulsation in torque” is selected, the filter switching unit 7a4 outputs iq2* that is an output of the secondary system transfer function filter 7a5 as a q-axis current command values iq*. Even using this example in which two secondary system transfer function filters are prepared and are switched between the low-speed region and the high-speed region, stable, highly accurate, and highly efficient control characteristics can be realized as in Example 1.
In this example, a low-power region is switched to the “control mode for suppressing pulsation in a speed” and a high-power region is switched to the “control mode for suppressing pulsation in torque” by using a power value. The magnet motor 1 to the speed control computation unit 6, and the phase error estimation computation unit 10 to the coordinate conversion unit 12 in
In a vector control computation unit 9a, in addition to the computations of (Mathematical Formula 6) to (Mathematical Formula 8), effective power Pa, which is an inner product of the power of the magnet motor 1, is computed in accordance with (mathematical Formula 17) by using voltage command values vdc and vqc**, and current detection values idc and iqc of the dc-axis and qc-axis.
A control mode selection unit 7b3 determines a low-power region and a high-power region in accordance with the magnitude of the effective power Pa. The control mode selection unit 7b3 selects the “control mode for suppressing pulsation in a speed” when it is determined as the low-power region, and the “control mode for suppressing pulsation in torque” when it is determined as the high-power region. Even when using this example in which switching is performed between the low-power region and the high-power region, highly stable, highly accurate, and highly efficient control characteristics can be realized as in Example 1.
Here, in the secondary system transfer function filter 7b1 of the secondary system transfer function filter computation unit 7b, the attenuation ratios ζa and ζb are switched between the low-power region and the high-power region, but two secondary system filters may be provided as in Example 2 for the low-power region and the high-power region, respectively.
In this example, state quantities of the control (for example, the voltage command value, the current detection value, the phase error, and the first speed estimation value or the second speed estimation value) are fed back to an IOT controller 31 as a higher-level device. Then, the attenuation ratios ζa and ζb, the progress time constant Ta, and the switching speed ωchg learned by the IOT controller 31 through machine learning are reset to the secondary system transfer function filter computation unit 7 and the primary system transfer function filter computation unit 8 in this example.
A magnet motor 1 to a coordinate conversion unit 12 in
Even when using this example, highly stable, highly accurate, and highly efficient control characteristics can be realized as in Example 1 without adjustment.
The magnet motor 1 that is a constituent element in
The control unit 20a is constituted by a semiconductor integrated circuit (arithmetic control unit) such as a microcomputer and a digital signal processor (DSP). Some or all of functions of the control unit 20a can be constituted by hardware such as an application specific integrated circuit (ASIC) and a field programmable gate array (FPGA). A central processing unit (CPU) of the control unit 20a reads out a program stored in a recording device such as a memory, and executes processing of each unit such as the coordinate conversion unit 5.
The power converter 2 to the current detector 4 in
In a case where this example is applied to a compressor system driven by a magnet motor, it is possible to realize highly stable, highly accurate, and highly efficient control characteristics even in position sensorless vector control. In addition, “Ψchg that is the low-speed/high-speed region switching speed 26”, and “ζa and ζb which are attenuation ratios of the low-speed/high-speed region, and time constant Ta” 27 may be set on a field bus such as a programmable logic controller, a local area network connected to a computer, and an IOT controller.
Furthermore, although this example has been disclosed by using Example 1, Example 2 to Example 4 may also be used.
According to this example, parameters to be set in the primary system transfer function filter computation unit and the secondary system transfer function filter computation unit can be set and changed from an outer side.
In Example 1 to Example 5, voltage correction values Δvdc and Δvqc are created from the current command values id* and iq* and the current detection values idc, iqc, and the voltage correction values are added to the voltage reference value for vector control as computation shown in (Mathematical Formula 8). Without limitation thereto, intermediate current command values id** and iq** shown in (Mathematical Formula 18) used for vector control computation may be created from the current command values id* and iq* and the current detection values idc, iqc, and vector control computation shown in (Mathematical Formula 19) may be performed by using the speed estimation value ωdc and the electric circuit parameters of the magnet motor 1.
Here,
Alternatively, from the current command values id* and iq* and the current detection values idc and iqc, a voltage correction value Δvd_p* of a dc-axis proportional computation component, a voltage correction value Δvd_i*of do-axis integral computation component, a voltage correction value Δvq_p* of qc-axis proportional computation component, and a voltage correction value Δvq_i* of qc-axis integral calculation component, which are used for vector control computation, are created in accordance with (Mathematical Formula 20). Then, vector control computation shown in (Mathematical Formula 21) may be performed by using the speed estimation value ωdc and the electrical circuit parameters of the magnet motor 1.
Here,
In addition, vector control computation shown in (Mathematical Formula 22) may be performed by using a primary delay signal iqctd of the dc-axis current command value id* and a qc-axis current detection value iqc, the speed estimation value ωdc, and the electrical circuit parameters of the magnet motor 1.
Note that, in Example 1 to Example 5, switching elements constituting the power converter 2 may be silicon (Si) semiconductor elements or wide band gap semiconductor elements such as silicon carbide (SiC) and gallium nitride (GaN).
| Number | Date | Country | Kind |
|---|---|---|---|
| 2022-185834 | Nov 2022 | JP | national |
| Filing Document | Filing Date | Country | Kind |
|---|---|---|---|
| PCT/JP2023/023716 | 6/27/2023 | WO |