CN109063408A - A kind of electromagnetism electromechanical transient simulation algorithm based on straight algorithm - Google Patents
A kind of electromagnetism electromechanical transient simulation algorithm based on straight algorithm Download PDFInfo
- Publication number
- CN109063408A CN109063408A CN201811285344.6A CN201811285344A CN109063408A CN 109063408 A CN109063408 A CN 109063408A CN 201811285344 A CN201811285344 A CN 201811285344A CN 109063408 A CN109063408 A CN 109063408A
- Authority
- CN
- China
- Prior art keywords
- electro
- previous frame
- equation
- transient
- magnetic transient
- 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
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
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)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Supply And Distribution Of Alternating Current (AREA)
- Complex Calculations (AREA)
Abstract
The invention discloses a kind of electromagnetism electromechanical transient simulation algorithm based on straight algorithm, belong to Simulating technique in Electric Power System field, it include generator model, circuit model, load module and transformer model in power system network, it calculates separately the electro-magnetic transient matrix of generator model, circuit model, load module and transformer model comprising following steps: (1) corresponding electro-magnetic transient equation being solved according to different models respectively；(2) the electro-magnetic transient equation in step (1) is organized into difierence equation；(3) difierence equation is arranged and deforms to obtain electro-magnetic transient matrix；And substitute into electro-magnetic transient matrix in the Electrical-Machanical Transient Simulation of Power System method based on straight algorithm, to carry out the emulation of electro-magnetic transient effect, it can be realized electromagnetism, electromechanical transient in the simulation calculation of one to reach.
Description
Technical field
The invention belongs to Simulating technique in Electric Power System fields, electromechanical in particular to a kind of electromagnetism based on straight algorithm
Transient emulation algorithm.
Background technique
In electric system electromechanical transient process and electromagnetic transient be two with different mathematics characterization, have difference
The physical process of time constant.In traditional Power System Analysis tool, it is imitative that usually the two processes are carried out with number respectively
Very.Relative to electromechanical transient simulation, power equipment mathematical model is more complicated under electro-magnetic transient, and dynamic time constant of interest is more
It is small, it considerably increases the calculation amount of emulation and calculates the time, realize that the real-timedigital simulation of certain scale system wants much more difficult.
Electromagnetic transient Digital Simulation be with numerical computation method in electric system from several microseconds between the several seconds
Electromagnetic transient carries out analogue simulation.Transmission Line Distributed Parameter characteristic and parameter must be taken into consideration in electromagnetic transient emulation
Frequency characteristic, the electromagnetism of generator and electromechanical transient process and a series of elements (arrester, transformer, reactor etc.) it is non-
Linear characteristic.Therefore, the mathematical model of electromagnetic transient simulation must be set up the algebra or differential, partial differential of these elements and system
Equation.The numerical integration method generally used is implicit integration algorithm.
Since electromagnetic transient simulation does not require nothing more than to the dynamic element of electric system using detailed nonlinear model, also want
The transient process of meter and network, also needs to use differential equation, so that the simulation scale of electromagnetic transients program receives
Limitation.When general progress electromagnetic transient simulation, equivalent simplification will be carried out to electric system.
The emulation of electromechanical transient process, it is main to study electric system by the transient stability after large disturbances and by microvariations
Steady stability performance afterwards.Wherein transient stability analysis is research electric system by such as short trouble, cuts off route, power generation
Machine, load, generator lose under the effect of the large disturbances such as excitation or impact load, and the dynamic behaviour and holding of electric system are same
Walk the ability of stable operation.
The algorithm of Electrical-Machanical Transient Simulation of Power System is simultaneous solution electric system differential equation group and Algebraic Equation set, with
Obtain the time solution of physical quantity.The method for solving of differential equation group mainly has hiding-trapezium integral method, improves outstanding daraf(reciprocal of farad), Long Ge-
Ku Tafa etc., wherein hiding-trapezium integral method is more and more applied due to better numerical value stability.Algebraic Equation set
Method for solving is mainly using the Newton method for being suitable for solving Groebner Basis.According to asking for the differential equation and algebraic equation
Solution sequence can be divided into alternately solution and simultaneous solution method.
In Patent No. ZL201410142938.7, denomination of invention are as follows: the more power supplys of three-phase symmetrical of a word chain and attachment formula
The non-looped network of the more power supplys of three-phase symmetrical that a kind of operation result is accurate, arithmetic speed is fast is provided in the non-straight algorithm of looped network electric system
The straight algorithm of trend；
Application No. is CN201610783305.3, denominations of invention are as follows: a kind of straight calculation method based on looped network electric system,
The error calculated that the patent mainly solves the iterative method of existing power flow of ring net algorithm application is big, does not restrain, arithmetic speed
The problems such as slow；
And application No. is CN201810219284.1, denominations of invention are as follows: a kind of electric system based on straight algorithm is electromechanical
Transient emulation method, which overcome the defect of traditional simulation method, no iteration, calculating speed are fast, precision is high and error is small, truly
The variation characteristic of power grid is reflected, such as each route, load and the impedance of transformer change, power grid as the frequency of power grid changes
In the frequency of each generator also can dynamically change by respective rule；
For it is above-mentioned in the prior art, although calculated result can be made more based on straight algorithm in traditional calculation method
To be accurate, improving operational speed still, can not be directly real in the technical solution documented by above three patent document
The simulation calculation of existing electro-magnetic transient.
Summary of the invention
In view of this, in order to solve the above problems existing in the present technology, the purpose of the present invention is to provide one kind to be based on
The electromagnetism electromechanical transient simulation algorithm of straight algorithm can be realized electro-magnetic transient in the simulation calculation of one to reach.
The technical scheme adopted by the invention is as follows: a kind of electromagnetism electromechanical transient simulation algorithm based on straight algorithm, in electric power
Include generator model, circuit model, load module and transformer model in system emulation, calculates separately generator model, line
The electro-magnetic transient matrix of road model, load module and transformer model comprising following steps:
(1) corresponding electro-magnetic transient equation is solved according to different models respectively；
(2) the electro-magnetic transient equation in step (1) is organized into difierence equation；
(3) difierence equation is arranged to obtain electro-magnetic transient matrix；And electro-magnetic transient matrix is substituted into and is calculated based on straight
In the Electrical-Machanical Transient Simulation of Power System method of method, carry out electromagnetism electromechanics transient effect emulation, with reach can be realized electromagnetism,
Electromechanical transient is in the simulation calculation of one.
Further, it is the series circuit of resistance R and inductance L in the circuit model, enables the current vector of the circuit beVoltage along the direction both ends of current vector is respectively
If:
Solve electro-magnetic transient equation:
It is organized into difierence equation:
Wherein: Im0For the amplitude of previous frame electric current, Im0·ejθFor previous frame electric currentPredicted value, i.e.,Time interval of the Δ T between adjacent two frame；
It arranges:
Deformation obtains electro-magnetic transient matrix:
Further, it is the series circuit of resistance R and capacitor C in the circuit model, enables the current vector of the circuit beVoltage along the direction both ends of current vector is respectively
If:
Solve electro-magnetic transient equation:
Both sides differential:
It substitutes intoWithArrangement obtains:
It is organized into difierence equation:
Wherein: Im0For the amplitude of previous frame electric current, Im0·ejθFor previous frame electric currentPredicted value, i.e.,Time interval of the Δ T between adjacent two frame；
U1m0For the amplitude of previous frame left end voltage,For previous frame left end voltagePredicted value, i.e.,
U2m0For the amplitude of previous frame right end voltage,For previous frame right end voltagePredicted value,
I.e.
It deforms:
It arranges:
Finally:
Deformation obtains electro-magnetic transient matrix:
Further, the load circuit after connecting for resistance R with inductance L in the load module, enables the load circuit
Current vector isThe voltage at load circuit both ends is respectively
If:
Solve electro-magnetic transient equation:
Arrange is difierence equation:
Wherein: Im0For the amplitude of previous frame electric current, Im0·ejθFor previous frame electric currentPredicted value, i.e.,Time interval of the Δ T between adjacent two frame；
It arranges:
Finally:
Deformation obtains electro-magnetic transient matrix:
Further, the load circuit after connecting for resistance R with capacitor C in the load module, enables the load circuit
Current vector isThe voltage at load circuit both ends is respectively
If:
Solve electro-magnetic transient equation:
After the differential of both sides:
Arrange is difierence equation:
Wherein: Im0For the amplitude of previous frame electric current,For previous frame electric currentPredicted value, i.e.,
Um0For the amplitude of previous frame voltage,For previous frame voltagePredicted value, i.e.,
Time interval of the Δ T between adjacent two frame；
It arranges:
Finally:
Deformation obtains electro-magnetic transient matrix:
Further, string has resistance R and inductance L on the series circuit of generator G in the generator model, enables power generation
The electromotive force vector of machine G isThe current vector of the series circuit isThe voltage at series circuit both ends is respectively
If:
Solve electro-magnetic transient equation:
Arrange is difierence equation:
Wherein: Im0For the amplitude of previous frame electric current, Im0·ejθFor previous frame electric currentPredicted value, i.e.,Time interval of the Δ T between adjacent two frame；
It arranges:
Finally:
Deformation obtains electro-magnetic transient matrix:
Further, in the transformer model, if transformer is transformer, the mutual inductance of transformer is M, transformer
In the electrodynamic capacities of two coils be respectively L1 and L2, and electrodynamic capacity is the current vector and voltage vector point in the coil of L1
It is notWithAnd electrodynamic capacity be L2 coil in current vector and voltage vector be respectivelyWith
If:Wherein 0 <
K≤1 solves electro-magnetic transient equation:
Arrange is difierence equation:
Wherein: Im10For the amplitude of left end previous frame electric current,For left end previous frame electric currentPrediction
Value, i.e.,
Im20For the amplitude of right end previous frame electric current,For right end previous frame electric currentPredicted value, i.e.,
Time interval of the Δ T between adjacent two frame；
It is deformed into electro-magnetic transient matrix:
The invention has the benefit that
1. in the present invention by generator model present in power system network, circuit model, load module and
Transformer model carry out respectively solve electro-magnetic transient equation, to electro-magnetic transient equation carry out it is a series of solve, arrange with
And deformation, electro-magnetic transient matrix is finally obtained, in a kind of prior art " Electrical-Machanical Transient Simulation of Power System side based on straight algorithm
In the calculating process of method ", it has been all made of the matrix of the matrix of route, the matrix of load, the matrix of generator and transformer, with
The electro-magnetic transient matrix of different models corresponds to above-mentioned matrix and is replaced, with realize electromagnetism, electromechanical transient in one emulation, with
Guarantee that final calculation result is consistent with the actual conditions of power grid, and can more really reflect the variation of power grid.
Detailed description of the invention
Fig. 1 is that the first circuit model shows in the electromagnetism electromechanical transient simulation algorithm provided by the invention based on straight algorithm
It is intended to；
Fig. 2 is that second circuit model shows in the electromagnetism electromechanical transient simulation algorithm provided by the invention based on straight algorithm
It is intended to；
Fig. 3 is that the first load module shows in the electromagnetism electromechanical transient simulation algorithm provided by the invention based on straight algorithm
It is intended to；
Fig. 4 is that second load module shows in the electromagnetism electromechanical transient simulation algorithm provided by the invention based on straight algorithm
It is intended to；
Fig. 5 is the signal of generator model in the electromagnetism electromechanical transient simulation algorithm provided by the invention based on straight algorithm
Figure；
Fig. 6 is the signal of transformer model in the electromagnetism electromechanical transient simulation algorithm provided by the invention based on straight algorithm
Figure；
Fig. 7 is route Π type equivalent circuit schematic diagram provided by the invention；
Fig. 8 is route T-type equivalent circuit schematic diagram provided by the invention.
Specific embodiment
In order to make the object, technical scheme and advantages of the embodiment of the invention clearer, below in conjunction with the embodiment of the present invention
In attached drawing, technical scheme in the embodiment of the invention is clearly and completely described, it is clear that described embodiment is
A part of the embodiment of the present invention, instead of all the embodiments.The present invention being usually described and illustrated herein in the accompanying drawings is implemented
The component of example can be arranged and be designed with a variety of different configurations.
Therefore, the detailed description of the embodiment of the present invention provided in the accompanying drawings is not intended to limit below claimed
The scope of the present invention, but be merely representative of selected embodiment of the invention.Based on the embodiments of the present invention, this field is common
Technical staff's every other embodiment obtained without creative efforts belongs to the model that the present invention protects
It encloses.
It should be noted that in the absence of conflict, the feature in embodiment and embodiment in the present invention can phase
Mutually combination.
The present invention provides a kind of electromagnetism electromechanical transient simulation algorithm based on straight algorithm, the electricity solved by this method
Magnetic transient state matrix can replace being used in the prior art " a kind of Electrical-Machanical Transient Simulation of Power System method based on straight algorithm "
Matrix, in power system network, mainly by one or any of several in generator, route, load and transformer model
It is composed, now in conjunction with attached drawing, a variety of different models is solved respectively, specific as follows:
As shown in Figure 1, being the series circuit of resistance R and inductance L in the circuit model, the current vector of the circuit is enabled to beVoltage along the direction both ends of current vector is respectivelyCurrent vector isDirection beExtremelyThat is the electricity
The voltage at road middle left and right both ends is respectively
If:Wherein,The π of ω=2 f is angular speed, and f is frequency, and t is
Time,For initial phase.
Solve electro-magnetic transient equation:
It is organized into difierence equation:
Wherein: Im0For the amplitude of previous frame electric current, Im0·ejθFor previous frame electric currentPredicted value, i.e.,Time interval of the Δ T between adjacent two frame；
It arranges:
Deformation obtains electro-magnetic transient matrix:
Above-mentioned deformation method can refer to Patent No. ZL201410142938.7, denomination of invention are as follows: a word chain and attachment formula
The non-straight algorithm of looped network electric system of the more power supplys of three-phase symmetrical in the calculation method used.
As shown in Fig. 2, being the series circuit of resistance R and capacitor C in the circuit model, the current vector of the circuit is enabled to be
Voltage along the direction both ends of current vector is respectivelyCurrent vector isDirection beExtremelyThat is the circuit
The voltage at middle left and right both ends is respectively
If:
Solve electro-magnetic transient equation:
Differential process is carried out to both sides:
It substitutes intoWithArrangement obtains:
It is organized into difierence equation:
Wherein: Im0For the amplitude of previous frame electric current, Im0·ejθFor previous frame electric currentPredicted value, i.e.,Time interval of the Δ T between adjacent two frame；
U1m0For the amplitude of previous frame left end voltage,For previous frame left end voltagePredicted value, i.e.,
U2m0For the amplitude of previous frame right end voltage,For previous frame right end voltagePredicted value,
I.e.
It arranges:
It arranges:
Finally:
Deformation obtains electro-magnetic transient matrix:
Above-mentioned deformation method can refer to Patent No. ZL201410142938.7, denomination of invention are as follows: a word chain and attachment formula
The non-straight algorithm of looped network electric system of the more power supplys of three-phase symmetrical in the calculation method used.
As shown in figure 3, the load circuit after connecting for resistance R with inductance L in the load module, enables the load circuit
Current vector isThe voltage at load circuit both ends is respectivelyThat is the voltage at the circuit middle left and right both ends is respectively
If:
Solve electro-magnetic transient equation:
Arrange is difierence equation:
Wherein: Im0For the amplitude of previous frame electric current, Im0·ejθFor previous frame electric currentPredicted value, i.e.,Time interval of the Δ T between adjacent two frame；
It arranges:
Finally:
Deformation obtains electro-magnetic transient matrix:
Above-mentioned deformation method can refer to Patent No. ZL201410142938.7, denomination of invention are as follows: a word chain and attachment formula
The non-straight algorithm of looped network electric system of the more power supplys of three-phase symmetrical in the calculation method used.
As shown in figure 4, the load circuit after connecting for resistance R with capacitor C in the load module, enables the load circuit
Current vector isThe voltage at load circuit both ends is respectivelyThe as voltage difference at the circuit middle left and right both ends
ForWherein,
If:
Solve electro-magnetic transient equation:
After the differential of both sides:
Arrange is difierence equation:
Wherein: Im0For the amplitude of previous frame electric current,For previous frame electric currentPredicted value, i.e.,
Um0For the amplitude of previous frame voltage,For previous frame voltagePredicted value, i.e.,
Time interval of the Δ T between adjacent two frame；
It arranges:
Finally:
Deformation obtains electro-magnetic transient matrix:
Above-mentioned deformation method can refer to Patent No. ZL201410142938.7, denomination of invention are as follows: a word chain and attachment formula
The non-straight algorithm of looped network electric system of the more power supplys of three-phase symmetrical in the calculation method used.
As shown in figure 5, string has resistance R and inductance L on the series circuit of generator G in the generator model, hair is enabled
The electromotive force vector of motor G isThe current vector of the series circuit isThe voltage at series circuit both ends is respectivelyThe voltage at the as circuit middle left and right both ends is respectively
If:
Solve electro-magnetic transient equation:
Arrange is difierence equation:
Wherein: Im0For the amplitude of previous frame electric current, Im0·ejθFor previous frame electric currentPredicted value, i.e.,Time interval of the Δ T between adjacent two frame；
It arranges:
Finally:
Deformation obtains electro-magnetic transient matrix:
Above-mentioned deformation method can refer to Patent No. ZL201410142938.7, denomination of invention are as follows: a word chain and attachment formula
The non-straight algorithm of looped network electric system of the more power supplys of three-phase symmetrical in the calculation method used.
As shown in fig. 6, in the transformer model, if the mutual inductance of transformer is M, in transformer two coils from
Sense coefficient is respectively L1 and L2, and electrodynamic capacity is that current vector and voltage vector in the coil of L1 is respectivelyWithAnd
Electrodynamic capacity is that the current vector and voltage vector in the coil of L2 are respectivelyWith
If:Wherein 0 <
K≤1 solves electro-magnetic transient equation:
Arrange is difierence equation:
Wherein: Im10For the amplitude of left end previous frame electric current,For left end previous frame electric currentPrediction
Value, i.e.,
Im20For the amplitude of right end previous frame electric current,For right end previous frame electric currentPredicted value, i.e.,
Time interval of the Δ T between adjacent two frame；
It is deformed into electro-magnetic transient matrix:
Above-mentioned deformation method can refer to Patent No. ZL201410142938.7, denomination of invention are as follows: a word chain and attachment formula
The non-straight algorithm of looped network electric system of the more power supplys of three-phase symmetrical in the calculation method used.
It to sum up solves, arrange and deforms, for the electro-magnetic transient matrix of route in the first circuit model are as follows:
For the electro-magnetic transient matrix of route in second of circuit model are as follows:
For the electro-magnetic transient matrix loaded in the first load module are as follows:
For the electro-magnetic transient matrix loaded in second of load module are as follows:
For the electro-magnetic transient matrix of generator model are as follows:
For the electro-magnetic transient matrix of transformer model are as follows:
(1), application No. is CN201810219284.1, denominations of invention are as follows: a kind of power system machine based on straight algorithm
In electric transient emulation method, Electrical-Machanical Transient Simulation of Power System is carried out by step A- step I, is initialized in step
When parameters of electric power system, online system frequency initial value is set in step A4 as fW,0=50Hz, according to fW,0Determine each node
Reactance and susceptance, finally calculate the initial matrix of all nodes；
Step A4 specifically includes following procedure: setting in electric system, in original frequency fW,0When=50Hz, the resistance of load
For Ri,0It is X with reactancei,0；Every kilometer of resistance of route is ri,0, every kilometer of reactance be xi,0, every kilometer of conductance be gi,0, every kilometer of electricity
It receives as bi,0It is l with line lengthi；The conductance of transformer is Gti,0, susceptance Bti,0, resistance Rti,0, reactance Xti,0, former
Side the number of turns is ni,1It is n with secondary side the number of turnsi,2；The internal resistance of generator is ri',0It is x with reactancei',0；Then each node initial matrix is such as
Under:
Load initial matrix are as follows:
Route initial matrix are as follows:
Wherein, zi,0=ri,0+jxi,0, yi,0=gi,0+jbi,0,
Transformer initial matrix are as follows:
Generator initial matrix are as follows:
Similarly, in specifically to solve in the present embodiment, load initial matrix is replaced with electro-magnetic transient matrix in load module；
Generator initial matrix is replaced with electro-magnetic transient matrix in generator model；It replaces becoming with electro-magnetic transient matrix in transformer model
Depressor initial matrix, then the trend based on electric system directly calculate step A- step I provided by method and carry out continuing to calculate, with into
Row Electrical-Machanical Transient Simulation of Power System, result calculated just can be realized electromagnetism, electromechanical transient and work as Δ in the emulation of one
T is got over hour, and electro-magnetic transient effect is more obvious, and when Δ T is bigger, electro-magnetic transient effect is weaker.
Wherein, for route initial matrix:It can be equivalent at route Π type
Equivalent circuit or route T-type equivalent circuit；
Wherein, as shown in fig. 7, when the equivalent type equivalent circuit at Π,
Then
If
Then the electro-magnetic transient matrix of route distributed parameter model can be obtained by three following electro-magnetic transient matrix multiples:
Wherein:
As shown in figure 8, when the equivalent equivalent circuit at T-type,
Then
If
Then the electro-magnetic transient matrix of route distributed parameter model can be obtained by three following electro-magnetic transient matrix multiples:
Wherein:
It, will application No. is CN201810219284.1, denominations of invention according to above-mentioned are as follows: a kind of power train based on straight algorithm
Load initial matrix, route initial matrix, transformer initial matrix and the initial square of generator in electromechanical transient simulation method of uniting
Battle array carries out after replacing accordingly, referring to its disclosed in " a kind of Electrical-Machanical Transient Simulation of Power System method based on straight algorithm "
He continues operation at step, and electromagnetism electromechanical transient can be realized in the simulation calculation of one.
The present invention is not limited to above-mentioned optional embodiment, anyone can show that other are various under the inspiration of the present invention
The product of form, however, make any variation in its shape or structure, it is all to fall into the claims in the present invention confining spectrum
Technical solution, be within the scope of the present invention.
Claims (7)
1. a kind of electromagnetism electromechanical transient simulation algorithm based on straight algorithm includes generator model, line in electric system simulation
Road model, load module and transformer model, which is characterized in that calculate separately generator model, circuit model, load module and
The electro-magnetic transient matrix of transformer model comprising following steps:
(1) corresponding electro-magnetic transient equation is solved according to different models respectively；
(2) the electro-magnetic transient equation in step (1) is organized into difierence equation；
(3) difierence equation is arranged to obtain electro-magnetic transient matrix；And electro-magnetic transient matrix is substituted into based on straight algorithm
In Electrical-Machanical Transient Simulation of Power System method, the emulation of electro-magnetic transient effect is carried out.
2. the electromagnetism electromechanical transient simulation algorithm according to claim 1 based on straight algorithm, which is characterized in that the route
It is the series circuit of resistance R and inductance L in model, enables the current vector of the circuit beElectricity along the direction both ends of current vector
Pressure is respectively
If:
Solve electro-magnetic transient equation:
It is organized into difierence equation:
Wherein: Im0For the amplitude of previous frame electric current, Im0·ejθFor previous frame electric currentPredicted value, i.e.,
Time interval of the Δ T between adjacent two frame；
It arranges:
Deformation obtains electro-magnetic transient matrix:
3. the electromagnetism electromechanical transient simulation algorithm according to claim 1 based on straight algorithm, which is characterized in that the route
It is the series circuit of resistance R and capacitor C in model, enables the current vector of the circuit beElectricity along the direction both ends of current vector
Pressure is respectively
If:
Solve electro-magnetic transient equation:
Both sides differential:
It substitutes intoWithArrangement obtains:
It is organized into difierence equation:
Wherein: Im0For the amplitude of previous frame electric current, Im0·ejθFor previous frame electric currentPredicted value, i.e.,
U1m0For the amplitude of previous frame left end voltage,For previous frame left end voltagePredicted value, i.e.,
U2m0For the amplitude of previous frame right end voltage,For previous frame right end voltagePredicted value, i.e.,
Time interval of the Δ T between adjacent two frame；
It deforms:
It arranges:
Finally:
Deformation obtains electro-magnetic transient matrix:
4. the electromagnetism electromechanical transient simulation algorithm according to claim 1 based on straight algorithm, which is characterized in that the load
Load circuit after connecting for resistance R with inductance L in model, enables the current vector of the load circuit beLoad circuit both ends
Voltage be respectively
If:
Solve electro-magnetic transient equation:
Arrange is difierence equation:
Wherein: Im0For the amplitude of previous frame electric current, Im0·ejθFor previous frame electric currentPredicted value, i.e.,
Time interval of the Δ T between adjacent two frame；
It arranges:
Finally:
Deformation obtains electro-magnetic transient matrix:
5. the electromagnetism electromechanical transient simulation algorithm according to claim 1 based on straight algorithm, which is characterized in that the load
Load circuit after connecting for resistance R with capacitor C in model, enables the current vector of the load circuit beLoad circuit both ends
Voltage be respectively
If:
Solve electro-magnetic transient equation:
After the differential of both sides:
Arrange is difierence equation:
:
Wherein: Im0For the amplitude of previous frame electric current,For previous frame electric currentPredicted value, i.e.,
Um0For the amplitude of previous frame voltage,For previous frame voltagePredicted value, i.e.,
Time interval of the Δ T between adjacent two frame；
It arranges:
Finally:
Deformation obtains electro-magnetic transient matrix:
6. the electromagnetism electromechanical transient simulation algorithm according to claim 1 based on straight algorithm, which is characterized in that the power generation
String has resistance R and inductance L on the series circuit of generator G in machine model, enables the electromotive force vector of generator G beThe string
Connection circuit current vector beThe voltage at series circuit both ends is respectively
If:
Solve electro-magnetic transient equation:
Arrange is difierence equation:
Wherein: Im0For the amplitude of previous frame electric current, Im0·ejθFor previous frame electric currentPredicted value, i.e.,
Time interval of the Δ T between adjacent two frame；
It arranges:
Finally:
Deformation obtains electro-magnetic transient matrix:
7. the electromagnetism electromechanical transient simulation algorithm according to claim 1 based on straight algorithm, which is characterized in that the transformation
In device model, if the mutual inductance of transformer is M, the electrodynamic capacity of two coils is respectively L1 and L2 in transformer, and self-induction
Coefficient is that the current vector and voltage vector in the coil of L1 are respectivelyWithAnd electrodynamic capacity is the electricity in the coil of L2
Flow vector and voltage vector are respectivelyWith
If: Wherein 0 < k≤1 is asked
Solve electro-magnetic transient equation:
Arrange is difierence equation:
Wherein: Im10For the amplitude of left end previous frame electric current,For left end previous frame electric currentPredicted value, i.e.,
Im20For the amplitude of right end previous frame electric current,For right end previous frame electric currentPredicted value, i.e.,
Time interval of the Δ T between adjacent two frame；
It is deformed into electro-magnetic transient matrix:
Priority Applications (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811285344.6A CN109063408B (en) | 2018-10-31 | 2018-10-31 | Electromagnetic electromechanical transient simulation method based on direct algorithm |
| PCT/CN2019/111295 WO2020088245A1 (en) | 2018-10-31 | 2019-10-15 | Direct algorithm-based electromagnetic and electromechanical transient simulation algorithm |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811285344.6A CN109063408B (en) | 2018-10-31 | 2018-10-31 | Electromagnetic electromechanical transient simulation method based on direct algorithm |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109063408A true CN109063408A (en) | 2018-12-21 |
| CN109063408B CN109063408B (en) | 2022-04-01 |
Family
ID=64789096
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811285344.6A Active CN109063408B (en) | 2018-10-31 | 2018-10-31 | Electromagnetic electromechanical transient simulation method based on direct algorithm |
Country Status (2)
| Country | Link |
|---|---|
| CN (1) | CN109063408B (en) |
| WO (1) | WO2020088245A1 (en) |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110061506A (en) * | 2019-04-23 | 2019-07-26 | 邓宏伟 | A kind of circuit-breaker switching on-off simulation algorithm based on straight algorithm |
| CN110442895A (en) * | 2019-05-24 | 2019-11-12 | 华北电力大学 | A kind of single-phase double winding high frequency transformer electro-magnetic transient equivalent modeling method considering capacity effect |
| WO2020088245A1 (en) * | 2018-10-31 | 2020-05-07 | 邓宏伟 | Direct algorithm-based electromagnetic and electromechanical transient simulation algorithm |
| WO2023240844A1 (en) * | 2022-06-15 | 2023-12-21 | 南方电网科学研究院有限责任公司 | Electromagnetic transient modeling method and system for high-efficiency synchronous machine, and device |
Families Citing this family (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112084624A (en) * | 2020-07-31 | 2020-12-15 | 清华大学 | Electromagnetic-electromechanical hybrid simulation electromagnetic transient side interface power calculation method and device |
| CN112234598B (en) * | 2020-08-28 | 2024-04-23 | 国网天津市电力公司电力科学研究院 | Electromagnetic transient simulation initialization method |
| CN113076675B (en) * | 2021-04-12 | 2022-11-08 | 中国电子科技集团公司第三十三研究所 | Electromagnetic environment effect simulation design method for air cushion landing boat |
| CN113536589B (en) * | 2021-07-30 | 2023-08-04 | 广东电网有限责任公司广州供电局 | Electromagnetic transient modeling method and system for alternating current power grid |
| WO2024082163A1 (en) * | 2022-10-19 | 2024-04-25 | 云南电网有限责任公司电力科学研究院 | Method for electromagnetic transient simulation of transformer, computer device, and storage medium |
| CN116722563A (en) * | 2023-05-30 | 2023-09-08 | 杭州盛星能源技术有限公司 | Electromagnetic transient simulation frequency domain expansion method and device based on dynamic phasors |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20110257943A1 (en) * | 2010-04-16 | 2011-10-20 | Texas Instruments Incorporated | Node-based transient acceleration method for simulating circuits with latency |
| CN103077268A (en) * | 2012-12-27 | 2013-05-01 | 天津大学 | State space automatic modeling method orienting electromagnetic transient simulation of power system |
| CN103956741A (en) * | 2014-04-10 | 2014-07-30 | 邓宏伟 | Horizontal chain type three-phase symmetry multi-source non-looped network electric power system straight algorithm and branched chain type three-phase symmetry multi-source non-looped network electric power system straight algorithm |
| CN105224754A (en) * | 2015-10-14 | 2016-01-06 | 清华大学 | A kind of simulation of power electronic method based on Interpolation compensation current switch model |
| CN106372339A (en) * | 2016-09-05 | 2017-02-01 | 清华大学 | Multi-rate simulation method and device for power-electronizing power system |
| CN108365629A (en) * | 2018-03-16 | 2018-08-03 | 邓宏伟 | A kind of Electrical-Machanical Transient Simulation of Power System method based on straight algorithm |
Family Cites Families (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7089163B2 (en) * | 2002-08-27 | 2006-08-08 | Synopsys, Inc. | Smooth operators in optimization of structures |
| CN106451456B (en) * | 2016-08-30 | 2019-01-29 | 邓宏伟 | A kind of straight calculation method based on looped network electric system |
| CN108536925B (en) * | 2018-03-21 | 2021-03-16 | 武汉大学 | Isolated dynamic whole-process real-time hybrid simulation interface system |
| CN109063408B (en) * | 2018-10-31 | 2022-04-01 | 邓宏伟 | Electromagnetic electromechanical transient simulation method based on direct algorithm |
-
2018
- 2018-10-31 CN CN201811285344.6A patent/CN109063408B/en active Active
-
2019
- 2019-10-15 WO PCT/CN2019/111295 patent/WO2020088245A1/en active Application Filing
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20110257943A1 (en) * | 2010-04-16 | 2011-10-20 | Texas Instruments Incorporated | Node-based transient acceleration method for simulating circuits with latency |
| CN103077268A (en) * | 2012-12-27 | 2013-05-01 | 天津大学 | State space automatic modeling method orienting electromagnetic transient simulation of power system |
| CN103956741A (en) * | 2014-04-10 | 2014-07-30 | 邓宏伟 | Horizontal chain type three-phase symmetry multi-source non-looped network electric power system straight algorithm and branched chain type three-phase symmetry multi-source non-looped network electric power system straight algorithm |
| CN105224754A (en) * | 2015-10-14 | 2016-01-06 | 清华大学 | A kind of simulation of power electronic method based on Interpolation compensation current switch model |
| CN106372339A (en) * | 2016-09-05 | 2017-02-01 | 清华大学 | Multi-rate simulation method and device for power-electronizing power system |
| CN108365629A (en) * | 2018-03-16 | 2018-08-03 | 邓宏伟 | A kind of Electrical-Machanical Transient Simulation of Power System method based on straight algorithm |
Non-Patent Citations (3)
| Title |
|---|
| MAURO CALABRIA等: "《Modeling power inverter interactions in a low voltage grid》", 《2014 IEEE 15TH WORKSHOP ON CONTROL AND MODELING FOR POWER ELECTRONICS(COMPEL)》 * |
| TIAN FANG等: "Realization of Electromechanical Transient and", 《2005 IEEE/PES TRANSMISSION AND DISTRIBUTION》 * |
| 李伟等: "电磁–机电暂态混合仿真精度改善技术研究", 《广东电力》 * |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2020088245A1 (en) * | 2018-10-31 | 2020-05-07 | 邓宏伟 | Direct algorithm-based electromagnetic and electromechanical transient simulation algorithm |
| CN110061506A (en) * | 2019-04-23 | 2019-07-26 | 邓宏伟 | A kind of circuit-breaker switching on-off simulation algorithm based on straight algorithm |
| CN110061506B (en) * | 2019-04-23 | 2022-11-11 | 邓宏伟 | Circuit breaker opening and closing simulation method based on direct algorithm |
| CN110442895A (en) * | 2019-05-24 | 2019-11-12 | 华北电力大学 | A kind of single-phase double winding high frequency transformer electro-magnetic transient equivalent modeling method considering capacity effect |
| CN110442895B (en) * | 2019-05-24 | 2021-05-14 | 华北电力大学 | High-frequency transformer electromagnetic transient equivalent modeling method considering capacitance effect |
| WO2023240844A1 (en) * | 2022-06-15 | 2023-12-21 | 南方电网科学研究院有限责任公司 | Electromagnetic transient modeling method and system for high-efficiency synchronous machine, and device |
Also Published As
| Publication number | Publication date |
|---|---|
| WO2020088245A1 (en) | 2020-05-07 |
| CN109063408B (en) | 2022-04-01 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109063408A (en) | A kind of electromagnetism electromechanical transient simulation algorithm based on straight algorithm | |
| CN108072847B (en) | A kind of estimation method of dynamic lithium battery identification of Model Parameters and remaining capacity | |
| CN106021682B (en) | Sub-synchronous oscillation analysis method and device based on impedance network model | |
| CN106451456B (en) | A kind of straight calculation method based on looped network electric system | |
| CN110221237A (en) | A kind of capacitance type potential transformer harmonic voltage measurement method and apparatus | |
| CN106896267A (en) | A kind of Distribution Network Harmonics resonance improves modal analysis method | |
| CN103984822A (en) | Three phase flow implementation method based on advanced digital power system simulator | |
| CN105529711A (en) | BPA data based alternating current system harmonic impedance scanning method | |
| CN105160058B (en) | A kind of network boundary equivalence method calculated for electric network electromagnet transient | |
| CN105243197A (en) | Single-core power cable resonance analysis method based on finite element simulation and equivalent circuit | |
| CN108054757A (en) | A kind of embedded idle and voltage N-1 Close loop security check methods | |
| CN104156542A (en) | Implicit-projection-based method for simulating stability of active power distribution system | |
| CN104979840B (en) | A kind of active distribution network three-phase idle work optimization method | |
| Bush | Structural analysis by electric circuit analogies | |
| CN106056479A (en) | Three-phase modeling method for distribution network | |
| CN107271941A (en) | A kind of zero-sequence mutual inductance impedance computation method on multiple-loop line transmission line of electricity | |
| Parashar | Continuum modeling of power networks | |
| CN106294966A (en) | A kind of list telegram in reply cable core transient state temperature rise acquisition methods being independent of skin temperature | |
| CN108763790A (en) | A kind of power system electromagnetic transient simulation method based on extension critical damping adjusting method | |
| CN110808593B (en) | Method and device for analyzing damping of electromagnetic oscillation caused by multiple controllers in weak-connection power grid | |
| CN104361260A (en) | Electric power system node dynamic frequency response characteristic analysis method based on curved surface fitting | |
| CN104331620A (en) | Method for analyzing subsynchronous oscillation of near-region dense thermal power multi-direct current transmission system | |
| Dimitrovski et al. | Analyzing subsynchronous resonance using component connection method | |
| Ourari et al. | Coherency approach for dynamic equivalents of large power systems | |
| CA2661753A1 (en) | Method of loadflow computation for electrical power system |
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 |