CN113139264A - Electromechanical-electromagnetic hybrid simulation stability detection method and system - Google Patents
Electromechanical-electromagnetic hybrid simulation stability detection method and system Download PDFInfo
- Publication number
- CN113139264A CN113139264A CN202010052403.6A CN202010052403A CN113139264A CN 113139264 A CN113139264 A CN 113139264A CN 202010052403 A CN202010052403 A CN 202010052403A CN 113139264 A CN113139264 A CN 113139264A
- Authority
- CN
- China
- Prior art keywords
- electromechanical
- interface
- matrix
- simulation
- hybrid simulation
- 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
- 238000004088 simulation Methods 0.000 title claims abstract description 133
- 238000001514 detection method Methods 0.000 title claims description 7
- 239000011159 matrix material Substances 0.000 claims abstract description 168
- 230000005540 biological transmission Effects 0.000 claims abstract description 52
- 238000000034 method Methods 0.000 claims abstract description 28
- 238000004458 analytical method Methods 0.000 claims abstract description 15
- 238000004364 calculation method Methods 0.000 claims abstract description 11
- 230000001052 transient effect Effects 0.000 claims description 23
- 238000006243 chemical reaction Methods 0.000 claims description 16
- 238000012545 processing Methods 0.000 claims description 16
- 230000001131 transforming effect Effects 0.000 claims description 10
- 238000000605 extraction Methods 0.000 claims description 7
- 238000010586 diagram Methods 0.000 description 8
- 239000003795 chemical substances by application Substances 0.000 description 7
- 238000004590 computer program Methods 0.000 description 7
- 230000006870 function Effects 0.000 description 6
- 230000008569 process Effects 0.000 description 5
- 230000006378 damage Effects 0.000 description 4
- 238000004422 calculation algorithm Methods 0.000 description 3
- 230000001939 inductive effect Effects 0.000 description 2
- 239000000243 solution Substances 0.000 description 2
- 238000013459 approach Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000005538 encapsulation Methods 0.000 description 1
- 238000002347 injection Methods 0.000 description 1
- 239000007924 injection Substances 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
Images
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention provides a method for detecting the stability of electromechanical-electromagnetic hybrid simulation, which comprises the following steps: extracting an interface matrix of electromechanical-electromagnetic hybrid simulation, and solving a network interface matrix; judging whether the interface on the electromechanical side has the problem of interface non-convergence or not by analyzing positive definite characteristics of the network interface matrix on different frequencies to obtain the passivity characteristics of the interface matrix on different frequencies and the range of the envelope curve of the difference between the state variable and the transmission current in the electromechanical side simulation; judging the stability of the electromechanical-electromagnetic hybrid simulation based on the passivity characteristic and the interface convergence condition of the electromechanical side; the stability analysis of the electromechanical-electromagnetic hybrid simulation case is carried out, the instability reason is determined, support can be provided for example adjustment, and the stability and reliability of hybrid simulation calculation can be improved.
Description
Technical Field
The invention belongs to a simulation method of an electric power system, and relates to a method and a system for detecting the stability of electromechanical-electromagnetic hybrid simulation
Background
In the electromechanical-electromagnetic transient hybrid simulation analysis of the alternating current and direct current power grid, the electromechanical side simulation power grid needs to send an equivalent network of the local side power grid at a hybrid simulation interface to an electromagnetic transient program and participate in calculation. The equivalent prototype network is physically passive, but when the port equivalent network is solved, the matrix representing the prototype network is developed to the port based on fundamental frequency in the process of contracting the port, and the equivalent network obtained in some non-fundamental frequency bands may present active characteristics, namely passive destruction. The passive damage is not allowed in the electromagnetic transient simulation, and rich harmonic waves in the electromagnetic transient simulation can be unreasonably amplified even bring the computation instability of the electromechanical-electromagnetic transient hybrid simulation interface through the equivalent circuit of the passive damage. Before the hybrid simulation starts to calculate, stability of hybrid simulation calculation can be judged in advance by detecting stability of the interface equivalent circuit, hybrid simulation interface schemes with calculation instability risks are adjusted in advance, and stability and reliability of hybrid simulation calculation can be improved.
The traditional modeling simulation tool provides a closed whole set of simulation solution based on the traditional operating system. If the user has the third-party function, the user needs to link at the compiling level through a static linking tool. The external library mode causes that a user terminal needs to install a compiling environment required by the simulation tool, and simultaneously, the function encapsulation of the user terminal also needs to adopt a compiling and linking means completely consistent with the simulation tool. In addition, the disadvantage of this library-packaged and linked emulation approach is that it is difficult for user functionality to support the use of multiple instances. These limitations greatly restrict the development of user-side models and increase the complexity of modeling simulation.
Disclosure of Invention
Aiming at the bottleneck problem of the existing analysis method of the logical association relationship of the anti-misoperation lockout, the invention provides a method and a system for detecting the stability of electromechanical and electromagnetic hybrid simulation,
an electromechanical-electromagnetic hybrid simulation stability detection method comprises the following steps:
extracting an interface matrix of electromechanical-electromagnetic hybrid simulation, and solving a network interface matrix;
judging whether the interface on the electromechanical side has the problem of interface non-convergence or not by analyzing positive definite characteristics of the network interface matrix on different frequencies to obtain the passivity characteristics of the interface matrix on different frequencies and the range of the envelope curve of the difference between the state variable and the transmission current in the electromechanical side simulation;
and judging the stability of the electromechanical-electromagnetic hybrid simulation based on the passivity characteristic and the interface convergence condition of the electromechanical side.
Preferably, the extracting of the interface matrix of the electromechanical-electromagnetic hybrid simulation includes:
converting the equivalent matrix of the electromechanical network into a real part matrix and an imaginary part matrix;
transforming the real part matrix into a single-phase real part matrix of an ABC phase, and simultaneously transforming the imaginary part matrix into a single-phase imaginary part matrix of the ABC phase;
for the real part matrix YABCInversion processing to obtain ARInterface matrix, inverting imaginary matrix to obtain ALAn interface matrix.
Preferably, the obtaining the network interface matrix includes:
based on the ALCarrying out frequency w conversion on the interface matrix to select a proper frequency point;
and obtaining the network interface matrix under different frequencies through calculation.
Preferably, said is based on said ALThe interface matrix performs frequency w conversion to select a suitable frequency point, which includes:
the frequency is increased from 0.1Hz close to direct current to the range of 10KHz, and m points are taken as the points of frequency scanning according to an equidistant logarithm form, wherein m is a positive integer.
Preferably, the network interface matrix at different frequencies is calculated according to the following formula:
AG(w)=w*AL
in the formula, ALIs an interface matrix; w is the frequency; a. theGIs a network interface matrix.
Preferably, the obtaining of the passive characteristics of the interface matrix at different frequencies by analyzing the positive definite characteristics of the network interface matrix at different frequencies includes:
obtaining the current injected into the interface in the interface matrix simulation as a state variable;
acquiring current transmitted from the electromechanical side simulation;
and solving the difference between the state variable and the transmission current in the electromechanical side simulation, and judging the passivity characteristic.
Preferably, the determining whether the interface on the electromechanical side has the interface non-convergence problem in a range where the envelope of the difference between the state variable and the transmission current in the electromechanical side simulation is located includes:
if the envelope curve of the difference between the state variable and the transmission current in the electromechanical side simulation is between the electromechanical side transmission variable and the interface element state variable and is smaller than the electromechanical side transmission variable, the electromechanical transient interface in the electromechanical-electromagnetic hybrid simulation has the problem of non-convergence stability;
and if the envelope curve of the difference between the state variable and the transmission current in the electromechanical side simulation is between the electromechanical side transmission variable and the interface element state variable and is greater than the electromechanical side transmission variable, the electromechanical transient interface in the electromechanical-electromagnetic hybrid simulation is converged.
Preferably, the determining the stability of the electromechanical-electromagnetic hybrid simulation based on the passivity characteristic and the interface convergence condition on the electromechanical side includes:
if the different frequencies of the interface matrix do not have passivity characteristics and the interface on the electromechanical side is converged, the electromechanical-electromagnetic hybrid simulation is stable;
and if the different frequencies of the interface matrix have passivity characteristics or the interface on the electromechanical side does not converge, the electromechanical-electromagnetic hybrid simulation is unstable.
An electromechanical-electromagnetic hybrid simulation stability detection system, comprising:
an extraction module: the interface matrix is used for extracting the electromechanical-electromagnetic hybrid simulation and solving a network interface matrix;
an analysis module: the method is used for judging whether the interface on the electromechanical side has the problem of interface non-convergence or not by analyzing positive definite characteristics of the network interface matrix on different frequencies and obtaining the passivity characteristics of the interface matrix on different frequencies and the range of the envelope curve of the difference between the state variable and the transmission current in the electromechanical side simulation;
a judging module: and the stability of the electromechanical-electromagnetic hybrid simulation is judged based on the passivity characteristic and the interface convergence condition of the electromechanical side.
Preferably, the extraction module includes: a processing submodule and a solving submodule;
the processing submodule is used for:
converting the equivalent matrix of the electromechanical network into a real part matrix and an imaginary part matrix;
transforming the real part matrix into a single-phase real part matrix of an ABC phase, and simultaneously transforming the imaginary part matrix into a single-phase imaginary part matrix of the ABC phase;
for the real part matrix YABCInversion processing to obtain ARInterface matrix, inverting imaginary matrix to obtain ALAn interface matrix;
the obtaining submodule is used for:
based on the ALCarrying out frequency w conversion on the interface matrix to select a proper frequency point;
obtaining network interface matrixes under different frequencies through calculation;
wherein, A isLThe interface matrix performs frequency w conversion to select a suitable frequency point, which includes:
the frequency is increased to the range of 10KHz from 0.1Hz which is close to direct current, and m points are taken as frequency scanning points according to an equidistant logarithm form, wherein m is a positive integer;
the network interface matrix under different frequencies is calculated according to the following formula:
AG(w)=w*AL
in the formula, ALIs an interface matrix; w is the frequency; a. theGIs a network interface matrix.
Preferably, the analysis module includes: acquiring a submodule and a judging submodule;
the acquisition submodule is used for:
obtaining the current injected into the interface in the interface matrix simulation as a state variable;
acquiring current transmitted from the electromechanical side simulation;
solving the difference between the state variable and the transmission current in the electromechanical side simulation, and judging the passivity characteristic;
the judgment submodule is used for: judging whether an electromechanical transient interface in electromechanical-electromagnetic hybrid simulation converges;
if the envelope curve of the difference between the state variable and the transmission current in the electromechanical side simulation is between the electromechanical side transmission variable and the interface element state variable and is smaller than the electromechanical side transmission variable, the electromechanical transient interface in the electromechanical-electromagnetic hybrid simulation has the problem of non-convergence stability;
and if the envelope curve of the difference between the state variable and the transmission current in the electromechanical side simulation is between the electromechanical side transmission variable and the interface element state variable and is greater than the electromechanical side transmission variable, the electromechanical transient interface in the electromechanical-electromagnetic hybrid simulation is converged.
Has the advantages that:
1. the invention provides a method for detecting the stability of electromechanical-electromagnetic hybrid simulation, which comprises the following steps: extracting an interface matrix of electromechanical-electromagnetic hybrid simulation, and solving a network interface matrix; judging whether the interface on the electromechanical side has the problem of interface non-convergence or not by analyzing positive definite characteristics of the network interface matrix on different frequencies to obtain the passivity characteristics of the interface matrix on different frequencies and the range of the envelope curve of the difference between the state variable and the transmission current in the electromechanical side simulation; judging the stability of the electromechanical-electromagnetic hybrid simulation based on the passivity characteristic and the interface convergence condition of the electromechanical side; the stability analysis of the electromechanical-electromagnetic hybrid simulation case is carried out, the instability reason is determined, support can be provided for example adjustment, and the stability and reliability of hybrid simulation calculation can be improved.
Drawings
FIG. 1 is a flow chart of a method for detecting stability of an electromechanical-electromagnetic hybrid simulation according to the present invention;
FIG. 2 is a schematic flow chart of an embodiment of the present invention;
Detailed Description
Example 1
The invention provides a method for detecting the stability of electromechanical-electromagnetic hybrid simulation, which comprises the following steps:
s1, extracting an interface matrix of the electromechanical-electromagnetic hybrid simulation, and solving a network interface matrix;
s2, judging whether the interface on the electromechanical side has the problem of interface non-convergence by analyzing the positive definite characteristic of the network interface matrix on different frequencies and obtaining the passive characteristic of the interface matrix on different frequencies and the range where the envelope curve of the difference between the state variable and the transmission current in the electromechanical side simulation is located;
and S3, judging the stability of the electromechanical-electromagnetic hybrid simulation based on the passivity characteristic and the interface convergence condition of the electromechanical side.
S1, extracting an interface matrix of the electromechanical-electromagnetic hybrid simulation, and solving a network interface matrix;
1. interface matrix A for extracting electromechanical-electromagnetic hybrid simulationLAnd ARThe interface matrix is not the traditional transmission matrix for transmitting the electromagnetic side from the electromechanical side program, but the A for reserving the characteristics of the original interface matrixLAn inductor and ARInterface matrix conversion of the resistive elements.
Inductance interface matrix ALAnd a resistor interface matrix ARThe acquisition method is the key point of the patent, and the equivalent matrix G of the electromechanical network needs to be acquired012Conversion to real matrix Y012And an imaginary matrix B012Two parts, then to the real part matrix Y012Transforming into single-phase real part matrix Y of ABC phaseABCWhile simultaneously aligning the imaginary matrix B012Single-phase imaginary matrix B transformed into ABC phaseABCThen, the real part matrix Y is processed by inversionABCInversion processing to obtain ARInverting the imaginary matrix to obtain AL. As shown in the following formula:
wherein: al is an inductance interface matrix, AR is a capacitance interface demonstration, Tr012 is a conversion matrix from 012 to ABC,
im () is the imaginary part of the matrix, Re is the real part of the matrix, and G012 is the interface on the electromechanical side.
2. To ALThe inductive element interface matrix is subjected to frequency conversion to be converted into an admittance interface matrix A on each frequency point wG(w)。
ALThe frequency w conversion of the inductance element interface matrix needs to reasonably select a proper frequency point, w is increased from 0.1Hz close to direct current to a range of 10KHz, and 250 points are taken as frequency scanning points in the frequency range according to an equidistant logarithm mode. Then, the user can use the device to perform the operation,AG(w)=w*ALthe network interface matrix under different frequencies can be obtained by the formula.
S2, judging whether the interface on the electromechanical side has the problem of interface non-convergence by analyzing the positive definite characteristic of the network interface matrix on different frequencies and obtaining the passive characteristic of the interface matrix on different frequencies and the range where the envelope curve of the difference between the state variable and the transmission current in the electromechanical side simulation is located;
3. analysis AR+AGAnd (w) and the eigenvalues of the sum matrix, and finding whether negative characteristic roots exist.
Analysis A by Standard characteristic root analysis methodR+AG(w) recording the position of the feature root. In the process of feature root analysis, the real number feature root of the feature root is mainly concerned.
4. And confirming the sensitive frequency w of the electromechanical and electromagnetic hybrid simulation which can be subjected to the passive destruction according to the analysis condition.
5. Performing electromechanical-electromagnetic hybrid simulation, extracting state variables of the interface element through an algorithm, analyzing the transient process and trend of the interface element, and determining whether the system has active problems or not
In the process of electromechanical-electromagnetic transient hybrid simulation, the state variable of the interface element is recorded at each step. The state variable is the current injected into the interface in the previous simulation of the interface matrix. The difference between the state variable and the current transmitted from the electromechanical side simulation is determined.
6. And comparing the difference between the state variable of the interface element and the transmission variable of the electromechanical side, and judging whether the interface of the electromechanical side has the problem of interface non-convergence.
The envelope curve of the difference between the state variable of the interface element and the state variable of the electromechanical side is smaller than the state variable of the electromechanical side between the state variable of the electromechanical side and the state variable of the interface element, and the stability problem of non-convergence of the electromechanical transient interface in the electromechanical-electromagnetic hybrid simulation exists; if the envelope of the difference between the state variable of the interface element and the state variable of the electromechanical side is larger than the state variable of the electromechanical side and is between the state variable of the electromechanical side and the state variable of the interface element, the instability phenomenon of the electromechanical-electromagnetic hybrid simulation is not S3 caused by the problem that the electromechanical transient interface is not converged, and the stability of the electromechanical-electromagnetic hybrid simulation is judged based on the passivity characteristic and the interface convergence condition of the electromechanical side
7. The stability and feasibility of this hybrid simulation were evaluated systematically in accordance with the conclusions of 5 and 6.
Example 2
The invention relates to a method for detecting simulation stability in power grid simulation, which takes common simulation instability of a system containing three direct current machines as an implementation example.
Step 1, opening and establishing a hybrid simulation example, wherein simulation instability phenomenon occurs in normal simulation.
Step 2, obtaining an inductance interface matrix A of the hybrid simulation interface through an interface matrix extraction algorithm integrated into the simulation functionLAnd a resistor interface matrix AR。
Step 3, for ALThe inductive element interface matrix is subjected to frequency conversion, the frequency range is selected to be 0.1 to 2500Hz, 500 w are selected at equal intervals of logarithm, and the frequency is converted into an admittance interface matrix A on each frequency point wG(w)。
Step 4, analyzing 500AR+AGAnd (w) and the eigenvalues of the sum matrix, and finding whether negative characteristic roots exist.
And 5, performing electromechanical-electromagnetic hybrid simulation, loading an algorithm, analyzing transient change of an extraction variable of an interface injection current source, and discovering an active problem.
And 6, comparing the difference between the state variable of the interface element and the electromechanical side transmission variable, and finding out that the envelope curve of the difference between the state variable of the interface element and the electromechanical side transmission variable is between the electromechanical side transmission variable and the state variable of the interface element and is larger than the electromechanical side transmission variable.
Example 3
The invention also provides an electromechanical-electromagnetic hybrid simulation stability detection system based on the same concept, which comprises:
an extraction module: the interface matrix is used for extracting the electromechanical-electromagnetic hybrid simulation and solving a network interface matrix;
an analysis module: the method is used for judging whether the interface on the electromechanical side has the problem of interface non-convergence or not by analyzing positive definite characteristics of the network interface matrix on different frequencies and obtaining the passivity characteristics of the interface matrix on different frequencies and the range of the envelope curve of the difference between the state variable and the transmission current in the electromechanical side simulation;
a judging module: and the stability of the electromechanical-electromagnetic hybrid simulation is judged based on the passivity characteristic and the interface convergence condition of the electromechanical side.
The extraction module comprises: a processing submodule and a solving submodule;
the processing submodule is used for:
converting the equivalent matrix of the electromechanical network into a real part matrix and an imaginary part matrix;
transforming the real part matrix into a single-phase real part matrix of an ABC phase, and simultaneously transforming the imaginary part matrix into a single-phase imaginary part matrix of the ABC phase;
for the real part matrix YABCInversion processing to obtain ARInterface matrix, inverting imaginary matrix to obtain ALAn interface matrix;
the obtaining submodule is used for:
based on the ALCarrying out frequency w conversion on the interface matrix to select a proper frequency point;
obtaining network interface matrixes under different frequencies through calculation;
wherein, A isLThe interface matrix performs frequency w conversion to select a suitable frequency point, which includes:
the frequency is increased to the range of 10KHz from 0.1Hz which is close to direct current, and m points are taken as frequency scanning points according to an equidistant logarithm form, wherein m is a positive integer;
the network interface matrix under different frequencies is calculated according to the following formula:
AG(w)=w*AL
in the formula, ALIs an interface matrix; w is the frequency; a. theGIs a network interface matrix.
The analysis module comprises: acquiring a submodule and a judging submodule;
the acquisition submodule is used for:
obtaining the current injected into the interface in the interface matrix simulation as a state variable;
acquiring current transmitted from the electromechanical side simulation;
solving the difference between the state variable and the transmission current in the electromechanical side simulation, and judging the passivity characteristic;
the judgment submodule is used for: judging whether an electromechanical transient interface in electromechanical-electromagnetic hybrid simulation converges;
if the envelope curve of the difference between the state variable and the transmission current in the electromechanical side simulation is between the electromechanical side transmission variable and the interface element state variable and is smaller than the electromechanical side transmission variable, the electromechanical transient interface in the electromechanical-electromagnetic hybrid simulation has the problem of non-convergence stability;
and if the envelope curve of the difference between the state variable and the transmission current in the electromechanical side simulation is between the electromechanical side transmission variable and the interface element state variable and is greater than the electromechanical side transmission variable, the electromechanical transient interface in the electromechanical-electromagnetic hybrid simulation is converged.
As will be appreciated by one skilled in the art, embodiments of the present application may be provided as a method, system, or computer program product. Accordingly, the present application may take the form of an entirely hardware embodiment, an entirely software embodiment or an embodiment combining software and hardware aspects. Furthermore, the present application may take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, and the like) having computer-usable program code embodied therein.
The present application is described with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of the application. It will be understood that each flow and/or block of the flow diagrams and/or block diagrams, and combinations of flows and/or blocks in the flow diagrams and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to a processor of a general purpose computer, special purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory produce an article of manufacture including instruction means which implement the function specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
Finally, it should be noted that: the above embodiments are only for illustrating the technical solutions of the present invention and not for limiting the same, and although the present invention is described in detail with reference to the above embodiments, those of ordinary skill in the art should understand that: modifications and equivalents may be made to the embodiments of the invention without departing from the spirit and scope of the invention, which is to be covered by the claims.
Claims (11)
1. A method for detecting the stability of electromechanical-electromagnetic hybrid simulation is characterized by comprising the following steps:
extracting an interface matrix of electromechanical-electromagnetic hybrid simulation, and solving a network interface matrix;
judging whether the interface on the electromechanical side has the problem of interface non-convergence or not by analyzing positive definite characteristics of the network interface matrix on different frequencies to obtain the passivity characteristics of the interface matrix on different frequencies and the range of the envelope curve of the difference between the state variable and the transmission current in the electromechanical side simulation;
and judging the stability of the electromechanical-electromagnetic hybrid simulation based on the passivity characteristic and the interface convergence condition of the electromechanical side.
2. The method for detecting the stability of the electromechanical-electromagnetic hybrid simulation according to claim 1, wherein the extracting the interface matrix of the electromechanical-electromagnetic hybrid simulation comprises:
converting the equivalent matrix of the electromechanical network into a real part matrix and an imaginary part matrix;
transforming the real part matrix into a single-phase real part matrix of an ABC phase, and simultaneously transforming the imaginary part matrix into a single-phase imaginary part matrix of the ABC phase;
for the real part matrix YABCInversion processing to obtain ARInterface matrix, inverting imaginary matrix to obtain ALAn interface matrix.
3. The method for detecting stability of electromechanical-electromagnetic hybrid simulation according to claim 2, wherein the obtaining the network interface matrix comprises:
based on the ALCarrying out frequency w conversion on the interface matrix to select a proper frequency point;
and obtaining the network interface matrix under different frequencies through calculation.
4. The method according to claim 3, wherein the stability detection method based on the electromechanical-electromagnetic hybrid simulation is based on the ALThe interface matrix performs frequency w conversion to select a suitable frequency point, which includes:
the frequency is increased from 0.1Hz close to direct current to the range of 10KHz, and m points are taken as the points of frequency scanning according to an equidistant logarithm form, wherein m is a positive integer.
5. The method according to claim 3, wherein the network interface matrix at different frequencies is calculated according to the following formula:
AG(w)=w*AL
in the formula, ALIs an interface matrix; w is the frequency; a. theGIs a network interface matrix.
6. The method for detecting stability of electromechanical-electromagnetic hybrid simulation according to claim 3, wherein the obtaining of the passivity characteristics of the interface matrix at different frequencies by analyzing the positive definite characteristics of the network interface matrix at different frequencies comprises:
obtaining the current injected into the interface in the interface matrix simulation as a state variable;
acquiring current transmitted from the electromechanical side simulation;
and solving the difference between the state variable and the transmission current in the electromechanical side simulation, and judging the passivity characteristic.
7. The method for detecting stability of a multi-point electromagnetic hybrid simulation according to claim 6, wherein the range of the envelope of the difference between the state variable and the transmission current in the simulation at the electromechanical side determines whether the interface at the electromechanical side has an interface non-convergence problem, and the method comprises the following steps:
if the envelope curve of the difference between the state variable and the transmission current in the electromechanical side simulation is between the electromechanical side transmission variable and the interface element state variable and is smaller than the electromechanical side transmission variable, the electromechanical transient interface in the electromechanical-electromagnetic hybrid simulation has the problem of non-convergence stability;
and if the envelope curve of the difference between the state variable and the transmission current in the electromechanical side simulation is between the electromechanical side transmission variable and the interface element state variable and is greater than the electromechanical side transmission variable, the electromechanical transient interface in the electromechanical-electromagnetic hybrid simulation is converged.
8. The method for detecting stability of a multi-point electromagnetic hybrid simulation according to claim 7, wherein the determining the stability of the electromechanical-electromagnetic hybrid simulation based on the passivity characteristic and the interface convergence condition of the electromechanical side comprises:
if the different frequencies of the interface matrix do not have passivity characteristics and the interface on the electromechanical side is converged, the electromechanical-electromagnetic hybrid simulation is stable;
and if the different frequencies of the interface matrix have passivity characteristics or the interface on the electromechanical side does not converge, the electromechanical-electromagnetic hybrid simulation is unstable.
9. An electromechanical-electromagnetic hybrid simulation stability detection system, comprising:
an extraction module: the interface matrix is used for extracting the electromechanical-electromagnetic hybrid simulation and solving a network interface matrix;
an analysis module: the method is used for judging whether the interface on the electromechanical side has the problem of interface non-convergence or not by analyzing positive definite characteristics of the network interface matrix on different frequencies and obtaining the passivity characteristics of the interface matrix on different frequencies and the range of the envelope curve of the difference between the state variable and the transmission current in the electromechanical side simulation;
a judging module: and the stability of the electromechanical-electromagnetic hybrid simulation is judged based on the passivity characteristic and the interface convergence condition of the electromechanical side.
10. The system for detecting stability of electro-mechanical-electro-magnetic hybrid simulation according to claim 9, wherein the extracting module comprises: a processing submodule and a solving submodule;
the processing submodule is used for:
converting the equivalent matrix of the electromechanical network into a real part matrix and an imaginary part matrix;
transforming the real part matrix into a single-phase real part matrix of an ABC phase, and simultaneously transforming the imaginary part matrix into a single-phase imaginary part matrix of the ABC phase;
for the real part matrix YABCInversion processing to obtain ARInterface matrix, inverting imaginary matrix to obtain ALAn interface matrix;
the obtaining submodule is used for:
based on the ALInterface matrix for frequency w conversion selectionA resultant frequency point;
obtaining network interface matrixes under different frequencies through calculation;
wherein, A isLThe interface matrix performs frequency w conversion to select a suitable frequency point, which includes:
the frequency is increased to the range of 10KHz from 0.1Hz which is close to direct current, and m points are taken as frequency scanning points according to an equidistant logarithm form, wherein m is a positive integer;
the network interface matrix under different frequencies is calculated according to the following formula:
AG(w)=w*AL
in the formula, ALIs an interface matrix; w is the frequency; a. theGIs a network interface matrix.
11. The system according to claim 9, wherein the analysis module comprises: acquiring a submodule and a judging submodule;
the acquisition submodule is used for:
obtaining the current injected into the interface in the interface matrix simulation as a state variable;
acquiring current transmitted from the electromechanical side simulation;
solving the difference between the state variable and the transmission current in the electromechanical side simulation, and judging the passivity characteristic;
the judgment submodule is used for: judging whether an electromechanical transient interface in electromechanical-electromagnetic hybrid simulation converges;
if the envelope curve of the difference between the state variable and the transmission current in the electromechanical side simulation is between the electromechanical side transmission variable and the interface element state variable and is smaller than the electromechanical side transmission variable, the electromechanical transient interface in the electromechanical-electromagnetic hybrid simulation has the problem of non-convergence stability;
and if the envelope curve of the difference between the state variable and the transmission current in the electromechanical side simulation is between the electromechanical side transmission variable and the interface element state variable and is greater than the electromechanical side transmission variable, the electromechanical transient interface in the electromechanical-electromagnetic hybrid simulation is converged.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010052403.6A CN113139264A (en) | 2020-01-17 | 2020-01-17 | Electromechanical-electromagnetic hybrid simulation stability detection method and system |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010052403.6A CN113139264A (en) | 2020-01-17 | 2020-01-17 | Electromechanical-electromagnetic hybrid simulation stability detection method and system |
Publications (1)
Publication Number | Publication Date |
---|---|
CN113139264A true CN113139264A (en) | 2021-07-20 |
Family
ID=76808471
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010052403.6A Pending CN113139264A (en) | 2020-01-17 | 2020-01-17 | Electromechanical-electromagnetic hybrid simulation stability detection method and system |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113139264A (en) |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102521452A (en) * | 2011-12-14 | 2012-06-27 | 中国电力科学研究院 | Computing system of large power grid closed loop |
CN102708250A (en) * | 2012-05-10 | 2012-10-03 | 天津大学 | Power system digital hybrid simulation method based on tree hierarchical bidirectional iteration |
US20130124181A1 (en) * | 2011-11-11 | 2013-05-16 | Sigrity, Inc. | Methods, Systems, and Computer-readable Media for Improving Accuracy of Network Parameter in Electromagnetic Simulation and Modeling |
CN105095590A (en) * | 2015-08-11 | 2015-11-25 | 华北电力大学 | Method for modeling of electromechanical transient simulation system based on three-sequence equivalent impedance |
-
2020
- 2020-01-17 CN CN202010052403.6A patent/CN113139264A/en active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20130124181A1 (en) * | 2011-11-11 | 2013-05-16 | Sigrity, Inc. | Methods, Systems, and Computer-readable Media for Improving Accuracy of Network Parameter in Electromagnetic Simulation and Modeling |
CN102521452A (en) * | 2011-12-14 | 2012-06-27 | 中国电力科学研究院 | Computing system of large power grid closed loop |
CN102708250A (en) * | 2012-05-10 | 2012-10-03 | 天津大学 | Power system digital hybrid simulation method based on tree hierarchical bidirectional iteration |
CN105095590A (en) * | 2015-08-11 | 2015-11-25 | 华北电力大学 | Method for modeling of electromechanical transient simulation system based on three-sequence equivalent impedance |
Non-Patent Citations (3)
Title |
---|
H.T. SU 等: "Parallel interaction protocol for electromagnetic and electromechanical hybrid simulation", IEE PROCEEDINGS - GENERATION, vol. 152, no. 3, pages 406 - 414, XP006023994, DOI: 10.1049/ip-gtd:20049019 * |
张怡 等: "电磁–机电暂态混合仿真中的频率相关网络等值", 中国电机工程学报, vol. 32, no. 13, pages 61 - 68 * |
蔡普成等: "基于ADPSS的含背靠背MMC-HVDC***的交直流电网机电-电磁混合仿真研究", 电网技术, vol. 42, no. 12, pages 3888 - 3894 * |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US5243538A (en) | Comparison and verification system for logic circuits and method thereof | |
Henrion et al. | Solving polynomial static output feedback problems with PENBMI | |
Margaria et al. | Efficient test-based model generation for legacy reactive systems | |
Caplain | Finding invariant assertions for proving programs | |
CN104360920A (en) | Automatic testing method and device for interface | |
EP3912076A1 (en) | Lookup table optimization for programming languages that target synchronous digital circuits | |
Yu et al. | Generating vulnerability signatures for string manipulating programs using automata-based forward and backward symbolic analyses | |
CN110543407A (en) | Static analysis method for performance of identity intelligent contract | |
Xudong | Decentralized adaptive regulation with unknown high-frequency-gain signs | |
Rapaić et al. | Stability regions of fractional systems in the space of perturbed orders | |
CN114153422A (en) | Intelligent contract code design and generation method and system based on formalized model | |
CN113139264A (en) | Electromechanical-electromagnetic hybrid simulation stability detection method and system | |
US20210357563A1 (en) | Signal integrity simulation method for encryption hybrid model | |
EP2548329A1 (en) | Data expansion using an approximate method | |
Zhang et al. | Smart contract vulnerability detection method based on bi-lstm neural network | |
CN106649081B (en) | A kind of method and device that commissioning test wraps offline in client application | |
CN115150143B (en) | Network access authentication method, device, equipment and storage medium for industrial control equipment | |
Broutin et al. | The dual tree of a recursive triangulation of the disk | |
Cao et al. | An exponential transform based dual-layer Krylov subspace method for computing unstable modes of power systems | |
CN114090011A (en) | Software development method convenient for developer to use | |
CN110347968B (en) | FPGA-based FFT optimization algorithm and device | |
CN110908685A (en) | Device-level constant value batch processing method and system based on custom script | |
CN113704126A (en) | Verification method and device, computer storage medium and processor | |
CN112543104A (en) | Application program packaging method and device, electronic equipment and storage medium | |
KR102541065B1 (en) | Apparatus and method for producing s-parameter |
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 |