A method and apparatus for initializing an electromagnetic transient simulation circuit

Abstract

The application relates to a method and device for initializing an electromagnetic transient simulation circuit, which comprises the following steps: obtaining the historical time equivalent current of a node to be initialized in the electromagnetic transient simulation circuit; and determining the initialization voltage of the node to be initialized in the electromagnetic transient simulation circuit according to the historical time equivalent current of the node to be initialized in the electromagnetic transient simulation circuit. According to the application, the initialization voltage of the node in the electromagnetic transient simulation circuit is obtained according to the historical time equivalent current of the node to be simulated in the electromagnetic transient simulation circuit, numerical oscillation in simulation is avoided, the method is simple, simulation time is saved, and the method is easy to implement.

Description

Technical Field

[0001] This invention relates to the field of electromagnetic transient simulation technology, specifically to an electromagnetic transient simulation circuit initialization method and apparatus. Background Technology

[0002] With the continuous expansion and interconnection of power systems and the increasing number of power electronic devices, electromagnetic transient analysis technology has become a key technology for the safe and stable operation of power grids. However, when performing electromagnetic transient simulations of power systems, a reasonable initial operating point that satisfies the KVL and KCL constraints of the electrical system is required to avoid numerical oscillations in the simulation.

[0003] The commonly used initialization method for electromagnetic transient simulation in the prior art is zero-state simulation initialization, which is widely used in EMTP-type software. However, its disadvantages include slow program startup speed and the tendency for numerical oscillations to occur in ill-conditioned networks such as pure capacitor and power supply loops and pure inductor and current source cutsets, requiring the addition of numerical oscillation identification and suppression programs, thus increasing the complexity of the simulation program. In addition, EMTP-type software only provides the function of assigning initial values ​​to capacitor voltages, but not to inductor currents, meaning the system cannot start the simulation completely from a given state, resulting in significant simulation time consumption. Therefore, there is an urgent need in this field for a simple, easy-to-implement, and time-saving initialization method. Summary of the Invention

[0004] To address the shortcomings of existing technologies, the present invention aims to provide an electromagnetic transient simulation circuit initialization method and apparatus. This method obtains the node initialization voltage in the electromagnetic transient simulation circuit based on the equivalent current of the node to be simulated at historical moments, thus avoiding numerical oscillations during simulation. Furthermore, the method is simple, saves simulation time, and is easy to implement.

[0005] The objective of this invention is achieved through the following technical solution:

[0006] This invention provides an initialization method for electromagnetic transient simulation circuits, the improvement of which is that the method includes:

[0007] Obtain the historical equivalent current of the node to be initialized in the electromagnetic transient simulation circuit;

[0008] The initialization voltage of the node to be initialized in the electromagnetic transient simulation circuit is determined based on the equivalent current of the node to be initialized at a historical moment in the electromagnetic transient simulation circuit.

[0009] Preferably, obtaining the historical equivalent current of the node to be initialized in the electromagnetic transient simulation circuit includes:

[0010] The equivalent current at the historical moment of the i-th node to be initialized in the electromagnetic transient simulation circuit is determined by the following formula.

[0011]

[0012] In the formula, Let be the equivalent current of all branch capacitors of the i-th node to be initialized in the electromagnetic transient simulation circuit at all historical moments. Let t0 be the equivalent current of the branch inductance of the i-th node to be initialized in the electromagnetic transient simulation circuit at a historical moment, where t0 is the historical moment, i∈[S+1,N], N+1 is the total number of nodes in the electromagnetic transient simulation circuit, and S+1 is the total number of nodes with known voltage in the electromagnetic transient simulation circuit.

[0013] Furthermore, the current of the branch capacitor of the i-th node to be initialized in the electromagnetic transient simulation circuit at a historical moment is determined by the following formula.

[0014]

[0015] In the formula, C i,j Let Δt be the capacitance value of the j-th branch capacitor of the i-th node to be initialized in the electromagnetic transient simulation circuit, and Δt be the simulation step size. Let be the initial voltage of the j-th branch capacitor of the i-th node to be initialized in the electromagnetic transient simulation circuit;

[0016] The equivalent current of the branch inductance of the i-th node to be initialized in the electromagnetic transient simulation circuit at a historical moment is determined by the following formula:

[0017]

[0018] In the formula, Let be the initial current of the inductor of the j-th branch of the i-th node to be initialized in the electromagnetic transient simulation circuit.

[0019] Preferably, determining the initialization voltage of the node to be initialized in the electromagnetic transient simulation circuit based on the historical equivalent current of the node to be initialized in the electromagnetic transient simulation circuit includes:

[0020] The initial voltage u of the node to be initialized in the electromagnetic transient simulation circuit is determined by the following formula. X :

[0021] u X =(G XX ) -1 (II t0 -G XS u S )

[0022] In the formula, I is the column vector of branch current source current values ​​of the node to be initialized in the electromagnetic transient simulation circuit, I = [I S+1 … I i … I N ] T Let i ∈ [S+1, N], where N+1 is the total number of nodes in the electromagnetic transient simulation circuit, S+1 is the total number of nodes with known voltage in the electromagnetic transient simulation circuit, and X is the total number of nodes to be initialized in the electromagnetic transient simulation circuit, S+X=N, I i Let t0 be the initial value of the branch current source of the i-th node to be initialized in the electromagnetic transient simulation circuit, and t0 be the historical time. This represents the equivalent current column vector of the node to be initialized at historical moments in the electromagnetic transient simulation circuit. Let u be the equivalent current at a historical moment of the i-th node to be initialized in the electromagnetic transient simulation circuit. S For the initial voltage of a node with known voltage in an electromagnetic transient simulation circuit, G XX Let G be the admittance matrix formed by rows S+1 to N and columns S+1 to N of the admittance matrix of all nodes in the electromagnetic transient simulation circuit. XS Let G be the admittance matrix formed by rows S+1 to N and columns 0 to S of the admittance matrix of all nodes in the electromagnetic transient simulation circuit.

[0023] The admittance matrix of all nodes in the electromagnetic transient simulation circuit is obtained by discretizing the electromagnetic transient simulation circuit using the backward Euler method. The nodes with known voltage include grounding nodes and nodes connected to the positive terminal of the grounding voltage source.

[0024] Furthermore, the admittance matrix G of all nodes in the electromagnetic transient simulation circuit is determined by the following formula:

[0025]

[0026] In the formula, n∈{0,…,S,S+1,…,N}, n'∈{0,…,S,S+1,…,N}, when n≠n', G nn' G is the mutual admittance of the nth node in the electromagnetic transient simulation circuit. When n = n', G nn' Let {S} be the self-admittance of the nth node in the electromagnetic transient simulation circuit, {0,…,S} be nodes with known voltages, and {S+1,…,N} be nodes to be initialized.

[0027] This invention provides an electromagnetic transient simulation circuit initialization device, the improvement of which is that the device includes:

[0028] The acquisition unit is used to acquire the historical equivalent current of the node to be initialized in the electromagnetic transient simulation circuit;

[0029] An initialization unit is used to determine the initialization voltage of the node to be initialized in the electromagnetic transient simulation circuit based on the equivalent current of the node to be initialized at a historical moment in the electromagnetic transient simulation circuit.

[0030] Preferably, the acquisition unit is specifically used for

[0031] The equivalent current at the historical moment of the i-th node to be initialized in the electromagnetic transient simulation circuit is determined by the following formula.

[0032]

[0033] In the formula, Let be the equivalent current of all branch capacitors of the i-th node to be initialized in the electromagnetic transient simulation circuit at all historical moments. Let t0 be the equivalent current of the branch inductance of the i-th node to be initialized in the electromagnetic transient simulation circuit at a historical moment, where t0 is the historical moment, i∈[S+1,N], N+1 is the total number of nodes in the electromagnetic transient simulation circuit, and S+1 is the total number of nodes with known voltage in the electromagnetic transient simulation circuit.

[0034] Furthermore, the current of the branch capacitor of the i-th node to be initialized in the electromagnetic transient simulation circuit at a historical moment is determined by the following formula.

[0035]

[0036] In the formula, C i,j Let Δt be the capacitance value of the j-th branch capacitor of the i-th node to be initialized in the electromagnetic transient simulation circuit, and Δt be the simulation step size. Let be the initial voltage of the j-th branch capacitor of the i-th node to be initialized in the electromagnetic transient simulation circuit;

[0037] The equivalent current of the branch inductance of the i-th node to be initialized in the electromagnetic transient simulation circuit at a historical moment is determined by the following formula:

[0038]

[0039] In the formula, Let be the initial current of the inductor of the j-th branch of the i-th node to be initialized in the electromagnetic transient simulation circuit.

[0040] Furthermore, the initialization unit is specifically used for:

[0041] The initial voltage u of the node to be initialized in the electromagnetic transient simulation circuit is determined by the following formula. X :

[0042] u X =(GXX ) -1 (II t0 -G XS u S )

[0043] In the formula, I is the column vector of branch current sources of the node to be initialized in the electromagnetic transient simulation circuit, I = [I S+1 … I i … I N ] T Let i ∈ [S+1, N], where N+1 is the total number of nodes in the electromagnetic transient simulation circuit, S+1 is the total number of nodes with known voltage in the electromagnetic transient simulation circuit, and X is the total number of nodes to be initialized in the electromagnetic transient simulation circuit, S+X=N, I i Let t0 be the initial value of the branch current source of the i-th node to be initialized in the electromagnetic transient simulation circuit, and t0 be the historical time. This represents the equivalent current column vector of the node to be initialized at historical moments in the electromagnetic transient simulation circuit. Let u be the equivalent current at a historical moment of the i-th node to be initialized in the electromagnetic transient simulation circuit. S For the initial voltage of a node with known voltage in an electromagnetic transient simulation circuit, G XX Let G be the admittance matrix of row S+1 to row N and column S+1 to column N of the admittance matrix of all nodes in the electromagnetic transient simulation circuit. XS Let G be the admittance matrix from row S+1 to row N and column 0 to column S in the admittance matrix G of all nodes in the electromagnetic transient simulation circuit.

[0044] The admittance matrix of all nodes in the electromagnetic transient simulation circuit is obtained by discretizing the electromagnetic transient simulation circuit using the backward Euler method. The nodes with known voltage include grounding nodes and nodes connected to the positive terminal of the grounding voltage source.

[0045] Furthermore, the admittance matrix G of all nodes in the electromagnetic transient simulation circuit is determined by the following formula:

[0046]

[0047] In the formula, n∈{0,…,S,S+1,…,N}, n'∈{0,…,S,S+1,…,N}, when n≠n', G nn' G is the mutual admittance of the nth node in the electromagnetic transient simulation circuit. When n = n', G nn' Let {S} be the self-admittance of the nth node in the electromagnetic transient simulation circuit, {0,…,S} be nodes with known voltages, and {S+1,…,N} be nodes to be initialized.

[0048] Compared with the closest existing technology, the present invention has the following advantages:

[0049] This invention provides an initialization method and apparatus for an electromagnetic transient simulation circuit, comprising: acquiring the historical equivalent current of a node to be initialized in the electromagnetic transient simulation circuit; and determining the initialization voltage of the node to be initialized in the electromagnetic transient simulation circuit based on the historical equivalent current of the node to be initialized in the electromagnetic transient simulation circuit. This invention obtains the node initialization voltage in the electromagnetic transient simulation circuit based on the historical equivalent current of the node to be simulated in the electromagnetic transient simulation circuit, which not only avoids numerical oscillations in the simulation but also simplifies the method, saves simulation time, and is easy to implement. Attached Figure Description

[0050] Figure 1 This is a flowchart of the electromagnetic transient simulation circuit initialization method of the present invention;

[0051] Figure 2 This is an electromagnetic transient simulation circuit diagram in an embodiment of the present invention;

[0052] Figure 3 This is a structural diagram of the electromagnetic transient simulation circuit initialization device of the present invention. Detailed Implementation

[0053] The specific embodiments of the present invention will be further described in detail below with reference to the accompanying drawings.

[0054] To make the objectives, technical solutions, and advantages of the embodiments of the present invention clearer, the technical solutions of the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some embodiments of the present invention, not all embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those skilled in the art without creative effort are within the scope of protection of the present invention.

[0055] This invention provides an initialization method for electromagnetic transient simulation circuits, such as... Figure 1 As shown, the method includes:

[0056] Obtain the historical equivalent current of the node to be initialized in the electromagnetic transient simulation circuit;

[0057] The initialization voltage of the node to be initialized in the electromagnetic transient simulation circuit is determined based on the equivalent current of the node to be initialized at a historical moment in the electromagnetic transient simulation circuit.

[0058] The above method obtains the historical equivalent current of the node to be initialized in the electromagnetic transient simulation circuit, including:

[0059] The equivalent current at the historical moment of the i-th node to be initialized in the electromagnetic transient simulation circuit is determined by the following formula.

[0060]

[0061] In the formula, Let be the equivalent current of all branch capacitors of the i-th node to be initialized in the electromagnetic transient simulation circuit at all historical moments. Let t0 be the equivalent current of the branch inductance of the i-th node to be initialized in the electromagnetic transient simulation circuit at a historical moment, where t0 is the historical moment, i∈[S+1,N], N+1 is the total number of nodes in the electromagnetic transient simulation circuit, and S+1 is the total number of nodes with known voltage in the electromagnetic transient simulation circuit.

[0062] Furthermore, the current of the branch capacitor of the i-th node to be initialized in the electromagnetic transient simulation circuit at a historical moment is determined by the following formula.

[0063]

[0064] In the formula, C i,j Let Δt be the capacitance value of the j-th branch capacitor of the i-th node to be initialized in the electromagnetic transient simulation circuit, and Δt be the simulation step size. Let be the initial voltage of the j-th branch capacitor of the i-th node to be initialized in the electromagnetic transient simulation circuit;

[0065] The equivalent current of the branch inductance of the i-th node to be initialized in the electromagnetic transient simulation circuit at a historical moment is determined by the following formula:

[0066]

[0067] In the formula, Let be the initial current of the inductor of the j-th branch of the i-th node to be initialized in the electromagnetic transient simulation circuit.

[0068] The method described above, which determines the initialization voltage of the node to be initialized in the electromagnetic transient simulation circuit based on the historical equivalent current of the node to be initialized in the electromagnetic transient simulation circuit, includes:

[0069] The initial voltage u of the node to be initialized in the electromagnetic transient simulation circuit is determined by the following formula. X :

[0070] u X =(G XX ) -1 (II t0 -G XS u S )

[0071] In the formula, I is the column vector of branch current sources of the node to be initialized in the electromagnetic transient simulation circuit, I = [I S+1 … I i … I N ] T Let i ∈ [S+1, N], where N+1 is the total number of nodes in the electromagnetic transient simulation circuit, S+1 is the total number of nodes with known voltage in the electromagnetic transient simulation circuit, and X is the total number of nodes to be initialized in the electromagnetic transient simulation circuit, S+X=N, I i Let t0 be the initial value of the branch current source of the i-th node to be initialized in the electromagnetic transient simulation circuit, and t0 be the historical time. This represents the equivalent current column vector of the node to be initialized at historical moments in the electromagnetic transient simulation circuit. Let u be the equivalent current at a historical moment of the i-th node to be initialized in the electromagnetic transient simulation circuit. S For the initial voltage of a node with known voltage in an electromagnetic transient simulation circuit, G XX Let G be the admittance matrix of row S+1 to row N and column S+1 to column N of the admittance matrix of all nodes in the electromagnetic transient simulation circuit. XS Let G be the admittance matrix from row S+1 to row N and column 0 to column S in the admittance matrix G of all nodes in the electromagnetic transient simulation circuit.

[0072] The admittance matrix of all nodes in the electromagnetic transient simulation circuit is obtained by discretizing the electromagnetic transient simulation circuit using the backward Euler method. The nodes with known voltage include grounding nodes and nodes connected to the positive terminal of the grounding voltage source.

[0073] Furthermore, the admittance matrix G of all nodes in the electromagnetic transient simulation circuit is determined by the following formula:

[0074]

[0075] In the formula, n∈{0,…,S,S+1,…,N}, n'∈{0,…,S,S+1,…,N}, when n≠n', G nn' G is the mutual admittance of the nth node in the electromagnetic transient simulation circuit. When n = n', G nn' Let {S} be the self-admittance of the nth node in the electromagnetic transient simulation circuit, {0,…,S} be nodes with known voltages, and {S+1,…,N} be nodes to be initialized.

[0076] In an embodiment of the present invention, the backward Euler method is used to calculate the step size Δt = 1 μs, such as... Figure 2As shown, the electromagnetic transient simulation circuit includes 6 nodes. Nodes 0 and 1 are nodes with known voltages, totaling 2. Nodes 3 to 6 are nodes to be initialized, totaling 4. The given electromagnetic transient simulation circuit includes: voltage source U. s1 The amplitude is 2V, the frequency is 1Hz, the capacitors are C1=1F, C2=2F, C3=0.5F, L1=2H, L2=1H, L3=2H, R1=1Ω, R2=2Ω. At a given historical time t0=t-Δt, the initial voltages of capacitors C1, C2, and C3 are 0V, 2V, and 8V respectively, and the initial currents of inductors L1, L2, and L3 are 4A, 4A, and 6A respectively.

[0077] Among them, I t0 =[0 4000000 10 -10] T I = [0 0 0 0] T u S =[u0 u1] T =[0 2] T ;

[0078] Discretizing it using the backward Euler method includes:

[0079] Discretize all nodes in the electromagnetic transient simulation circuit using the following method:

[0080] 1) The circuit equation for the resistor is:

[0081] u m1 -u n1 =R·i (1)

[0082] Among them, u m1 u n1 Let be the voltage across the nodes of the resistor element, and i be the current flowing through the resistor element.

[0083] Discretize equation (1) using the backward Euler method:

[0084]

[0085] 2) The circuit equation for the capacitor is:

[0086]

[0087] Among them, u m2 u n2 Let u be the node voltage across the capacitor element. c =u m2 -u n2 .

[0088] Discretize equation (3) using the backward Euler method:

[0089]

[0090] The capacitor C is equivalent to a resistance of... The resistor and the current source at the historical time t0 = t - Δt are connected in parallel.

[0091] 3) The circuit equation for an inductor is:

[0092]

[0093] Among them, u m3 u n3 Let u be the node voltage across the inductor element. L =u m3 -u n3 .

[0094] Discretize equation (5) using the backward Euler method:

[0095]

[0096] The inductance L is equivalent to a resistance of... The resistor and the current source at the historical time t0 = t - Δt are connected in parallel.

[0097] Its admittance matrix is ​​obtained as follows:

[0098]

[0099] in,

[0100]

[0101] Initial voltage u of the node to be initialized X =[2 10 10 -20] T Therefore, the node voltage of the electromagnetic transient simulation circuit

[0102] Let u = [0 2 2 10 10 -20] T

[0103] This invention provides an electromagnetic transient simulation circuit initialization device, such as... Figure 3 As shown, the device includes:

[0104] The acquisition unit is used to acquire the historical equivalent current of the node to be initialized in the electromagnetic transient simulation circuit;

[0105] An initialization unit is used to determine the initialization voltage of the node to be initialized in the electromagnetic transient simulation circuit based on the equivalent current of the node to be initialized at a historical moment in the electromagnetic transient simulation circuit.

[0106] Preferably, the acquisition unit is specifically used for

[0107] The equivalent current at the historical moment of the i-th node to be initialized in the electromagnetic transient simulation circuit is determined by the following formula.

[0108]

[0109] In the formula, Let be the equivalent current of all branch capacitors of the i-th node to be initialized in the electromagnetic transient simulation circuit at all historical moments. Let t0 be the equivalent current of the branch inductance of the i-th node to be initialized in the electromagnetic transient simulation circuit at a historical moment, where t0 is the historical moment, i∈[S+1,N], N+1 is the total number of nodes in the electromagnetic transient simulation circuit, and S+1 is the total number of nodes with known voltage in the electromagnetic transient simulation circuit.

[0110] Furthermore, the current of the branch capacitor of the i-th node to be initialized in the electromagnetic transient simulation circuit at a historical moment is determined by the following formula.

[0111]

[0112] In the formula, C i,j Let Δt be the capacitance value of the j-th branch capacitor of the i-th node to be initialized in the electromagnetic transient simulation circuit, and Δt be the simulation step size. Let be the initial voltage of the j-th branch capacitor of the i-th node to be initialized in the electromagnetic transient simulation circuit;

[0113] The equivalent current of the branch inductance of the i-th node to be initialized in the electromagnetic transient simulation circuit at a historical moment is determined by the following formula:

[0114]

[0115] In the formula, Let be the initial current of the inductor of the j-th branch of the i-th node to be initialized in the electromagnetic transient simulation circuit.

[0116] Furthermore, the initialization unit is specifically used for:

[0117] The initial voltage u of the node to be initialized in the electromagnetic transient simulation circuit is determined by the following formula. X :

[0118] u X =(G XX ) -1 (II t0 -G XS u S )

[0119] In the formula, I is the column vector of branch current sources of the node to be initialized in the electromagnetic transient simulation circuit, I = [I S+1 … I i … I N ] T Let i ∈ [S+1, N], where N+1 is the total number of nodes in the electromagnetic transient simulation circuit, S+1 is the total number of nodes with known voltage in the electromagnetic transient simulation circuit, and X is the total number of nodes to be initialized in the electromagnetic transient simulation circuit, S+X=N, I i Let t0 be the initial value of the branch current source of the i-th node to be initialized in the electromagnetic transient simulation circuit, and t0 be the historical time. This represents the equivalent current column vector of the node to be initialized at historical moments in the electromagnetic transient simulation circuit. Let u be the equivalent current at a historical moment of the i-th node to be initialized in the electromagnetic transient simulation circuit. S For the initial voltage of a node with known voltage in an electromagnetic transient simulation circuit, G XX Let G be the admittance matrix of row S+1 to row N and column S+1 to column N of the admittance matrix of all nodes in the electromagnetic transient simulation circuit. XS Let G be the admittance matrix from row S+1 to row N and column 0 to column S in the admittance matrix G of all nodes in the electromagnetic transient simulation circuit.

[0120] The admittance matrix of all nodes in the electromagnetic transient simulation circuit is obtained by discretizing the electromagnetic transient simulation circuit using the backward Euler method. The nodes with known voltage include grounding nodes and nodes connected to the positive terminal of the grounding voltage source.

[0121] Furthermore, the admittance matrix G of all nodes in the electromagnetic transient simulation circuit is determined by the following formula:

[0122]

[0123] In the formula, n∈{0,…,S,S+1,…,N}, n'∈{0,…,S,S+1,…,N}, when n≠n', G nn' G is the mutual admittance of the nth node in the electromagnetic transient simulation circuit. When n = n', G nn' Let {S} be the self-admittance of the nth node in the electromagnetic transient simulation circuit, {0,…,S} be nodes with known voltages, and {S+1,…,N} be nodes to be initialized.

[0124] Compared with the closest existing technology, the present invention has the following advantages:

[0125] This invention provides an initialization method and apparatus for an electromagnetic transient simulation circuit, comprising: acquiring the historical equivalent current of a node to be initialized in the electromagnetic transient simulation circuit; and determining the initialization voltage of the node to be initialized in the electromagnetic transient simulation circuit based on the historical equivalent current of the node to be initialized in the electromagnetic transient simulation circuit. This invention obtains the node initialization voltage in the electromagnetic transient simulation circuit based on the historical equivalent current of the node to be simulated in the electromagnetic transient simulation circuit, which not only avoids numerical oscillations in the simulation but also simplifies the method, saves simulation time, and is easy to implement.

[0126] Those skilled in the art will understand that embodiments of this application can be provided as methods, systems, or computer program products. Therefore, this application can take the form of a completely hardware embodiment, a completely software embodiment, or an embodiment combining software and hardware aspects. Furthermore, this application can 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, etc.) containing computer-usable program code.

[0127] This application is described with reference to flowchart illustrations and / or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of this application. It will be understood that each block of the flowchart illustrations and / or block diagrams, and combinations of blocks in the flowchart illustrations and / or block diagrams, can be implemented by computer program instructions. These computer program instructions can 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, generate instructions for implementing the flowchart... Figure 1 One or more processes and / or boxes Figure 1 A device that provides the functions specified in one or more boxes.

[0128] These computer program instructions may also be stored in a computer-readable storage medium that can direct a computer or other programmable data processing device to function in a particular manner, such that the instructions stored in the computer-readable storage medium produce an article of manufacture including instruction means, which are implemented in a process Figure 1 One or more processes and / or boxes Figure 1 The function specified in one or more boxes.

[0129] These computer program instructions may also be loaded onto a computer or other programmable data processing equipment to cause a series of operational steps to be performed on the computer or other programmable equipment to produce a computer-implemented process, thereby providing instructions that execute on the computer or other programmable equipment for implementing the process. Figure 1 One or more processes and / or boxes Figure 1The steps of the function specified in one or more boxes.

[0130] Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention and not to limit it. Although the present invention has been described in detail with reference to the above embodiments, those skilled in the art should understand that modifications or equivalent substitutions can still be made to the specific implementation of the present invention. Any modifications or equivalent substitutions that do not depart from the spirit and scope of the present invention should be covered within the protection scope of the claims of the present invention.

Claims

1. A method for initializing an electromagnetic transient simulation circuit, characterized in that, The method includes: Obtain the historical equivalent current of the node to be initialized in the electromagnetic transient simulation circuit; The initial voltage of the node to be initialized in the electromagnetic transient simulation circuit is determined based on the equivalent current of the node to be initialized at a historical moment in the electromagnetic transient simulation circuit. The step of determining the initialization voltage of the node to be initialized in the electromagnetic transient simulation circuit based on the equivalent current of the node to be initialized at a historical moment in the electromagnetic transient simulation circuit includes: The initial voltage of the node to be initialized in the electromagnetic transient simulation circuit is determined by the following formula. : In the formula, This is a column vector of branch current source current values ​​for the node to be initialized in the electromagnetic transient simulation circuit. , N is the number of nodes in the electromagnetic transient simulation circuit. This represents the total number of nodes with known voltages in the electromagnetic transient simulation circuit. This represents the total number of nodes to be initialized in the electromagnetic transient simulation circuit. , The first in the electromagnetic transient simulation circuit Initialization values ​​of branch current sources for each node to be initialized. For a historic moment, , This represents the equivalent current column vector of the node to be initialized at historical moments in the electromagnetic transient simulation circuit. Electromagnetic In transient simulation circuits, the first The equivalent current at historical moments of each node to be initialized. The initial voltage of the node with known voltage in the electromagnetic transient simulation circuit. Let G be the admittance matrix formed by rows S+1 to N and columns S+1 to N of the admittance matrix of all nodes in the electromagnetic transient simulation circuit. Let G be the admittance matrix formed by rows S+1 to N and columns 0 to S of the admittance matrix of all nodes in the electromagnetic transient simulation circuit. The admittance matrix of all nodes in the electromagnetic transient simulation circuit is obtained by discretizing the electromagnetic transient simulation circuit using the backward Euler method. The nodes with known voltage include grounding nodes and nodes connected to the positive terminal of the grounding voltage source.

2. The method as described in claim 1, characterized in that, The acquisition of the historical equivalent current of the node to be initialized in the electromagnetic transient simulation circuit includes: The first step in the electromagnetic transient simulation circuit is determined by the following formula. Historical equivalent current of each node to be initialized : In the formula, The first in the electromagnetic transient simulation circuit The equivalent current of all branch capacitors of the node to be initialized at any historical moment. The first in the electromagnetic transient simulation circuit The equivalent current of the branch inductance of each node to be initialized at any given time in its historical data. For a historic moment, denoted as the total number of nodes with known voltage in the electromagnetic transient simulation circuit, and N is the number of nodes in the electromagnetic transient simulation circuit.

3. The method as described in claim 2, characterized in that, The equivalent current of the branch capacitor at the historical moment of the i-th node to be initialized in the electromagnetic transient simulation circuit is determined by the following formula: In the formula, The first in the electromagnetic transient simulation circuit The node to be initialized. The capacitance value of each branch capacitor. For simulating step size, The first in the electromagnetic transient simulation circuit The node to be initialized. The initial voltage of each branch capacitor; The electromagnetic transient simulation circuit is determined by the following formula: The equivalent current of the branch inductance at a given historical moment for each node to be initialized. : In the formula, The first in the electromagnetic transient simulation circuit The node to be initialized. The initial current of each branch inductor.

4. The method as described in claim 1, characterized in that, The admittance matrix G of all nodes in the electromagnetic transient simulation circuit is determined by the following formula: In the formula, , ,when hour, The first in the electromagnetic transient simulation circuit The mutual admittance of each node, when hour, The first in the electromagnetic transient simulation circuit The self-admittance of each node, For nodes with known voltages, This is the node to be initialized.

5. An electromagnetic transient simulation circuit initialization device, characterized in that, The device includes: The acquisition unit is used to acquire the historical equivalent current of the node to be initialized in the electromagnetic transient simulation circuit; An initialization unit is used to determine the initialization voltage of the node to be initialized in the electromagnetic transient simulation circuit based on the equivalent current of the node to be initialized at a historical moment in the electromagnetic transient simulation circuit. The initialization unit is specifically used for: The initial voltage of the node to be initialized in the electromagnetic transient simulation circuit is determined by the following formula. : In the formula, This represents the column vector of branch current sources for the node to be initialized in the electromagnetic transient simulation circuit. , N is the number of nodes in the electromagnetic transient simulation circuit. This represents the total number of nodes with known voltages in the electromagnetic transient simulation circuit. This represents the total number of nodes to be initialized in the electromagnetic transient simulation circuit. , The first in the electromagnetic transient simulation circuit Initialization values ​​of branch current sources for each node to be initialized. For a historic moment, , This represents the equivalent current column vector of the node to be initialized at historical moments in the electromagnetic transient simulation circuit. The first in the electromagnetic transient simulation circuit The equivalent current at historical moments of each node to be initialized. The initial voltage of the node with known voltage in the electromagnetic transient simulation circuit. Let G be the admittance matrix formed by rows S+1 to N and columns S+1 to N of the admittance matrix of all nodes in the electromagnetic transient simulation circuit. Let G be the admittance matrix formed by rows S+1 to N and columns 0 to S of the admittance matrix of all nodes in the electromagnetic transient simulation circuit. The admittance matrix of all nodes in the electromagnetic transient simulation circuit is obtained by discretizing the electromagnetic transient simulation circuit using the backward Euler method. The nodes with known voltage include grounding nodes and nodes connected to the positive terminal of the grounding voltage source.

6. The apparatus as claimed in claim 5, characterized in that, The acquisition unit is specifically used for The first step in the electromagnetic transient simulation circuit is determined by the following formula. Historical equivalent current of each node to be initialized : In the formula, The first in the electromagnetic transient simulation circuit The equivalent current of all branch capacitors of the node to be initialized at any historical moment. The first in the electromagnetic transient simulation circuit The equivalent current of the branch inductance of each node to be initialized at any given time in its historical data. For a historic moment, , denoted as the total number of nodes with known voltage in the electromagnetic transient simulation circuit, and N is the number of nodes in the electromagnetic transient simulation circuit.

7. The apparatus as claimed in claim 6, characterized in that, The current of the branch capacitor of the i-th node to be initialized in the electromagnetic transient simulation circuit at a historical moment is determined by the following formula: In the formula, The first in the electromagnetic transient simulation circuit The node to be initialized. The capacitance value of each branch capacitor. For simulating step size, The first in the electromagnetic transient simulation circuit The node to be initialized. The initial voltage of each branch capacitor; The electromagnetic transient simulation circuit is determined by the following formula: The equivalent current of the branch inductance at a given historical moment for each node to be initialized. : In the formula, The first in the electromagnetic transient simulation circuit The node to be initialized. The initial current of each branch inductor.

8. The apparatus as claimed in claim 5, characterized in that, The admittance matrix G of all nodes in the electromagnetic transient simulation circuit is determined by the following formula: In the formula, , ,when hour, The first in the electromagnetic transient simulation circuit The mutual admittance of each node, when hour, The first in the electromagnetic transient simulation circuit The self-admittance of each node, For nodes with known voltages, This is the node to be initialized.

Images