CN108595798A - A kind of load model parameters computational methods suitable for electromechanical transient simulation - Google Patents
A kind of load model parameters computational methods suitable for electromechanical transient simulation Download PDFInfo
- Publication number
- CN108595798A CN108595798A CN201810322844.6A CN201810322844A CN108595798A CN 108595798 A CN108595798 A CN 108595798A CN 201810322844 A CN201810322844 A CN 201810322844A CN 108595798 A CN108595798 A CN 108595798A
- Authority
- CN
- China
- Prior art keywords
- formula
- induction conductivity
- model
- simulation
- rated
- 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.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/30—Circuit design
- G06F30/36—Circuit design at the analogue level
- G06F30/367—Design verification, e.g. using simulation, simulation program with integrated circuit emphasis [SPICE], direct methods or relaxation methods
Landscapes
- Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Microelectronics & Electronic Packaging (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Control Of Ac Motors In General (AREA)
Abstract
The present invention relates to the technical fields of power system transient simulation, more particularly, to a kind of load model parameters computational methods suitable for electromechanical transient simulation, including:Based on induction conductivity manufacture nameplate data, the mathematical model that the optimization of induction conductivity simulation parameters calculates is established;The induction conductivity simulation model is the Type Equivalent Circuit Model for ignoring excitatory resistance, the mathematical model is equal to corresponding name plate rating as equality constraint using stator current, input reactive power and maximum electromagnetic torque, using computational efficiency and the relative deviation minimum of name plate rating as object function;Model is solved using sequential quadratic programming SQP algorithms, obtains three rank model parameters of electromechanical transient state simulation of induction conductivity.The model parameter that the present invention obtains induction conductivity is more accurate, calculates simply, and has preferable analysis precision.
Description
Technical field
The present invention relates to the technical fields of power system transient simulation, being suitable for electromechanical transient more particularly, to one kind
The load model parameters computational methods of emulation.
Background technology
With the development of distant-range high voltage direct current transmission construction, multi-infeed DC large size receiving-end system has been gradually formed.By
Ratio of powering outside end system area is big, and the motor property dynamic load such as industry, air-conditioning ratio is high so that voltage stabilization sex chromosome mosaicism day
Benefit is prominent.It is a dynamic problem in voltage stability question essence, many dynamic factors play an important role in system.
Wherein, dynamic load characteristic has it significant impact, and the accuracy of load model and parameter has become restriction electric system
The key of Voltage stability analysis precision.Induction conductivity occupies larger specific gravity in power system load.It is counted in power train
It calculates in analysis software package, synthetic load is all made of the model of induction conductivity and static load parallel connection.Determine that induction conductivity is negative
Lotus model parameter is necessary to power system load modeling and transient stability simulation calculation.
Currently, several allusion quotations that the induction motor model parameter used is directly recommended with reference to IEEE load modelings group 1987
Shape parameter.However, these typical model parameters to be applied to the simulation analysis of current alternating current-direct current bulk power grid, adaptability and accurate
Property can be restricted.The model parameter that acquisition induction conductivity is surveyed based on experiment is the most accurate, but motor in electric system
Type, the quantity of load are very more, and measurement method is also restrained in the application.
Invention content
It is an object of the invention to overcome the deficiencies of the prior art and provide a kind of load moulds suitable for electromechanical transient simulation
Shape parameter computational methods are established the optimization of induction conductivity simulation parameters and are calculated based on induction conductivity manufacture nameplate data
Mathematical model;Model is solved using sequential quadratic programming SQP algorithms, it is imitative to obtain three rank electromechanical transient of induction conductivity
True mode parameter.
In order to solve the above technical problems, the technical solution adopted by the present invention is:
A kind of load model parameters computational methods suitable for electromechanical transient simulation, which is characterized in that include the following steps:
S1. induction conductivity manufacture nameplate data are based on, the number that the optimization of induction conductivity simulation parameters calculates is established
Learn model;The induction conductivity simulation model is the Type Equivalent Circuit Model for ignoring excitatory resistance, and the mathematical model is with stator
It is equality constraint that electric current, input reactive power and maximum electromagnetic torque, which are equal to corresponding name plate rating, with computational efficiency and inscription
The relative deviation minimum of board value is as object function;
S2. sequential quadratic programming SQP algorithms are based on and calculate three rank electromechanical transient simulation analog parameter of induction conductivity.
Preferably, equivalent circuit described in step S1 includes the stator resistance R being connected in series withs, stator reactance XsAnd it is excitatory
Reactance Xm, the excitatory reactance XmBoth ends be connected in parallel to single cage model rotor reactance XrAnd single cage model rotor resistance Rr。
Preferably, nameplate data of dispatching from the factory described in step S1 include specified mechanical output Pn, rated voltage Un, rated current
In, rated speed nn, rated efficiency ηn, rated power factor cos θn, torque capacity multiple Km, stall or starting torque multiple Kst
And stall or starting current times KIst。
Preferably, Km=Tem/Tn, TnFor nominal torque, TemFor maximum or critical electromagnetic torque;Kst=Tst/Tn, TstTo rise
Dynamic torque;KIst=Ist/In, IstFor starting current, InFor rated current.
Preferably, step S1 includes the following steps:
S11. the stator resistance R of induction conductivity is calculateds;
S12. induction motor load simulation parameters X=[X are calculateds, Xm, Xr, Rr]。
Preferably, stator resistance R in step S11sComputational methods it is as follows:
1) the rated slip s of induction conductivity is calculatedn:
In formula (1), nnFor induction conductivity rated speed, ns=60f/P is induction conductivity synchronous rotational speed;F is system frequency
Rate, P are number of pole-pairs;
2) the specified active input power P of induction conductivity is calculated1With nominal reactive input power Q1:
In formula (2), UnFor rated voltage, InFor rated current, cos θnFor rated power factor;
3) stator winding resistance R is calculated according to stator copper losss:
In formula (3)~formula (4), P1For specified active input power, InFor rated current, PenFor specified electromagnetic power, PnFor
Specified mechanical output, snFor the rated slip of induction conductivity.
Preferably, the computational methods of induction motor load simulation parameters are as follows in step S12:
1) corresponding name plate rating is equal to as equality constraint item using stator current, input reactive power and maximum electromagnetic torque
Part is established using computational efficiency and the relative deviation minimum of name plate rating as object function and calculates induction motor load emulation mould
Shape parameter X=[Xs, Xm, Xr, Rr] mathematical model:
In formula (5)~(6), induction electric engine efficiency is η (sn)=Pn(sn)/P1(sn);
2) according to the nameplate rated power P of induction conductivityn, slippage snWith maximum electromagnetic torque multiple Km, calculate induced electricity
The specified electromagnetic torque T of motivationnWith maximum electromagnetic torque Tem:
In formula (7)~(8), Ωs=2 π f/P are synchronous angular velocity;
3) stator impedance of induction conductivity single cage model and excitatory reaction component are subjected to Thevenin's equivalence, equivalent impedance
ZthFor:
In formula (9), ZthFor the Thevenin's equivalence impedance of motor stator and excitatory reactance, RthFor ZthReal part, XthFor
ZthImaginary part, j is imaginary number;
4) according to Thevenin's equivalence circuit, breakdown slip s is calculatedm:
5) equivalent impedance Zm(s) it is about the calculation expression of s:
In formula (11), Zr(s)=Rr/s+jXr, it is induction conductivity rotor impedance;
6) stator current of induction conductivityAnd rotor currentCalculation expression about s is:
In formula (12),For stator terminal phase voltage.
7) the electromagnetic torque T of induction conductivitye(s) it is about the calculation expression of s:
8) the active input power P calculated according to induction motor simulation model1(s), idle input power Q1(s) and machinery is defeated
Go out power Pn(s) it is:
9) the rated slip s that formula (1) is calculatednSubstitution formula (12) and formula (14), are calculated corresponding to snStator
Electric current and rotor current, input power P1(sn) and Q1(sn), output mechanical power Pn(sn);Formula (10) is calculated critical
Slippage smIt substitutes into (13) and obtains maximum electromagnetic torque Te(sm)。
Preferably, the calculating of step S2 includes the following steps:
S21.SQP algorithms solve Solution of Nonlinear Optimal Problem general mathematics model be:
min F(X) (15)
In formula (15), F (X) is object function, is a nonlinear function;In formula (15), Hi(X)=0 it is non-linear etc.
Formula constraints, Gj(X)≤0 it is nonlinear complementary problem condition, LB≤X≤UB is variable bound constraints;
S22. in iteration point Xk, quadratic programming subproblem corresponding with formula (17) and formula (18) is expressed as form:
In formula (17), ▽ F (Xk) it is F (Xk) in XkThe derivative at place;In formula (18), ▽ H (Xk) it is H (Xk) in xkPlace is led
Number, ▽ G (Xk) it is G (Xk) in XkThe derivative at place;dkAnd λ *k*、vk *For the optimal solution of formula (17), respectively kth time iteration point XkPlace
The direction of search and Lagrange multiplier correction amount;B(Xk,λk) be Lagrangian second dervative true Hessian squares
Battle array, wherein Lagrangian is:
L (X, λ, ν)=F (X)+λTH(X)-νTG(X) (19)
The true Hessian matrixes of second dervative are:
S23. BFGS is used to correct, when correcting the non-timing of curvature, using damping BFGS modification methods:
In formula (21), yk=▽XL(Xk+1,λk,νk)-▽XL(Xk,λk,νk), sk=Xk+1-Xk;
S24. next iteration point X is obtained by formula (22) based on SQP algorithmsk+1:
S25. the parameter X=[X of induction motor simulation model are givens, Xm, Rr, Xr] assign initial value:
S26. with rated voltage UnWith specified mechanical output PnOn the basis of be worth, by the famous value reduction of parameter be perunit value, calculate
Motor inertia time constant H:
The load model parameters computational methods suitable for electromechanical transient simulation of the present invention, based on induction conductivity manufacture inscription
Board data establish the mathematical model that the optimization of induction conductivity simulation parameters calculates;Using sequential quadratic programming SQP algorithms pair
Model is solved, and three rank model parameters of electromechanical transient state simulation of induction conductivity is obtained.The present invention obtains the mould of induction conductivity
Shape parameter is more accurate, and has preferable analysis precision.
Description of the drawings
Fig. 1 is the structure of the equivalent circuit of the load model parameters computational methods suitable for electromechanical transient simulation of the present invention
Schematic diagram.
Fig. 2 is the flow chart that sequential quadratic programming algorithm calculates induction conductivity parameter.
Fig. 3 is 15 induction conductivities manufacture nameplate data in embodiment two.
Fig. 4 is 15 induction conductivity simulation parameters being calculated.
Fig. 5 is the comparison I of induction conductivity simulation model calculation of performance indicators and the name plate rating that dispatches from the factory.
Fig. 6 is the comparison II of induction conductivity simulation model calculation of performance indicators and the name plate rating that dispatches from the factory.
Fig. 7 is the comparison III of induction conductivity simulation model calculation of performance indicators and the name plate rating that dispatches from the factory.
Specific implementation mode
The present invention is further illustrated With reference to embodiment.Wherein, attached drawing only for illustration,
What is indicated is only schematic diagram rather than pictorial diagram, should not be understood as the limitation to this patent;Reality in order to better illustrate the present invention
Example is applied, the certain components of attached drawing have omission, zoom in or out, and do not represent the size of actual product;To those skilled in the art
For, the omitting of some known structures and their instructions in the attached drawings are understandable.
The same or similar label correspond to the same or similar components in the attached drawing of the embodiment of the present invention;In retouching for the present invention
In stating, it is to be understood that if it is based on attached drawing to have the orientation or positional relationship of the instructions such as term "upper", "lower", "left", "right"
Shown in orientation or positional relationship, be merely for convenience of description of the present invention and simplification of the description, do not indicate or imply the indicated
Device or element must have a particular orientation, with specific azimuth configuration and operation, therefore position relationship described in attached drawing
Term only for illustration, should not be understood as the limitation to this patent, for the ordinary skill in the art, can
To understand the concrete meaning of above-mentioned term as the case may be.
Embodiment 1
It is as shown in Figure 1 to Figure 2 the of the load model parameters computational methods suitable for electromechanical transient simulation of the present invention
One embodiment, includes the following steps:
S1. induction conductivity manufacture nameplate data are based on, the number that the optimization of induction conductivity simulation parameters calculates is established
Learn model;The induction conductivity simulation model is the Type Equivalent Circuit Model for ignoring excitatory resistance, and the mathematical model is with stator
It is equality constraint that electric current, input reactive power and maximum electromagnetic torque, which are equal to corresponding name plate rating, with computational efficiency and inscription
The relative deviation minimum of board value is as object function;
S2. sequential quadratic programming SQP algorithms are based on and calculate three rank electromechanical transient simulation analog parameter of induction conductivity.
Wherein, as shown in Figure 1, equivalent circuit described in step S1 includes the stator resistance R being connected in series withs, stator reactance Xs
And excitatory reactance Xm, the excitatory reactance XmBoth ends be connected in parallel to single cage model rotor reactance XrAnd single cage model turns
Sub- resistance Rr.Nameplate data of dispatching from the factory include specified mechanical output Pn, rated voltage Un, rated current In, rated speed nn, specified effect
Rate ηn, rated power factor cos θn, torque capacity multiple Km, stall or starting torque multiple KstAnd stall or starting current are again
Number KIst.Wherein, Km=Tem/Tn, TnFor nominal torque, TemFor maximum or critical electromagnetic torque;Kst=Tst/Tn, TstTo play turn
Square;KIst=Ist/In, IstFor starting current, InFor rated current.
Step S1 includes the following steps:
S11. the stator resistance R of induction conductivity is calculateds;
S12. induction motor load simulation parameters X=[X are calculateds, Xm, Xr, Rr]。
Wherein, stator resistance R in step S11sComputational methods it is as follows:
1) the rated slip s of induction conductivity is calculatedn:
In formula (1), nnFor induction conductivity rated speed, ns=60f/P is induction conductivity synchronous rotational speed;F is system frequency
Rate, P are number of pole-pairs;
2) the specified active input power P of induction conductivity is calculated1With nominal reactive input power Q1:
In formula (2), UnFor rated voltage, InFor rated current, cos θnFor rated power factor;
3) stator winding resistance R is calculated according to stator copper losss:
In formula (3)~formula (4), P1For specified active input power, InFor rated current, PenFor specified electromagnetic power, PnFor
Specified mechanical output, snFor the rated slip of induction conductivity.
Preferably, as shown in Fig. 2, the computational methods of induction motor load simulation parameters are as follows in step S12:
1) corresponding name plate rating is equal to as equality constraint item using stator current, input reactive power and maximum electromagnetic torque
Part is established using computational efficiency and the relative deviation minimum of name plate rating as object function and calculates induction motor load emulation mould
Shape parameter X=[Xs, Xm, Xr, Rr] mathematical model:
In formula (5)~(6), induction electric engine efficiency is η (sn)=Pn(sn)/P1(sn);
2) according to the nameplate rated power P of induction conductivityn, slippage snWith maximum electromagnetic torque multiple Km, calculate induced electricity
The specified electromagnetic torque T of motivationnWith maximum electromagnetic torque Tem:
In formula (7)~(8), Ωs=2 π f/P are synchronous angular velocity;
3) stator impedance of induction conductivity single cage model and excitatory reaction component are subjected to Thevenin's equivalence, equivalent impedance
ZthFor:
In formula (9), ZthFor the Thevenin's equivalence impedance of motor stator and excitatory reactance, RthFor ZthReal part, XthFor
ZthImaginary part, j is imaginary number;
4) according to Thevenin's equivalence circuit, breakdown slip s is calculatedm:
5) equivalent impedance Zm(s) it is about the calculation expression of s:
In formula (11), Zr(s)=Rr/s+jXr, it is induction conductivity rotor impedance;
6) stator current of induction conductivityAnd rotor currentCalculation expression about s is:
In formula (12),For stator terminal phase voltage.
7) the electromagnetic torque T of induction conductivitye(s) it is about the calculation expression of s:
8) the active input power P calculated according to induction motor simulation model1(s), idle input power Q1(s) and machinery is defeated
Go out power Pn(s) it is:
9) the rated slip s that formula (1) is calculatednSubstitution formula (12) and formula (14), are calculated corresponding to snStator
Electric current and rotor current, input power P1(sn) and Q1(sn), output mechanical power Pn(sn);Formula (10) is calculated critical
Slippage smIt substitutes into (13) and obtains maximum electromagnetic torque Te(sm)。
The calculating of step S2 includes the following steps:
S21.SQP algorithms solve Solution of Nonlinear Optimal Problem general mathematics model be:
min F(X) (15)
In formula (15), F (X) is object function, is a nonlinear function;In formula (15), Hi(X)=0 it is non-linear etc.
Formula constraints, Gj(X)≤0 it is nonlinear complementary problem condition, LB≤X≤UB is variable bound constraints;
S22. in iteration point Xk, quadratic programming subproblem corresponding with formula (17) and formula (18) is expressed as form:
In formula (17), ▽ F (Xk) it is F (Xk) in XkThe derivative at place;In formula (18), ▽ H (Xk) it is H (Xk) in xkPlace is led
Number, ▽ G (Xk) it is G (Xk) in XkThe derivative at place;dkAnd λ *k*、vk *For the optimal solution of formula (17), respectively kth time iteration point XkPlace
The direction of search and Lagrange multiplier correction amount;B(Xk,λk) be Lagrangian second dervative true Hessian squares
Battle array, wherein Lagrangian is:
L (X, λ, ν)=F (X)+λTH(X)-νTG(X) (19)
The true Hessian matrixes of second dervative are:
S23. BFGS is used to correct, when correcting the non-timing of curvature, using damping BFGS modification methods:
In formula (21), yk=▽XL(Xk+1,λk,νk)-▽XL(Xk,λk,νk), sk=Xk+1-Xk;
S24. next iteration point X is obtained by formula (22) based on SQP algorithmsk+1:
S25. the parameter X=[X of induction motor simulation model are givens, Xm, Rr, Xr] assign initial value:
S26. with rated voltage UnWith specified mechanical output PnOn the basis of be worth, by the famous value reduction of parameter be perunit value, calculate
Motor inertia time constant H:
Embodiment two
The present embodiment calculates 15 induction conductivity simulation parameters, 15 induced electricities by the method for embodiment one
The manufacture nameplate data of motivation are as shown in figure 3, induction conductivity changed power ranging from 22~1400kW, the variation of rotary inertia
Ranging from 0.76~673kg.m2, rated voltage is respectively 380V, 6kV and 10kV, rated operation frequency 50Hz.Using implementation
The computational methods that example one is proposed calculate the corresponding simulation parameters of each induction conductivity, result of calculation such as Fig. 4 institutes in Fig. 3
Show.
The induction conductivity simulation parameters being calculated according to Fig. 4, calculate the rated active power of induction conductivity
P1(kW), rated reactive power Q1(kvar), rated current In(A), specified mechanical output Pn(kW), rated efficiency ηn, specified work(
Rate factor cos θn, torque capacity multiple Km, starting torque multiple Kst, starting current times KIstAnd with manufacture nameplate data into
Row comparison, is as a result shown in shown in Fig. 5 to Fig. 7.Shown in Fig. 5 to Fig. 7, induced electricity that the method that is provided by embodiment one is calculated
The parameter for the induction conductivity that motivation simulation parameters are calculated corresponds to the difference very little of parameter on nameplate, and the present invention obtains
The accuracy of the model parameter of the induction conductivity taken is verified.
Obviously, the above embodiment of the present invention be only to clearly illustrate example of the present invention, and not be pair
The restriction of embodiments of the present invention.For those of ordinary skill in the art, may be used also on the basis of the above description
To make other variations or changes in different ways.There is no necessity and possibility to exhaust all the enbodiments.It is all this
All any modification, equivalent and improvement etc., should be included in the claims in the present invention made by within the spirit and principle of invention
Protection domain within.
Claims (8)
1. a kind of load model parameters computational methods suitable for electromechanical transient simulation, which is characterized in that include the following steps:
S1. induction conductivity manufacture nameplate data are based on, the mathematical modulo that the optimization of induction conductivity simulation parameters calculates is established
Type;The induction conductivity simulation model is the Type Equivalent Circuit Model for ignoring excitatory resistance, the mathematical model with stator current,
Input reactive power and maximum electromagnetic torque to be equal to corresponding name plate rating be equality constraint, with computational efficiency and name plate rating
Relative deviation minimum is as object function;
S2. sequential quadratic programming SQP algorithms are based on and calculate three rank electromechanical transient simulation analog parameter of induction conductivity.
2. the load model parameters computational methods according to claim 1 suitable for electromechanical transient simulation, which is characterized in that
Equivalent circuit described in step S1 includes the stator resistance R being connected in series withs, stator reactance XsAnd excitatory reactance Xm, described excitatory
Reactance XmBoth ends be connected in parallel to single cage model rotor reactance XrAnd single cage model rotor resistance Rr。
3. the load model parameters computational methods according to claim 1 suitable for electromechanical transient simulation, which is characterized in that
Nameplate data of dispatching from the factory described in step S1 include specified mechanical output Pn, rated voltage Un, rated current In, rated speed nn, volume
Determine efficiency etan, rated power factor cos θn, torque capacity multiple Km, stall or starting torque multiple KstAnd stall or start electricity
Flow multiple KIst。
4. the load model parameters computational methods according to claim 3 suitable for electromechanical transient simulation, which is characterized in that
Km=Tem/Tn, TnFor nominal torque, TemFor maximum or critical electromagnetic torque;Kst=Tst/Tn, TstFor starting torque;KIst=
Ist/In, IstFor starting current, InFor rated current.
5. the load model parameters computational methods according to claim 1 suitable for electromechanical transient simulation, which is characterized in that
Step S1 includes the following steps:
S11. the stator resistance R of induction conductivity is calculateds;
S12. induction motor load simulation parameters X=[X are calculateds, Xm, Xr, Rr]。
6. the load model parameters computational methods according to claim 5 suitable for electromechanical transient simulation, which is characterized in that
Stator resistance R in step S11sComputational methods it is as follows:
1) the rated slip s of induction conductivity is calculatedn:
In formula (1), nnFor induction conductivity rated speed, ns=60f/P is induction conductivity synchronous rotational speed;F is system frequency, P
For number of pole-pairs;
2) the specified active input power P of induction conductivity is calculated1With nominal reactive input power Q1:
In formula (2), UnFor rated voltage, InFor rated current, cos θnFor rated power factor;
3) stator winding resistance R is calculated according to stator copper losss:
In formula (3)~formula (4), P1For specified active input power, InFor rated current, PenFor specified electromagnetic power, PnIt is specified
Mechanical output, snFor the rated slip of induction conductivity.
7. the load model parameters computational methods according to claim 1 suitable for electromechanical transient simulation, which is characterized in that
The computational methods of induction motor load simulation parameters are as follows in step S12:
1) corresponding name plate rating is equal to as equality constraint using stator current, input reactive power and maximum electromagnetic torque,
Using computational efficiency and the relative deviation minimum of name plate rating as object function, establishes and calculate induction motor load simulation model ginseng
Number X=[Xs, Xm, Xr, Rr] mathematical model:
In formula (5)~(6), induction electric engine efficiency is η (sn)=Pn(sn)/P1(sn);
2) according to the nameplate rated power P of induction conductivityn, slippage snWith maximum electromagnetic torque multiple Km, calculate induction conductivity
Specified electromagnetic torque TnWith maximum electromagnetic torque Tem:
In formula (7)~(8), Ωs=2 π f/P are synchronous angular velocity;
3) stator impedance of induction conductivity single cage model and excitatory reaction component are subjected to Thevenin's equivalence, equivalent impedance Zth
For:
In formula (9), ZthFor the Thevenin's equivalence impedance of motor stator and excitatory reactance, RthFor ZthReal part, XthFor ZthVoid
Portion, j are imaginary number;
4) according to Thevenin's equivalence circuit, breakdown slip s is calculatedm:
5) equivalent impedance Zm(s) it is about the calculation expression of s:
In formula (11), Zr(s)=Rr/s+jXr, it is induction conductivity rotor impedance;
6) stator current of induction conductivityAnd rotor currentCalculation expression about s is:
In formula (12),For stator terminal phase voltage.
7) the electromagnetic torque T of induction conductivitye(s) it is about the calculation expression of s:
8) the active input power P calculated according to induction motor simulation model1(s), idle input power Q1(s) and mechanical output work
Rate Pn(s) it is:
9) the rated slip s that formula (1) is calculatednSubstitution formula (12) and formula (14), are calculated corresponding to snStator current
With rotor current, input power P1(sn) and Q1(sn), output mechanical power Pn(sn);The critical slippage that formula (10) is calculated
smIt substitutes into (13) and obtains maximum electromagnetic torque Te(sm)。
8. the load model parameters computational methods according to claim 1 suitable for electromechanical transient simulation, which is characterized in that
The calculating of step S2 includes the following steps:
S21.SQP algorithms solve Solution of Nonlinear Optimal Problem general mathematics model be:
min F(X) (15)
In formula (15), F (X) is object function, is a nonlinear function;In formula (15), Hi(X)=0 it is Nonlinear Equality Constrained
Condition, Gj(X)≤0 it is nonlinear complementary problem condition, LB≤X≤UB is variable bound constraints;
S22. in iteration point Xk, quadratic programming subproblem corresponding with formula (17) and formula (18) is expressed as form:
In formula (17), ▽ F (Xk) it is F (Xk) in XkThe derivative at place;In formula (18), ▽ H (Xk) it is H (Xk) in xkThe derivative at place, ▽ G
(Xk) it is G (Xk) in XkThe derivative at place;dkAnd λ *k*、vk *For the optimal solution of formula (17), respectively kth time iteration point XkThe search at place
Direction and Lagrange multiplier correction amount;B(Xk,λk) be Lagrangian second dervative true Hessian matrixes, wherein
Lagrangian is:
L (X, λ, ν)=F (X)+λTH(X)-νTG(X) (19)
The true Hessian matrixes of second dervative are:
S23. BFGS is used to correct, when correcting the non-timing of curvature, using damping BFGS modification methods:
In formula (21), yk=▽XL(Xk+1,λk,νk)-▽XL(Xk,λk,νk), sk=Xk+1-Xk;
S24. next iteration point X is obtained by formula (22) based on SQP algorithmsk+1:
S25. the parameter X=[X of induction motor simulation model are givens, Xm, Rr, Xr] assign initial value:
S26. with rated voltage UnWith specified mechanical output PnOn the basis of be worth, by the famous value reduction of parameter be perunit value, calculate it is electronic
Machine inertia time constant H:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810322844.6A CN108595798A (en) | 2018-04-11 | 2018-04-11 | A kind of load model parameters computational methods suitable for electromechanical transient simulation |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810322844.6A CN108595798A (en) | 2018-04-11 | 2018-04-11 | A kind of load model parameters computational methods suitable for electromechanical transient simulation |
Publications (1)
Publication Number | Publication Date |
---|---|
CN108595798A true CN108595798A (en) | 2018-09-28 |
Family
ID=63621984
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810322844.6A Pending CN108595798A (en) | 2018-04-11 | 2018-04-11 | A kind of load model parameters computational methods suitable for electromechanical transient simulation |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108595798A (en) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109638830A (en) * | 2019-01-18 | 2019-04-16 | 广东电网有限责任公司 | A kind of electric load model building method, device and equipment |
CN110501614A (en) * | 2019-08-31 | 2019-11-26 | 广东电网有限责任公司 | Induction-motor load transient process calculation method suitable for electromagnetic transient simulation |
CN110580371A (en) * | 2019-05-10 | 2019-12-17 | 广东电网有限责任公司 | Motor model parameter conversion calculation method suitable for electromagnetic transient simulation program |
CN114677532A (en) * | 2022-04-20 | 2022-06-28 | 重庆邮电大学 | Induction motor parameter identification method by utilizing factory data and characteristic curve |
CN115441798A (en) * | 2022-09-01 | 2022-12-06 | 东南大学成贤学院 | Induction motor parameter identification method based on factory data and maximum torque formula |
KR20230003300A (en) * | 2020-05-08 | 2023-01-05 | 지멘스 악티엔게젤샤프트 | Method and system for providing simulation models of electric rotating machines |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101727522A (en) * | 2009-12-11 | 2010-06-09 | 中国电力科学研究院 | Acquisition method of model parameters of electromechanical transient state simulation of motor based on factory data |
-
2018
- 2018-04-11 CN CN201810322844.6A patent/CN108595798A/en active Pending
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101727522A (en) * | 2009-12-11 | 2010-06-09 | 中国电力科学研究院 | Acquisition method of model parameters of electromechanical transient state simulation of motor based on factory data |
Non-Patent Citations (1)
Title |
---|
毛晓明等: "基于铭牌数据和序列二次规划算法的感应电动机", 《电力自动化设备》 * |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109638830A (en) * | 2019-01-18 | 2019-04-16 | 广东电网有限责任公司 | A kind of electric load model building method, device and equipment |
CN109638830B (en) * | 2019-01-18 | 2022-03-22 | 广东电网有限责任公司 | Power load model construction method, device and equipment |
CN110580371A (en) * | 2019-05-10 | 2019-12-17 | 广东电网有限责任公司 | Motor model parameter conversion calculation method suitable for electromagnetic transient simulation program |
CN110580371B (en) * | 2019-05-10 | 2022-04-19 | 广东电网有限责任公司 | Motor model parameter conversion calculation method suitable for electromagnetic transient simulation program |
CN110501614A (en) * | 2019-08-31 | 2019-11-26 | 广东电网有限责任公司 | Induction-motor load transient process calculation method suitable for electromagnetic transient simulation |
WO2021036279A1 (en) * | 2019-08-31 | 2021-03-04 | 广东电网有限责任公司 | Motor load transient process calculation method suitable for electromagnetic transient simulation |
KR102559570B1 (en) | 2020-05-08 | 2023-07-25 | 지멘스 악티엔게젤샤프트 | Method and system for providing simulation models of electric rotating machines |
KR20230003300A (en) * | 2020-05-08 | 2023-01-05 | 지멘스 악티엔게젤샤프트 | Method and system for providing simulation models of electric rotating machines |
US11790135B2 (en) | 2020-05-08 | 2023-10-17 | Siemens Aktiengesellschaft | Method and systems for provision of a simulation model of an electric rotating machine |
CN114677532A (en) * | 2022-04-20 | 2022-06-28 | 重庆邮电大学 | Induction motor parameter identification method by utilizing factory data and characteristic curve |
CN114677532B (en) * | 2022-04-20 | 2024-05-28 | 重庆邮电大学 | Induction motor parameter identification method utilizing factory data and characteristic curve |
CN115441798A (en) * | 2022-09-01 | 2022-12-06 | 东南大学成贤学院 | Induction motor parameter identification method based on factory data and maximum torque formula |
CN115441798B (en) * | 2022-09-01 | 2023-09-15 | 东南大学成贤学院 | Induction motor parameter identification method based on factory data and maximum torque formula |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108595798A (en) | A kind of load model parameters computational methods suitable for electromechanical transient simulation | |
CN106357073B (en) | High winding coefficient permanent magnetic brushless and its design and fault tolerant control method | |
CN107590323A (en) | Electromagnetism, temperature rise, machinery based on the phase claw pole electric generator of composite excitation six, the modeling method of electrical coupling system | |
CN103853891A (en) | Finite element analysis-based variable-element permanent magnet synchronous motor modeling method | |
CN103995734A (en) | Electric power system hybrid real-time simulation system and simulation method based on RTDS | |
CN102609575A (en) | Power system transient stability simulating method based on implicit numerical integration | |
CN105720594B (en) | The compensation method of pure active three-phase unbalance load compensation capacity | |
CN105488315B (en) | The grid-connected analysis method for generating sub-synchronous oscillation of wind-driven generator | |
WO2014181367A1 (en) | Method and relevant system for converting mechanical energy from a generator actuated by a turbine into electric energy | |
CN103500245B (en) | A kind of field road transient state-transient state coupling simulation method based on m ulti-loop m ethod | |
CN102928672B (en) | Method for realizing resistance measurement of asynchronous motor fixed rotor | |
CN108923464A (en) | The building of doubly fed induction generator three-phase imbalance steady-state model, emulation mode and device | |
CN107645253A (en) | The three-phase simulation device of current-responsive type permagnetic synchronous motor and its drive system | |
CN105259462A (en) | Three-phase short circuit current calculating method for double-fed motor | |
CN102611380A (en) | Online identification method for parameters of double-fed motor | |
CN110501614A (en) | Induction-motor load transient process calculation method suitable for electromagnetic transient simulation | |
CN102545263B (en) | Power system transient stability simulation method based on explicit numerical integration | |
CN115021638A (en) | Electromagnetic transient modeling method, system and equipment of efficient synchronous machine | |
CN102918406A (en) | AC electric charge measurement device, and AC electric charge measurement method | |
Elizondo et al. | Model reduction, validation, and calibration of wind power plants for dynamic studies | |
CN108573094A (en) | The method for building up and system of the VBR electromagnetic transient simulation models of synchronous generator | |
CN105184418A (en) | Optimal power flow calculating method for equivalent interconnected power network on the basis of consistency of power flow, sensitivity and constraint | |
CN109359266A (en) | Method for solving, system, equipment and the storage medium of induction conductivity transient response | |
CN103208795A (en) | Method for optimizing parameters of supplementary sub-synchronous damping controller for high-voltage direct-current transmission | |
CN110580371B (en) | Motor model parameter conversion calculation method suitable for electromagnetic transient simulation program |
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 | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20180928 |