CN109782625B - Real-time simulation method and system for circuit model - Google Patents
Real-time simulation method and system for circuit model Download PDFInfo
- Publication number
- CN109782625B CN109782625B CN201811562834.6A CN201811562834A CN109782625B CN 109782625 B CN109782625 B CN 109782625B CN 201811562834 A CN201811562834 A CN 201811562834A CN 109782625 B CN109782625 B CN 109782625B
- Authority
- CN
- China
- Prior art keywords
- sub
- circuit
- simulation
- decoupling
- real
- 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.)
- Active
Links
- 238000004088 simulation Methods 0.000 title claims abstract description 106
- 238000000034 method Methods 0.000 title claims abstract description 38
- 238000004364 calculation method Methods 0.000 claims abstract description 42
- 238000009826 distribution Methods 0.000 claims abstract description 12
- 230000005540 biological transmission Effects 0.000 claims description 38
- 239000011159 matrix material Substances 0.000 claims description 33
- 238000004422 calculation algorithm Methods 0.000 claims description 25
- 230000010354 integration Effects 0.000 claims description 12
- 238000012545 processing Methods 0.000 claims description 8
- 230000008054 signal transmission Effects 0.000 claims description 4
- 238000005070 sampling Methods 0.000 claims description 3
- 230000001939 inductive effect Effects 0.000 claims description 2
- 238000010586 diagram Methods 0.000 description 17
- 238000004590 computer program Methods 0.000 description 7
- 230000000694 effects Effects 0.000 description 6
- 230000006870 function Effects 0.000 description 5
- 230000008859 change Effects 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 238000000926 separation method Methods 0.000 description 3
- 238000003860 storage Methods 0.000 description 3
- 238000010521 absorption reaction Methods 0.000 description 2
- 230000002457 bidirectional effect Effects 0.000 description 2
- 230000033228 biological regulation Effects 0.000 description 2
- 238000011217 control strategy Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000005457 optimization Methods 0.000 description 2
- 230000010355 oscillation Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000003990 capacitor Substances 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000005520 cutting process Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000001914 filtration Methods 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
- 230000006798 recombination Effects 0.000 description 1
- 238000005215 recombination Methods 0.000 description 1
- 230000006641 stabilisation Effects 0.000 description 1
- 238000011105 stabilization Methods 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
The application provides a real-time simulation method and a system of a circuit model, wherein the method comprises the following steps: decoupling the whole system into a plurality of sub-networks according to the distribution condition of switching value in the circuit system; based on a set state space equation, calculating each sub-network by adopting a parallel calculation mode to obtain an equivalent circuit corresponding to each sub-network; based on the equivalent circuits of all the subnetworks, the real-time simulation of the circuit system is realized by joint simulation, after the circuit is decoupled, the state space equation of the circuit is divided into a plurality of small state space equations, so that the use of memory space and calculation amount are reduced, and when the state of a circuit switch changes, the state space equation is not required to be solved again, thereby greatly accelerating the simulation speed.
Description
Technical Field
The application relates to the field of direct current converters, in particular to a real-time simulation method and system of a circuit model.
Background
The actual circuit is usually formed by connecting electrical components according to a circuit principle, and the problem solving through simulation becomes a necessary approach, but some circuit models have complicated connection relations, and real-time simulation is not easy to realize, for example, a direct current converter is introduced in the circuit for solving the system efficiency, and the direct current converter is a power electronic device for converting direct current electric energy into direct current electric energy with controllable voltage or current required by a load. The method comprises the steps of cutting constant direct current voltage into a series of pulse voltages through quick on-off control of a power electronic device, changing the pulse width of the pulse series through control of duty ratio change so as to realize the regulation of the average value of output voltage, filtering by an output filter, and obtaining direct current electric energy with controllable current or voltage on a controlled load.
The full-bridge direct-current converter is a direct-current converter, is suitable for the field of high-voltage high-power direct-current transmission, can fully utilize the variable-voltage parameters, can ensure the working state of a soft switch in a full-load range, has constant-voltage regulation capability, and has higher efficiency. Along with the wide development of the high-voltage direct-current transmission technology, the direct-current converter is the key point in the research of the high-voltage direct-current transmission technology, because the direct-current converter is used as a key link of a power electronic transformer, is the key point for realizing the stable power balance of the whole network voltage and is difficult to control strategies, when the system is connected with a grid, the whole system can only be used as a state space equation to calculate the system due to the fact that the switching value is too large, the method needs larger memory space and calculated amount, real-time simulation of a model is difficult to realize, and inconvenience is brought to the research of the direct-current converter.
Disclosure of Invention
In order to solve the defects that the whole system can only be used as a state space equation to calculate the system due to overlarge switching value in the prior art, larger memory space and calculation amount are needed, and real-time simulation of a model is difficult to realize, the application provides a simulation method and a simulation system of a circuit model. The method adopts a split-core decoupling mode to carry out switching value recombination aiming at the topological structures of the full-bridge direct-current converter under different working conditions in the model, and performs split-core calculation, so that the problem that the model with too large switching value cannot be real-time is solved, and the calculation efficiency of the simulation model can be increased.
The technical scheme provided by the application is as follows: a real-time simulation method of a circuit model, comprising:
decoupling the whole system into a plurality of sub-networks according to the distribution condition of switching value in the circuit system;
based on a set state space equation, calculating each sub-network by adopting a parallel calculation mode to obtain an equivalent circuit corresponding to each sub-network;
and carrying out joint simulation based on the equivalent circuits of all the subnetworks to realize real-time simulation of the circuit system.
Preferably, the calculating each sub-network by adopting a parallel calculation mode based on the set state space equation to obtain an equivalent circuit corresponding to each sub-network includes:
the operation parameters of each sub-network are put into a set state space equation to obtain an equivalent circuit equation corresponding to each sub-network;
and carrying out parallel calculation on the equivalent circuit equation by adopting a high-order integration algorithm based on a multiprocessor to obtain an equivalent circuit of each sub-network.
Preferably, the state space equation is as follows:
y＝C n x+D n u
wherein:inputting a vector; y: outputting a vector; x: state variables in the circuit model; u: a current or voltage variable; a is that n : a matrix A; b (B) n : a matrix B; c (C) n : a matrix C; d (D) n : a matrix D; n: the order of the coefficient matrix.
Preferably, the equivalent circuit equation is as follows:
y t+Δt ＝y h +W n u t+Δt
wherein y is t+Δt : an equivalent circuit equation; y is h : a state space group history output; w (W) n : an admittance matrix; u (u) t+Δt : the voltage variable at the next moment; Δt: step length of discrete sampling;
wherein the state space group history output y h Admittance matrix W n Calculated as follows:
wherein:a discretized matrix A; x is x t : the current time value of the intermediate variable; />A discretized matrix B; u (u) t : the voltage variable at the current moment.
Preferably, the decoupling the whole system into a plurality of sub-networks includes:
processing the circuit system according to a set decoupling mode;
and decoupling the processed circuit system into a plurality of sub-networks according to the distribution condition of the switching value.
Preferably, the decoupling mode includes: a first decoupling mode and a second decoupling mode;
the first decoupling mode is to separate the switches by adopting a series direct-current voltage source and a parallel controllable alternating-current source;
the second decoupling mode is to make two IGBTs connected in series equivalent to an S function algorithm, so that the number of switches is reduced.
Preferably, the performing the joint simulation based on the equivalent circuit of each sub-network to realize the real-time simulation of the circuit system includes:
summarizing each sub-network into a plurality of sub-systems based on the topology of the circuitry;
in a simulation time step, based on the equivalent circuit of each sub-network in each sub-system, carrying out joint simulation by adopting a post Euler discrete integration algorithm to obtain an electric transmission signal of the sub-system;
and the subsystems transmit electrical transmission signals based on the transmission lines, so that the real-time simulation of the circuit system is realized. Preferably, the subsystem transmits an electrical transmission signal based on a transmission line, so as to realize real-time simulation of the circuit system, and the method comprises the following steps:
and after delaying the signal transmission among the subsystems by one simulation step, simulating based on the received electrical signals until the simulation is completed to obtain the real-time simulation of the whole circuit.
Preferably, the subsystem transmits an electrical transmission signal based on a transmission line, and after implementing the real-time simulation of the circuit system, the method further includes:
and verifying the topological structure of the circuit model and simulating in real time through an external controller.
Based on the same inventive concept, the application also provides a real-time simulation system of the circuit model, comprising:
the decoupling module is used for decoupling the whole system into a plurality of sub-networks according to the distribution condition of the switching value in the circuit system;
the equivalent circuit obtaining module is used for respectively calculating all the sub-networks in a parallel calculation mode based on a set state space equation to obtain equivalent circuits corresponding to all the sub-networks;
and the circuit simulation module is used for carrying out joint simulation based on the equivalent circuits of all the subnetworks to realize real-time simulation of the circuit system.
Preferably, the obtaining an equivalent circuit module includes:
the equivalent circuit equation unit is used for bringing the operation parameters of each sub-network into a set state space equation to obtain an equivalent circuit equation corresponding to each sub-network;
and the equivalent circuit unit is used for carrying out parallel calculation on the equivalent circuit equation by adopting a high-order integration algorithm based on the multiprocessor to obtain an equivalent circuit of each sub-network.
Compared with the closest prior art, the technical scheme provided by the application has the following beneficial effects:
1. according to the technical scheme provided by the application, aiming at the distribution condition of switching value in a circuit system, the whole system is decoupled into a plurality of sub-networks; based on a set state space equation, calculating each sub-network by adopting a parallel calculation mode to obtain an equivalent circuit corresponding to each sub-network; based on the equivalent circuits of all the subnetworks, the joint simulation is carried out to realize the real-time simulation of the circuit system, after the circuits are decoupled, the state space equation of the whole circuit is divided into a plurality of small state space equations, so that the use of memory space and calculation amount are reduced, and when the state of a circuit switch changes, the state space equation is not required to be solved again, thereby greatly accelerating the simulation speed.
2. The split-core decoupling algorithm provided by the application can solve the problem of real-time model of a large-scale full-bridge direct-current converter, can be used for decoupling any large-scale power electronic switching device, solves the problem that real-time simulation of the model cannot be realized due to too large switching value of the model, and improves the simulation efficiency.
3. According to the technical scheme provided by the application, the full-bridge direct-current converter model can verify the rationality of the topological structure of the model and the effect of the controller through the external controller in real time, and a foundation is laid for the application of the follow-up theory to practice.
4. By adopting the technical scheme provided by the application, the simulation calculation of two parts is realized in the same simulation step length, and in the initial stage of calculation, the calculation among all the subnetworks is mutually independent, and the calculation efficiency is improved by parallel calculation; and in the later stage of calculation, the node network of the whole circuit is combined with simulation so as to avoid artificial delay caused by decoupling and improve the simulation precision and numerical stability.
5. According to the technical scheme provided by the application, the circuit is decoupled according to the switching value to generate a plurality of sub-networks, and each sub-network can be distributed to a plurality of cores for parallel calculation, so that the simulation speed of the model is greatly increased.
6. According to the technical scheme provided by the application, different discrete integral algorithms can be adopted for solving in two stages, a high-order integral algorithm can be used in the initial calculation stage, and when the equivalent circuits of all the sub-networks are combined for solving, a discrete integral algorithm such as a rear Euler algorithm can be used for solving, so that the numerical oscillation caused by abrupt change of the input variables of the system can be effectively eliminated by utilizing interpolation and iteration algorithms commonly used in classical node method solving.
7. The technical scheme provided by the application is used for the situation that a large number of switch power electronic devices in the model are difficult to decouple in real time and the verification of a model control strategy is influenced.
Drawings
FIG. 1 is a flow chart of a simulation method provided by the application;
FIG. 2 is a schematic diagram of a first decoupling mode according to the present application;
FIG. 3 is a schematic diagram of a second decoupling mode according to the present application;
fig. 4 is a schematic diagram of a circuit decoupling module according to the present application;
FIG. 5 is a schematic diagram of a circuit model according to an embodiment of the present application;
FIG. 6 is a topology diagram of a bi-directional full-bridge DC converter according to an embodiment of the present application;
FIG. 7 is a diagram showing the effect of converting a bi-directional full-bridge DC converter according to an embodiment of the present application;
FIG. 8 is a topology diagram of a unidirectional full-bridge DC converter in an embodiment of the application;
fig. 9 is an effect diagram of the unidirectional full-bridge dc converter converted according to the embodiment of the present application;
fig. 10 is an effect diagram of the split core decoupling method according to the embodiment of the present application.
Detailed Description
For a better understanding of the present application, reference is made to the following description, drawings and examples.
Example 1
FIG. 1 is a flow chart of a real-time simulation method of a circuit model provided by the application, which comprises the following steps:
step S1, decoupling the whole system into a plurality of sub-networks according to the distribution condition of switching value in a circuit system;
step S2, based on a set state space equation, calculating each sub-network in a parallel calculation mode to obtain an equivalent circuit corresponding to each sub-network;
and step S3, carrying out joint simulation based on the equivalent circuits of all the sub-networks to realize real-time simulation of the circuit system.
The basic principle of the split-core decoupling algorithm provided by the application is that a circuit model is divided into a plurality of subsystems, each subsystem is decoupled into a plurality of sub-networks according to the distribution condition of switching value, each sub-network is respectively modeled and calculated by using a state space method based on MATLAB software, and Dai Weining and Norton equivalent circuits of each sub-network are solved; then, establishing a solved network node equation between all sub-networks by using a classical node method; and finally, carrying out joint simulation on the equivalent circuits of all the subsystems to obtain the calculation result of the whole circuit.
The decoupling algorithm realizes the simulation calculation of two parts in the same simulation step length, the calculation among all the subnetworks is independent in the initial stage of calculation, and a plurality of processor cores can be adopted for parallel calculation, namely core-division independent calculation, according to the configuration condition of the simulator, so that the calculation efficiency is improved. And in the later calculation stage, the integrated simulation in the whole subsystem is completed in the same simulation step length, so that the delay caused by decoupling is avoided, the final simulation result is influenced, and the simulation precision and the numerical stability are improved.
Step S2, based on a set state space equation, calculating each sub-network by adopting a parallel calculation mode to obtain an equivalent circuit corresponding to each sub-network, wherein the step comprises the following steps:
each subnetwork can be represented using state space equation (1):
wherein: x is the state variable of the system (the number of state quantities is determined by the conditions of the capacitor, the inductance and the transformer in the network topology), u is the input vector (i.e. the current or voltage variable), A n ，B n ，C n ，D n Is the coefficientA matrix. n represents the order of the coefficient matrix.
Discretizing the formula (1) to obtain:
wherein: Δt is the step size of the discrete samples,and->Is discretized by matrix A and matrix B. u (u) t+Δt U is the voltage variable at the next moment t Is the voltage variable at the current moment. X is x t For the current time value of the intermediate variable, x t+Δt Is the next time value of the intermediate variable.
The formula (2) can be further simplified and combined to obtain:
i.e. y t+Δt Consists of variables based on the known result of the previous step and the variables to be solved for the current step.
Formula (3) can be expressed as:
y t+Δt ＝y h +W n u t+Δt (4)
in the method, in the process of the application,W n as admittance matrix, u t+Δt And y is t+Δt Forming a classical Norton equivalent circuit equation; y is h : state space group history output.
In the two resolving stages of the network decoupling method of the embodiment, different discrete integration algorithms can be adopted for solving, a high-order integration algorithm can be used in parallel computation, and a post-Euler discrete integration algorithm can be used in joint solution of equivalent circuits of all sub-networks, so that interpolation and iteration algorithms commonly used in classical node method resolving can be utilized, and numerical oscillation caused by abrupt change of input variables of a system can be effectively eliminated.
As shown in fig. 2 and fig. 3, according to the decoupling principle, two modes suitable for decoupling the model are set up in this embodiment, and fig. 2 is a mode of adopting a series direct-current voltage source and a parallel controllable alternating-current source to realize separation of each switch, so as to reduce the scale of a state space equation and achieve the decoupling effect. In the fig. 3, the two IGBTs which are originally connected in series are expressed in the form of an sfunction algorithm in an equivalent way, so that the number of switches is reduced, and the state space equation is simpler, thereby achieving the decoupling effect.
Step S3, based on the equivalent circuit of each sub-network, carrying out joint simulation to realize real-time simulation of the circuit system, wherein the step comprises the following steps:
the decoupling principle of the sub-network only decomposes the original state space equation into several state space equations, but the signal transmission between the sub-systems needs to realize the transmission of the electrical signal by means of the form of line decoupling.
As shown in fig. 4, the method is implemented by adopting simulation N distributed line parameters, and delaying the electrical characteristic of line transmission by one simulation step. Because the transmission of the electrical signals at the two ends of the line, i.e. the signals at the input end are transmitted to the output end after passing through the line, the process only has a delay of a simulation step length, and the electrical characteristics are not changed.
Example 2
Taking the topological diagram of the whole circuit model structure in fig. 5 as an example, the embodiment of the application is used for explaining the simulation method provided by the application, the system is composed of 7 direct current converters, the whole voltage stabilization and the power balance control of the system after the grid connection of the full-bridge direct current converters are mainly realized, 3 of the simulation method are unidirectional direct current converter models (P3, P4 and P5), and the fixed absorption power is 5MW by adopting the power control; 2 unidirectional DC converter models (P1, P2) are arranged, so that power control and voltage control can be realized, and the fixed power is 5MW;1 is a bidirectional direct current converter model (P6), and fixed absorption or power emission of 5MW can be realized according to different parallel current source parameters; and 1 unidirectional direct current converter model (P7) which can be used for power control and voltage control, wherein the power control sends out 5MW, and the voltage control depends on specific working conditions. In addition, there is 1 ac balancing node (P8). The circuit model totally comprises 122 power electronic switching devices, if a split-core decoupling mode is not adopted, the model cannot run in a real-time simulation machine, the calculation result of the model is slow, the overtime condition occurs, and the simulation working condition of an external controller cannot be limited.
Firstly, a circuit model needs to be subjected to kernel division optimization before decoupling, and in the embodiment, the circuit model is divided into three subsystems according to a network topology structure of the model, wherein the first subsystem is P1, P2 and P8, the second subsystem is P5 and P7, and the third subsystem is P3, P4 and P6; the transmission between the cores adopts a line decoupling mode, the core-splitting mode has reasonable resource occupancy rate, and the overtime phenomenon can not be caused by the transmission between the cores.
Then, the topological structure of the full-bridge direct current converter aiming at different working conditions is divided into a plurality of sub-networks in each sub-system, and in-core decoupling analysis is carried out in each sub-network. P6 in FIG. 5 is a bidirectional DC converter model, the structure totally comprises 2*8 IGBT switches, 2*7 switch states exist, the structure is composed of two topological structures in FIG. 6, if a split core decoupling mode is not adopted, the simulation calculation is performed according to a state space matrix of network topology, the method is extended into the whole network, a state space equation is huge, the resource occupancy rate is high, the simulation calculation speed is slow, and real-time simulation of the model cannot be realized. And taking P6 as a sub-network to perform decoupling processing.
As shown in fig. 7, in order to reduce simulation time and complete real-time simulation of the model and ensure calculation accuracy, the topology in fig. 6 is implemented in a form of a decoupling module 1 to implement conversion of the model topology. According to the method, the original 2*8 IGBT switches are separated by a series voltage source or a parallel controllable alternating current source, so that the whole topological state space equation is divided into a plurality of small state space equations, the state space equations are fixed when the model is repeatedly calculated each time, and the state space equations are not required to be solved additionally.
For the unidirectional full-bridge dc converter shown in fig. 8, decoupling needs to be performed by combining the two decoupling modes of fig. 2 and 3, and a specific topological diagram after decoupling is shown in fig. 9. The switch close to the parallel network point is decoupled in a second decoupling mode, the other end of the switch is decoupled in a first decoupling mode, and only 8 topological decoupling actual calculation of the original 16 switches are performed after decoupling, so that the state space equation is simpler.
As shown in fig. 10, after the decoupling of the circuit model is completed, the model needs to be further subjected to core separation optimization, and transmission among three subsystems adopts a line decoupling mode, so that the core separation mode has reasonable resource occupancy rate and the transmission among cores cannot cause a timeout phenomenon.
Specifically, after the simulation of each subsystem is completed, the obtained electrical signal is transmitted to other subsystems according to the connected transmission line in the next simulation step, and each sub-network in the other subsystems continues to simulate according to the received electrical signal until the simulation is completed, so as to obtain the simulation result of the whole circuit model.
Example 3
Based on the same inventive concept, the present embodiment further provides a real-time simulation system of a circuit model, including:
the decoupling module is used for decoupling the whole system into a plurality of sub-networks according to the distribution condition of the switching value in the circuit system;
the equivalent circuit obtaining module is used for respectively calculating all the sub-networks in a parallel calculation mode based on a set state space equation to obtain equivalent circuits corresponding to all the sub-networks;
and the circuit simulation module is used for carrying out joint simulation based on the equivalent circuits of all the subnetworks to realize real-time simulation of the circuit system.
In an embodiment, the acquiring an equivalent circuit module includes:
the equivalent circuit equation unit is used for bringing the operation parameters of each sub-network into a set state space equation to obtain an equivalent circuit equation corresponding to each sub-network;
and the equivalent circuit unit is used for carrying out parallel calculation on the equivalent circuit equation by adopting a high-order integration algorithm based on the multiprocessor to obtain an equivalent circuit of each sub-network.
In an embodiment, the circuit simulation module includes:
the circuit decoupling unit is used for inducing each sub-network into a plurality of sub-systems based on the topological structure of the circuit system;
the electric transmission signal unit is used for carrying out joint simulation by adopting a post Euler discrete integration algorithm based on the equivalent circuit of each sub-network in each sub-system in one simulation time step to obtain electric transmission signals of the sub-systems;
the circuit simulation unit is used for transmitting electric transmission signals based on the transmission line among the subsystems and realizing real-time simulation of the circuit system.
It will be appreciated by those skilled in the art that 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 flowchart illustrations and/or block diagrams, and combinations of flows and/or blocks in the flowchart illustrations 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.
The foregoing is illustrative of the present application and is not to be construed as limiting thereof, but rather as providing for the use of additional embodiments and advantages of all such modifications, equivalents, improvements and similar to the present application are intended to be included within the scope of the present application as defined by the appended claims.
Claims (10)
1. A method for real-time simulation of a circuit model, comprising:
decoupling the whole system into a plurality of sub-networks according to the distribution condition of switching value in the circuit system;
based on a set state space equation, calculating each sub-network by adopting a parallel calculation mode to obtain an equivalent circuit corresponding to each sub-network;
based on the equivalent circuits of all the sub-networks, carrying out joint simulation to realize real-time simulation of the circuit system;
the method for calculating the sub-networks by adopting the parallel calculation mode based on the set state space equation to obtain the equivalent circuit corresponding to each sub-network comprises the following steps:
the operation parameters of each sub-network are put into a set state space equation to obtain an equivalent circuit equation corresponding to each sub-network;
performing parallel calculation on the equivalent circuit equation by adopting a high-order integration algorithm based on a multiprocessor to obtain an equivalent circuit of each sub-network;
the decoupling the entire system into a plurality of sub-networks includes:
processing the circuit system according to a set decoupling mode;
decoupling the processed circuit system into a plurality of sub-networks according to the distribution condition of the switching value;
the decoupling mode comprises the following steps: a first decoupling mode and a second decoupling mode;
the first decoupling mode is to separate the switches by adopting a series direct-current voltage source and a parallel controllable alternating-current source;
the second decoupling mode is to make two IGBTs connected in series equivalent to an S function algorithm, so that the number of switches is reduced;
the equivalent circuit based on each sub-network carries out joint simulation to realize real-time simulation of the circuit system, and the method comprises the following steps:
summarizing each sub-network into a plurality of sub-systems based on the topology of the circuitry;
in a simulation time step, based on the equivalent circuit of each sub-network in each sub-system, carrying out joint simulation by adopting a post Euler discrete integration algorithm to obtain an electric transmission signal of the sub-system;
and the subsystems transmit electrical transmission signals based on the transmission lines, so that the real-time simulation of the circuit system is realized.
2. The method of claim 1, wherein the state space equation is as follows:
y＝C n x+D n u
wherein:inputting a vector; y: outputting a vector; x: state variables in the circuit model; u: a current or voltage variable; a is that n : a matrix A; b (B) n : a matrix B; c (C) n : a matrix C; d (D) n : a matrix D; n: the order of the coefficient matrix.
3. The method of claim 2, wherein the equivalent circuit equation is represented by the formula:
y t+Δt ＝y h +W n u t+Δt
wherein y is t+Δt : an equivalent circuit equation; y is h : a state space group history output; w (W) n : an admittance matrix; u (u) t+Δt : the voltage variable at the next moment; Δt: step length of discrete sampling;
wherein the state space group history output y h Admittance matrix W n Calculated as follows:
wherein:a discretized matrix A; x is x t : the current time value of the intermediate variable; />A discretized matrix B; u (u) t : the voltage variable at the current moment.
4. The method of claim 1, wherein the transmitting electrical transmission signals between the subsystems based on the transmission line, implementing real-time simulation of the circuitry, comprises:
and after delaying the signal transmission among the subsystems by one simulation step, simulating based on the received electrical signals until the simulation is completed to obtain the real-time simulation of the whole circuit.
5. The method of claim 1, wherein the transmitting electrical transmission signals between the subsystems based on the transmission line, after implementing the real-time simulation of the circuit system, further comprises:
and verifying the topological structure of the circuit model and simulating in real time through an external controller.
6. A real-time simulation system of a circuit model, comprising:
the decoupling module is used for decoupling the whole system into a plurality of sub-networks according to the distribution condition of the switching value in the circuit system;
the equivalent circuit obtaining module is used for respectively calculating all the sub-networks in a parallel calculation mode based on a set state space equation to obtain equivalent circuits corresponding to all the sub-networks;
the circuit simulation module is used for carrying out joint simulation based on the equivalent circuits of all the sub-networks to realize real-time simulation of the circuit system;
the decoupling module decouples the whole system into a plurality of sub-networks, including:
processing the circuit system according to a set decoupling mode;
decoupling the processed circuit system into a plurality of sub-networks according to the distribution condition of the switching value;
the decoupling mode comprises the following steps: a first decoupling mode and a second decoupling mode;
the first decoupling mode is to separate the switches by adopting a series direct-current voltage source and a parallel controllable alternating-current source;
the second decoupling mode is to make two IGBTs connected in series equivalent to an S function algorithm, so that the number of switches is reduced;
the obtaining equivalent circuit module includes:
the equivalent circuit equation unit is used for bringing the operation parameters of each sub-network into a set state space equation to obtain an equivalent circuit equation corresponding to each sub-network;
the equivalent circuit unit is used for carrying out parallel calculation on the equivalent circuit equation by adopting a high-order integration algorithm based on a multiprocessor to obtain an equivalent circuit of each sub-network;
the circuit simulation module comprises:
the circuit decoupling unit is used for inducing each sub-network into a plurality of sub-systems based on the topological structure of the circuit system;
the electric transmission signal unit is used for carrying out joint simulation by adopting a post Euler discrete integration algorithm based on the equivalent circuit of each sub-network in each sub-system in one simulation time step to obtain electric transmission signals of the sub-systems;
the circuit simulation unit is used for transmitting electric transmission signals based on the transmission line among the subsystems and realizing real-time simulation of the circuit system.
7. The system of claim 6, wherein the state space equation is as follows:
y＝C n x+D n u
wherein:inputting a vector; y: outputting a vector; x: state variables in the circuit model; u: a current or voltage variable; a is that n : a matrix A; b (B) n : a matrix B; c (C) n : a matrix C; d (D) n : a matrix D; n: the order of the coefficient matrix.
8. The system of claim 7, wherein the equivalent circuit equation is represented by the formula:
y t+Δt ＝y h +W n u t+Δt
wherein y is t+Δt : an equivalent circuit equation; y is h : a state space group history output; w (W) n : an admittance matrix; u (u) t+Δt : the voltage variable at the next moment; Δt: step length of discrete sampling;
wherein the state space group history output y h Admittance matrix W n Calculated as follows:
wherein:a discretized matrix A; x is x t : the current time value of the intermediate variable; />A discretized matrix B; u (u) t : the voltage variable at the current moment.
9. The system of claim 6, wherein the subsystem performs real-time simulation of the circuit system based on transmission of electrical transmission signals by the transmission line, and the system comprises:
and after delaying the signal transmission among the subsystems by one simulation step, simulating based on the received electrical signals until the simulation is completed to obtain the real-time simulation of the whole circuit.
10. The system of claim 6, wherein the subsystems are configured to transmit electrical transmission signals based on a transmission line, and wherein the real-time simulation of the circuitry is further configured to:
and verifying the topological structure of the circuit model and simulating in real time through an external controller.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811562834.6A CN109782625B (en) | 2018-12-20 | 2018-12-20 | Real-time simulation method and system for circuit model |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811562834.6A CN109782625B (en) | 2018-12-20 | 2018-12-20 | Real-time simulation method and system for circuit model |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109782625A CN109782625A (en) | 2019-05-21 |
| CN109782625B true CN109782625B (en) | 2023-11-14 |
Family
ID=66497465
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811562834.6A Active CN109782625B (en) | 2018-12-20 | 2018-12-20 | Real-time simulation method and system for circuit model |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109782625B (en) |
Families Citing this family (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111209675B (en) * | 2020-01-10 | 2023-04-28 | 南方电网科学研究院有限责任公司 | Simulation method and device of power electronic device, terminal equipment and storage medium |
| CN111464400B (en) * | 2020-04-23 | 2021-09-24 | 中国人民解放军国防科技大学 | Double-end multi-state network reliability evaluation method based on state space decomposition |
| CN111695318A (en) * | 2020-05-29 | 2020-09-22 | 苏州复鹄电子科技有限公司 | Method for calculating semiconductor device model of circuit simulation program |
| CN111783279A (en) * | 2020-06-08 | 2020-10-16 | 中国电力科学研究院有限公司 | Joint simulation method and device of system, electronic equipment and readable storage medium |
| CN111709107B (en) * | 2020-06-16 | 2023-03-24 | 全球能源互联网研究院有限公司 | Power network simulation method and system containing large-scale power electronic equipment |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH06282596A (en) * | 1992-03-02 | 1994-10-07 | Nec Corp | Parallel circuit simulation system |
| CN102346442A (en) * | 2011-05-12 | 2012-02-08 | 天津大学 | Method for carrying out parallel transient simulation computing on microgrid suitable for nodal analysis frameworks |
| CN103678900A (en) * | 2013-12-06 | 2014-03-26 | 云南电力试验研究院（集团）有限公司电力研究院 | Network decoupling calculation method used for regional power distribution network real-time simulation |
| CN106372339A (en) * | 2016-09-05 | 2017-02-01 | 清华大学 | Multi-rate simulation method and device for power-electronizing power system |
| CN106909702A (en) * | 2015-12-23 | 2017-06-30 | 加拿大欧泊实时技术有限公司北京代表处 | A kind of modular multilevel rectifier and transformer station, Modeling and simulation platform, method |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6694044B1 (en) * | 1999-09-16 | 2004-02-17 | Hewlett-Packard Development Company, L.P. | Method for motion classification using switching linear dynamic system models |
-
2018
- 2018-12-20 CN CN201811562834.6A patent/CN109782625B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH06282596A (en) * | 1992-03-02 | 1994-10-07 | Nec Corp | Parallel circuit simulation system |
| CN102346442A (en) * | 2011-05-12 | 2012-02-08 | 天津大学 | Method for carrying out parallel transient simulation computing on microgrid suitable for nodal analysis frameworks |
| CN103678900A (en) * | 2013-12-06 | 2014-03-26 | 云南电力试验研究院（集团）有限公司电力研究院 | Network decoupling calculation method used for regional power distribution network real-time simulation |
| CN106909702A (en) * | 2015-12-23 | 2017-06-30 | 加拿大欧泊实时技术有限公司北京代表处 | A kind of modular multilevel rectifier and transformer station, Modeling and simulation platform, method |
| CN106372339A (en) * | 2016-09-05 | 2017-02-01 | 清华大学 | Multi-rate simulation method and device for power-electronizing power system |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109782625A (en) | 2019-05-21 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109782625B (en) | Real-time simulation method and system for circuit model | |
| CN111709107B (en) | Power network simulation method and system containing large-scale power electronic equipment | |
| CN101719182B (en) | Parallel partition electromagnetic transient digital simulation method of AC and DC power system | |
| CN102436534B (en) | Simulating method for electric power system | |
| CN103678900B (en) | A kind of network decoupling calculation method for regional distribution network real-time simulation | |
| CN103810646B (en) | Improved projection integral algorithm based active power distribution system dynamic simulation method | |
| CN105938578A (en) | Large-scale photovoltaic power station equivalent modeling method based on clustering analysis | |
| CN108959671B (en) | Real-time simulation modeling method for half-bridge and full-bridge modular multilevel converter | |
| CN103020385A (en) | RTDS-based 500-Kv major network modeling simulation systems for power network | |
| CN103913996A (en) | Electromagnetic electromechanical sub-network real-time simulation method based on broadband equivalence large-scale power system | |
| CN106067792A (en) | High-power fractional order electric capacity and the control method thereof that a kind of exponent number is more than 1 | |
| CN103049617A (en) | Method for simplifying passive electromagnetic transient-state simulation model of large-scale power distribution network | |
| CN109638931A (en) | The power converter system model predictive control method and system of more DC-DC parallel connections | |
| CN113865895B (en) | Locomotive and motor train unit control algorithm test system and application method thereof | |
| CN105449705B (en) | The emulation mode and system of tap transformer real time load pressure regulation in electric system | |
| CN103729502A (en) | Method for increasing electromagnetic transient simulation speed of power system | |
| CN109149566B (en) | Modeling method of simulation model for predicting lowest frequency point under high-power deficiency | |
| CN110705035A (en) | Discrete state event driven decoupling simulation method of large-capacity power electronic system | |
| CN117269838A (en) | Method and system for determining short-circuit current of network-structured power electronic equipment | |
| CN103279626B (en) | A kind of Building of Simulation Model method and system | |
| CN111079268A (en) | Online simulation method and system for LCC direct current transmission electromechanical electromagnetic hybrid system | |
| CN106961223A (en) | The current-sharing control method and device of a kind of switching rectifier | |
| CN115173467A (en) | Wind power plant simulation system and method through flexible direct current grid connection and storage medium | |
| CN111725818B (en) | Grid-connected simulation method and simulation terminal for three-phase weak current network converter group | |
| CN108763788B (en) | Simulation method and device of information physical system |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |