JP2919336B2 - Power system calculator - Google Patents

Description

DETAILED DESCRIPTION OF THE INVENTION

[0001]

BACKGROUND OF THE INVENTION 1. Field of the Invention The present invention relates to a power system calculation device, and more particularly to a power system calculation device for performing an analysis calculation on a power system including a plurality of generators and a plurality of loads and a transmission network interconnecting them.

[0002]

2. Description of the Related Art The electric power system is the core of the electric power business, and its contents include power generation facilities (hereinafter, referred to as power generation facilities).
A power transmission facility that distributes the power generated by the generator, a substation facility that sends out load power (hereinafter called a load), and a demand facility that consumes the power. It is necessary to operate the power system efficiently without causing a frequency change or a local voltage change between these facilities.

[0003] The analysis calculation of the electric power system is performed by expanding the system scale,
Demands for high-speed processing are increasing due to demands for immediate response.

On the other hand, recent advances in microprocessors have been remarkable. By performing power system calculations using a parallel computer system in which a number of these microprocessors are coupled, it is expected that a power system calculation device with much better cost performance than conventional systems can be constructed.

[0005] Generally, in the analysis calculation of the power system,
To solve the simultaneous ordinary differential equations for the internal states of multiple generators and multiple loads, calculate the data at each time of the simultaneous ordinary differential equations for the internal states of the generators and loads, and establish them between them. Repeated grid calculations.

[0006] The percentage of time for these calculations depends on the algorithm used. Therefore, it is necessary to efficiently perform both the calculation of the simultaneous ordinary differential equations relating to the internal states of the generator and the load, and the calculation of the transmission network.

[0007] In the latter transmission network calculation, the main content is the solution of simultaneous equations using LU decomposition. In general, in the case of analytical calculation of a power system, the number of couplings between a generator and a load is small. Therefore, when a simultaneous equation to be solved is represented as a matrix, the non-zero element rate of this matrix becomes very small.

When such simultaneous equations of a matrix having a very small non-zero element ratio are solved in parallel, the granularity of the parallel processing is fine, and a parallel memory system of a distributed memory type requiring explicit communication between processors is required. Is more efficient in a bus-coupled shared-memory parallel computer system in which a plurality of processors perform operations while referencing / updating data arranged on a shared memory connected to a bus (Hironori Kasahara:
System software for parallel processing, information processing Vo
l. 34, no. 9, pp. 1134-1138 (19
93)).

[0009]

On the other hand, in the calculation of the simultaneous ordinary differential equations relating to the internal state of the generator and the load, the data is completely independent for each generator and load. When such a calculation is performed by a parallel computer system of the bus-coupled shared memory system, when each processor refers to the shared memory via the bus at the same time, even though the data is completely independent, one processor is simultaneously used. Since only the processor can refer to the shared memory, there is a problem that other processors cannot proceed with the calculation.

[0010] Furthermore, in the parallel computer system of the bus-coupled shared memory system, there is a problem that a wait state occurs when referring not only to data but also to a program (instruction group executed by a processor).

[0011] It is an object of the present invention to secure the superiority of a parallel computer system of a bus-coupled shared memory system in power transmission network calculation, while maintaining the independence of operation in calculating simultaneous ordinary differential equations relating to the internal state of a generator and a load. An object of the present invention is to provide a power system calculation device that makes effective use of the power system.

[0012]

SUMMARY OF THE INVENTION A power system calculation apparatus according to the present invention performs analysis calculation for a power system composed of a plurality of generators and a plurality of loads and a transmission network interconnecting them, every hour step. In a power system calculation device that calculates a simultaneous ordinary differential equation relating to the internal state of a generator and a load by repeatedly solving a simultaneous equation for a transmission network calculation established between them, a plurality of hardware A bus that enables the connection of hardware,
A central processing unit connected to the bus and performing calculations;
A plurality of processors each comprising a program executed by the central processing unit and a local memory, and holding data relating to the internal state of the generator and the load assigned to the processor among the data relating to the internal state of the generator and the load; And a shared memory connected to the bus and holding data necessary for the transmission network calculation.In the step of calculating the simultaneous ordinary differential equations relating to the internal states of the generator and the load, each of the processors In the step of performing a calculation using data on the internal state of the generator and the load on the local memory and calculating the power grid calculation by solving simultaneous equations, the processors transmit the power transmission data on the shared memory. For net calculation
It is characterized in that the allocated transmission network is calculated using necessary data .

Further, in the power system calculating apparatus according to the present invention, in the step of calculating the simultaneous ordinary differential equations relating to the internal states of the generator and the load, each of the processors includes:
After reading the connection terminal voltage data of the power grid regarding the generator and the load assigned to the processor from the shared memory, the internal state of the generator and the load assigned to the processor and existing in the local memory is stored. data using the independently calculated based on the Runge-Kutta method, to transfer the internal induced voltage of the resulting generator, the Adomittan <br/> scan of the admittance and the load of the generator in the shared memory, wherein the transmission net calculate the coalition side
In the step of calculating by solving the equation, each of the processors is connected to the inside of the generator on the shared memory.
Induced voltage, generator admittance and load admittance
Perform grid calculations assigned to the processor while referring to and updating the stance .

[0014] The power system calculation device according to the present invention includes:
In the process of calculating a simultaneous ordinary differential equation relating to the internal state of the generator and the load, a correction equation of a variable indicating the internal state of the generator and the load based on the trapezoidal formula, and a variable indicating the state of the power transmission network. In the case where each state variable solution is obtained by iterative correction using the Newton-Raphson method by simultaneously linking the correction equation and the system, after reading the connection terminal voltage data of the power grid regarding the generator and the load from the shared memory, the generator and The residuals of the equations relating to the internal state variables of the load and the elements of the modified equation coefficient matrix (Jacobi matrix) are independently calculated using the data relating to the internal state of the generator and the load assigned to the processor, and stored in the local memory. The residuals of the equations relating to the state variables of the power grid and the elements of the modified equation coefficient matrix (Jacobi matrix) The calculated using the data necessary to the grid computing assigned to the processor with data necessary for stored was power grid calculated local memory and the shared memory, stored the result in the shared memory,
In the step of calculating the correction equation, for the correction equation coefficient matrix relating to the generator and the load, refer to the residual and the Jacobian matrix element of the local memory for the LU factorization and forward elimination of the elements allocated to the processor. If the data of the portion relating to the power transmission network needs to be updated in accordance with the calculation of the generator and the load by the LU decomposition and the forward elimination operation, the data on the shared memory is updated. , For each of the corrected equation coefficient matrix and the residual of the power transmission network, each processor executes an LU factorization and forward elimination operation of elements assigned to the processor while referring / updating data on the shared memory, Each processor uses the data on the shared memory to determine the amount of correction related to the state variable of the power grid. Determined by performing the backward substitution calculation on the residual subjected to forward elimination, the amount of correction of the state variables related to the generator and the load, each processor from the shared memory the amount of correction data of the state variable of the power grid from the shared memory After the reading, the Jacobi matrix on the local memory and the residual data subjected to forward elimination and back substitution are determined by back substitution calculation to correct all state variables, and the correction amount of all state variables is reduced. Repeat calculations until it is.

[0015]

Next, the present invention will be described in detail with reference to the drawings.

[0016]

FIG. 1 is a block diagram showing the configuration of a power system calculation apparatus 1 according to one embodiment of the present invention. The power system calculation device 1 of the present embodiment has a configuration in which a plurality of processors 10 and a shared memory 16 are connected by a bus 15.

Each processor 10 includes a central processing unit 11
And the local memory 12.

The local memory 12 stores a program 13 to be executed by the processor 10 and data (generator /
14 (abbreviated as load data).

The shared memory 16 retains data (abbreviated as power transmission network calculation data) 17 necessary for power transmission network calculation.

FIG. 2 shows the overall processing flow when the transient stability calculation of the power system calculation is executed by the Runge-Kutta method in the power system calculation device 1 of the present embodiment. ing. Referring to FIG.
This processing is basically performed in each time step in the calculation step 21 relating to the generator / load and the transmission network calculation step 2.
Step 2 is repeated four times.

The transient stability calculations described in the power system computing device 1 of the present embodiment is basically a M _g-number of the generator, with respect to the M _l number of loads, simultaneous ordinary differential including transmission network calculation This is the problem of solving equations.

[0022] The this is a method for solving in parallel, if there is n (a positive integer) base processor 10, and M _g-number generator, and M _l number of loads, equally assigned to each processor 10.

The transmission network calculation step 22 is assigned to each processor 10 for each row.

In general, for one generator, the input is a voltage

[Outside 1] _g and current

[Outside 2] _g , a simultaneous ordinary differential equation exists, and as a result, the internal induced voltage of the generator

[Outside 3] _g and internal impedance 1 /

[Outside 4] _g is obtained.

On the other hand, for one load, the terminal voltage

[Outside 5] Admittance with ₁ as input

[Outside 6] ₁ is obtained.

Since these generators and loads are connected to each other through a power transmission network, the relationship expressed by the following equation is established.

[0027]

(Equation 1)

Here,

[Outside 7] Is the terminal voltage, for generators

[Outside 8] _g for load

[Outside 9] _Is one. Also,

[Outside 10] Is the grid admittance matrix. Solving Equation 1 gives the terminal voltage

[Outside 11] Is obtained.

The internal induced voltage of the generator

[Outside 12] _g and the terminal voltage of the generator

[Outside 13] _g and current

[Outside 14] _Since the relationship expressed by Equation 2 holds between _g and _g , the generator current

[Outside 15] _g is obtained.

[0030]

(Equation 2)

FIG. 3 is a flowchart for explaining in more detail the calculation step 21 relating to the generator / load in FIG. Each processor 10 performs the generator operation process 30 by the number of assigned generators.

The following ordinary differential equation is established for the operation step 30 of the generator (see Yasuji Sekine: Transient analysis of power system, Ohmsha (1984), pp. 377).

Equations (3) and (4) hold for the equation of rotor motion.

[0034]

(Equation 3)

[0035]

(Equation 4)

Equations (5), (6), (7) and (8) hold for the change equation of the linkage flux.

[0037]

(Equation 5)

[0038]

(Equation 6)

[0039]

(Equation 7)

[0040]

(Equation 8)

With respect to the characteristic formula of the automatic voltage regulator, Equation 9
Holds.

[0042]

(Equation 9)

Equation (10) holds for the governor characteristic equation.

[0044]

(Equation 10)

Each processor 10 calculates the terminal voltage of the generator obtained from the shared memory 16 in the immediately preceding transmission network calculation step 22.

[Outside 16] _g , the current of the generator

[Outside 17] _g and the internal variables in the generator are calculated based on the Runge-Kutta method, and finally the admittance

[Outside 18] _g and generator internal induced voltage

[Outside 19] The coefficient value on the left side and the value on the right side of Equation 1 calculated from _g are written in the shared memory 16.

Thereafter, each processor 10 executes the load calculation step 31 by the number of assigned loads.

Regarding the load calculating step 31, the following equation (11) is used.
Equations (12) and (13) hold (Masato Yamamoto: Dynamic Load Model for System Voltage Stability Analysis, Transactions of the Institute of Electrical Engineers of Japan, B, Vol. 114-B, No.
12, pp. 1273-1279 (1994)).

[0048]

[Equation 11]

[0049]

(Equation 12)

[0050]

(Equation 13)

Each processor 10 reads the voltage of the load obtained from the shared memory 16 in the immediately preceding transmission network calculation step 22.

[Outside 20] _{Based on 1} , the value of the internal variable in the load is calculated based on the Runge-Kutta method, and finally the admittance is calculated.

[Outside 21] Writing the number 1 of the left side of coefficients calculated from ₁ to the shared memory 16.

The execution time required for the operation step 34 based on the Runge-Kutta method is the same as that of the shared memory read step 32.
And the time required for the shared memory writing step 33 is much longer.

In the operation process 34 based on the Runge-Kutta method, each processor 10 can execute the operation completely independently.

FIG. 4 is a flowchart for explaining the transmission network calculation step 22 in FIG. 2 in more detail.

In the transmission network calculation step 22, the equation 1 is obtained by calculating the internal induced voltage of the generator obtained in the calculation step 34 based on the Runge-Kutta method.

[Outside 22] _g , generator admittance

[Outside 23] _g and load admittance

[Outside 24] Solve based on ₁ . That is, from equation 1, the terminal voltage

[Outside 25] Ask for.

The solution of the simultaneous equations is occupied by an LU decomposition / forward elimination step 40 and a backward substitution step 41. Among these processes, the LU decomposition process has the longest processing time. Therefore, the parallel execution of LU decomposition will be described in more detail with reference to FIG.

In general, LU decomposition is performed in the following steps. In the following, an n * n matrix to be subjected to LU decomposition is defined as (n is the size of the matrix).

[Outside 26] And

[Outside 27] Let the element at the i-th row and j-th column be represented by a _{i, j} . At this time, the processing of the outer product form of the LU decomposition is represented as follows (Hario Tsuda: Numerical Processing Programming, pp. 1)
21-123, Iwanami Shoten).

The following is repeated for the k-th row (k = 1 to n). -A _{k, k} = 1 / a _{k, k} -For the i-th row (i = k + 1 to n), the following is repeated. * A _{i, k} = a _{i, k} * _{ak, k} * The following is repeated for the j-th column (j = k + 1 to n).・ A _{i, j} = a _{i, j} −a _{i, k} * a _{K, J}

In the LU decomposition, each processor 10 is in charge of processing for each row. The process itself is performed in a process (50) in which the k-th row to be decomposed is repeated from k = 1 to n.

In this process, the processor 10 in charge of the k-th row notifies other processors 1 that the row has been reached.
Notify 0. Other processors 10 wait for processing until a notification is received.

Thereafter, each processor 10 executes the (k +
1) With respect to the i-th row assigned to the processor 10 out of the rows from the n-th row, a new value of each element of the assigned row is calculated based on the value of the k-th row (step 51). That is, after updating the value of the element a _{i, k} in the i-th row and the k-th column, the calculation of each element in the j-th column from the k-th column to the n-th column in the i-th row is performed (step 52). _Find the value of a _{i, j} . In the actual calculation, since the matrix is a canonical matrix, the operation is performed only on elements having non-zero values.

In the LU factorization / forward erasure step 40, the values of each element a _{i, j} of the Jacobian matrix are all arranged on the shared memory 16. In this part, the inter-processor communication that occurs when reading the value of the pivot row is performed by the shared memory 16.
It is acting by using.

In the LU disassembly / forward erasure step 40 in which many processes share data, the superiority of the bus-coupled shared memory type parallel computer system is secured.

FIG. 6 shows the flow of the entire processing when the transient stability calculation in the power system calculation device 1 of the present embodiment is executed by the trapezoidal rule.

Basically, in each time step, the processing 6 of the Newton-Raphson method is performed until the convergence condition is satisfied.
0, that is, a residual calculation step 61 relating to a state equation of a generator, a load, and a transmission network; a Jacobian matrix element calculation step 62 of a correction equation for obtaining a correction amount of a state variable of the generator, load, and the transmission network; , Load and transmission network correction equation coefficient matrix (hereinafter, Jacobi matrix), LU decomposition, forward elimination and backward substitution step 63 are repeated.

Here, the solution by the trapezoidal rule will be described. Ordinary differential equations for the generator (Equations 3 to 10)
And an equation obtained by time-differentiating the ordinary differential equation (Equation 11) regarding the load by using the trapezoidal rule. When this equation, the algebraic equation relating to the generator and the load (equation 2 and a similar equation omitted from the description), and the algebraic equation relating to the power transmission network are simultaneously solved, this equation is represented by equation 14.

[0067]

[Equation 14]

Here,

[Outside 28] Indicates all the state variables of the generator and load and the grid.

The simultaneous equations in Equation 14 are nonlinear equations, and when the Newton-Raphson method is used as a numerical solution, the correction equation becomes Equation 15.

[0070]

(Equation 15)

Here, 1, 2,..., I,..., N are identification numbers assigned to each equation and each state variable of the simultaneous equations, and the matrix on the left side is a coefficient matrix of the simultaneous equations, which is called a Jacobian matrix. is there. {Δx _i } indicates a correction amount vector to be added to each state variable provisional solution, and {r _i } indicates a residual vector when each state variable provisional solution is substituted into each equation.

FIG. 7 shows the processing 60 of the Newton-Raphson method of FIG. 6 in more detail.

Each processor 10 is mechanically assigned a generator and a load. As for the power grid,
The connection terminals of the transmission network to the generators or loads are mechanically allocated to the processors 10 to which the generators and loads are assigned, by the number of nodes.

FIG. 8 shows the structure of the Jacobian matrix. Blocks in which non-zero matrix elements are present in blocks surrounded by dotted lines are hatched, and blocks having all matrix elements are zero. Is not hatched. However, even in the hatched blocks, not all matrix elements are non-zero.

Blocks that are not hatched do not need to be operated in all cases. The calculation of the Newton-Raphson method uses the structural features shown in FIG. 8 to arrange data in a memory so that the time required for reading and writing data described below is reduced. That is, the generator / load data 14 is stored in the local memory 12, and the transmission network calculation data 17 is stored in the shared memory 16.

Each processor 10 has a terminal voltage of the generator with respect to the generator and load in charge.

[Outside 29] _g and load terminal voltage

[Outside 30] After reading ₁ from the shared memory 16 (step 32), using the local memory 12 of each processor 10, the residual calculation step 61 relating to the state equation of the generator, load and transmission network in charge, and the generator, load and A Jacobi matrix element calculation step 62 of a correction equation relating to the power transmission network and an LU decomposition and forward elimination calculation step 71 of the Jacobi matrix relating to the generator / load are performed. And the resulting generator
A step 33 of storing update data of the Jacobian matrix element at the point of contact of the load with the power transmission network in the shared memory 16 is performed.

It should be noted that a residual calculation step 61 relating to the state equation of the generator, load and transmission network in charge, a Jacobian matrix element calculation step 62 for correction equations relating to the generator, load and transmission network, and a Jacobian matrix relating to the generator and load The execution time required for the LU decomposition and forward erasure calculation step 71 is much longer than the time required for the shared memory read step 32 and the shared memory write step 33.

Further, a residual calculation step 61 relating to a state equation of a generator, a load and a transmission network in charge, a Jacobian matrix element calculation step 62 of a correction equation relating to a generator, a load and a transmission network, and a Jacobian matrix relating to a generator / load In the LU decomposition and forward elimination calculation step 71, each processor 10 can execute an operation completely independently.

Thereafter, each processor 10 refers to the shared memory 16 and frequently performs the LU decomposition and forward elimination calculation step 72 of the Jacobian matrix relating to the power transmission network (while communicating by updating the data in the shared memory 16). And a backward substitution calculation step 73 for the power transmission network.

Steps 72 and 73 are described in LU
Processing is performed in the same manner as in the disassembling / forward erasing step 40 and the backward substitution step 41.

Then, after reading the amount of correction of the state variable of the power grid in the shared memory 16, each processor 10 independently performs the backward substitution step 74 regarding the generator and load in the local memory 12.

Thereafter, each processor 10 executes the correction step 75 with reference to the local memory 12 for all of the state variables assigned to each processor 10.

As a result of the backward substitution regarding the power transmission network, the correction amount of the provisional solution of the state variable regarding the power transmission network is obtained and stored in the shared memory 16.

Subsequently, each processor 10 reads, from the shared memory 16, the generator and the amount of correction of the provisional solution of the state variable at the point of contact between the load and the transmission network, which are shared by the respective processors 10, and uses them to generate the generator Then, a backward substitution calculation step 74 relating to the load and the load is performed to obtain a correction amount of the provisional solution of the state variable relating to the generator and the load, and the correction amount is stored in each local memory 12.

Thereafter, each processor 10 executes a provisional solution of the state variables relating to the generator and load and the grid.

[Outside 31] 'To the correction amount Δ

[Outside 32] The state variable solution is updated in a correction calculation step 75 according to the equation (16).

[0086]

(Equation 16)

This series of operations is repeated until the correction amount of all state variables becomes sufficiently small, and the provisional solution at the end of the repetition is regarded as a solution.

[0088]

As described above, according to the present invention,
As for the transmission network calculation process, there is an effect that the superiority of the parallel computer system of the bus-coupled shared memory system is secured.

On the other hand, in the process of calculating the simultaneous ordinary differential equations relating to the internal state of the generator and the load, each processor performs an operation using the local memory of the processor, so that it is not necessary to refer to the shared memory on the bus. Therefore, there is an advantage that performance degradation due to bus reference conflict does not occur.

Further, in both the transmission network calculation step and the generator / load calculation step, each processor reads the instruction of the processor from the local memory, so that there is an advantage that the performance of reading the instruction does not decrease.

[Brief description of the drawings]

FIG. 1 is a block diagram illustrating a configuration of a power system calculation device according to an embodiment of the present invention.

FIG. 2 is a diagram showing a runge and a power system in the power system calculation apparatus according to the embodiment;
4 is a flowchart illustrating an example of a power system calculation based on the Kutta method.

FIG. 3 is a flowchart illustrating in more detail a process of calculating a simultaneous ordinary differential equation relating to the internal state of the generator / load in FIG. 2;

FIG. 4 is a flowchart illustrating a transmission network calculation step in FIG. 2 in more detail;

FIG. 5 is a flowchart illustrating parallel execution of an LU disassembly / forward erasure process in FIG. 4;

FIG. 6 is a flowchart illustrating an example of a power system calculation based on the trapezoidal rule in the power system calculation device of the present embodiment.

FIG. 7 is a flowchart for explaining the processing of the Newton-Raphson method in FIG. 6 in more detail;

FIG. 8 is a diagram showing a structure of a Jacobi matrix.

[Explanation of symbols]

REFERENCE SIGNS LIST 1 power system calculation device 10 processor 11 central processing unit 12 local memory 13 program 14 generator / load data 15 bus 16 shared memory 17 data for power grid calculation 21 calculation process for generator / load 22 power grid calculation process 30 generator Calculation step 31 Load calculation step 32 Shared memory reading step 33 Shared memory writing step 34 Calculation step based on Runge-Kutta method 40 LU disassembly / forward erasing step 41 Backward substitution step 42 Forward erasing step 50 Pivot row (k) from 1 Steps to be repeated up to n 51 Steps of updating the element of the assigned row (i) from the (k + 1) th row to the nth row 52 Updating the elements of the (k + 1) th column to the nth column of the ith row Process 60 Processing process of Newton-Raphson method 61 For generator, load and transmission network state variables Calculating the residuals of the equations to be performed 62 generating the Jacobian matrix elements of the generator, load and transmission network 63 LU decomposition, forward elimination and backward substitution of the Jacobian matrix of the generator, load and transmission network 64 correcting the state variables and time Update Step 71 LU Decomposition of Jacobian Matrix and Forward Elimination Calculation Step for Generator / Load 72 LU Decomposition of Jacobian Matrix and Forward Elimination Calculation Step for Transmission Network 73 Backward Substitution Calculation Step for Transmission Network 74 Backward Substitution Calculation Step for Generator / Load 75 Correction process of all state variables

Continuation of the front page (56) References JP-A-58-33933 (JP, A) JP-A-6-261453 (JP, A) JP-A-3-256525 (JP, A) (58) Fields investigated (Int) .Cl. ⁶ , DB name) H02J 3/00-5/00

Claims (3)

(57) [Claims] 1. An analytical calculation for a power system composed of a plurality of generators and a plurality of loads and a transmission network interconnecting the plurality of generators, and a simultaneous ordinary differential equation relating to the internal state of the generator and the load for each time step. A power system calculation device that calculates the transmission network calculated between them by repeatedly solving simultaneous equations; a bus that enables connection of a plurality of hardware; and a connection to the bus. A central processing unit that performs calculations,
A plurality of processors, each of which includes a program executed by the central processing unit and a local memory, and retains data relating to the internal state of the generator and the load assigned to the processor among data relating to the internal state of the generator and the load; And a shared memory connected to the bus and holding data necessary for the transmission network calculation.In the step of calculating the simultaneous ordinary differential equations relating to the internal state of the generator and the load, each of the processors ,
The generator and load assigned to the processor.
Power supply connection terminal voltage data from the shared memory
After reading, the process that exists in the local memory
Internal state of the generator and load assigned to the generator
Calculated independently using data on
Generator induced internal voltage, generator admittance and
Network and load admittance to the shared memory
In the step of transferring as data necessary for calculation, and calculating the power grid calculation by solving simultaneous equations,
Each of the processors is configured to operate the generator on the shared memory.
Internal induced voltage, generator admittance and load
The processor while referring to and updating dommitance
A power system calculation device for executing a transmission network calculation assigned to a power system. 2. In the step of calculating the simultaneous ordinary differential equations relating to the internal state of the generator and the load, each of the processors converts the connection terminal voltage data of the transmission network relating to the generator and the load assigned to the processor to the processor. After reading from the shared memory, it is independently calculated based on the Runge-Kutta method using data on the internal state of the generator and load assigned to the processor, which is present in the local memory, and the resulting internal Transferring the data of the induced voltage and the admittance to the shared memory, in the power grid calculation step, each of the processors,
The power system calculation device according to claim 1, wherein the power system calculation device executes the power network calculation assigned to the processor while referring to and updating data necessary for the power network calculation on the shared memory. 3. A step of calculating a simultaneous ordinary differential equation relating to the internal state of the generator and the load, wherein the equation is a correction equation of a variable indicating the internal state of the generator and the load based on a trapezoidal rule; In the case where the state variable correction equations are simultaneously obtained and the state variable solutions are obtained by iterative correction using the Newton-Raphson method, the connection terminal voltage data of the power grid regarding the generator and the load are read from the shared memory. Later, the residuals of the equations for the internal state variables of the generator and the load and the elements of the modified equation coefficient matrix (Jacobi matrix) are independently calculated using the data on the internal state of the generator and the load assigned to the processor. In the local memory, and the residuals of the equations relating to the state variables of the power grid and the modified equation coefficient matrix (Yako Matrix) is calculated using the data required for the power grid calculation assigned to the processor using the data required for the power grid calculation stored in the local memory and the shared memory, and the result is calculated as In the correction equation calculation step, for the correction equation coefficient matrix relating to the generator and the load, the LU decomposition and forward elimination of the elements allocated to the processor are performed by the residual and Jacobian of the local memory. When it is necessary to independently execute the matrix while referring to / updating the elements of the matrix, and further to update the data of the portion relating to the power transmission network with the calculation of the generator and the load by the LU decomposition and forward elimination, Is updated, and the corrected equation coefficient matrix and the residual of the power grid are assigned to the processor. The respective processors execute the LU factorization and forward elimination operation of the assigned elements while referring to / updating data on the shared memory, and each processor determines the amount of correction related to the state variable of the power grid on the shared memory. Is determined by performing a backward substitution calculation on the residual subjected to forward elimination using the data of the above, the amount of correction of the state variables related to the generator and load, the respective processors, from the shared memory the state variable of the power grid After reading the correction amount data of, the Jacobi matrix on the local memory and the residual data subjected to forward elimination and back substitution are determined by back substitution calculation, and all state variables are corrected. The power system calculation device according to claim 1, wherein the calculation is repeatedly performed until the correction amount of the power system becomes small.