CN114611061B - Method and device for calculating closing transient current of reactor and electronic equipment - Google Patents
Method and device for calculating closing transient current of reactor and electronic equipment Download PDFInfo
- Publication number
- CN114611061B CN114611061B CN202210411344.6A CN202210411344A CN114611061B CN 114611061 B CN114611061 B CN 114611061B CN 202210411344 A CN202210411344 A CN 202210411344A CN 114611061 B CN114611061 B CN 114611061B
- Authority
- CN
- China
- Prior art keywords
- magnetic
- reactor
- target
- target reactor
- density
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 230000001052 transient effect Effects 0.000 title claims abstract description 56
- 238000000034 method Methods 0.000 title claims abstract description 37
- XEEYBQQBJWHFJM-UHFFFAOYSA-N Iron Chemical group [Fe] XEEYBQQBJWHFJM-UHFFFAOYSA-N 0.000 claims abstract description 123
- 230000004907 flux Effects 0.000 claims abstract description 69
- 230000035699 permeability Effects 0.000 claims abstract description 47
- 238000004364 calculation method Methods 0.000 claims abstract description 35
- 229910000976 Electrical steel Inorganic materials 0.000 claims abstract description 14
- 239000000463 material Substances 0.000 claims abstract description 12
- 230000006870 function Effects 0.000 claims description 54
- 239000011162 core material Substances 0.000 description 13
- 238000004891 communication Methods 0.000 description 12
- 230000008569 process Effects 0.000 description 6
- 230000009471 action Effects 0.000 description 3
- 230000008859 change Effects 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 229920006395 saturated elastomer Polymers 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- General Physics & Mathematics (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Theoretical Computer Science (AREA)
- Operations Research (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Regulation Of General Use Transformers (AREA)
Abstract
The invention provides a method and a device for calculating a closing transient current of a reactor and electronic equipment, wherein target parameters of a target reactor are obtained; constructing a voltage balance equation of the target reactor based on the target parameters, and solving the voltage balance equation to obtain a function relation of magnetic flux of the target reactor along with time variation; based on the material characteristics of the silicon steel sheet of the target reactor, obtaining a relation curve between the relative magnetic permeability and the magnetic density of the target reactor; solving a magnetic circuit equivalent equation of the target reactor to obtain a function relation of the magnetic flux leakage density of the target reactor changing along with time; and solving an objective function based on a relation curve between the relative magnetic permeability and the magnetic density of the objective reactor to obtain the magnetic densities of the iron cores at different moments, and then obtaining the closing transient current values at different moments corresponding to the magnetic densities of the iron cores at different moments, thereby realizing the accurate calculation of the closing transient current of the reactor.
Description
Technical Field
The invention relates to the field of power transformer design, in particular to a method and a device for calculating a reactor closing transient current and electronic equipment.
Background
When the reactor is switched on and put into a power grid, due to saturation of magnetic flux of an iron core of the reactor and nonlinear characteristics of iron core materials, switching-on transient current with quite large amplitude can be generated at the moment of switching on, and therefore differential protection misoperation of the reactor can be caused. For the reactor product itself, because the current of the inductance coil cannot be suddenly changed and the surge current is required to bear the surge current, if the surge current protection fixed value of the reactor is selected to be too low, the surge current of the reactor may exceed the surge current protection fixed value, so that problems are caused to the reactor, such as a magnetic field established by the surge current may generate larger force on the coil, and meanwhile, the surge current causes problems such as impedance reduction of the reactor, and the reactor cannot normally operate. In addition, the switching transient current contains a plurality of harmonic components and direct current components, which can reduce the power supply quality of the power system, and higher harmonics in the inrush current have extremely strong damage to sensitive power electronic devices connected to the power system.
At present, a method for calculating the closing transient current of a transformer is more studied, an estimation method is mainly adopted for calculating the closing transient current of a reactor with an air gap core, and the method does not have more accurate calculation. Unlike transformer core structure, the core column of the reactor with air gap is formed by the core cake and the air gap which are staggered along the axial direction of the core column, the magnetic flux path passing through the core column is along the direction vertical to the section of the core cake, and part of the magnetic flux flows out of the outer surface of the core, bypasses the air gap, flows to the outer surface of the core, and then enters the core, so that the magnetic circuit distribution is complex, the magnetic circuit distribution phase difference is larger after the core column is saturated, and the accurate calculation of the transient current of the switch is very difficult.
Disclosure of Invention
In view of the above, the embodiment of the invention provides a method and a device for calculating a reactor switching-on transient current and electronic equipment, so as to realize the calculation of the reactor switching-on transient current.
In order to achieve the above object, the embodiment of the present invention provides the following technical solutions:
a method for calculating the closing transient current of a reactor comprises the following steps:
acquiring target parameters of a target reactor;
constructing a voltage balance equation of the target reactor based on the target parameters, and solving the voltage balance equation to obtain a function relation of magnetic flux of the target reactor along with time variation;
based on the material characteristics of the silicon steel sheet of the target reactor, obtaining a relation curve between the relative magnetic permeability and the magnetic density of the target reactor;
Solving a magnetic circuit equivalent equation of the target reactor to obtain a function relation of the magnetic flux leakage density of the target reactor changing along with time;
solving an objective function based on a relation curve between relative magnetic permeability and magnetic density of a target reactor to obtain magnetic densities of iron cores at different moments, wherein the objective function is that Wherein, B c (t) is the magnetic flux density of the iron core at time t, S c is the effective cross-sectional area of the iron core of the target reactor, Σδ is the total length of the air gap of the iron core of the target reactor, l is the window height of the target reactor, μ r (t) is the relative magnetic permeability of the target reactor at time t, S 0 is the effective cross-sectional area of the magnetic leakage of the target reactor, and Φ (t) is the magnetic flux of the target reactor at time t;
and acquiring closing transient current values at different moments corresponding to the magnetic densities of the iron cores at different moments.
Optionally, in the foregoing solution, the target parameter includes:
the direct current resistance, the main magnetic flux effective area and the leakage magnetic flux effective area of the target reactor.
Optionally, in the above scheme, the functional relationship of the magnetic flux over time is:
Wherein, phi is the closing angle of the target reactor, R is the direct current resistance of the target reactor, L is the inductance of the target reactor, phi m is the average magnetic flux peak value of the target reactor, phi (t) is the magnetic flux peak value of the target reactor at the moment t, and omega is the angular frequency of the voltage applied to the target reactor
Optionally, in the above solution, the relationship curve between the relative permeability and the magnetic density is:
Mu r＝f(BC), wherein mu r is the relative magnetic permeability of the target reactor, and B c is the magnetic density of the iron core of the target reactor.
Optionally, in the above solution, the functional relationship of the magnetic flux leakage density over time is:
The reactor window height of the target reactor is l, the flux leakage density of the target reactor at the moment t is B 0 (t), the sigma delta is the total length of an iron core air gap of the target reactor, and the relative magnetic permeability of the target reactor at the moment t is mu r (t).
Optionally, in the above solution, the solving the objective function based on a relation curve between relative permeability and magnetic density of the objective reactor to obtain magnetic densities of the iron cores at different moments includes:
configuring magnetic density of a primary iron core corresponding to each moment;
Substituting the magnetic densities of the primary selected iron cores corresponding to the moments into the objective function, and judging whether the value of F (t) is in a preset range or not;
When the value of F (t) is not in the preset range, the primary iron core is more magnetic dense until the value of F (t) corresponding to the primary iron core magnetic dense is in the preset range;
when the value of F (t) is in a preset range, taking the magnetic density of the primary selected iron core as the magnetic density of the iron core at the corresponding moment;
and obtaining the magnetic density of the iron core at each moment.
Optionally, in the above solution, the obtaining closing transient current values at different moments corresponding to the magnetic densities of the iron core at different moments includes:
substituting magnetic densities of iron cores at different moments into formula And calculating to obtain closing transient current values I (t) at different moments corresponding to the magnetic densities of the iron cores at different moments, wherein Sigma sigma refers to the total length of an air gap of the target reactor, and N is the number of turns of the reactor of the target reactor.
A reactor closing transient current calculation device, comprising:
a parameter acquisition unit for acquiring a target parameter of a target reactor;
The first calculation unit is used for constructing a voltage balance equation of the target reactor based on the target parameters, and solving the voltage balance equation to obtain a function relation of magnetic flux of the target reactor along with time variation;
The second calculation unit is used for obtaining a relation curve between the relative magnetic permeability and the magnetic density of the target reactor based on the material characteristics of the silicon steel sheet of the target reactor;
The third calculation unit is used for solving a magnetic circuit equivalent equation of the target reactor to obtain a function relation of the magnetic flux leakage density of the target reactor changing along with time;
a fourth calculation unit, configured to solve an objective function based on a relationship curve between relative permeability and magnetic density of the objective reactor, to obtain magnetic densities of the iron cores at different times, where the objective function is Wherein, B c (t) is the magnetic flux density of the iron core at time t, S c is the effective cross-sectional area of the iron core of the target reactor, Σδ is the total length of the air gap of the iron core of the target reactor, l is the window height of the target reactor, μ r (t) is the relative magnetic permeability of the target reactor at time t, S 0 is the effective cross-sectional area of the magnetic leakage of the target reactor, and Φ (t) is the magnetic flux of the target reactor at time t;
And the fifth calculation unit is used for obtaining closing transient current values at different moments corresponding to the magnetic densities of the iron cores at different moments.
An electronic device, comprising:
a memory and a processor; the memory stores a program adapted for execution by the processor, the program for:
acquiring target parameters of a target reactor;
constructing a voltage balance equation of the target reactor based on the target parameters, and solving the voltage balance equation to obtain a function relation of magnetic flux of the target reactor along with time variation;
based on the material characteristics of the silicon steel sheet of the target reactor, obtaining a relation curve between the relative magnetic permeability and the magnetic density of the target reactor;
Solving a magnetic circuit equivalent equation of the target reactor to obtain a function relation of the magnetic flux leakage density of the target reactor changing along with time;
solving an objective function based on a relation curve between relative magnetic permeability and magnetic density of a target reactor to obtain magnetic densities of iron cores at different moments, wherein the objective function is that Wherein, B c (t) is the magnetic flux density of the iron core at time t, S c is the effective cross-sectional area of the iron core of the target reactor, Σδ is the total length of the air gap of the iron core of the target reactor, l is the window height of the target reactor, μ r (t) is the relative magnetic permeability of the target reactor at time t, S 0 is the effective cross-sectional area of the magnetic leakage of the target reactor, and Φ (t) is the magnetic flux of the target reactor at time t;
and acquiring closing transient current values at different moments corresponding to the magnetic densities of the iron cores at different moments.
Optionally, in the electronic device, the electronic device is a computer or a PC.
Based on the technical scheme, the target parameters of the target reactor are obtained through the scheme provided by the embodiment of the invention; constructing a voltage balance equation of the target reactor based on the target parameters, and solving the voltage balance equation to obtain a function relation of magnetic flux of the target reactor along with time variation; based on the material characteristics of the silicon steel sheet of the target reactor, obtaining a relation curve between the relative magnetic permeability and the magnetic density of the target reactor; solving a magnetic circuit equivalent equation of the target reactor to obtain a function relation of the magnetic flux leakage density of the target reactor changing along with time; and solving an objective function based on a relation curve between the relative magnetic permeability and the magnetic density of the objective reactor to obtain the magnetic densities of the iron cores at different moments, and then obtaining the closing transient current values at different moments corresponding to the magnetic densities of the iron cores at different moments, thereby realizing the accurate calculation of the closing transient current of the reactor.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings that are required to be used in the embodiments or the description of the prior art will be briefly described below, and it is obvious that the drawings in the following description are only embodiments of the present invention, and that other drawings can be obtained according to the provided drawings without inventive effort for a person skilled in the art.
Fig. 1 is a schematic flow chart of a method for calculating a closing transient current of a reactor according to an embodiment of the present application;
Fig. 2 is a schematic flow chart of a method for calculating a closing transient current of a reactor according to another embodiment of the present application;
fig. 3 is a schematic structural diagram of a reactor switching-on transient current calculating device according to an embodiment of the present application;
fig. 4 is a schematic structural diagram of an electronic device according to an embodiment of the present application.
Detailed Description
The following description of the embodiments of the present invention will be made clearly and completely with reference to the accompanying drawings, in which it is apparent that the embodiments described are only some embodiments of the present invention, but not all embodiments. All other embodiments, which can be made by those skilled in the art based on the embodiments of the invention without making any inventive effort, are intended to be within the scope of the invention.
The application discloses a calculation scheme of a reactor closing transient current, which is used for obtaining a functional relation of magnetic flux changing along with time by constructing a magnetic circuit transient equation and a voltage balance equation. And then according to the structural characteristics of the reactor, based on a preset silicon steel sheet B-u function relation, an objective function is established, the calculation of the closing transient current of the reactor at different moments is carried out by using solution methods such as loop iteration, linear interpolation, curve fitting and the like, and the attenuation process of the closing transient current is simulated.
Specifically, referring to fig. 1, the application discloses a method for calculating a reactor closing transient current, which comprises the following steps:
Step S101: and obtaining target parameters of the target reactor.
In this solution, the target parameters of the target reactor include, but are not limited to, a direct current resistance R of the target reactor, a main magnetic flux effective area of the target reactor, a leakage magnetic flux effective area S 0 of the target reactor, and the like.
Step S102: and constructing a voltage balance equation of the target reactor based on the target parameters, and solving the voltage balance equation to obtain a function relation of the magnetic flux of the target reactor along with the change of time.
In this step, a voltage balance equation is first constructedSolving the voltage balance equation to obtain the functional relation of the magnetic flux change along with time: wherein, psi is the closing angle of the target reactor, R is the direct current resistance of the target reactor, L is the inductance of the target reactor, phi m is the average magnetic flux peak value of the target reactor, t is the moment, phi (t) is the magnetic flux peak value of the target reactor at the moment t, and omega is the voltage angular frequency applied on the target reactor.
Step S103: and obtaining a relation curve between the relative magnetic permeability and the magnetic density of the target reactor based on the material characteristics of the silicon steel sheet of the target reactor.
In this step, it is necessary to obtain a relationship curve of the relative permeability and the magnetic flux density B-u according to the characteristics of the silicon steel sheet material of the target reactor. Specifically, the functional relationship between the relative magnetic permeability and the magnetic flux density may be represented as μ r＝f(BC, where μ r is the relative magnetic permeability of the target reactor, and B c is the magnetic flux density of the core of the target reactor, and the unit is T.
Step S104: and solving a magnetic circuit equivalent equation of the target reactor to obtain a function relation of the magnetic flux leakage density of the target reactor changing along with time.
The sum of magnetic sections consisting of the iron core and the air gap is equal to the magnetic circuit of the main air gap, and a functional relation of the magnetic leakage flux density changing along with time is obtained:
the sum delta is the total length of an iron core air gap of the target reactor, the unit is m, and the unit is m, wherein l is the reactor window height of the target reactor.
Step S105: solving an objective function based on a relation curve between relative magnetic permeability and magnetic density of a target reactor to obtain magnetic densities of iron cores at different moments, wherein the objective function is thatWherein, B c (t) is the magnetic flux density of the iron core at the time t, S c is the effective cross section area of the iron core of the target reactor, sigma delta is the total length of the air gap of the iron core of the target reactor, l is the window height of the target reactor, mu r (t) is the relative magnetic permeability of the target reactor at the time t, S 0 is the effective cross section area of the magnetic flux leakage of the target reactor at the time t, and phi (t) is the magnetic flux of the target reactor at the time t.
Step S106: and acquiring closing transient current values at different moments corresponding to the magnetic densities of the iron cores at different moments.
In this step, after the magnetic densities of the carved iron cores of the target reactor at different moments are obtained, the closing transient currents of the target reactor at different moments are calculated based on the magnetic densities of the carved iron cores at different moments, so that the attenuation process of the closing transient currents can be simulated, and specifically, the magnetic densities of the carved iron cores at different moments are substituted into a formulaBased on the formula, the closing transient currents of the target reactor at different moments can be calculated, B c (T) is the magnetic density of an iron core of the target reactor at different moments, the unit is T, sigma delta is the total length of an air gap of the iron core of the target reactor, m, mu r (T) is the relative magnetic permeability of the target reactor at different moments, l is the window height of the reactor of the target reactor, the unit is m, N is the number of turns of the reactor of the target reactor, and Sigma delta refers to the total length of the air gap of the target reactor.
In the technical solution disclosed in another embodiment of the present application, when solving the objective function based on the relation curve between the relative permeability and the magnetic density of the objective reactor to obtain the magnetic densities of the iron cores at different moments, the magnetic densities of the iron cores at different moments are obtained by calculation in a continuous iterative manner, specifically, referring to fig. 2, the process may include:
Step S201: and configuring the magnetic density of the primary iron core corresponding to each moment.
The magnetic density value of the primary iron core of the target reactor corresponding to each moment can be set according to the user demand.
Step S202: and substituting the magnetic densities of the primary selected iron cores corresponding to the moments into the objective function, and judging whether the value of F (t) (the calculation result of the objective function) is in a preset range or not.
Step S203: when the value of F (t) is not in the preset range, the primary iron core is more magnetic dense until the value of F (t) corresponding to the primary iron core magnetic dense is in the preset range;
step S204: when the value of F (t) is in a preset range, taking the magnetic density of the primary selected iron core as the magnetic density of the iron core at the corresponding moment;
step S205: and obtaining the magnetic density of the iron core at each moment.
In the scheme, firstly, a primary iron core magnetic density is configured at a moment t, then, according to a given silicon steel sheet B-u relation curve, the relative magnetic permeability matched with the primary iron core magnetic density is obtained through linear interpolation, whether the value of an objective function F (t) is in a preset range is checked, when the value of F (t) is in the preset range, the fact that F (t) =0 is met is indicated, if the value of F (t) is not in the preset range, the fact that F (t) =0 is not met is not indicated, the primary iron core magnetic density is revised again, the revised primary iron core magnetic density is substituted into the objective formula to be recalculated, and therefore, the objective formula is calculated in a continuous loop iteration mode until the value of F (t) is in the preset range, and the iron core magnetic density at the moment is taken as the iron core magnetic density corresponding to the moment. Then solving the magnetic density of the iron core at the next moment.
In this embodiment, corresponding to the above method, a device for calculating a transient current for closing a reactor is also disclosed, and specific working contents of each unit in the device are referred to in the content of the above method embodiment.
The following describes a reactor switching-on transient current calculation device provided by the embodiment of the present invention, and the reactor switching-on transient current calculation device described below and the reactor switching-on transient current calculation method described above may be referred to correspondingly.
Referring to fig. 3, a reactor closing transient current calculation apparatus disclosed in an embodiment of the present application may include,
A parameter acquisition unit a, a first calculation unit B, a second calculation unit C, a third calculation unit D, a fourth calculation unit E, and a fifth calculation unit F.
Corresponding to the method, the parameter obtaining unit A is used for obtaining the target parameter of the target reactor;
Corresponding to the method, the first calculating unit B is configured to construct a voltage balance equation of the target reactor based on the target parameter, and solve the voltage balance equation to obtain a functional relationship of magnetic flux of the target reactor over time;
corresponding to the method, the second calculating unit C is used for obtaining a relation curve between the relative magnetic permeability and the magnetic density of the target reactor based on the material characteristics of the silicon steel sheet of the target reactor;
corresponding to the method, the third calculation unit D is configured to solve a magnetic circuit equivalent equation of the target reactor, so as to obtain a functional relationship of the magnetic flux leakage density of the target reactor changing with time;
Corresponding to the method, the fourth calculating unit E is configured to solve an objective function based on a relationship curve between relative permeability and magnetic density of the objective reactor, so as to obtain magnetic densities of the iron cores at different moments, where the objective function is Wherein, B c (t) is the magnetic flux density of the iron core at time t, S c is the effective cross-sectional area of the iron core of the target reactor, Σδ is the total length of the air gap of the iron core of the target reactor, l is the window height of the target reactor, μ r (t) is the relative magnetic permeability of the target reactor at time t, S 0 is the effective cross-sectional area of the magnetic leakage of the target reactor, and Φ (t) is the magnetic flux of the target reactor at time t;
Corresponding to the method, the fifth calculating unit F is configured to obtain closing transient current values at different moments corresponding to the magnetic densities of the iron cores at different moments.
The specific working contents of the above units are described in the method embodiments, and are not described here.
Corresponding to the method, the application also discloses an electronic device, referring to fig. 4, the device comprises: at least one processor 100, at least one communication interface 200, at least one memory 300, and at least one communication bus 400;
in the embodiment of the present invention, the number of the processor 100, the communication interface 200, the memory 300 and the communication bus 400 is at least one, and the processor 100, the communication interface 200 and the memory 300 complete the communication with each other through the communication bus 400; it will be apparent that the communication connection schematic shown in the processor 100, the communication interface 200, the memory 300 and the communication bus 400 shown in fig. 4 is only optional;
alternatively, the communication interface 200 may be an interface of a communication module, such as an interface of a GSM module;
Processor 100 may be a central processing unit CPU, or an Application-specific integrated Circuit ASIC (Application SPECIFIC INTEGRATED Circuit), or one or more integrated circuits configured to implement embodiments of the present invention.
Memory 300 may comprise high-speed RAM memory or may further comprise non-volatile memory (non-volatile memory), such as at least one disk memory.
The processor 100 is specifically configured to:
acquiring target parameters of a target reactor;
constructing a voltage balance equation of the target reactor based on the target parameters, and solving the voltage balance equation to obtain a function relation of magnetic flux of the target reactor along with time variation;
based on the material characteristics of the silicon steel sheet of the target reactor, obtaining a relation curve between the relative magnetic permeability and the magnetic density of the target reactor;
Solving a magnetic circuit equivalent equation of the target reactor to obtain a function relation of the magnetic flux leakage density of the target reactor changing along with time;
solving an objective function based on a relation curve between relative magnetic permeability and magnetic density of a target reactor to obtain magnetic densities of iron cores at different moments, wherein the objective function is that Wherein, B c (t) is the magnetic flux density of the iron core at time t, S c is the effective cross-sectional area of the iron core of the target reactor, Σδ is the total length of the air gap of the iron core of the target reactor, l is the window height of the target reactor, μ r (t) is the relative magnetic permeability of the target reactor at time t, S 0 is the effective cross-sectional area of the magnetic leakage of the target reactor, and Φ (t) is the magnetic flux of the target reactor at time t;
and acquiring closing transient current values at different moments corresponding to the magnetic densities of the iron cores at different moments.
In this scheme, the electronic device is a computer or a PC.
For convenience of description, the above system is described as being functionally divided into various modules, respectively. Of course, the functions of each module may be implemented in the same piece or pieces of software and/or hardware when implementing the present invention.
In this specification, each embodiment is described in a progressive manner, and identical and similar parts of each embodiment are all referred to each other, and each embodiment mainly describes differences from other embodiments. In particular, for a system or system embodiment, since it is substantially similar to a method embodiment, the description is relatively simple, with reference to the description of the method embodiment being made in part. The systems and system embodiments described above are merely illustrative, wherein the elements illustrated as separate elements may or may not be physically separate, and the elements shown as elements may or may not be physical elements, may be located in one place, or may be distributed over a plurality of network elements. Some or all of the modules may be selected according to actual needs to achieve the purpose of the solution of this embodiment. Those of ordinary skill in the art will understand and implement the present invention without undue burden.
Those of skill would further appreciate that the various illustrative elements and algorithm steps described in connection with the embodiments disclosed herein may be implemented as electronic hardware, computer software, or combinations of both, and that the various illustrative elements and steps are described above generally in terms of functionality in order to clearly illustrate the interchangeability of hardware and software. Whether such functionality is implemented as hardware or software depends upon the particular application and design constraints imposed on the solution. Skilled artisans may implement the described functionality in varying ways for each particular application, but such implementation decisions should not be interpreted as causing a departure from the scope of the present invention.
The steps of a method or algorithm described in connection with the embodiments disclosed herein may be embodied directly in hardware, in a software module executed by a processor, or in a combination of the two. The software modules may be disposed in Random Access Memory (RAM), memory, read Only Memory (ROM), electrically programmable ROM, electrically erasable programmable ROM, registers, hard disk, a removable disk, a CD-ROM, or any other form of storage medium known in the art.
It is further noted that relational terms such as first and second, and the like are used solely to distinguish one entity or action from another entity or action without necessarily requiring or implying any actual such relationship or order between such entities or actions. Moreover, the terms "comprises," "comprising," or any other variation thereof, are intended to cover a non-exclusive inclusion, such that a process, method, article, or apparatus that comprises a list of elements does not include only those elements but may include other elements not expressly listed or inherent to such process, method, article, or apparatus. Without further limitation, an element defined by the phrase "comprising one … …" does not exclude the presence of other like elements in a process, method, article, or apparatus that comprises the element.
The previous description of the disclosed embodiments is provided to enable any person skilled in the art to make or use the present invention. Various modifications to these embodiments will be readily apparent to those skilled in the art, and the generic principles defined herein may be applied to other embodiments without departing from the spirit or scope of the invention. Thus, the present invention is not intended to be limited to the embodiments shown herein but is to be accorded the widest scope consistent with the principles and novel features disclosed herein.
Claims (8)
1. The method for calculating the closing transient current of the reactor is characterized by comprising the following steps of:
acquiring target parameters of a target reactor;
constructing a voltage balance equation of the target reactor based on the target parameters, and solving the voltage balance equation to obtain a function relation of magnetic flux of the target reactor along with time variation;
based on the material characteristics of the silicon steel sheet of the target reactor, obtaining a relation curve between the relative magnetic permeability and the magnetic density of the target reactor;
Solving a magnetic circuit equivalent equation of the target reactor to obtain a function relation of the magnetic flux leakage density of the target reactor changing along with time;
solving an objective function based on a relation curve between relative magnetic permeability and magnetic density of a target reactor to obtain magnetic densities of iron cores at different moments, wherein the objective function is that Wherein, B c (t) is the magnetic flux density of the iron core at time t, S c is the effective cross-sectional area of the iron core of the target reactor, Σδ is the total length of the air gap of the iron core of the target reactor, l is the window height of the target reactor, μ r (t) is the relative magnetic permeability of the target reactor at time t, S 0 is the effective cross-sectional area of the magnetic flux leakage of the target reactor at time t, and Φ (t) is the magnetic flux of the target reactor at time t;
acquiring closing transient current values at different moments corresponding to the magnetic densities of the iron cores at different moments;
Solving the objective function based on a relation curve between relative magnetic permeability and magnetic density of the objective reactor to obtain magnetic densities of the iron cores at different moments, comprising:
configuring magnetic density of a primary iron core corresponding to each moment;
Substituting the magnetic densities of the primary selected iron cores corresponding to the moments into the objective function, and judging whether the value of F (t) is in a preset range or not;
when the value of F (t) is not in the preset range, changing the magnetic density of the primary iron core until the value of F (t) corresponding to the magnetic density of the primary iron core is in the preset range;
when the value of F (t) is in a preset range, taking the magnetic density of the primary selected iron core as the magnetic density of the iron core at the corresponding moment;
Acquiring the magnetic density of the iron core at each moment;
The obtaining the closing transient current values at different moments corresponding to the magnetic densities of the iron cores at different moments comprises the following steps:
substituting magnetic densities of iron cores at different moments into formula And calculating to obtain closing transient current values I (t) at different moments corresponding to the magnetic densities of the iron cores at different moments, wherein Sigma sigma refers to the total length of an air gap of the target reactor, and N is the number of turns of the reactor of the target reactor.
2. The reactor closing transient current calculation method of claim 1, wherein the target parameters include:
the direct current resistance, the main magnetic flux effective area and the leakage magnetic flux effective area of the target reactor.
3. The reactor closing transient current calculation method of claim 2, wherein the magnetic flux as a function of time is:
wherein, psi is the closing angle of the target reactor, R is the direct current resistance of the target reactor, L is the inductance of the target reactor, phi m is the average magnetic flux peak value of the target reactor, t is the moment, phi (t) is the magnetic flux peak value of the target reactor at the moment t, and omega is the voltage angular frequency applied on the target reactor.
4. The reactor closing transient current calculation method of claim 3, wherein the relationship between the relative permeability and the magnetic density is:
Mu r＝f(BC), wherein mu r is the relative magnetic permeability of the target reactor, and B c is the magnetic density of the iron core of the target reactor.
5. The method for calculating a closing transient current of a reactor according to claim 4, wherein the leakage flux density is a function of time:
The reactor window height of the target reactor is l, the flux leakage density of the target reactor at the moment t is B 0 (t), the sigma delta is the total length of an iron core air gap of the target reactor, and the relative magnetic permeability of the target reactor at the moment t is mu r (t).
6. A reactor closing transient current calculation apparatus, comprising:
a parameter acquisition unit for acquiring a target parameter of a target reactor;
The first calculation unit is used for constructing a voltage balance equation of the target reactor based on the target parameters, and solving the voltage balance equation to obtain a function relation of magnetic flux of the target reactor along with time variation;
The second calculation unit is used for obtaining a relation curve between the relative magnetic permeability and the magnetic density of the target reactor based on the material characteristics of the silicon steel sheet of the target reactor;
The third calculation unit is used for solving a magnetic circuit equivalent equation of the target reactor to obtain a function relation of the magnetic flux leakage density of the target reactor changing along with time;
a fourth calculation unit, configured to solve an objective function based on a relationship curve between relative permeability and magnetic density of the objective reactor, to obtain magnetic densities of the iron cores at different times, where the objective function is Wherein, B c (t) is the magnetic flux density of the iron core at time t, S c is the effective cross-sectional area of the iron core of the target reactor, Σδ is the total length of the air gap of the iron core of the target reactor, l is the window height of the target reactor, μ r (t) is the relative magnetic permeability of the target reactor at time t, S 0 is the effective cross-sectional area of the magnetic leakage of the target reactor, and Φ (t) is the magnetic flux of the target reactor at time t;
A fifth calculation unit, configured to obtain closing transient current values at different moments corresponding to magnetic densities of the iron cores at different moments;
The fourth calculation unit solves the objective function based on a relation curve between relative permeability and magnetic density of the objective reactor to obtain magnetic densities of the iron cores at different moments, and comprises:
configuring magnetic density of a primary iron core corresponding to each moment;
Substituting the magnetic densities of the primary selected iron cores corresponding to the moments into the objective function, and judging whether the value of F (t) is in a preset range or not;
when the value of F (t) is not in the preset range, changing the magnetic density of the primary iron core until the value of F (t) corresponding to the magnetic density of the primary iron core is in the preset range;
when the value of F (t) is in a preset range, taking the magnetic density of the primary selected iron core as the magnetic density of the iron core at the corresponding moment;
Acquiring the magnetic density of the iron core at each moment;
the fifth calculating unit obtains closing transient current values at different moments corresponding to the magnetic densities of the iron cores at different moments, and the fifth calculating unit comprises:
substituting magnetic densities of iron cores at different moments into formula And calculating to obtain closing transient current values I (t) at different moments corresponding to the magnetic densities of the iron cores at different moments, wherein Sigma sigma refers to the total length of an air gap of the target reactor, and N is the number of turns of the reactor of the target reactor.
7. An electronic device, comprising:
a memory and a processor; the memory stores a program adapted for execution by the processor, the program for:
acquiring target parameters of a target reactor;
constructing a voltage balance equation of the target reactor based on the target parameters, and solving the voltage balance equation to obtain a function relation of magnetic flux of the target reactor along with time variation;
based on the material characteristics of the silicon steel sheet of the target reactor, obtaining a relation curve between the relative magnetic permeability and the magnetic density of the target reactor;
Solving a magnetic circuit equivalent equation of the target reactor to obtain a function relation of the magnetic flux leakage density of the target reactor changing along with time;
solving an objective function based on a relation curve between relative magnetic permeability and magnetic density of a target reactor to obtain magnetic densities of iron cores at different moments, wherein the objective function is that Wherein, B c (t) is the magnetic flux density of the iron core at time t, S c is the effective cross-sectional area of the iron core of the target reactor, Σδ is the total length of the air gap of the iron core of the target reactor, l is the window height of the target reactor, μ r (t) is the relative magnetic permeability of the target reactor at time t, S 0 is the effective cross-sectional area of the magnetic leakage of the target reactor, and Φ (t) is the magnetic flux of the target reactor at time t;
acquiring closing transient current values at different moments corresponding to the magnetic densities of the iron cores at different moments;
Solving the objective function based on a relation curve between relative magnetic permeability and magnetic density of the objective reactor to obtain magnetic densities of the iron cores at different moments, comprising:
configuring magnetic density of a primary iron core corresponding to each moment;
Substituting the magnetic densities of the primary selected iron cores corresponding to the moments into the objective function, and judging whether the value of F (t) is in a preset range or not;
when the value of F (t) is not in the preset range, changing the magnetic density of the primary iron core until the value of F (t) corresponding to the magnetic density of the primary iron core is in the preset range;
when the value of F (t) is in a preset range, taking the magnetic density of the primary selected iron core as the magnetic density of the iron core at the corresponding moment;
Acquiring the magnetic density of the iron core at each moment;
The obtaining the closing transient current values at different moments corresponding to the magnetic densities of the iron cores at different moments comprises the following steps:
substituting magnetic densities of iron cores at different moments into formula And calculating to obtain closing transient current values I (t) at different moments corresponding to the magnetic densities of the iron cores at different moments, wherein Sigma sigma refers to the total length of an air gap of the target reactor, and N is the number of turns of the reactor of the target reactor.
8. The electronic device of claim 7, wherein the electronic device is a computer or a PC.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210411344.6A CN114611061B (en) | 2022-04-19 | 2022-04-19 | Method and device for calculating closing transient current of reactor and electronic equipment |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210411344.6A CN114611061B (en) | 2022-04-19 | 2022-04-19 | Method and device for calculating closing transient current of reactor and electronic equipment |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN114611061A CN114611061A (en) | 2022-06-10 |
| CN114611061B true CN114611061B (en) | 2024-06-04 |
Family
ID=81869205
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202210411344.6A Active CN114611061B (en) | 2022-04-19 | 2022-04-19 | Method and device for calculating closing transient current of reactor and electronic equipment |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114611061B (en) |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2008122210A (en) * | 2006-11-10 | 2008-05-29 | Tokyo Metropolitan Univ | Device for measuring reactor loss and its measurement method |
| CN108962557A (en) * | 2018-06-04 | 2018-12-07 | 国网浙江省电力有限公司绍兴供电公司 | A kind of design method of the high pressure iron core filter reactor based on equivalent flux density method |
| CN112487677A (en) * | 2020-11-20 | 2021-03-12 | 西安西电变压器有限责任公司 | Method and device for calculating iron core loss of reactor |
| CN112710910A (en) * | 2020-11-26 | 2021-04-27 | 广西电网有限责任公司电力科学研究院 | Method for calculating excitation inrush current of core type shunt reactor and evaluating safety |
| CN113343538A (en) * | 2021-07-06 | 2021-09-03 | 华北电力大学（保定） | Reactor structure optimization method based on three-dimensional fixed-point harmonic balance finite element |
-
2022
- 2022-04-19 CN CN202210411344.6A patent/CN114611061B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2008122210A (en) * | 2006-11-10 | 2008-05-29 | Tokyo Metropolitan Univ | Device for measuring reactor loss and its measurement method |
| CN108962557A (en) * | 2018-06-04 | 2018-12-07 | 国网浙江省电力有限公司绍兴供电公司 | A kind of design method of the high pressure iron core filter reactor based on equivalent flux density method |
| CN112487677A (en) * | 2020-11-20 | 2021-03-12 | 西安西电变压器有限责任公司 | Method and device for calculating iron core loss of reactor |
| CN112710910A (en) * | 2020-11-26 | 2021-04-27 | 广西电网有限责任公司电力科学研究院 | Method for calculating excitation inrush current of core type shunt reactor and evaluating safety |
| CN113343538A (en) * | 2021-07-06 | 2021-09-03 | 华北电力大学（保定） | Reactor structure optimization method based on three-dimensional fixed-point harmonic balance finite element |
Also Published As
| Publication number | Publication date |
|---|---|
| CN114611061A (en) | 2022-06-10 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Wang et al. | Principle and design of a single-phase inverter-based grounding system for neutral-to-ground voltage compensation in distribution networks | |
| Piasecki et al. | Mitigating ferroresonance in voltage transformers in ungrounded MV networks | |
| CN109119971B (en) | Zero sequence overcurrent protection method, system and device and computer readable storage medium | |
| Rico et al. | The study of inrush current phenomenon using operational matrices | |
| CN105609300B (en) | A kind of Transformer shielding layer method for designing of inverse-excitation type switch power-supply | |
| Akbari et al. | Air gap inductance calculation for transformer transient model | |
| CN112035967B (en) | Simulation method and device for no-load exciting current of transformer under direct-current magnetic bias working condition | |
| Niyomsatian et al. | Time‐domain homogenization of multiturn windings based on RL Cauer ladder networks | |
| Singh et al. | Hyperbolic tangent function‐based least mean‐square control algorithm for distribution static compensator | |
| CN114611061B (en) | Method and device for calculating closing transient current of reactor and electronic equipment | |
| Tokic et al. | Numerical calculations of three-phase transformer transients | |
| Andreev et al. | Settings determination for numerical transformer differential protection via its detailed mathematical model | |
| Sahoo et al. | Power system harmonics estimation using adaptive filters | |
| Xu et al. | Harmonic suppression analysis of a harmonic filtering distribution transformer with integrated inductors based on field–circuit coupling simulation | |
| Preis et al. | Transient behavior of large transformer windings taking capacitances and eddy currents into account | |
| Dugan et al. | An example of 3-phase transformer modeling for distribution system analysis | |
| Vasilj et al. | Modeling of current-limiting air-core series reactor for transient recovery voltage studies | |
| Abbasi et al. | Voltage transformer ferroresonance analysis using multiple scales method and chaos theory | |
| Lesniewska et al. | Behaviour of measuring current transformers with cores composed from different magnetic materials at non‐rated loads and overcurrents | |
| Klempka et al. | Bank harmonic filters operation in power supply system–cases studies | |
| Riu et al. | Turbine generator modeling by non-integer order systems | |
| Zadeh et al. | A compensation scheme for CVT transient effects using artificial neural network | |
| Chen et al. | An improved numerical distance relay based on CCVT transient characteristic matching | |
| Li et al. | Core form transformer topological duality based transient model validation for exciting current calculation within DC bias | |
| CN104967107B (en) | A kind of method for inhibiting intermittent arc grounding overvoltage |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |