Electromagnetic transient simulation method and system

Abstract

The invention provides an electromagnetic transient simulation method and system. An initialization module is used for caching element parameters and impedance matrix parameters; when an element updating enabling instruction transmitted by the main control module is obtained, the element updating module is used for obtaining historical current of each element in a simulation time period according to element parameters and a node voltage vector of a previous simulation time period; and when a solving instruction transmitted by the main control module is obtained, the solving module is used for carrying out solving according to the historical current and impedance matrix parameters of each element in the simulation time period so as to obtain a node current vector and / or a node voltage vector in the simulation time period. The historical current of each element is quickly and accurately obtained through the element updating module, so that the solving module can quickly determine the node current vector and / or the node voltage vector, and the microsecond-level or nanosecond-level simulation step length requirement is met.

Description

Technical Field

[0001] The present invention relates to the field of power systems, and in particular to an electromagnetic transient simulation method and system. Background Art

[0002] With the continuous increase in energy demand and the rapid development of renewable energy power generation, the power grid system presents the "double high" characteristics of high proportion of renewable energy and high proportion of power electronic equipment access. A large number of renewable energy power generation systems such as photovoltaic and wind power are connected to the power grid, which affects the security and stability of the power grid. In order to ensure the safe and stable operation of the power grid, power grid practitioners and renewable energy manufacturers are seeking more accurate and reliable solutions to analyze the impact of equipment access on the stability of the power grid. Therefore, electromagnetic transient simulation plays an increasingly important role in system planning, design and operation.

[0003] Hardware-in-the-loop simulation is currently a commonly used means to assess the safety risks of equipment access to the network. Different from large-scale system-level full electromagnetic transient simulation, the scale of hardware-in-the-loop simulation is often limited. It requires real-time simulation of the device under test and the digital model while ensuring accuracy, that is, the simulation results must be synchronized with the real-time clock to achieve real-time data interaction between the real physical device and the digital model. In real-time simulation, the simulation step size is a key factor in reproducing high-frequency transient characteristics. For power electronic systems, there are a large number of power electronic devices that open and close quickly, and the switching frequency is generally at the kHz level. According to the sampling law, in order to accurately sample and simulate the pulse width modulation (PWM) signal of such systems, the simulation cycle needs to be 1 / 50 or 1 / 100 of the PWM cycle, that is, at least 50 to 100 samples are required in one PWM cycle to better simulate the impact of the control signal on the power electronic system. Therefore, a simulation step size of microseconds or nanoseconds is required to meet the requirements of electromagnetic transient simulation of high-frequency power electronic devices. Summary of the invention

[0004] The object of the present invention is to provide an electromagnetic transient simulation method and system to improve the above-mentioned problem.

[0005] In order to achieve the above purpose, the technical solution adopted by the embodiment of the present invention is as follows:

[0006] In a first aspect, an embodiment of the present invention provides an electromagnetic transient simulation system, the electromagnetic transient simulation system comprising: a main control module, an initialization module, a component update module and a solution module, the main control module is respectively connected to the component update module, the solution module and the initialization module, the component update module is connected to the solution module, the initialization module is respectively connected to the component update module and the solution module, and the initialization module is used to cache component parameters and impedance matrix parameters;

[0007] When the component update enable instruction transmitted by the main control module is obtained, the component update module is used to obtain the historical current of each component in the simulation period according to the component parameters and the node voltage vector of the previous simulation period;

[0008] After obtaining the solution instruction transmitted by the main control module, the solution module is used to perform a solution according to the historical current of each component in the simulation period and the impedance matrix parameters to obtain the node current vector and / or node voltage vector in the simulation period;

[0009] The node current vector includes the node current of each node in the simulation period, and the node voltage vector includes the node voltage of each node in the simulation period.

[0010] Optionally, the component update module includes any one or more of an RLC component update module, a power supply type component update module, a switch type component update module, a four-port component update module and a six-port component update module;

[0011] The RLC element updating module is used to obtain the historical current of each RLC element in the simulation period;

[0012] The power supply type component update module is used to obtain the historical current of each power supply type component during the simulation period;

[0013] The switch type element updating module is used to obtain the historical current of each switch type element in the simulation period;

[0014] The four-port element updating module is used to obtain the historical current of each four-port element in the simulation period;

[0015] The six-port element updating module is used to obtain the historical current of each six-port element in the simulation period.

[0016] Optionally, the calculation formula of the RLC element update module is:

[0017]

[0018] I k (t) = G k ×U k (t)-I h，k (t-dT) (2)

[0019] I h，k (t) = P k ×U k (t)+Q k ×I k (t) (3)

[0020] Among them, U k (t) represents the branch voltage of the kth element, represents the voltage value of the first node of the kth element in the previous simulation period, represents the voltage value of the end node of the kth element in the previous simulation period, I k (t) represents the branch current corresponding to the kth element, G k represents the equivalent admittance coefficient of the kth element, I h,k (t-dT) represents the historical current of the kth element in the previous simulation period, I h,k (t) represents the historical current of the kth component in the current simulation period, P k represents the branch voltage term coefficient corresponding to the kth element, Q k Represents the branch current term coefficient corresponding to the kth element.

[0021] Optionally, the calculation formula of the power supply type component update module is:

[0022]

[0023] I k (t) = G k ×U k (t)-I h，k (t-dT) (5)

[0024] I h，k (t) = P k ×Ctrl k (t) (6)

[0025] Among them, U k (t) represents the branch voltage of the kth element, represents the voltage value of the first node of the kth element in the previous simulation period, represents the voltage value of the end node of the kth element in the previous simulation period, I k (t) represents the branch current corresponding to the kth element, G k represents the equivalent admittance coefficient of the kth element, I h,k (t-dT) represents the historical current of the kth element in the previous simulation period, I h,k (t) represents the historical current of the kth component in the current simulation period, P k Indicates the controlled signal coefficient corresponding to the kth element, Ctrl k (t) represents the controlled signal corresponding to the kth element.

[0026] Optionally, the calculation formula of the switch element update module is:

[0027]

[0028] I k (t) = G k ×U k (t)-I h，k (t-dT) (8)

[0029] When the switch control signal SW k When the switch corresponding to the kth element is turned off, we have:

[0030] I h，k (t) = P 1k ×U k (t)+Q 1k ×I k (t) (9)

[0031] When the switch control signal SW k When the switch corresponding to the kth element is turned on, we have:

[0032] I h，k (t) = P 2k ×U k (t)+Q 2k ×I k (t) (10)

[0033] Among them, U k (t) represents the branch voltage of the kth element, represents the voltage value of the first node of the kth element in the previous simulation period, represents the voltage value of the end node of the kth element in the previous simulation period, I k (t) represents the branch current corresponding to the kth element, G k represents the equivalent admittance coefficient of the kth element, I h,k (t-dT) represents the historical current of the kth element in the previous simulation period, I h,k (t) represents the historical current of the kth component in the current simulation period, P 1k It represents the branch voltage coefficient corresponding to the kth element when the switch is turned off, Q 1k It represents the branch current coefficient corresponding to the kth element when the switch is turned off, P 2k It represents the branch voltage coefficient corresponding to the kth element when the switch is turned on, Q 2k It represents the branch current coefficient corresponding to the kth element when the switch is turned on.

[0034] Optionally, the calculation formula of the four-port element update module is:

[0035]

[0036] Among them, U 1k (t) represents the first branch voltage of the kth element, U 2k (t) represents the second branch voltage of the kth element, represents the voltage value of the first node of the first branch of the kth element in the previous simulation period, represents the voltage value of the first node of the second branch of the kth element in the previous simulation period, represents the voltage value of the end node of the first branch of the kth element in the previous simulation period, represents the voltage value of the end node of the second branch of the kth element in the previous simulation period, I 1k (t) represents the first branch current of the kth element, I 2k (t) represents the second branch current of the kth element, G k represents the 2*2 dimensional equivalent admittance matrix corresponding to the kth element, I h1,k (t-dT) represents the historical current of the first branch of the kth element in the previous simulation period, I h2,k (t-dT) represents the historical current of the second branch of the kth element in the previous simulation period, I h1,k (t) represents the historical current of the first branch of the kth element in the current simulation period, I h2,k (t) represents the historical current of the second branch of the kth element in the current simulation period, P k represents the 2*2-dimensional branch voltage coefficient matrix corresponding to the kth element, Q k Represents the 2*2 dimensional branch current term coefficient matrix corresponding to the kth element.

[0037] Optionally, the calculation formula of the six-port element update module is:

[0038]

[0039] Among them, U ik (t) represents the voltage of the i-th branch of the k-th element, represents the voltage value of the first node of the i-th branch of the k-th element in the previous simulation period, represents the voltage value of the end node of the i-th branch of the k-th element in the previous simulation period, I ik (t) represents the i-th branch current of the k-th element, G k represents the 3*3 dimensional equivalent admittance matrix corresponding to the kth element, I hi,k (t-dT) represents the historical current of the i-th branch of the k-th element in the previous simulation period, Ihi,k (t) represents the historical current of the i-th branch of the k-th element in the current simulation period, P k represents the 3*3-dimensional branch voltage term coefficient matrix corresponding to the kth element, Q k Represents the 3*3 dimensional branch current term coefficient matrix corresponding to the kth element, 1≤i≤3.

[0040] Optionally, the initialization module is used to communicate with an external device to receive initialization parameters after acquiring the initialization start instruction transmitted by the main control module, wherein the initialization parameters include simulation parameters, component parameters and impedance matrix parameters;

[0041] The simulation parameters include the total number of nodes in the electromagnetic transient simulation system, the simulation integration step, the total simulation time and the data interaction interval;

[0042] The component parameters include type information of each component, the number of components of each type, node information, and input parameter information, where the input parameter information represents the management relationship between the component and the corresponding component update module.

[0043] Optionally, the main control module is also used to determine whether the simulation is ended based on the simulation duration after the solution module completes the solution; if not ended, determine the start time of the next simulation period based on the start time of the current simulation period and the simulation integration step, and transmit the component update enable instruction to the component update module and the solution instruction to the solution module at the start time of the next simulation period.

[0044] In a second aspect, an embodiment of the present invention provides an electromagnetic transient simulation method, which is applied to the above-mentioned electromagnetic transient simulation system, and the method includes:

[0045] When the component update enable instruction transmitted by the main control module is obtained, the component update module obtains the historical current of each component in the simulation period according to the component parameters and the node voltage vector of the previous simulation period;

[0046] After obtaining the solution instruction transmitted by the main control module, the solution module performs a solution according to the historical current of each component in the simulation period and the impedance matrix parameters to obtain the node current vector and / or node voltage vector in the simulation period;

[0047] The node current vector includes the node current of each node in the simulation period, and the node voltage vector includes the node voltage of each node in the simulation period.

[0048] Compared with the prior art, an electromagnetic transient simulation method and system provided by an embodiment of the present invention includes a main control module, an initialization module, a component update module and a solution module, wherein the main control module is connected to the component update module, the solution module and the initialization module respectively, the component update module is connected to the solution module, the initialization module is connected to the component update module and the solution module respectively, and the initialization module is used to cache component parameters and impedance matrix parameters; when obtaining the component update enable instruction transmitted by the main control module, the component update module is used to obtain the historical current of each component in the simulation period according to the component parameters and the node voltage vector of the previous simulation period; when obtaining the solution instruction transmitted by the main control module, the solution module is used to solve according to the historical current and impedance matrix parameters of each component in the simulation period to obtain the node current vector and / or node voltage vector in the simulation period; wherein the node current vector includes the node current of each node in the simulation period, and the node voltage vector includes the node voltage of each node in the simulation period. The historical current of each component is quickly and accurately obtained through the component update module, so that the solution module can quickly determine the node current vector and / or node voltage vector, thereby meeting the simulation step requirements of microsecond or nanosecond level.

[0049] In order to make the above-mentioned objects, features and advantages of the present invention more obvious and easy to understand, preferred embodiments are given below and described in detail with reference to the accompanying drawings. BRIEF DESCRIPTION OF THE DRAWINGS

[0050] In order to more clearly illustrate the technical solutions of the embodiments of the present invention, the drawings required for use in the embodiments will be briefly introduced below. It should be understood that the following drawings only show certain embodiments of the present invention and therefore should not be regarded as limiting the scope. For ordinary technicians in this field, other related drawings can be obtained based on these drawings without paying creative work.

[0051] Figure 1 A schematic diagram of the structure of an electromagnetic transient simulation system provided in an embodiment of the present invention.

[0052] Figure 2 A structural diagram of an RLC element update module provided by an embodiment of the present invention.

[0053] Figure 3 A structural diagram of a power supply type component update module provided in an embodiment of the present invention.

[0054] Figure 4 A structural diagram of a switch-type element update module provided in an embodiment of the present invention.

[0055] Figure 5 A structural diagram of a four-port element update module provided in an embodiment of the present invention.

[0056] Figure 6 A structural diagram of a six-port element update module provided in an embodiment of the present invention.

[0057] Figure 7 A schematic flow chart of an electromagnetic transient simulation method provided in an embodiment of the present invention.

[0058] In the figure: 10 - main control module; 20 - initialization module; 30 - component update module; 40 - solution module. DETAILED DESCRIPTION

[0059] In order to make the purpose, technical solutions and advantages of the embodiments of the present invention clearer, the technical solutions in the embodiments of the present invention will be clearly and completely described below in conjunction with the drawings in the embodiments of the present invention. Obviously, the described embodiments are part of the embodiments of the present invention, not all of the embodiments. Generally, the components of the embodiments of the present invention described and shown in the drawings here can be arranged and designed in various different configurations.

[0060] Therefore, the following detailed description of the embodiments of the present invention provided in the accompanying drawings is not intended to limit the scope of the invention claimed for protection, but merely represents selected embodiments of the present invention. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of the present invention.

[0061] It should be noted that similar reference numerals and letters represent similar items in the following drawings, so once an item is defined in one drawing, it does not need to be further defined and explained in the subsequent drawings. At the same time, in the description of the present invention, the terms "first", "second", etc. are only used to distinguish the description and cannot be understood as indicating or implying relative importance.

[0062] It should be noted that, in this article, relational terms such as first and second, etc. are only used to distinguish one entity or operation from another entity or operation, and do not necessarily require or imply any such actual relationship or order between these entities or operations. Moreover, the terms "include", "comprise" or any other variants thereof are intended to cover non-exclusive inclusion, so that a process, method, article or device including a series of elements includes not only those elements, but also other elements not explicitly listed, or also includes elements inherent to such process, method, article or device. In the absence of further restrictions, the elements defined by the sentence "comprise a ..." do not exclude the existence of other identical elements in the process, method, article or device including the elements.

[0063] In the description of the present invention, it should be noted that the terms "upper", "lower", "inside", "outside", etc. indicate directions or positional relationships based on the directions or positional relationships shown in the accompanying drawings, or are directions or positional relationships in which the product of the invention is usually placed when in use. They are only for the convenience of describing the present invention and simplifying the description, and do not indicate or imply that the device or element referred to must have a specific direction, be constructed and operated in a specific direction, and therefore should not be understood as a limitation on the present invention.

[0064] In the description of the present invention, it is also necessary to explain that, unless otherwise clearly specified and limited, the terms "disposed" and "connected" should be understood in a broad sense, for example, it can be a fixed connection, a detachable connection, or an integral connection; it can be a mechanical connection or an electrical connection; it can be a direct connection, or it can be an indirect connection through an intermediate medium, or it can be the internal communication of two elements. For ordinary technicians in this field, the specific meanings of the above terms in the present invention can be understood according to specific circumstances.

[0065] Some embodiments of the present invention are described in detail below in conjunction with the accompanying drawings. In the absence of conflict, the following embodiments and features in the embodiments can be combined with each other.

[0066] The embodiment of the present invention provides an electromagnetic transient simulation system based on a field-programmable gate array (FPGA). Figure 1 , Figure 1 A schematic diagram of the structure of an electromagnetic transient simulation system provided in an embodiment of the present invention.

[0067] The electromagnetic transient simulation system includes: a main control module 10, an initialization module 20, a component update module 30 and a solution module 40. The main control module 10 is connected to the component update module 30, the solution module 40 and the initialization module 20 respectively. The component update module 30 is connected to the solution module 40. The initialization module 20 is connected to the component update module 30 and the solution module 40 respectively. The initialization module 20 is used to cache component parameters and impedance matrix parameters.

[0068] When the component update enable instruction transmitted by the main control module 10 is obtained, the component update module 30 is used to obtain the historical current of each component in the simulation period according to the component parameters and the node voltage vector of the previous simulation period.

[0069] It should be noted that a simulation cycle includes multiple simulation time periods, and each simulation time period includes multiple clock cycles.

[0070] After obtaining the solution instruction transmitted by the main control module 10, the solution module 40 is used to solve according to the historical current and impedance matrix parameters of each component in the simulation period to obtain the node current vector and / or node voltage vector in the simulation period.

[0071] The node current vector includes the node current of each node in the simulation period, and the node voltage vector includes the node voltage of each node in the simulation period.

[0072] In the electromagnetic transient simulation system provided in an embodiment of the present invention, the historical current of each component is quickly and accurately obtained through the component update module, so that the solution module can quickly determine the node current vector and / or node voltage vector, thereby meeting the simulation step requirements of microseconds or nanoseconds.

[0073] Because there are various types of components in the power grid system, including RLC components, power type components, switch type components, four-port components, six-port components, etc. Among them, RLC components represent resistance, inductance and capacitance components. In order to accurately obtain the historical current of different types of components, the embodiment of the present invention also provides an optional implementation method, please refer to the following.

[0074] The component update module 30 includes any one or more of an RLC component update module, a power supply component update module, a switch component update module, a four-port component update module, and a six-port component update module. It should be noted that the number of each type of component update module included is not limited to 1, and can be greater than or equal to 1.

[0075] The RLC component update module is used to obtain the historical current of each RLC component during the simulation period.

[0076] The power supply component update module is used to obtain the historical current of each power supply component during the simulation period.

[0077] The switch element update module is used to obtain the historical current of each switch element during the simulation period.

[0078] The four-port component update module is used to obtain the historical current of each four-port component during the simulation period.

[0079] The six-port component update module is used to obtain the historical current of each six-port component during the simulation period.

[0080] It should be noted that each component update module 30 can calculate the historical current of multiple components during the simulation period. Based on the parallel computing characteristics of FPGA, the parameters of each component type are input using a pipeline. For example, assuming that 26 FPGA clocks (clk) are required to calculate the historical current of the RLC component, and the RLC component update module contains 5 RLC components. Then the first component parameter is input in the first clk, and the second component parameter is input in the second clk, until the fifth clk completes the parameter input of all components. After 26 clks after the first component parameter is input, that is, the 27th clk completes the historical current calculation of the first RLC component, and the 31st clk completes the historical current calculation of all RLC components.

[0081] Based on the foregoing, the component update modules of each component type are described in detail as follows.

[0082] Please refer to Figure 2 , Figure 2 The structure diagram of the RLC element update module provided in the embodiment of the present invention. The node current merging and network equation solving module in the figure corresponds to the solving module 40 in the present application. The electromagnetic transient simulation system also includes a node voltage storage unit for storing the node voltage vector U of the previous simulation period fed back by the solving module. n (t-dT), the component update module can obtain the node voltage vector U stored in the node voltage storage unit based on the node label n (t-dT) Gets the voltage value of the node in the previous simulation period.

[0083] Optionally, the calculation formula of the RLC element update module is:

[0084]

[0085] I k (t) = G k ×U k (t)-I h，k (t-dT) (2)

[0086] I h，k (t) = P k ×U k (t)+Q k ×I k (t) (3)

[0087] Among them, U k (t) represents the branch voltage of the kth element, represents the voltage value of the first node of the kth element in the previous simulation period, represents the voltage value of the end node of the kth element in the previous simulation period, I k(t) represents the branch current corresponding to the kth element, G k represents the equivalent admittance coefficient of the kth element, I h,k (t-dT) represents the historical current of the kth element in the previous simulation period, I h,k (t) represents the historical current of the kth component in the current simulation period, P k represents the branch voltage term coefficient corresponding to the kth element, Q k Represents the branch current term coefficient corresponding to the kth element.

[0088] Please refer to Figure 3 , Figure 3 A structural diagram of a power supply type component update module provided in an embodiment of the present invention.

[0089] Optionally, the calculation formula of the power supply type component update module is:

[0090]

[0091] I k (t) = G k ×U k (t)-I h，k (t-dT) (5)

[0092] I h，k (t) = P k ×Ctrl k (t) (6)

[0093] Among them, U k (t) represents the branch voltage of the kth element, represents the voltage value of the first node of the kth element in the previous simulation period, represents the voltage value of the end node of the kth element in the previous simulation period, I k (t) represents the branch current corresponding to the kth element, G k represents the equivalent admittance coefficient of the kth element, I h,k (t-dT) represents the historical current of the kth element in the previous simulation period, I h,k (t) represents the historical current of the kth component in the current simulation period, P k Indicates the controlled signal coefficient corresponding to the kth element, Ctrl k (t) represents the controlled signal corresponding to the kth element (also called the power controlled signal).

[0094] Please refer to Figure 4 , Figure 4 A structural diagram of a switch-type element update module provided in an embodiment of the present invention.

[0095] Optionally, the calculation formula of the switch type element update module is:

[0096]

[0097] I k (t) = G k ×U k (t)-I h，k (t-dT) (8)

[0098] When the switch control signal SW k When the switch corresponding to the kth element is turned off, we have:

[0099] I h，k (t) = P 1k ×U k (t)+Q 1k ×I k (t) (9)

[0100] When the switch control signal SW k When the switch corresponding to the kth element is turned on, we have:

[0101] I h，k (t) = P 2k ×U k (t)+Q 2k ×I k (t) (10

[0102] Among them, U k (t) represents the branch voltage of the kth element, represents the voltage value of the first node of the kth element in the previous simulation period, represents the voltage value of the end node of the kth element in the previous simulation period, I k (t) represents the branch current corresponding to the kth element, G k represents the equivalent admittance coefficient of the kth element, I h,k (t-dT) represents the historical current of the kth element in the previous simulation period, I h,k (t) represents the historical current of the kth component in the current simulation period, P 1k It represents the branch voltage coefficient corresponding to the kth element when the switch is turned off, Q 1k It represents the branch current coefficient corresponding to the kth element when the switch is turned off, P 2k It represents the branch voltage coefficient corresponding to the kth element when the switch is turned on, Q 2k It represents the branch current coefficient corresponding to the kth element when the switch is turned on.

[0103] It should be noted that the main control module 10 can also send a data interaction enable instruction to the component update module 30, so that the component update module 30 interacts with the external device, thereby obtaining the switch control signal SW corresponding to the kth component. k The controlled signal Ctrl corresponding to the kth element k (t).

[0104] Please refer to Figure 5 , Figure 5 A structural diagram of a four-port element update module provided in an embodiment of the present invention.

[0105] Optionally, the calculation formula of the four-port element update module is:

[0106]

[0107] Among them, U 1k (t) represents the first branch voltage of the kth element, U 2k (t) represents the second branch voltage of the kth element, represents the voltage value of the first node of the first branch of the kth element in the previous simulation period, represents the voltage value of the first node of the second branch of the kth element in the previous simulation period, represents the voltage value of the end node of the first branch of the kth element in the previous simulation period, represents the voltage value of the end node of the second branch of the kth element in the previous simulation period, I 1k (t) represents the first branch current of the kth element, I 2k (t) represents the second branch current of the kth element, G k represents the 2*2 dimensional equivalent admittance matrix corresponding to the kth element, I h1,k (t-dT) represents the historical current of the first branch of the kth element in the previous simulation period, I h2,k (t-dT) represents the historical current of the second branch of the kth element in the previous simulation period, I h1,k (t) represents the historical current of the first branch of the kth element in the current simulation period, I h2,k (t) represents the historical current of the second branch of the kth element in the current simulation period, P k represents the 2*2-dimensional branch voltage coefficient matrix corresponding to the kth element, Q k Represents the 2*2 dimensional branch current term coefficient matrix corresponding to the kth element.

[0108] Please refer to Figure 6 , Figure 6 A structural diagram of a six-port element update module provided in an embodiment of the present invention.

[0109] Optionally, the calculation formula of the six-port element update module is:

[0110]

[0111] Among them, U ik (t) represents the voltage of the i-th branch of the k-th element, represents the voltage value of the first node of the i-th branch of the k-th element in the previous simulation period, represents the voltage value of the end node of the i-th branch of the k-th element in the previous simulation period, I ik (t) represents the i-th branch current of the k-th element, G k represents the 3*3 dimensional equivalent admittance matrix corresponding to the kth element, I hi,k (t-dT) represents the historical current of the i-th branch of the k-th element in the previous simulation period, I hi,k (t) represents the historical current of the i-th branch of the k-th element in the current simulation period, P k represents the 3*3-dimensional branch voltage term coefficient matrix corresponding to the kth element, Q k Represents the 3*3 dimensional branch current term coefficient matrix corresponding to the kth element, 1≤i≤3.

[0112] Please continue to refer to Figure 1 In an optional implementation, the initialization module 20 is used to communicate with an external device after acquiring the initialization start instruction transmitted by the main control module 10 to receive initialization parameters, wherein the initialization parameters include simulation parameters, component parameters, and impedance matrix parameters.

[0113] The simulation parameters include the total number of nodes in the electromagnetic transient simulation system, the simulation integration step, the total simulation time, and the data interaction interval (including the interaction interval of the controlled signal and the switch control signal).

[0114] The component parameters include the type information of each component, the number of components of each type, node information, and input parameter information, which indicates the management relationship between the component and the corresponding component update module 30 (including the above-mentioned G k , P k and Q k and other parameters).

[0115] Optionally, the main control module 10 is also used to determine whether the simulation is ended based on the simulation duration after the solution module 40 completes the solution. If not ended, the start time of the next simulation period is determined based on the start time of the current simulation period and the simulation integration step, and the component update enable instruction is transmitted to the component update module 30 and the solution instruction is transmitted to the solution module 40 at the start time of the next simulation period.

[0116] The embodiment of the present invention further provides an electromagnetic transient simulation method, which is applied to the above-mentioned electromagnetic transient simulation system. Figure 7 , electromagnetic transient simulation methods include

[0117] S310, when the component update enable instruction transmitted by the main control module is obtained, the component update module obtains the historical current of each component in the simulation period according to the component parameters and the node voltage vector of the previous simulation period.

[0118] S410, after obtaining the solution instruction transmitted by the main control module, the solution module performs a solution according to the historical current and impedance matrix parameters of each component in the simulation period to obtain the node current vector and / or node voltage vector in the simulation period.

[0119] The node current vector includes the node current of each node during the simulation period, and the node voltage vector includes the node voltage of each node during the simulation period.

[0120] It should be noted that the electromagnetic transient simulation method provided in this embodiment can perform the functions and purposes shown in the above electromagnetic transient simulation system embodiment to achieve the corresponding technical effects. For the sake of brief description, for parts not mentioned in this embodiment, reference can be made to the corresponding contents in the above embodiments.

[0121] In summary, an electromagnetic transient simulation method and system provided by an embodiment of the present invention includes a main control module, an initialization module, a component update module and a solution module, wherein the main control module is connected to the component update module, the solution module and the initialization module respectively, the component update module is connected to the solution module, the initialization module is connected to the component update module and the solution module respectively, and the initialization module is used to cache component parameters and impedance matrix parameters; when the component update enable instruction transmitted by the main control module is obtained, the component update module is used to obtain the historical current of each component in the simulation period according to the component parameters and the node voltage vector of the previous simulation period; when the solution instruction transmitted by the main control module is obtained, the solution module is used to solve according to the historical current and impedance matrix parameters of each component in the simulation period to obtain the node current vector and / or node voltage vector in the simulation period; wherein the node current vector includes the node current of each node in the simulation period, and the node voltage vector includes the node voltage of each node in the simulation period. The historical current of each component is quickly and accurately obtained through the component update module, so that the solution module can quickly determine the node current vector and / or node voltage vector, thereby meeting the simulation step requirements of microsecond or nanosecond level.

[0122] The above description is only a preferred embodiment of the present invention and is not intended to limit the present invention. For those skilled in the art, the present invention may have various modifications and variations. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the present invention shall be included in the protection scope of the present invention.

[0123] It will be apparent to those skilled in the art that the invention is not limited to the details of the exemplary embodiments described above and that the invention can be implemented in other specific forms without departing from the spirit or essential features of the invention. Therefore, the embodiments should be considered exemplary and non-limiting in all respects, and the scope of the invention is defined by the appended claims rather than the foregoing description, and it is intended that all variations falling within the meaning and scope of the equivalent elements of the claims be included in the invention. Any reference numeral in a claim should not be considered as limiting the claim to which it relates.

Claims

1. An electromagnetic transient simulation system, characterized in that: The electromagnetic transient simulation system comprises: a main control module, an initialization module, a component update module and a solution module, wherein the main control module is respectively connected to the component update module, the solution module and the initialization module, the component update module is connected to the solution module, the initialization module is respectively connected to the component update module and the solution module, and the initialization module is used to cache component parameters and impedance matrix parameters; When the component update enable instruction transmitted by the main control module is obtained, the component update module is used to obtain the historical current of each component in the simulation period according to the component parameters and the node voltage vector of the previous simulation period; After obtaining the solution instruction transmitted by the main control module, the solution module is used to perform a solution according to the historical current of each component in the simulation period and the impedance matrix parameters to obtain the node current vector and / or node voltage vector in the simulation period; The node current vector includes the node current of each node in the simulation period, and the node voltage vector includes the node voltage of each node in the simulation period.

2. The electromagnetic transient simulation system according to claim 1, characterized in that: The component update module includes any one or more of an RLC component update module, a power supply type component update module, a switch type component update module, a four-port component update module and a six-port component update module; The RLC element updating module is used to obtain the historical current of each RLC element in the simulation period; The power supply type component update module is used to obtain the historical current of each power supply type component during the simulation period; The switch type element updating module is used to obtain the historical current of each switch type element in the simulation period; The four-port element updating module is used to obtain the historical current of each four-port element in the simulation period; The six-port element updating module is used to obtain the historical current of each six-port element in the simulation period.

3. The electromagnetic transient simulation system according to claim 2, characterized in that: The calculation formula of the RLC element update module is: I k (t)=G k ×U k (t)-I h，k (t-dT) (2) I h，k (t)=P k ×U k (t)+Q k ×I k (t) (3) Among them, U k (t) represents the branch voltage of the kth element, represents the voltage value of the first node of the kth element in the previous simulation period, represents the voltage value of the end node of the kth element in the previous simulation period, I k (t) represents the branch current corresponding to the kth element, G k represents the equivalent admittance coefficient of the kth element, I h,k (t-dT) represents the historical current of the kth element in the previous simulation period, I h,k (t) represents the historical current of the kth component in the current simulation period, P k represents the branch voltage term coefficient corresponding to the kth element, Q k Represents the branch current term coefficient corresponding to the kth element.

4. The electromagnetic transient simulation system according to claim 2, characterized in that: The calculation formula of the power supply type component update module is: I k (t)=G k ×U k (t)-I h，k (t-dT) (5) I h，k (t)=P k ×Ctrl k (t) (6) Among them, U k (t) represents the branch voltage of the kth element, represents the voltage value of the first node of the kth element in the previous simulation period, represents the voltage value of the end node of the kth element in the previous simulation period, I k (t) represents the branch current corresponding to the kth element, G k represents the equivalent admittance coefficient of the kth element, I h,k (t-dT) represents the historical current of the kth element in the previous simulation period, I h,k (t) represents the historical current of the kth component in the current simulation period, P k Indicates the controlled signal coefficient corresponding to the kth element, Ctrl k (t) represents the controlled signal corresponding to the kth element.

5. The electromagnetic transient simulation system according to claim 2, characterized in that: The calculation formula of the switch type element update module is: I k (t)=G k ×U k (t)-I h，k (t-dT) (8) When the switch control signal SW k When the switch corresponding to the kth element is turned off, we have: I h，k (t)=P 1k ×U k (t)+Q 1k ×I k (t) (9) When the switch control signal SW k When the switch corresponding to the kth element is turned on, we have: I h，k (t)=P 2k ×U k (t)+Q 2k ×I k (t) (10) Among them, U k (t) represents the branch voltage of the kth element, represents the voltage value of the first node of the kth element in the previous simulation period, represents the voltage value of the end node of the kth element in the previous simulation period, I k (t) represents the branch current corresponding to the kth element, G k represents the equivalent admittance coefficient of the kth element, I h,k (t-dT) represents the historical current of the kth element in the previous simulation period, I h,k (t) represents the historical current of the kth component in the current simulation period, P 1k It represents the branch voltage coefficient corresponding to the kth element when the switch is turned off, Q 1k It represents the branch current coefficient corresponding to the kth element when the switch is turned off, P 2k It represents the branch voltage coefficient corresponding to the kth element when the switch is turned on, Q 2k It represents the branch current coefficient corresponding to the kth element when the switch is turned on.

6. The electromagnetic transient simulation system according to claim 2, characterized in that: The calculation formula of the four-port element update module is: Among them, U 1k (t) represents the first branch voltage of the kth element, U 2k (t) represents the second branch voltage of the kth element, represents the voltage value of the first node of the first branch of the kth element in the previous simulation period, represents the voltage value of the first node of the second branch of the kth element in the previous simulation period, represents the voltage value of the end node of the first branch of the kth element in the previous simulation period, represents the voltage value of the end node of the second branch of the kth element in the previous simulation period, I 1k (t) represents the first branch current of the kth element, I 2k (t) represents the second branch current of the kth element, G k represents the 2*2 dimensional equivalent admittance matrix corresponding to the kth element, I h1,k (t-dT) represents the historical current of the first branch of the kth element in the previous simulation period, I h2,k (t-dT) represents the historical current of the second branch of the kth element in the previous simulation period, I h1,k (t) represents the historical current of the first branch of the kth element in the current simulation period, I h2,k (t) represents the historical current of the second branch of the kth element in the current simulation period, P k represents the 2*2-dimensional branch voltage coefficient matrix corresponding to the kth element, Q k Represents the 2*2 dimensional branch current term coefficient matrix corresponding to the kth element.

7. The electromagnetic transient simulation system according to claim 2, characterized in that: The calculation formula of the six-port element update module is: Among them, U ik (t) represents the voltage of the i-th branch of the k-th element, represents the voltage value of the first node of the i-th branch of the k-th element in the previous simulation period, represents the voltage value of the end node of the i-th branch of the k-th element in the previous simulation period, I ik (t) represents the i-th branch current of the k-th element, G k represents the 3*3 dimensional equivalent admittance matrix corresponding to the kth element, I hi,k (t-dT) represents the historical current of the i-th branch of the k-th element in the previous simulation period, I hi,k (t) represents the historical current of the i-th branch of the k-th element in the current simulation period, P k represents the 3*3 dimensional branch voltage term coefficient matrix corresponding to the kth element, Q k Represents the 3*3 dimensional branch current term coefficient matrix corresponding to the kth element, 1≤i≤3.

8. The electromagnetic transient simulation system according to claim 1, characterized in that: The initialization module is used to communicate with an external device to receive initialization parameters after acquiring the initialization start instruction transmitted by the main control module, wherein the initialization parameters include simulation parameters, component parameters and impedance matrix parameters; The simulation parameters include the total number of nodes in the electromagnetic transient simulation system, the simulation integration step, the total simulation time and the data interaction interval; The component parameters include type information of each component, the number of components of each type, node information, and input parameter information, where the input parameter information represents the management relationship between the component and the corresponding component update module.

9. The electromagnetic transient simulation system according to claim 8, characterized in that: The main control module is also used to determine whether the simulation is ended based on the simulation duration after the solution module completes the solution. If not ended, the start time of the next simulation period is determined based on the start time of the current simulation period and the simulation integration step, and at the start time of the next simulation period, a component update enable instruction is transmitted to the component update module and a solution instruction is transmitted to the solution module.

10. An electromagnetic transient simulation method, characterized in that: The electromagnetic transient simulation system applied to any one of claims 1 to 9, the method comprising: When the component update enable instruction transmitted by the main control module is obtained, the component update module obtains the historical current of each component in the simulation period according to the component parameters and the node voltage vector of the previous simulation period; After obtaining the solution instruction transmitted by the main control module, the solution module performs a solution according to the historical current of each component in the simulation period and the impedance matrix parameters to obtain the node current vector and / or node voltage vector in the simulation period; The node current vector includes the node current of each node in the simulation period, and the node voltage vector includes the node voltage of each node in the simulation period.

Images