CN111931343A - Stability analysis method and system for permanent magnet synchronous motor driving system - Google Patents

Stability analysis method and system for permanent magnet synchronous motor driving system Download PDF

Info

Publication number
CN111931343A
CN111931343A CN202010653858.3A CN202010653858A CN111931343A CN 111931343 A CN111931343 A CN 111931343A CN 202010653858 A CN202010653858 A CN 202010653858A CN 111931343 A CN111931343 A CN 111931343A
Authority
CN
China
Prior art keywords
permanent magnet
magnet synchronous
synchronous motor
driving system
motor driving
Prior art date
Legal status (The legal status 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 status listed.)
Granted
Application number
CN202010653858.3A
Other languages
Chinese (zh)
Other versions
CN111931343B (en
Inventor
林程
孙建侠
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beijing Institute of Technology BIT
Original Assignee
Beijing Institute of Technology BIT
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 Beijing Institute of Technology BIT filed Critical Beijing Institute of Technology BIT
Publication of CN111931343A publication Critical patent/CN111931343A/en
Application granted granted Critical
Publication of CN111931343B publication Critical patent/CN111931343B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/02Reliability analysis or reliability optimisation; Failure analysis, e.g. worst case scenario performance, failure mode and effects analysis [FMEA]

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Control Of Ac Motors In General (AREA)

Abstract

The invention provides a stability analysis method and a system for a permanent magnet synchronous motor driving system, wherein the method comprises the following steps of establishing a permanent magnet synchronous motor driving system model, and further comprises the following steps: and according to the permanent magnet synchronous motor driving system model, carrying out stability analysis on the permanent magnet synchronous motor driving system by adopting a Nyquist stability criterion to obtain a stable condition of the permanent magnet synchronous motor driving system. The stability condition provided by the invention limits parameters such as a filter, a motor, a controller and the like, and provides a basis for the design of a permanent magnet synchronous motor driving system.

Description

Stability analysis method and system for permanent magnet synchronous motor driving system
Technical Field
The invention relates to the field of stability analysis, in particular to a stability analysis method and system for a permanent magnet synchronous motor driving system.
Background
The permanent magnet synchronous motor is widely applied to the fields of electric automobiles, fully electrically propelled ships, multi-electric airplanes and the like due to the advantages of high efficiency, high response, high power density and the like. The driving system of the permanent magnet synchronous motor is generally composed of an inverter, a controller and the permanent magnet synchronous motor, and is a complex nonlinear system. The permanent magnet synchronous motor generally adopts strict double closed loop control, so that the load power is constant, and the current rises when the voltage drops, thereby presenting a negative impedance characteristic. Research shows that negative impedance characteristics cause system stability to be reduced, even oscillation or collapse of bus voltage in severe cases, so that misoperation of protection equipment or damage of electric equipment is caused, and negative influence is brought to stable operation of the system. In response to these problems, many scholars have conducted extensive research, which can be mainly classified into large signal stability analysis and small signal stability analysis. And the large signal stability analysis adopts an equivalent negative resistance or a current source for analysis, and the system stability is analyzed according to a Lypunov stability criterion. The method is suitable for system stability analysis during low-frequency large-signal disturbance, and analysis of high-frequency small-signal disturbance is lacked. The small signal stability analysis is suitable for the system stability analysis during high-frequency small signal disturbance, but the analysis process is complex and the calculation difficulty is high.
Disclosure of Invention
The invention provides a stability analysis method and a system for a permanent magnet synchronous motor driving system, wherein the method establishes a system integral mathematical model, defines the system stability condition of the permanent magnet synchronous motor driving system, can simply and quickly judge the system stability by judging whether the permanent magnet synchronous motor meets the stability condition, and provides a theoretical basis for the design of the permanent magnet synchronous motor driving system.
A first aspect of the present invention provides a stability analysis method for a permanent magnet synchronous motor driving system, the method including the steps of:
establishing a permanent magnet synchronous motor driving system model;
and according to the permanent magnet synchronous motor driving system model, carrying out stability analysis on the permanent magnet synchronous motor driving system by adopting a Nyquist stability criterion to obtain a stable condition of the permanent magnet synchronous motor driving system.
Preferably, the establishing of the permanent magnet synchronous motor driving system model comprises the following steps:
establishing a mathematical model and a small signal model of the permanent magnet synchronous motor;
establishing a mathematical model and a small signal model of the inverter;
establishing a mathematical model and a small signal model of the controller;
and establishing an admittance model of the permanent magnet synchronous motor driving system.
In any of the above schemes, preferably, the establishing a mathematical model and a small-signal model of the permanent magnet synchronous motor includes:
(11) determining a stator current equation of the permanent magnet synchronous motor under a dq coordinate system as follows:
Figure BDA0002575992360000021
wherein id、iqDq-axis stator currents, respectively; u. ofd、uqDq-axis stator voltages, respectively; rsIs stator resistance, Ld、LqRespectively, dq-axis inductance,. psifIs a permanent magnet flux linkage, omegaeThe electrical angular velocity of the motor;
(12) performing Laplace transformation on a stator current equation of the permanent magnet synchronous motor in a dq coordinate system to obtain id、iqAnd ud、uqTransfer function between:
Figure BDA0002575992360000022
(13) the dq-axis voltage current is represented in the form of X + X, where X is the steady-state operating point and X is the small-signal disturbance. According to a stator current equation of the permanent magnet synchronous motor, the steady-state working point of the stator voltage is as follows:
Figure BDA0002575992360000023
the corresponding small signal model is:
Figure BDA0002575992360000024
in any of the above aspects, preferably, the establishing a mathematical model and a small-signal model of the inverter includes:
(21) defining a switching function Sj,j=a,b,c;S j1 denotes upper arm closed, S j0 denotes that the lower arm is closed, thus obtaining:
Figure BDA0002575992360000025
wherein u isan、ubn、ucnThree-phase voltage of points a, b and c relative to the negative pole N of the power supply, udcIs a dc supply voltage.
(22) For symmetrical three-phase loads, there is uon=(uan+ubn+ucn) /3 wherein uonIs the voltage between the load neutral point O and the negative pole N of the power supply, and therefore has
Figure BDA0002575992360000031
Wherein u isao、ubo、ucoThe point a, b and c are phase voltages relative to the neutral point O of the load.
(23) Since the inverter switching frequency is much higher than the motor frequency, the switching function S is adjustedjApproximated by duty cycle dj. Performing Park conversion on the formula obtained in the step (22) to obtain
Figure BDA0002575992360000032
Wherein d isd、dqIs the duty cycle in dq coordinate system.
(24) Will dd、dqAnd udcExpressed in the form of X + X, the steady state operating point is, according to the formula obtained in step (23):
Figure BDA0002575992360000033
the corresponding small signal model is
Figure BDA0002575992360000034
Wherein D isd、Dq、UdcAre respectively dd、dqAnd udcSteady state operating point of dd、dqAnd udcCorresponding small signal perturbations.
In any of the above aspects, preferably, the establishing a mathematical model and a small-signal model of the controller includes:
(31) the transfer function of the current loop PI controller is Gc=Kp+KiS, wherein KpAnd KiRespectively a scaling factor and an integration factor.
(32) And after the current given and the current feedback are compared, the output signal passes through a PI controller and is used as dq axis voltage given after decoupling compensation. Thus, there are
Figure BDA0002575992360000035
The corresponding small signal model is
Figure BDA0002575992360000036
(33) The SVPWM modulation scheme can be regarded as an effect of superimposing a zero-sequence component on a corresponding SPWM modulation signal. According to the triangle similarity theorem, for j-phase, j is a, b, c, the duty cycle
Figure BDA0002575992360000037
Wherein T issIs one sampling period, TjFor the on-time of the j-phase leg in one sampling period,
Figure BDA0002575992360000038
given phase voltages of points a, b and c relative to the negative pole N of the power supply and given phase voltages of points a, b and c relative to the neutral point O of the load
Figure BDA0002575992360000039
Has a zero sequence component eoI.e. by
Figure BDA00025759923600000310
Thus, there are
Figure BDA0002575992360000041
(34) Performing Park transformation on the formula in (33) can obtain:
Figure BDA0002575992360000042
the corresponding small signal model is
Figure BDA0002575992360000043
In any of the above aspects, preferably, the establishing an admittance model of the permanent magnet synchronous motor drive system includes:
(41) for a PMSM drive system, there are times when power losses are neglected
Figure BDA0002575992360000044
Expressing the variables in the formula as X + X, and simplifying to obtain the system admittance Y:
Figure BDA0002575992360000045
(42) from the small-signal models established in (13), (24), (32), and (34), it can be found that:
Figure BDA0002575992360000046
(43) by substituting the formula in (42) into the formula for admittance Y in (41), Y can be obtainedd+YqWherein
Figure BDA0002575992360000047
(44) Let d-axis current be given as 0, U is known from the steady state modeld=-ωeLqIq,UdcIdc=3UqIq/2, mixing U withd=-ωeLqIqAnd UdcIdc=3UqIqThe formula of substituting/2 into (43) can be obtained:
Figure BDA0002575992360000051
in any of the above schemes, preferably, the performing, according to the permanent magnet synchronous motor driving system model, the stability analysis on the permanent magnet synchronous motor driving system by using the Nyquist stability criterion includes:
establishing a transfer function of a permanent magnet synchronous motor driving system;
obtaining a permanent magnet synchronous motor driving system stability condition by adopting a Nyquist stability criterion and combining a permanent magnet synchronous motor driving system transfer function;
and analyzing the stability of the permanent magnet synchronous motor driving system according to the stable condition.
In any of the above aspects, preferably, establishing the transfer function of the driving system of the permanent magnet synchronous motor comprises:
according to an impedance analysis method, equivalence is carried out on an admittance Y for a permanent magnet synchronous motor driving system, and impedance Z for a direct current end filterinIn an equivalent manner, the first and second electrodes,
Figure BDA0002575992360000052
Rfis a line resistance, LfIs a filter inductor and C is a support capacitor.
According to circuit principles, the transfer function of the established system is:
Figure BDA0002575992360000053
wherein u isinIs the equivalent supply voltage.
In any of the above schemes, preferably, the method for obtaining the stability condition of the driving system of the permanent magnet synchronous motor by using a Nyquist stability criterion in combination with the transfer function of the driving system of the permanent magnet synchronous motor comprises:
Figure BDA0002575992360000054
in any of the above aspects, preferably, analyzing the stability of the permanent magnet synchronous motor driving system according to the stable condition includes: if the permanent magnet synchronous motor driving system meets the stability condition, the system is stable, or the value range of one or more components of the permanent magnet synchronous motor driving system is judged according to the stability condition, so that the system is stable.
A second aspect of the present invention provides a stability analysis system for a driving system of a permanent magnet synchronous motor, which is configured to operate the stability analysis method for a driving system of a permanent magnet synchronous motor to analyze the stability of the driving system of a permanent magnet synchronous motor.
The stability analysis method and the system for the permanent magnet synchronous motor driving system establish an overall mathematical model of the permanent magnet synchronous motor driving system, and perform stability analysis on the system based on an impedance analysis method. The overall mathematical model of the permanent magnet synchronous motor driving system comprises an inverter model, a motor model and a controller model; establishing a small signal model of the system according to a small signal analysis method, and establishing an admittance model of a permanent magnet synchronous motor driving system; establishing a transfer function model of the system according to an impedance analysis method, and limiting an open-loop transfer function phase angle according to a Nyquist stability criterion to meet the stability requirement; and finally, deducing the stability condition of the system in the whole frequency domain according to the mentioned phase angle limiting condition, thereby meeting the global stability requirement. The stability conditions limit parameters such as a filter, a motor, a controller and the like, and provide a basis for the design of a permanent magnet synchronous motor driving system.
Drawings
Fig. 1 is a schematic flow diagram of a stability analysis method for a permanent magnet synchronous motor drive system according to the present invention.
Fig. 2 is a block diagram of a driving system of a permanent magnet synchronous motor.
Fig. 3 is a three-phase inverter topology.
Fig. 4 is a schematic diagram of the SVPWM modulation process.
Fig. 5 is an overall mathematical model of a permanent magnet synchronous motor driving system for a stability analysis method of the permanent magnet synchronous motor driving system according to the present invention.
Fig. 6 is a small signal model of a permanent magnet synchronous motor driving system for a stability analysis method of the permanent magnet synchronous motor driving system according to the present invention.
Fig. 7 is an equivalent circuit diagram of a permanent magnet synchronous motor driving system for a stability analysis method of the permanent magnet synchronous motor driving system according to the present invention.
Fig. 8 is an admittance phase-frequency diagram of a permanent magnet synchronous motor driving system for a stability analysis method of a permanent magnet synchronous motor driving system according to the present invention.
Detailed Description
For a better understanding of the present invention, reference will now be made in detail to the following examples.
Example 1
As shown in fig. 1, a stability analysis method for a permanent magnet synchronous motor driving system, the method comprising the steps of:
s1, establishing a permanent magnet synchronous motor driving system model;
and S2, according to the permanent magnet synchronous motor driving system model, adopting a Nyquist stability criterion to perform stability analysis on the permanent magnet synchronous motor driving system to obtain the stable condition of the permanent magnet synchronous motor driving system.
Fig. 2 is a block diagram of a driving system of a permanent magnet synchronous motor, in which a dc bus supplies power to an inverter after passing through a filter, and the inverter converts dc power into three-phase ac power and supplies power to the motor. The controller adopts a double closed loop PI control structure, a rotating speed loop provides a given signal for q-axis current, and d-axis current is given and usually set to be 0 or determined according to the field weakening degree of the motor. In order to improve the control effect of the dq-axis current, a voltage decoupling unit is usually added in a current loop, and a signal output by the current loop is subjected to decoupling compensation and then serves as a dq-axis voltage given signal. And after the given voltage signal is converted to an alpha beta coordinate system, the voltage type inverter is controlled by adopting an SVPWM (space vector pulse width modulation) technology to supply power to the permanent magnet synchronous motor, so that the motor drive control is realized. It can be seen that the pmsm driving system is a complex nonlinear system, including the motor, the inverter, the controller, etc., therefore, in step S1, when building the pmsm driving system model, the motor, the inverter, and the controller need to be modeled separately.
In step S1, the establishing of the driving system model of the permanent magnet synchronous motor includes the steps of:
s11, establishing a mathematical model and a small signal model of the permanent magnet synchronous motor;
s12, establishing a mathematical model and a small signal model of the inverter;
s13, establishing a mathematical model and a small signal model of the controller;
and S14, establishing an admittance model of the permanent magnet synchronous motor driving system.
In step S11, the permanent magnet synchronous motor model adopts a typical mathematical model in dq coordinate system, and establishes a corresponding motor small signal model according to a small signal analysis method, including:
(11) determining a stator current equation of the permanent magnet synchronous motor under a dq coordinate system as follows:
Figure BDA0002575992360000071
wherein id、iqDq-axis stator currents, respectively; u. ofd、uqDq-axis stator voltages, respectively; rsIs stator resistance, Ld、LqRespectively, dq-axis inductance,. psifIs a permanent magnet flux linkage, omegaeThe electrical angular velocity of the motor;
(12) performing Laplace transformation on a stator current equation of the permanent magnet synchronous motor in a dq coordinate system to obtain id、iqAnd ud、uqTransfer function between:
Figure BDA0002575992360000081
(13) will dThe q-axis voltage current is represented in the form of X + X, where X is the steady state operating point and X is the small signal disturbance. According to a stator current equation of the permanent magnet synchronous motor, the steady-state working point of the stator voltage is as follows:
Figure BDA0002575992360000082
the corresponding small signal model is:
Figure BDA0002575992360000083
in step S12, a mathematical model and a small-signal model of the inverter are established. The inverter circuit topology is as shown in fig. 3, the inverter adopts a mean value model, a mathematical model under a dq coordinate system is established through Park transformation, and a small signal model of the inverter is established according to a small signal analysis method. The establishment of the mathematical model and the small signal model of the inverter comprises the following steps:
(21) defining a switching function Sj,j=a,b,c;S j1 represents the upper arm closed, Sj 0 represents the lower arm closed, thus obtaining:
Figure BDA0002575992360000084
wherein u isan、ubn、ucnThree-phase voltage of points a, b and c relative to the negative pole N of the power supply, udcIs a dc supply voltage.
(22) For symmetrical three-phase loads, there is uon=(uan+ubn+ucn) /3 wherein uonIs the voltage between the load neutral point O and the negative pole N of the power supply, and therefore has
Figure BDA0002575992360000085
Wherein u isao、ubo、ucoThe point a, b and c are phase voltages relative to the neutral point O of the load.
(23) Since the inverter switching frequency is much higher than the motor frequency, the switching function S is adjustedjApproximated by duty cycle dj. Formulating the result of step (22)The line Park is transformed into
Figure BDA0002575992360000086
Wherein d isd、dqIs the duty cycle in dq coordinate system.
(24) Will dd、dqAnd udcExpressed in the form of X + X, the steady state operating point is, according to the formula obtained in step (23):
Figure BDA0002575992360000091
the corresponding small signal model is
Figure BDA0002575992360000092
Wherein D isd、Dq、UdcAre respectively dd、dqAnd udcSteady state operating point of dd、dqAnd udcCorresponding small signal perturbations.
In step S13, a mathematical model and a small-signal model of the controller are established. The controller comprises two links of PI control and SVPWM modulation. The current loop PI controller inputs signals for current giving and current feedback, the output signals are used as dq axis voltage giving after decoupling compensation, and the specific control process is shown in figure 2. A mathematical model of a PI control link and a corresponding small signal model can be established according to a control block diagram. The SVPWM modulation strategy modulates a given voltage signal output by the PI control link to complete waveform output, and the SVPWM modulation process is as shown in fig. 2. The SVPWM modulation scheme can be regarded as the effect of superimposing zero sequence component on the corresponding SPWM modulation signal, and the points a, b, and c can be regarded as the given phase voltage relative to the negative pole N of the power supply
Figure BDA0002575992360000093
Given phase voltages with points a, b, c relative to the load neutral point O
Figure BDA0002575992360000094
Has a zero sequence component eoTherefore, a mathematical relation between the duty ratio and the given phase voltage can be established, and the data can be obtained under the dq coordinate system after Park changeA mathematical model of the controller. Similarly, a corresponding small-signal model can be established according to the small-signal method. The establishment of the mathematical model and the small signal model of the controller comprises the following steps:
(31) the transfer function of the current loop PI controller is Gc=Kp+KiS, wherein KpAnd KiRespectively a scaling factor and an integration factor.
(32) And after the current given and the current feedback are compared, the output signal passes through a PI controller and is used as dq axis voltage given after decoupling compensation. Thus, there are
Figure BDA0002575992360000095
The corresponding small signal model is
Figure BDA0002575992360000096
(33) The SVPWM modulation scheme can be regarded as an effect of superimposing a zero-sequence component on a corresponding SPWM modulation signal. According to the triangle similarity theorem, for j-phase, j is a, b, c, the duty cycle
Figure BDA0002575992360000097
Wherein T issIs one sampling period, TjFor the on-time of the j-phase leg in one sampling period,
Figure BDA0002575992360000098
given phase voltages of points a, b and c relative to the negative pole N of the power supply and given phase voltages of points a, b and c relative to the neutral point O of the load
Figure BDA0002575992360000099
Has a zero sequence component eoI.e. by
Figure BDA00025759923600000910
Thus, there are
Figure BDA00025759923600000911
(34) Performing Park transformation on the formula in (33) can obtain:
Figure BDA0002575992360000101
the corresponding small signal model is
Figure BDA0002575992360000102
In step S14, an admittance model of the driving system of the permanent magnet synchronous motor is established. The overall mathematical model of the permanent magnet synchronous motor driving system is shown in fig. 5, and comprises a dq axis model of the motor, the inverter and the controller, and the corresponding system small signal model is shown in fig. 6. Because the time constant of the rotating speed loop is far larger than that of the current loop, the rotating speed loop can be ignored when establishing a small signal model of the controller, and only the influence of the current loop is considered. After the system model is built, respectively solving according to the small signal model of the system
Figure BDA0002575992360000103
Substituting into a calculation formula of the system admittance to obtain an expression of the system admittance Y. In particular, by giving the d-axis current 0, a further simplified admittance expression can be obtained. The establishment of the admittance model of the permanent magnet synchronous motor driving system comprises the following steps:
(41) for a PMSM drive system, there are times when power losses are neglected
Figure BDA0002575992360000104
Expressing the variables in the formula as X + X, and simplifying to obtain the system admittance Y:
Figure BDA0002575992360000105
(42) from the small-signal models established in (13), (24), (32), and (34), it can be found that:
Figure BDA0002575992360000106
(43) by substituting the formula in (42) into the formula for admittance Y in (41), Y can be obtainedd+YqWherein
Figure BDA0002575992360000111
(44) Let d-axis current be given as 0, U is known from the steady state modeld=-ωeLqIq,UdcIdc=3UqIq/2, mixing U withd=-ωeLqIqAnd UdcIdc=3UqIqThe formula of substituting/2 into (43) can be obtained:
Figure BDA0002575992360000112
it should be understood that the above step numbers for establishing the driving system model of the permanent magnet synchronous motor are only used to distinguish the steps, and not to limit the sequence of the steps, and those skilled in the art may make adaptive adjustments to the sequence of the steps according to actual situations, and these adjustments do not make the essence of the corresponding technical solution depart from the scope of the technical solution of the present invention, for example, the sequence of the steps S11-S13 may be arbitrarily adjusted.
The expression of the system admittance Y shows that Y is formed by YdAnd YqTwo parts, wherein Y isdIs positive admittance, and YqThen has a negative admittance characteristic, and in the low frequency band, the amplitude of Y is-Idc/UdcAnd the negative impedance characteristic can reduce the stability area of the system and even cause the instability of the system in severe cases, so that the stability analysis is carried out on the permanent magnet synchronous motor driving system by adopting a Nyquist stability criterion.
In step S2, according to the permanent magnet synchronous motor driving system model, performing stability analysis on the permanent magnet synchronous motor driving system by using a Nyquist stability criterion includes:
s21, establishing a transfer function of a permanent magnet synchronous motor driving system;
s22, obtaining a permanent magnet synchronous motor driving system stability condition by adopting a Nyquist stability criterion and combining a permanent magnet synchronous motor driving system transfer function;
and S23, analyzing the stability of the permanent magnet synchronous motor driving system according to the stable condition.
In step S21, a transfer function of the driving system of the permanent magnet synchronous motor is established.
According to the impedance analysis method, the driving system of the permanent magnet synchronous motor is equivalent to a circuit shown in fig. 7, the driving system of the permanent magnet synchronous motor is equivalent by admittance Y, thevenin is equivalent to the power supply side, and a direct current end filter is equivalent to impedance Zin
Figure BDA0002575992360000113
RfIs a line resistance, LfIs a filter inductor, C is a supporting capacitor, and the power supply voltage usIs correspondingly equivalent to uin. According to the circuit principle, the closed loop transfer function of the system is established as follows:
Figure BDA0002575992360000121
open loop transfer function of ZinY。
In step S22, obtaining the permanent magnet synchronous motor drive system stability condition by using the Nyquist stability criterion in combination with the transfer function of the permanent magnet synchronous motor drive system includes:
s221, according to a Nyquist stability criterion, ZinY may not encompass the (-1, j0) point, so ZinWhen the phase angle of Y is equal to +/-180 DEG, the amplitude value should be less than 1, namely ZinThe intersection of Y with the real axis is to the right of point (-1, j 0).
In the low frequency band having ZinYω→0=-IdcRf/UdcTo make ZinY amplitude is less than 1, condition A, R should be satisfiedf<Udc/Idc
In the high frequency band having ZinYω→∞0, located to the right of point (-1, j 0).
In the middle frequency bandinThe phase angle of Y is limited between (-180 degrees, 180 degrees), so that Z can be ensuredinY does not encompass the (-1, j0) point throughout the band.
Open loop transfer function Z for mid-rangeinPhase angle analysis was performed for Y: admittance Y of permanent magnet synchronous motor driving system is YdAnd YqTwo parts are composed of YqAnd YdIs limited to (0 DEG, 180 DEG), the phase angle of admittance Y must be YqAnd YdThe phase angle of (a) is shown in fig. 8 as the phase-frequency diagram of the admittance. According to YqAnd YdIs expressed byq/YdIs represented by-A1(s)A2(s) wherein
Figure BDA0002575992360000122
A is to be1(s)、A2The phase angles of(s) are all limited to the range of (-90 °,0 °).
If restriction A1The phase angle range of(s) is (-90 DEG, 0 DEG), and the condition B, K is satisfiedp>Uq/Iq
Substituting s-j ω into a2The expression for(s) can be found:
Figure BDA0002575992360000123
due to Ld<Lq,A2The imaginary part of(s) is determined to be negative, such that A2The phase angle of(s) is in the range of (-90 DEG, 0 DEG), provided that A is2The real part of(s) is positive, and therefore the conditions C and (R) should be satisfieds+Kp)2>Ki(Ld+Lq)。
The phase angle of admittance Y being at YqAnd YdAnd the upper limit of the angle is angle YqLower limit of less than YdThus, making ZinYdIs greater than-180 DEG, ZinYqThe phase angle of (A) is less than 180 DEG, that is, Z can be ensuredinThe phase angle of Y is between (-180 deg., 180 deg.). ZinHas a phase angle of (-90 DEG, 90 DEG), YdThe phase angle of (1) is in the range of (-90 °,90 °). Thus, ZinYdThe phase angle of (a) is greater than-180 °; y with s ═ j ω substitutedqAnd ZinExpression, expressed in complex form, of ZinYqIs positive, then Z may be madeinYqIs always less than 180 deg.. Therefore, the conditions D,
Figure BDA0002575992360000131
In summary, satisfying both conditions B, C and D in the IF range ensures an open-loop transfer function ZinThe phase angle of Y is in the range of (-180 DEG, 180 DEG), and there is no intersection with the real axis, and at the same time, because Z isinY is located right from the point of (-1, j0) at the intersection of the real axis at both the low and high frequency bands, and it is necessary to satisfy the condition A, and therefore, when the conditions A, B, C and D are satisfied in the entire frequency domain, Z isinThe Y curve does not surround the point (-1, j0), and the system is stable, namely the stable condition of the permanent magnet synchronous motor driving system is as follows:
Figure BDA0002575992360000132
in step S23, analyzing the stability of the driving system of the permanent magnet synchronous motor according to the stable condition includes: if the permanent magnet synchronous motor driving system meets the stability condition, the system is stable or the value range of one or more components of the permanent magnet synchronous motor driving system is judged according to the stability condition so as to keep the system stable.
Example 2
A stability analysis system for a permanent magnet synchronous motor driving system is used for operating the method for analyzing the stability of the permanent magnet synchronous motor driving system, and the stable condition of the permanent magnet synchronous motor driving system is determined. It should be understood that all or part of the steps of implementing the method for analyzing the stability of the driving system of the permanent magnet synchronous motor can be indicated by a computer program which can be stored in a non-volatile computer readable storage medium and can implement all or part of the steps of the method. If the stable condition of the permanent magnet synchronous motor driving system is stored in the nonvolatile computer readable storage medium, after the structural parameters of a certain permanent magnet synchronous motor driving system are input, the system automatically establishes a mathematical model for the permanent magnet synchronous motor driving system, and judges the value range of one or more components according to the comprehensive stable condition so as to keep the permanent magnet synchronous motor driving system stable.
It should be noted that the above embodiments are only used for illustrating the technical solution of the present invention, and not for limiting the same; although the foregoing embodiments illustrate the invention in detail, those skilled in the art will appreciate that: it is possible to modify the technical solutions described in the foregoing embodiments or to substitute some or all of the technical features thereof, without departing from the scope of the technical solutions of the present invention.

Claims (10)

1. A stability analysis method for a permanent magnet synchronous motor drive system includes the steps of: the method for establishing the permanent magnet synchronous motor driving system model is characterized by further comprising the following steps: and according to the permanent magnet synchronous motor driving system model, carrying out stability analysis on the permanent magnet synchronous motor driving system by adopting a Nyquist stability criterion to obtain a stable condition of the permanent magnet synchronous motor driving system.
2. The stability analysis method for a permanent magnet synchronous motor drive system according to claim 1, wherein establishing a permanent magnet synchronous motor drive system model comprises the steps of:
establishing a mathematical model and a small signal model of the permanent magnet synchronous motor;
establishing a mathematical model and a small signal model of the inverter;
establishing a mathematical model and a small signal model of the controller;
and establishing an admittance model of the permanent magnet synchronous motor driving system.
3. The stability analysis method for a permanent magnet synchronous motor driving system according to claim 2, characterized in that: the establishment of the mathematical model and the small signal model of the permanent magnet synchronous motor comprises the following steps:
(11) determining a stator current equation of the permanent magnet synchronous motor under a dq coordinate system as follows:
Figure FDA0002575992350000011
wherein id、iqDq-axis stator currents, respectively; u. ofd、uqDq-axis stator voltages, respectively; rsIs stator resistance, Ld、LqRespectively, dq-axis inductance,. psifIs a permanent magnet flux linkage, omegaeThe electrical angular velocity of the motor;
(12) performing Laplace transformation on a stator current equation of the permanent magnet synchronous motor in a dq coordinate system to obtain id、iqAnd ud、uqTransfer function between:
Figure FDA0002575992350000012
(13) the dq-axis voltage current is represented in the form of X + X, where X is the steady-state operating point and X is the small-signal disturbance. According to a stator current equation of the permanent magnet synchronous motor, the steady-state working point of the stator voltage is as follows:
Figure FDA0002575992350000013
the corresponding small signal model is:
Figure FDA0002575992350000014
4. the stability analysis method for a permanent magnet synchronous motor driving system according to claim 3, characterized in that: the establishment of the mathematical model and the small signal model of the inverter comprises the following steps:
(21) defining a switching function Sj,j=a,b,c;Sj1 representsUpper arm closed, Sj0 denotes that the lower arm is closed, thus obtaining:
Figure FDA0002575992350000021
wherein u isan、ubn、ucnThree-phase voltage of points a, b and c relative to the negative pole N of the power supply, udcIs a dc supply voltage.
(22) For symmetrical three-phase loads, there is uon=(uan+ubn+ucn) /3 wherein uonIs the voltage between the load neutral point O and the negative pole N of the power supply, and therefore has
Figure FDA0002575992350000022
Wherein u isao、ubo、ucoThe point a, b and c are phase voltages relative to the neutral point O of the load.
(23) Since the inverter switching frequency is much higher than the motor frequency, the switching function S is adjustedjApproximated by duty cycle dj. Performing Park conversion on the formula obtained in the step (22) to obtain
Figure FDA0002575992350000023
Wherein d isd、dqIs the duty cycle in dq coordinate system.
(24) Will dd、dqAnd udcExpressed in the form of X + X, the steady state operating point is, according to the formula obtained in step (23):
Figure FDA0002575992350000024
the corresponding small signal model is
Figure FDA0002575992350000025
Wherein D isd、Dq、UdcAre respectively dd、dqAnd udcSteady state operating point of dd、dqAnd udcCorresponding small signal perturbations.
5. The stability analysis method for a permanent magnet synchronous motor drive system according to claim 4, characterized in that: the establishment of the mathematical model and the small signal model of the controller comprises the following steps:
(31) the transfer function of the current loop PI controller is Gc=Kp+KiS, wherein KpAnd KiRespectively a scaling factor and an integration factor.
(32) And after the current given and the current feedback are compared, the output signal passes through a PI controller and is used as dq axis voltage given after decoupling compensation. Thus, there are
Figure FDA0002575992350000026
The corresponding small signal model is
Figure FDA0002575992350000027
(33) The SVPWM modulation scheme can be regarded as an effect of superimposing a zero-sequence component on a corresponding SPWM modulation signal. According to the triangle similarity theorem, for j-phase, j is a, b, c, the duty cycle
Figure FDA0002575992350000031
Wherein T issIs one sampling period, TjFor the on-time of the j-phase leg in one sampling period,
Figure FDA0002575992350000032
given phase voltages of points a, b and c relative to the negative pole N of the power supply and given phase voltages of points a, b and c relative to the neutral point O of the load
Figure FDA0002575992350000033
Has a zero sequence component eoI.e. by
Figure FDA0002575992350000034
Thus, there are
Figure FDA0002575992350000035
(34) Performing Park transformation on the formula in (33) can obtain:
Figure FDA0002575992350000036
the corresponding small signal model is
Figure FDA0002575992350000037
6. The stability analysis method for a permanent magnet synchronous motor driving system according to claim 5, characterized in that: the establishment of the admittance model of the permanent magnet synchronous motor driving system comprises the following steps:
(41) for a PMSM drive system, there are times when power losses are neglected
Figure FDA0002575992350000038
Expressing the variables in the formula as X + X, and simplifying to obtain the system admittance Y:
Figure FDA0002575992350000039
(42) from the small-signal models established in (13), (24), (32), and (34), it can be found that:
Figure FDA00025759923500000310
(43) by substituting the formula in (42) into the formula for admittance Y in (41), Y can be obtainedd+YqWherein
Figure FDA0002575992350000041
(44) Let d-axis current be given as 0, U is known from the steady state modeld=-ωeLqIq,UdcIdc=3UqIq/2, mixing U withd=-ωeLqIqAnd UdcIdc=3UqIqThe formula of substituting/2 into (43) can be obtained:
Figure FDA0002575992350000042
7. the stability analysis method for a permanent magnet synchronous motor drive system according to claim 6, characterized in that: according to the permanent magnet synchronous motor driving system model, the stability analysis of the permanent magnet synchronous motor driving system by adopting a Nyquist stability criterion comprises the following steps:
establishing a transfer function of a permanent magnet synchronous motor driving system;
obtaining a permanent magnet synchronous motor driving system stability condition by adopting a Nyquist stability criterion and combining a permanent magnet synchronous motor driving system transfer function;
and analyzing the stability of the permanent magnet synchronous motor driving system according to the stable condition.
8. The stability analysis method for a permanent magnet synchronous motor drive system according to claim 7, characterized in that: the method for establishing the transfer function of the permanent magnet synchronous motor driving system comprises the following steps:
according to an impedance analysis method, equivalence is carried out on an admittance Y for a permanent magnet synchronous motor driving system, and impedance Z for a direct current end filterinIn an equivalent manner, the first and second electrodes,
Figure FDA0002575992350000043
Rfas a lineResistance, LfIs a filter inductor and C is a support capacitor.
According to circuit principles, the transfer function of the established system is:
Figure FDA0002575992350000044
wherein u isinIs the equivalent supply voltage.
9. The stability analysis method for a permanent magnet synchronous motor drive system according to claim 8, characterized in that: and (3) obtaining the stability condition of the permanent magnet synchronous motor driving system by adopting a Nyquist stability criterion and combining the transfer function of the permanent magnet synchronous motor driving system as follows:
Figure FDA0002575992350000051
10. a stability analytic system for permanent magnet synchronous machine actuating system, its characterized in that: the stability analysis method for a permanent magnet synchronous motor drive system according to any one of claims 1 to 9 is operated to analyze the stability of the permanent magnet synchronous motor drive system.
CN202010653858.3A 2020-05-21 2020-07-08 Stability analysis method and system for permanent magnet synchronous motor driving system Active CN111931343B (en)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
CN202010438021 2020-05-21
CN2020104380217 2020-05-21

Publications (2)

Publication Number Publication Date
CN111931343A true CN111931343A (en) 2020-11-13
CN111931343B CN111931343B (en) 2023-03-24

Family

ID=73313785

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010653858.3A Active CN111931343B (en) 2020-05-21 2020-07-08 Stability analysis method and system for permanent magnet synchronous motor driving system

Country Status (1)

Country Link
CN (1) CN111931343B (en)

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9316701B1 (en) * 2014-05-06 2016-04-19 The Florida State University Research Foundation, Inc. Real-time small-signal stability assessment of power electronic-based components in contemporary power systems
CN106934174A (en) * 2017-03-27 2017-07-07 上海海事大学 A kind of modeling of Electrical Propulsion Ship experiment porch and method for analyzing stability
CN108988391A (en) * 2018-08-16 2018-12-11 西南交通大学 The method for analyzing stability of double-fed fan rotor side converter based on revolving speed control
CN110148967A (en) * 2019-06-18 2019-08-20 华北电力大学 A kind of research method based on the straight drive blower sub-synchronous oscillation characteristic of admittance analysis
CN110768299A (en) * 2019-12-02 2020-02-07 湖南大学 Sequence impedance modeling and stability analysis method of load virtual synchronous machine
CN110994668A (en) * 2019-10-29 2020-04-10 国网甘肃省电力公司电力科学研究院 Stability analysis method based on output impedance model of doubly-fed wind power plant grid-connected system

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9316701B1 (en) * 2014-05-06 2016-04-19 The Florida State University Research Foundation, Inc. Real-time small-signal stability assessment of power electronic-based components in contemporary power systems
CN106934174A (en) * 2017-03-27 2017-07-07 上海海事大学 A kind of modeling of Electrical Propulsion Ship experiment porch and method for analyzing stability
CN108988391A (en) * 2018-08-16 2018-12-11 西南交通大学 The method for analyzing stability of double-fed fan rotor side converter based on revolving speed control
CN110148967A (en) * 2019-06-18 2019-08-20 华北电力大学 A kind of research method based on the straight drive blower sub-synchronous oscillation characteristic of admittance analysis
CN110994668A (en) * 2019-10-29 2020-04-10 国网甘肃省电力公司电力科学研究院 Stability analysis method based on output impedance model of doubly-fed wind power plant grid-connected system
CN110768299A (en) * 2019-12-02 2020-02-07 湖南大学 Sequence impedance modeling and stability analysis method of load virtual synchronous machine

Also Published As

Publication number Publication date
CN111931343B (en) 2023-03-24

Similar Documents

Publication Publication Date Title
Corradini et al. Analysis of multisampled current control for active filters
CN111769591A (en) Double-mode combined control method for multi-inverter system based on double split transformers
WO2021169666A1 (en) Resonance control method for differentiated phase correction
CN110492774B (en) High-power fractional order impedance element implementation circuit and control method thereof
CN108418226B (en) Reactive compensation control method of open-winding double-inverter photovoltaic power generation system
Choi et al. Compensation of zero-current clamping effects in high-frequency-signal-injection-based sensorless PM motor drives
Al-Kaf et al. Hybrid current controller for permanent-magnet synchronous motors using robust switching techniques
Lin et al. A novel modulation-based current harmonic control strategy for PMSM considering current measurement error and asymmetric impedance
CN113422533A (en) Vector angle proportional-integral control method
CN113437895A (en) Matrix vector angle proportional resonance control method
Yepes et al. Torque ripple minimization in surface-mounted PM drives by means of PI+ multi-resonant controller in synchronous reference frame
CN111931343B (en) Stability analysis method and system for permanent magnet synchronous motor driving system
Zhang et al. Model predictive flux control of three-level inverter-fed induction motor drives based on space vector modulation
CN114301361B (en) Control method of electrolytic capacitor-free permanent magnet synchronous motor driving system based on bus current control
Lin et al. Input impedance characterization of a single-phase PFC in DQ frame
CN109150050B (en) Stator resistance identification method, motor controller and computer readable storage medium
Zhang et al. Adaptive PI parameter of flux-weakening controller based on voltage feedback for model predictive control of SPMSM
Lu et al. Improved sliding mode observer for position sensorless open-winding permanent magnet brushless motor drives
Abosh et al. Current control of permanent magnet synchronous machine with asymmetric phases
CN105633969A (en) Active power filter power source harmonic compensation system and method
CN113809967B (en) Permanent magnet synchronous motor modulation model prediction control method based on matrix converter
CN114512991A (en) Active harmonic and reactive power generation system and method
Zhu et al. Interior permanent magnet synchronous motor dead-time compensation combined with extended Kalman and neural network bandpass filter
Yang et al. Impedance Modeling and Stability Analysis of Three-phase Grid-connected Inverter with LCL Filter in dq-frame
Thapa et al. Torque ripple reduction in a traction ipmsm with resistance asymmetry using an adaptive pir current controller

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant