WO2024082163A1 - Procédé de simulation transitoire électromagnétique de transformateur, dispositif informatique et support de stockage - Google Patents

Procédé de simulation transitoire électromagnétique de transformateur, dispositif informatique et support de stockage Download PDF

Info

Publication number
WO2024082163A1
WO2024082163A1 PCT/CN2022/126082 CN2022126082W WO2024082163A1 WO 2024082163 A1 WO2024082163 A1 WO 2024082163A1 CN 2022126082 W CN2022126082 W CN 2022126082W WO 2024082163 A1 WO2024082163 A1 WO 2024082163A1
Authority
WO
WIPO (PCT)
Prior art keywords
field component
time step
current time
magnetic field
electric field
Prior art date
Application number
PCT/CN2022/126082
Other languages
English (en)
Chinese (zh)
Inventor
彭庆军
施勇
邹德旭
李国友
王欣
王浩州
Original Assignee
云南电网有限责任公司电力科学研究院
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by 云南电网有限责任公司电力科学研究院 filed Critical 云南电网有限责任公司电力科学研究院
Priority to PCT/CN2022/126082 priority Critical patent/WO2024082163A1/fr
Publication of WO2024082163A1 publication Critical patent/WO2024082163A1/fr

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation

Definitions

  • the present application relates to the technical field of transformers, and in particular to a transformer electromagnetic transient simulation method, a computer device, and a storage medium.
  • Power transformers are one of the most important equipment in the power system. Once a power transformer fails, it will cause great disturbances to the power system, and even interrupt the power supply of the power system, which will greatly affect industrial production and the normal life of residents. Moreover, the cost of large power transformers is very expensive. If the protection device cannot work correctly when the transformer fails, it may cause different degrees of damage to the transformer or even burn it, which will lead to very serious economic losses. Therefore, it is very necessary to study the transient process of power transformers. By studying the working principles of various transient processes of transformers, it is effective to avoid further expansion of power accidents and damage to transformers.
  • This solution applies the finite time-domain difference method to the simulation of the electromagnetic transient process of the power transformer, adopts the vector matching method to fit the frequency characteristic curve of the transformer equipment, and combines the time-domain convolution method to fully consider the frequency-varying characteristics of the transformer equipment in the time-domain transient simulation calculation, thereby greatly improving the time-domain simulation accuracy.
  • this patent has the advantages of high simulation efficiency, high accuracy, good stability, small memory space occupation, one-time calculation of broadband transient processes, explicit calculation without large matrix inversion, and convenient acceleration of calculation through high-performance graphics computing processors, etc., providing computing support for the manufacture of transformers and the analysis of various operation and fault conditions.
  • the present application provides a transformer electromagnetic transient simulation method, the method comprising:
  • the electric field component and the magnetic field component of the transformer model are updated to obtain the electric field component of the current time step and the magnetic field component of the current time step;
  • the method further includes: if the number of iterations corresponding to the current time step is greater than an iteration number threshold, the simulation ends.
  • establishing a transformer model based on the finite time-domain difference method includes: determining a target calculation area of the finite time-domain difference method according to the structural size of the transformer; and establishing the transformer model within the target calculation area.
  • setting the time step of the transformer model includes: setting the time step of the transformer model according to the following time step formula:
  • ⁇ x, ⁇ y, ⁇ z are the minimum grid sizes of the finite time-domain difference method grid in the three orthogonal directions of X, Y, and Z, c is the propagation speed of light in the corresponding medium, and ⁇ t is the maximum time step.
  • the electric field component and the magnetic field component of the transformer model are updated based on the frequency characteristic curve to obtain the electric field component of the current time step and the magnetic field component of the current time step, including: using a vector matching method to fit the frequency-varying function corresponding to the frequency characteristic curve to obtain a first-order rational function; based on the first-order rational function, the electric field component and the magnetic field component of the transformer model are updated to obtain the electric field component of the current time step and the magnetic field component of the current time step.
  • the electric field component of the transformer model is updated based on the first-order rational function to obtain the electric field component of the current time step and the magnetic field component of the current time step, including: inputting the convolution form of the first-order rational function into the electric field update equation of the finite time-domain difference method in the transformer model to determine the electric field component of the current time step.
  • updating the magnetic field component of the transformer model based on the frequency characteristic curve to obtain the magnetic field component of the current time step includes: inputting the convolution form of the first-order rational function into the magnetic field update equation of the finite time-domain difference method in the transformer model to determine the magnetic field component of the current time step.
  • the present application provides a computer device, including a memory and a processor, wherein the memory stores a computer program, and when the computer program is executed by the processor, the processor performs the following steps:
  • the electric field component and the magnetic field component of the transformer model are updated to obtain the electric field component of the current time step and the magnetic field component of the current time step;
  • the method further includes: if the number of iterations corresponding to the current time step is greater than an iteration number threshold, the simulation ends.
  • establishing a transformer model based on the finite time-domain difference method includes: determining a target calculation area of the finite time-domain difference method according to the structural size of the transformer; and establishing a transformer model within the target calculation area.
  • setting the time step of the transformer model includes: setting the time step of the transformer model according to the following time step formula:
  • ⁇ x, ⁇ y, ⁇ z are the minimum grid sizes of the finite time-domain difference method grid in the three orthogonal directions of X, Y, and Z, c is the propagation speed of light in the corresponding medium, and ⁇ t is the maximum time step.
  • the electric field component and the magnetic field component of the transformer model are updated based on the frequency characteristic curve to obtain the electric field component of the current time step and the magnetic field component of the current time step, including: using a vector matching method to fit the frequency-varying function corresponding to the frequency characteristic curve to obtain a first-order rational function; based on the first-order rational function, the electric field component and the magnetic field component of the transformer model are updated to obtain the electric field component of the current time step and the magnetic field component of the current time step.
  • the electric field component of the transformer model is updated based on the first-order rational function to obtain the electric field component of the current time step and the magnetic field component of the current time step, including: inputting the convolution form of the first-order rational function into the electric field update equation of the finite time-domain difference method in the transformer model to determine the electric field component of the current time step.
  • updating the magnetic field component of the transformer model based on the frequency characteristic curve to obtain the magnetic field component of the current time step includes: inputting the convolution form of the first-order rational function into the magnetic field update equation of the finite time-domain difference method in the transformer model to determine the magnetic field component of the current time step.
  • the present application provides a computer-readable storage medium storing a computer program, wherein when the computer program is executed by a processor, the processor performs the following steps:
  • the electric field component and the magnetic field component of the transformer model are updated to obtain the electric field component of the current time step and the magnetic field component of the current time step;
  • the method further includes: if the number of iterations corresponding to the current time step is greater than an iteration number threshold, the simulation ends.
  • establishing a transformer model based on the finite time-domain difference method includes: determining a target calculation area of the finite time-domain difference method according to the structural size of the transformer; and establishing a transformer model within the target calculation area.
  • setting the time step of the transformer model includes: setting the time step of the transformer model according to the following time step formula:
  • ⁇ x, ⁇ y, ⁇ z are the minimum grid sizes of the finite time-domain difference method grid in the three orthogonal directions of X, Y, and Z, c is the propagation speed of light in the corresponding medium, and ⁇ t is the maximum time step.
  • the electric field component and the magnetic field component of the transformer model are updated based on the frequency characteristic curve to obtain the electric field component of the current time step and the magnetic field component of the current time step, including: using a vector matching method to fit the frequency-varying function corresponding to the frequency characteristic curve to obtain a first-order rational function; based on the first-order rational function, the electric field component and the magnetic field component of the transformer model are updated to obtain the electric field component of the current time step and the magnetic field component of the current time step.
  • the electric field component of the transformer model is updated based on the first-order rational function to obtain the electric field component of the current time step and the magnetic field component of the current time step, including: inputting the convolution form of the first-order rational function into the electric field update equation of the finite time-domain difference method in the transformer model to determine the electric field component of the current time step.
  • updating the magnetic field component of the transformer model based on the frequency characteristic curve to obtain the magnetic field component of the current time step includes: inputting the convolution form of the first-order rational function into the magnetic field update equation of the finite time-domain difference method in the transformer model to determine the magnetic field component of the current time step.
  • the present invention provides a transformer electromagnetic transient simulation method, the method comprising: establishing a transformer model based on a finite time-domain difference method, and setting a time step for the transformer model; obtaining a frequency characteristic curve of the transformer; based on the frequency characteristic curve, updating the electric field component and the magnetic field component of the transformer model, and obtaining the electric field component of the current time step and the magnetic field component of the current time step; if the electric field component of the current time step or the magnetic field component of the current time step converges, the simulation ends.
  • the present application establishes a transformer model based on a finite time-domain difference method, does not require complex and tedious calculations and field tests, and can greatly improve the calculation efficiency; and the present application updates the electric field component and the magnetic field component of the transformer model based on the frequency characteristic curve, and can fully consider the frequency-varying transformer electromagnetic characteristic parameters in the time domain, greatly improving the simulation accuracy.
  • FIG1 is a schematic diagram of a flow chart of a transformer electromagnetic transient simulation method provided in an embodiment of the present application
  • FIG2 is a schematic diagram of an electric field vector in a grid provided in an embodiment of the present application.
  • FIG3 is a schematic diagram of a magnetic field vector in a grid provided in an embodiment of the present application.
  • FIG4 is a schematic diagram of an electromagnetic field vector in a grid provided in an embodiment of the present application.
  • FIG5 is a magnetic field frequency characteristic curve of a long straight conductor with a rectangular cross section and a vector matching method fitting curve thereof provided in an embodiment of the present application;
  • FIG6 is a current frequency characteristic curve of a long straight conductor with a rectangular cross section and a vector matching method fitting curve thereof provided in an embodiment of the present application;
  • FIG7 is a schematic structural diagram of a transformer electromagnetic transient simulation device provided in an embodiment of the present application.
  • FIG8 is a schematic diagram of the structure of a computer device provided in an embodiment of the present application.
  • This application establishes a transformer model based on the finite time-domain difference method, and sets the grid size, time step, and material properties of the transformer model; through experimental or theoretical analysis, the material parameter characteristic curve of the transformer that changes nonlinearly with frequency is obtained, and the characteristic curve is fitted into a set of rational functions based on frequency domain analysis; based on the rational function, the Laplace transform technology is applied to convert it into a time domain update equation, and replace the classic update equations of the electric field component and the magnetic field component that need to consider the frequency-varying material characteristics in the transformer model; after the iterative calculation starts, the electric field vector and the magnetic field vector in the transformer model are calculated based on the Maxwell equations or based on the rational function time domain update equation to obtain the electric field component of the current time step and the magnetic field component of the current time step; after multiple iterative calculations, if the electric field component of the current time step or the magnetic field component of the current time
  • the present application proposes a transformer electromagnetic transient simulation method, as shown in FIG1 , which is a flow chart of a transformer electromagnetic transient simulation method provided in an embodiment of the present application, and the method includes:
  • Step 101 Establish a transformer model based on a finite difference time domain method, and set a time step for the transformer model.
  • the finite-difference time-domain method is a global discrete time-domain simulation algorithm. Its calculation area includes not only all simulated models, but also the areas between simulated objects and their adjacent areas.
  • the transformer model is completed based on the model building software. Specifically, the transformer model is built in the calculation area of the finite time-domain difference method according to the transformer parameters.
  • the establishment of a transformer model based on the finite time-domain difference method includes: determining a target calculation area of the finite time-domain difference method according to the structure size of the transformer; and establishing the transformer model in the target calculation area.
  • the calculation area of the finite time difference method needs to be gridded first, that is, in this application, the calculation area of the finite time difference method is composed of multiple FDTD orthogonal grid arrangements.
  • the calculation domain of the finite time-domain difference method is discretized into a set of parallelepiped grids through FDTD orthogonal grids, and the electromagnetic field in each grid is set to be uniformly distributed.
  • Figures 2-3 Figure 2 is a schematic diagram of an electric field vector in a grid provided in an embodiment of the present application;
  • Figure 3 is a schematic diagram of a magnetic field vector in a grid provided in an embodiment of the present application.
  • the lower left vertex of each parallelepiped grid is taken as the origin, and the three edges connected to the origin are respectively defined to point to the three orthogonal directions of XYZ.
  • the electric field vector Ex, Ey, and Ez as shown in Figure 3, the lower left vertex of each parallelepiped grid is taken as the origin, and the three faces connected to the origin are respectively defined to point to the three orthogonal directions of XYZ. , Hy, Hz.
  • the electric field and magnetic field vectors in each direction need to set the corresponding material parameters according to the relative spatial position, including conductivity ⁇ , dielectric constant ⁇ , and magnetic permeability ⁇ .
  • Material parameters refer to the parameters of transformer electrical materials.
  • the grid size can be encrypted in areas where the electromagnetic field changes dramatically, such as the air-conductor interface, the conductor-dielectric interface, etc.; large-size grids can be used in areas where the electromagnetic field changes slowly, such as the inside of the dielectric.
  • the simulated objects that may contain fine structures can be discretized to obtain simulation accuracy that meets the requirements; on the other hand, it can also avoid unnecessary computational burdens caused by overly dense grids, reducing the overall computational implementation efficiency of the solution.
  • the electromagnetic field vectors are spaced apart by half a spatial step (i.e., the grid size), and the electromagnetic vectors surround and surround each other, that is, the electric field vector in a certain direction is surrounded by four magnetic field vectors, and vice versa.
  • the electromagnetic field vectors are also spaced apart by half a time step in time, that is, the overall electric field vector and the overall magnetic field vector are always 0.5 ⁇ t apart.
  • Figure 4 is a schematic diagram of electromagnetic field vectors in a grid provided in an embodiment of the present application. It should be noted that the above space-time characteristics meet the characteristics of solving Maxwell's discrete equations, because alternating step-by-step solutions of electric field vectors and magnetic field vectors can be achieved.
  • the target calculation area of the finite time-domain difference method is determined according to the structural size of the transformer, that is, the model construction area of the transformer model is determined.
  • the structural size of the transformer can be the transformer body and the adjacent area (which can be 50% of the transformer size), and the target calculation area of the finite time-domain difference method is not less than the structural size of the transformer.
  • setting the time step of the transformer model includes: setting the time step of the transformer model according to the following time step formula:
  • ⁇ x, ⁇ y, ⁇ z are the minimum grid sizes of the finite time-domain difference method grid in the three orthogonal directions of X, Y, and Z, c is the propagation speed of light in the corresponding medium, and ⁇ t is the maximum time step.
  • the selection of the maximum time step is determined by the minimum FDTD discrete grid size and should satisfy the Courant-Friedrich-Levy (CFL) criterion to prevent data divergence, oscillation, non-convergence and other problems that may occur in time domain calculations, that is,
  • ⁇ x, ⁇ y, ⁇ z are the minimum grid sizes of the FDTD grid in the three orthogonal directions of X, Y, and Z, and c is the propagation speed of light in the corresponding medium.
  • a transformer model is established based on the finite time-domain difference method, that is, the transformer body and the adjacent area (generally 50% of the transformer size) are modeled in the FDTD calculation area, and the transformer and the adjacent area are discretized using the FDTD grid.
  • the FDTD classical update equation is used to solve the electromagnetic field distribution in the time domain.
  • the frequency characteristic curve of the transformer is first obtained, and then the electromagnetic field distribution is solved based on the frequency characteristic curve of the transformer and the classical FDTD update equation. The scheme of nonlinear change of material parameters with frequency is described below.
  • Step 102 Obtain a frequency characteristic curve of the transformer.
  • the frequency-varying characteristic curve of the key structure of the power transformer is obtained through experiments, theoretical formula calculations or numerical algorithms, and the frequency range must cover all frequency bands of concern in the electromagnetic transient process.
  • the theoretical formula can be Bessel Function, etc.
  • the numerical algorithm can be Method of Moment (MoM), Partial Element Equivalent Circuit (PEEC), etc.
  • the frequency characteristic curve of the key components of the transformer (such as copper loss considering the frequency-varying characteristics) is not only related to the material, but also to the geometric structure size of the key components, so even for the same material, they need to be obtained separately. Since the internal structure of the transformer is relatively complex and the accuracy of the detection equipment is limited, in order to improve the accuracy of the frequency characteristic curve, this application uses theoretical formulas or numerical algorithms to obtain the corresponding characteristic curves.
  • the Bessel function can be used to obtain the current density distribution of the circular conductor, and then derive the frequency-varying loss.
  • the form of the Bessel function is as follows:
  • Z is impedance
  • R d mrd
  • m 2 j ⁇ ( ⁇ +j ⁇ )
  • r d is the conductor radius
  • ⁇ is the magnetic permeability
  • ⁇ is the dielectric constant
  • ⁇ is the conductivity
  • j represents a complex number
  • ⁇ is the angular frequency
  • I 0 and I 1 are the 0th and 1st order first kind modified form Bessel functions respectively.
  • the two-dimensional MoM algorithm can be used to discretize the two-dimensional cross section, obtain the current distribution on each discrete unit, and then obtain the overall frequency-dependent loss of the long straight conductor.
  • the three-dimensional PEEC algorithm can be used to discretize the conductor in space according to the accuracy requirements, and each discrete unit is equivalent to a circuit element, and then combined with the circuit solving algorithm to obtain its frequency characteristic curve.
  • Step 103 based on the frequency characteristic curve, update the electric field component and the magnetic field component of the transformer model to obtain the electric field component of the current time step and the magnetic field component of the current time step.
  • the electric field component and the magnetic field component of the transformer model are updated to obtain the electric field component of the current time step and the magnetic field component of the current time step, including: using a vector matching method to fit the frequency-varying function corresponding to the frequency characteristic curve to obtain a first-order rational function; based on the first-order rational function, the electric field component and the magnetic field component of the transformer model are updated to obtain the electric field component of the current time step and the magnetic field component of the current time step.
  • the corresponding frequency-variant function is fitted using the vector fitting technique (VFT), and the complex frequency-variant function is fitted into a first-order rational function, and then the rational function is substituted into the FDTD solution equation in the form of convolution to achieve the purpose of considering the frequency-variant parameters in the time domain calculation.
  • VFT vector fitting technique
  • the rational function is substituted into the FDTD solution equation in the form of convolution to achieve the purpose of considering the frequency-variant parameters in the time domain calculation.
  • the vector matching method has the advantages of high stability, few iterations, and fast convergence.
  • the frequency characteristic curve of the key components of the transformer can be approximately fitted in the frequency domain as a first-order rational function and.
  • the first-order rational function of the frequency characteristic curve fitted by the vector matching method is (the frequency-dependent impedance Z can be fitted in the complex frequency domain as):
  • s represents the complex frequency domain
  • d is the DC component
  • h is the inductive component
  • cm is the residue
  • am is the pole
  • N is the vector matching order. It has been verified that the relative error between the frequency characteristic curve fitted by the vector matching method and the frequency characteristic curve of the original transformer key component is less than 0.1%.
  • Figure 5 is a magnetic field frequency characteristic curve of a long straight conductor with a rectangular cross-section and its vector matching method fitting curve provided in an embodiment of the present application
  • Figure 6 is a current frequency characteristic curve of a long straight conductor with a rectangular cross-section and its vector matching method fitting curve provided in an embodiment of the present application. It should be noted that as many frequency points as possible should be uniformly selected within each frequency order of magnitude of the original characteristic frequency curve to improve the accuracy of curve fitting.
  • the updating of the electric field component of the transformer model based on the first-order rational function to obtain the electric field component of the current time step and the magnetic field component of the current time step includes: inputting the convolution form of the first-order rational function into the electric field update equation of the finite time-domain difference method in the transformer model to determine the electric field component of the current time step.
  • the updating of the magnetic field component of the transformer model based on the frequency characteristic curve to obtain the magnetic field component of the current time step includes: inputting the convolution form of the first-order rational function into the magnetic field update equation of the finite time-domain difference method in the transformer model to determine the magnetic field component of the current time step.
  • the classic FDTD update equation is applied to iteratively solve the electric field vector in the calculation area.
  • the first-order rational function needs to be substituted into the FDTD electric field update equation in the form of convolution, so as to fully consider the influence of the electric frequency-dependent characteristic parameters on the electromagnetic transient simulation results.
  • the iterative process involves the electric field vector value of the previous time step and the four magnetic field vectors surrounding the electric field vector.
  • the specific update equation is as follows:
  • Ey , and Ez are the electric field vectors in three orthogonal directions
  • i, j, and k are the position numbers of the electric field vector based on the FDTD grid numbering
  • q represents the number of time steps
  • ⁇ and ⁇ represent the equivalent conductivity and dielectric constant in the corresponding space, respectively
  • ⁇ t represents the FDTD time step.
  • the frequency characteristic curve of the impedance can be fitted into the form of a first-order rational function, and then the first-order rational function can be transformed into a time domain convolution form by an inverse Laplace transform, and then the electric field vector update equation in the time domain convolution form is replaced by the electric field vector update equation at the corresponding position in formula (4).
  • the update equation of the electric field vector E at the spatial position of the frequency-dependent characteristic material can be replaced by the following equation:
  • Kp d + h / ⁇ t
  • Kn -h / ⁇ t
  • Iq the current
  • d the DC component
  • h the inductive component
  • the time step can be appropriately reduced when combined with the time domain convolution calculation.
  • Each iteration needs to reduce the time step to between 0.8 and 0.9 times the time step of the previous iteration. Since such convolution calculations are performed in local areas, they do not affect the parallel computing capabilities of the FDTD program.
  • the electromagnetic transient time domain simulation of power transformers can still use high-performance graphics processing units (GPUs) to improve computing efficiency.
  • the classic FDTD update equation is used to iteratively solve the magnetic field vector in the calculation area.
  • the first-order rational function mentioned above needs to be substituted into the FDTD magnetic field update equation in the form of convolution, so as to fully consider the influence of the magnetic frequency-variable characteristic parameters on the electromagnetic transient simulation results.
  • the magnetic field vector of the entire calculation domain is iteratively calculated using the classic FDTD magnetic field vector update equation.
  • the specific equation is expressed as follows:
  • ⁇ and ⁇ m are magnetic permeability and permeability, and the permeability is generally set to 0.
  • Step 104 If the electric field component of the current time step or the magnetic field component of the current time step converges, the simulation ends.
  • the electric field component and the magnetic field component are iteratively calculated until the electric field component or the magnetic field component of the current time step converges, and the simulation ends.
  • the simulation ends.
  • the present invention provides a transformer electromagnetic transient simulation method, the method comprising: establishing a transformer model based on a finite time-domain difference method, and setting a time step for the transformer model; obtaining a frequency characteristic curve of the transformer; based on the frequency characteristic curve, updating the electric field component and the magnetic field component of the transformer model, and obtaining the electric field component of the current time step and the magnetic field component of the current time step; if the electric field component of the current time step or the magnetic field component of the current time step converges, the simulation ends.
  • the present application establishes a transformer model based on a finite time-domain difference method, does not require complex and tedious calculations and field tests, and can greatly improve the calculation efficiency; and the present application updates the electric field component and the magnetic field component of the transformer model based on the frequency characteristic curve, and can fully consider the frequency-varying transformer electromagnetic characteristic parameters in the time domain, greatly improving the simulation accuracy.
  • the present application provides a transformer electromagnetic transient simulation device, as shown in FIG7 , which is a schematic diagram of the structure of a transformer electromagnetic transient simulation device provided in an embodiment of the present application, the device comprising:
  • the model building module 701 is used to build a transformer model based on the finite time-domain difference method and set the time step of the transformer model.
  • the characteristic curve acquisition module 702 is used to acquire the frequency characteristic curve of the transformer.
  • the first simulation module 703 is used to update the electric field component and the magnetic field component of the transformer model based on the frequency characteristic curve to obtain the electric field component of the current time step and the magnetic field component of the current time step.
  • the second simulation module 704 is used to end the simulation if the electric field component of the current time step or the magnetic field component of the current time step converges.
  • the computer device may be a transformer electromagnetic transient simulation device, or a terminal or server connected to a transformer electromagnetic transient simulation device.
  • the computer device includes a processor, a memory, and a network interface connected via a system bus.
  • the memory includes a non-volatile storage medium and an internal memory.
  • the non-volatile storage medium of the computer device stores an operating system and may also store a computer program.
  • the processor may implement a transformer electromagnetic transient simulation method.
  • the internal memory may also store a computer program.
  • the processor may execute a transformer electromagnetic transient simulation method.
  • the network interface is used to communicate with an external device.
  • the network interface is used to communicate with an external device.
  • FIG8 is only a block diagram of a partial structure related to the present application scheme, and does not constitute a limitation on the computer device to which the present application scheme is applied.
  • the specific computer device may include more or fewer components than shown in the figure, or combine certain components, or have a different arrangement of components.
  • a transformer electromagnetic transient simulation method provided by the present application can be implemented in the form of a computer program, and the computer program can be run on a computer device as shown in FIG8.
  • the memory of the computer device can store various program templates constituting the transformer electromagnetic transient simulation device. For example, a model building module 701, a characteristic curve acquisition module 702, a first simulation module 703, and a second simulation module 704.
  • a computer device includes a memory and a processor, wherein the memory stores a computer program, and when the computer program is executed by the processor, the processor executes the following steps: establishing a transformer model based on a finite time-domain difference method, and setting a time step for the transformer model; obtaining a frequency characteristic curve of the transformer; updating an electric field component and a magnetic field component of the transformer model based on the frequency characteristic curve to obtain an electric field component of a current time step and a magnetic field component of a current time step; if the electric field component of the current time step or the magnetic field component of the current time step converges, the simulation ends.
  • the processor when the computer program is executed by the processor, the processor further executes the following steps: if the number of iterations corresponding to the current time step is greater than an iteration number threshold, the simulation ends.
  • the establishment of a transformer model based on the finite time-domain difference method includes: determining a target calculation area of the finite time-domain difference method according to the structure size of the transformer; and establishing the transformer model in the target calculation area.
  • setting the time step of the transformer model includes: setting the time step of the transformer model according to the following time step formula:
  • ⁇ x, ⁇ y, ⁇ z are the minimum grid sizes of the finite time-domain difference method grid in the three orthogonal directions of X, Y, and Z, c is the propagation speed of light in the corresponding medium, and ⁇ t is the maximum time step.
  • the updating of the electric field component and the magnetic field component of the transformer model based on the frequency characteristic curve to obtain the electric field component of the current time step and the magnetic field component of the current time step includes: fitting the frequency-varying function corresponding to the frequency characteristic curve using a vector matching method to obtain a first-order rational function; updating the electric field component and the magnetic field component of the transformer model based on the first-order rational function to obtain the electric field component of the current time step and the magnetic field component of the current time step.
  • the updating of the electric field component of the transformer model based on the first-order rational function to obtain the electric field component of the current time step and the magnetic field component of the current time step includes: inputting the convolution form of the first-order rational function into the electric field update equation of the finite time-domain difference method in the transformer model to determine the electric field component of the current time step.
  • updating the magnetic field component of the transformer model based on the frequency characteristic curve to obtain the magnetic field component of the current time step includes: inputting the convolution form of the first-order rational function into the magnetic field update equation of the finite time-domain difference method in the transformer model to determine the magnetic field component of the current time step.
  • a computer-readable storage medium stores a computer program.
  • the processor When the computer program is executed by a processor, the processor performs the following steps: establishing a transformer model based on a finite time-domain difference method, and setting a time step for the transformer model; obtaining a frequency characteristic curve of the transformer; updating an electric field component and a magnetic field component of the transformer model based on the frequency characteristic curve to obtain an electric field component of a current time step and a magnetic field component of a current time step; if the electric field component of the current time step or the magnetic field component of the current time step converges, the simulation ends.
  • the processor when the computer program is executed by the processor, the processor further executes the following steps: if the number of iterations corresponding to the current time step is greater than an iteration number threshold, the simulation ends.
  • the establishment of a transformer model based on the finite time-domain difference method includes: determining a target calculation area of the finite time-domain difference method according to the structure size of the transformer; and establishing the transformer model in the target calculation area.
  • setting the time step of the transformer model includes: setting the time step of the transformer model according to the following time step formula:
  • ⁇ x, ⁇ y, ⁇ z are the minimum grid sizes of the finite time-domain difference method grid in the three orthogonal directions of X, Y, and Z, c is the propagation speed of light in the corresponding medium, and ⁇ t is the maximum time step.
  • the updating of the electric field component and the magnetic field component of the transformer model based on the frequency characteristic curve to obtain the electric field component of the current time step and the magnetic field component of the current time step includes: fitting the frequency-varying function corresponding to the frequency characteristic curve using a vector matching method to obtain a first-order rational function; updating the electric field component and the magnetic field component of the transformer model based on the first-order rational function to obtain the electric field component of the current time step and the magnetic field component of the current time step.
  • the updating of the electric field component of the transformer model based on the first-order rational function to obtain the electric field component of the current time step and the magnetic field component of the current time step includes: inputting the convolution form of the first-order rational function into the electric field update equation of the finite time-domain difference method in the transformer model to determine the electric field component of the current time step.
  • updating the magnetic field component of the transformer model based on the frequency characteristic curve to obtain the magnetic field component of the current time step includes: inputting the convolution form of the first-order rational function into the magnetic field update equation of the finite time-domain difference method in the transformer model to determine the magnetic field component of the current time step.
  • the storage medium can be a disk, an optical disk, a read-only memory (ROM) or a random access memory (RAM), etc.

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • Supply And Distribution Of Alternating Current (AREA)

Abstract

La présente demande concerne un procédé de simulation transitoire électromagnétique d'un transformateur. Le procédé consiste à : construire un modèle de transformateur sur la base d'une méthode du domaine temporel aux différences finies, et définir l'étape temporelle du modèle de transformateur ; acquérir une courbe de réponse en fréquence du transformateur ; mettre à jour une composante de champ électrique et une composante de champ magnétique du modèle de transformateur sur la base de la courbe de réponse en fréquence pour obtenir une composante de champ électrique de l'étape temporelle actuelle et une composante de champ magnétique de l'étape temporelle actuelle ; et si la composante de champ électrique de l'étape temporelle actuelle ou la composante de champ magnétique de l'étape temporelle actuelle converge, mettre fin à la la simulation. Selon la présente demande, le modèle de transformateur est établi sur la base de la méthode du domaine temporel aux différences finies, sans qu'un calcul complexe et fastidieux et des tests sur site soient nécessaires, de sorte que l'efficacité de calcul peut être considérablement améliorée. De plus, la composante de champ électrique et la composante de champ magnétique du modèle de transformateur sont mises à jour sur la base de la courbe de réponse en fréquence. De cette manière, des paramètres de réponse électromagnétique d'un transformateur dépendant de la fréquence peuvent être pleinement pris en considération dans le domaine temporel, ce qui permet d'améliorer considérablement la précision de simulation.
PCT/CN2022/126082 2022-10-19 2022-10-19 Procédé de simulation transitoire électromagnétique de transformateur, dispositif informatique et support de stockage WO2024082163A1 (fr)

Priority Applications (1)

Application Number Priority Date Filing Date Title
PCT/CN2022/126082 WO2024082163A1 (fr) 2022-10-19 2022-10-19 Procédé de simulation transitoire électromagnétique de transformateur, dispositif informatique et support de stockage

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
PCT/CN2022/126082 WO2024082163A1 (fr) 2022-10-19 2022-10-19 Procédé de simulation transitoire électromagnétique de transformateur, dispositif informatique et support de stockage

Publications (1)

Publication Number Publication Date
WO2024082163A1 true WO2024082163A1 (fr) 2024-04-25

Family

ID=90736637

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/CN2022/126082 WO2024082163A1 (fr) 2022-10-19 2022-10-19 Procédé de simulation transitoire électromagnétique de transformateur, dispositif informatique et support de stockage

Country Status (1)

Country Link
WO (1) WO2024082163A1 (fr)

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102542111A (zh) * 2011-12-31 2012-07-04 天津大学 一种提高变压器电磁仿真效率的方法
CN104298822A (zh) * 2014-10-09 2015-01-21 天津大学 一种考虑非线性的变电站电磁暂态实时仿真方法
CN105205238A (zh) * 2015-09-14 2015-12-30 华南理工大学 基于pscad/emtdc软件的变压器电磁暂态仿真建模方法
CN107037312A (zh) * 2016-11-28 2017-08-11 国家电网公司 一种用于扫频阻抗法的变压器绕组仿真模型建立方法
WO2020088245A1 (fr) * 2018-10-31 2020-05-07 邓宏伟 Algorithme de simulation transitoire électromagnétique et électromécanique basé sur un algorithme direct
CN113283129A (zh) * 2021-04-13 2021-08-20 上海电力大学 基于时频域分析变压器局部放电信号接收器位置设计方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102542111A (zh) * 2011-12-31 2012-07-04 天津大学 一种提高变压器电磁仿真效率的方法
CN104298822A (zh) * 2014-10-09 2015-01-21 天津大学 一种考虑非线性的变电站电磁暂态实时仿真方法
CN105205238A (zh) * 2015-09-14 2015-12-30 华南理工大学 基于pscad/emtdc软件的变压器电磁暂态仿真建模方法
CN107037312A (zh) * 2016-11-28 2017-08-11 国家电网公司 一种用于扫频阻抗法的变压器绕组仿真模型建立方法
WO2020088245A1 (fr) * 2018-10-31 2020-05-07 邓宏伟 Algorithme de simulation transitoire électromagnétique et électromécanique basé sur un algorithme direct
CN113283129A (zh) * 2021-04-13 2021-08-20 上海电力大学 基于时频域分析变压器局部放电信号接收器位置设计方法

Similar Documents

Publication Publication Date Title
Edelvik A new technique for accurate and stable modeling of arbitrarily oriented thin wires in the FDTD method
Ren et al. A new 3-D nonspurious discontinuous Galerkin spectral element time-domain (DG-SETD) method for Maxwell’s equations
Edelvik et al. An unconditionally stable subcell model for arbitrarily oriented thin wires in the FETD method
Kettner et al. On the properties of the compound nodal admittance matrix of polyphase power systems
Yan et al. An unsymmetric FDTD subgridding algorithm with unconditional stability
Leone et al. Efficient broadband circuit-modeling approach for parallel-plane structures of arbitrary shape
CN116187258B (zh) 量子芯片版图的仿真方法、装置、计算设备及存储介质
Chen et al. A memory-efficient implementation of perfectly matched layer with smoothly varying coefficients in discontinuous Galerkin time-domain method
WO2024082163A1 (fr) Procédé de simulation transitoire électromagnétique de transformateur, dispositif informatique et support de stockage
Li et al. Transient analysis of lumped circuit networks-loaded thin wires by DGTD method
Sharma et al. A complete surface integral method for broadband modeling of 3D interconnects in stratified media
Hackl et al. Multifunction approach and specialized numerical integration algorithms for fast inductance evaluations in nonorthogonal PEEC systems
CN116151073A (zh) 一种基于混合有限元法的集成电路仿真方法及其装置
Schuhmann et al. The nonorthogonal finite integration technique applied to 2D-and 3D-eigenvalue problems
Law et al. A nonlinear relaxation formulation of the p-curl problem modelling high-temperature superconductors: A modified Yee's scheme
JP3703812B2 (ja) Ｆｄｔｄ法を用いた電磁界解析方法、電磁界解析における媒質表現方法、シミュレーション装置、及びプログラム
Nentchev Numerical analysis and simulation in microelectronics by vector finite elements
Aiello et al. A GMRES iterative solution of FEM‐BEM global systems in skin effect problems
Li et al. Quasi‐static numerical modeling of miniature RF circuits based on lumped equivalent circuits
Grylonakis et al. A note on solving the generalized Dirichlet to Neumann map on irregular polygons using Generic Factored Approximate Sparse Inverses
Zygiridis et al. Development of optimized operators based on spherical-harmonic expansions for 3D FDTD schemes
Hackl et al. Efficient simulation of magnetic components using the MagPEEC-method
Bokil et al. A dispersion minimized mimetic method for cold plasma
CN116882219B (zh) 一种基于dgtd与fdtd的场线耦合方法
Ijjeh et al. Frequency-Domain TLM Method with Cartesian Block Meshing