CN111104743A - Method for determining distribution of direct-current magnetic bias transient magnetic field and eddy current loss of transformer - Google Patents

Method for determining distribution of direct-current magnetic bias transient magnetic field and eddy current loss of transformer Download PDF

Info

Publication number
CN111104743A
CN111104743A CN201911312635.4A CN201911312635A CN111104743A CN 111104743 A CN111104743 A CN 111104743A CN 201911312635 A CN201911312635 A CN 201911312635A CN 111104743 A CN111104743 A CN 111104743A
Authority
CN
China
Prior art keywords
magnetic
transformer
eddy current
transient
magnetic field
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.)
Pending
Application number
CN201911312635.4A
Other languages
Chinese (zh)
Inventor
董霞
刘树林
俞晓冬
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Qilu University of Technology
Original Assignee
Qilu University of Technology
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 Qilu University of Technology filed Critical Qilu University of Technology
Priority to CN201911312635.4A priority Critical patent/CN111104743A/en
Publication of CN111104743A publication Critical patent/CN111104743A/en
Pending legal-status Critical Current

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00—Arrangements or instruments for measuring magnetic variables
    • G01R33/02—Measuring direction or magnitude of magnetic fields or magnetic flux

Landscapes

  • Physics & Mathematics (AREA)
  • Condensed Matter Physics & Semiconductors (AREA)
  • General Physics & Mathematics (AREA)
  • Measuring Magnetic Variables (AREA)

Abstract

The invention discloses a method for determining distribution of a transient magnetic field and eddy current loss of a direct current magnetic bias of a transformer, which adopts the technical scheme that: establishing a transient circuit model and a magnetic circuit model of the transformer under the direct current magnetic biasing, wherein in the circuit model, a direct current voltage source is added on a primary side winding of the transformer, and a circuit equation containing direct current voltage is listed; in the magnetic circuit model, different topological structures of a transformer core and the magnetic saturation characteristics of ferromagnetic materials are considered, magnetomotive force generated by eddy current is introduced into the model, and a magnetic circuit equation is listed; and coupling the magnetic circuit equation with the circuit equation to obtain a nonlinear differential-algebraic equation set, solving the equation set by using a Newton-Raphson method, and determining the transient eddy current loss and the magnetic field distribution of the transformer with different topological structures under the direct-current magnetic bias. The calculation process of the invention has fast convergence speed, and the error of the calculation result is smaller and more accurate, and the problems of local overheating, vibration, noise and the like of the transformer caused by direct current magnetic biasing can be further accurately analyzed by using the calculation result of the invention.

Description

Method for determining distribution of direct-current magnetic bias transient magnetic field and eddy current loss of transformer
Technical Field
The invention relates to the technical field of transformers in the professional directions of motors and electric appliances, in particular to a method for determining distribution of direct-current magnetic bias transient magnetic field and eddy current loss of a transformer.
Background
Dc magnetic biasing is an abnormal phenomenon in which dc current is generated in a winding of a transformer to generate dc magnetic flux in an iron core. The reason for the dc magnetic bias phenomenon is two: firstly, the solar magnetic storm generates an induced potential difference on the earth surface, so that a Geomagnetic Induction Current (GIC) is entered into a transformer winding with a grounded neutral point. The frequency of the geomagnetic induction current is very low, about 0.001-1 Hz, and can be approximated to a DC current. And secondly, the influence of direct current transmission engineering. When the high-voltage direct-current transmission project adopts a single-pole-earth return line operation mode or a double-pole unbalanced operation mode, larger direct current flows into the earth through the grounding electrode, so that direct current potential difference is generated on the surface of the earth, and direct current components enter a transformer winding with a grounded neutral point through the earth.
The dc bias causes the core of the transformer to quickly reach a saturated state, and the exciting current is distorted, thereby generating a series of electromagnetic effects, as shown in fig. 1. In addition, the dc magnetic bias may cause malfunction of a relay protection device of the power system, harmonic amplification of a capacitor, and unstable system voltage, which may affect normal operation of the power system and power equipment, and may even cause accidents such as capacitor explosion, system voltage breakdown, and blackout. Therefore, the dc magnetic bias phenomenon is not negligible, and is a hot point of research particularly for transient magnetic field and loss distribution of the transformer core.
A hysteresis model method is mostly adopted for calculating transient magnetic fields and losses of transformer cores at home and abroad, and the method utilizes hysteresis characteristic models such as a Preisach model or a Jiles-Atherton model and the like to estimate hysteresis loops of ferromagnetic materials during alternating current dynamic magnetization, expresses the relation between magnetic flux density B or magnetization M and magnetic field intensity H through a mathematical formula, is coupled with finite element numerical calculation to analyze the magnetic field distribution of the cores, and finally calculates the core losses through numerical integration of B and H by utilizing the Pongting theorem. The method relates to complex finite element numerical derivation, the programming calculation workload is large, the calculation result has larger error, and the engineering practicability is poor.
Disclosure of Invention
In order to overcome the defects of the existing method, the invention provides a method for determining the distribution of the transient magnetic field and the eddy current loss of the direct current magnetic bias of the transformer. And establishing a circuit-magnetic circuit coupling model of the transformer under the DC magnetic biasing and a relational expression of magnetic flux and iron core eddy current loss, and performing simulation calculation on the transient magnetic field and the eddy current loss distribution of the iron core by using a Newton-Raphson method. The method converts the solving process of the nonlinear equation into iterative solving of the corresponding linear equation, has the advantages of simple and practical calculating process, high convergence rate, more accurate calculating result and smaller error, and can be widely applied to the engineering field.
The invention adopts the following technical scheme:
a method for determining distribution of transient magnetic field and eddy current loss of direct current magnetic biasing of a transformer comprises the following steps:
establishing a transient circuit model and a magnetic circuit model of the transformer under the direct current magnetic biasing, wherein in the circuit model, a direct current voltage source is added on a primary side winding of the transformer, and a circuit equation containing direct current voltage is listed; in the magnetic circuit model, different topological structures of a transformer core and the magnetic saturation characteristics of ferromagnetic materials are considered, magnetomotive force generated by eddy current is introduced into the model, and a magnetic circuit equation is listed;
and coupling the magnetic circuit equation with the circuit equation to obtain a nonlinear differential-algebraic equation set, solving the equation set by using a Newton-Raphson method, and determining the transient magnetic field and eddy current loss distribution of the transformer with different topological structures under the direct-current magnetic bias.
Further, a no-load equivalent circuit model of the transformer under the direct-current magnetic biasing is established, and a corresponding circuit equation is obtained:
Figure BDA0002324952810000021
wherein u isjRepresenting an alternating voltage, U0Representing a direct voltage, rjpDenotes the winding resistance, LjpRepresenting the equivalent leakage inductance of the winding, ijpRepresents the current, ejpRepresenting induced electromotive force,. phijDenotes a magnetic flux, and p denotes a primary-side three-phase winding.
Further, in the magnetic circuit model, a magnetic circuit equation obtained based on the three-phase five-column transformer is as follows:
Figure BDA0002324952810000031
wherein j is 1 to 7, Fa、FbAnd FcMagnetomotive force, R, generated for primary side three-phase winding currentjIs the magnetic resistance of the iron core column, the magnet yoke and the side column, phijIs the magnetic flux of the core limb, the magnet yoke and the side limb, R8For magnetic leakage resistance,. phi8For leakage flux, FejMagnetomotive force generated by eddy currents;
for a three-phase three-column or single-phase three-column transformer, j is 1-5, R6For magnetic leakage resistance,. phi6Is the leakage flux.
Further, the magnetic saturation characteristic of the ferromagnetic material of the transformer core under the direct current magnetic bias is expressed by a B-H nonlinear equation:
B＝f(H)
where B is the magnetic flux density and H is the magnetic field strength.
Further, the nonlinear reluctance in the magnetic circuit is:
Figure BDA0002324952810000032
where A is the core cross-sectional area and l is the magnetic path length.
Further, R isjSubstituting into a magnetic circuit equation, then combining with the circuit equation to obtain a nonlinear differential-algebraic equation system, and performing iterative solution on the nonlinear differential-algebraic equation system by using a Newton-Raphson method to obtain instantaneous values of B, H phi.
Further, determining transient eddy current loss P of transformer coreeAnd Φ as:
pe＝f(Φ)；
solving to obtain the transient eddy current loss P of the iron coreeDistribution of (2).
Compared with the prior art, the invention has the beneficial effects that:
(1) the invention fully considers the different topological structures of the transformer iron core and the magnetic saturation characteristics of the ferromagnetic material, introduces the magnetomotive force generated by the eddy current into the magnetic circuit model, and then couples the magnetomotive force with the circuit model; the distribution of the transient magnetic field of the transformer under the direct current magnetic bias can be accurately and quickly obtained through the circuit-magnetic circuit model; and then the relation expression of the eddy current loss and the magnetic flux is utilized to quickly obtain the transient eddy current loss distribution of the iron core.
(2) Compared with the existing calculation method, the method is simple and clear, the calculation process is fast in convergence, the calculation result is more accurate, the error is smaller, and the method is suitable for transformers with different topological structures.
Drawings
The accompanying drawings, which form a part of the specification, are included to provide a further understanding of the application, and are incorporated in and constitute a part of this specification, illustrate embodiment(s) of the application and together with the description serve to explain the application and not to limit the application.
FIG. 1 is a schematic diagram of the effect of DC magnetic bias on a transformer;
FIG. 2 is a flowchart of a first embodiment of the present invention;
fig. 3 is a no-load circuit model of a three-phase transformer under dc magnetic biasing according to a first embodiment of the present invention;
fig. 4 is a magnetic circuit model of a three-phase five-limb transformer under dc magnetic bias according to a first embodiment of the present invention;
5(a) -5(b) are transient magnetic field distributions in the iron core of the three-phase five-limb transformer of the second embodiment of the invention when no DC current is provided;
6(a) -6(b) are transient magnetic field distributions in the iron core when the three-phase five-column transformer of the second embodiment of the present invention has DC;
fig. 7(a) -7(b) are graphs showing the variation of the magnetic field amplitude of the three-phase five-limb transformer according to the second embodiment of the present invention with dc current;
fig. 8 shows the distribution of the transient eddy current loss in the core of the three-phase five-limb transformer according to the second embodiment of the present invention without dc;
fig. 9 shows the transient eddy current loss distribution in the core when the three-phase five-limb transformer of the second embodiment of the invention has dc;
fig. 10 shows the eddy current loss amplitude of the three-phase five-limb transformer according to the second embodiment of the present invention as a function of dc current.
Detailed Description
It should be noted that the following detailed description is exemplary and is intended to provide further explanation of the disclosure. Unless defined otherwise, all technical and scientific terms used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this application belongs.
It is noted that the terminology used herein is for the purpose of describing particular embodiments only and is not intended to be limiting of example embodiments according to the present application. As used herein, the singular forms "a", "an", and/or "the" are intended to include the plural forms as well, and it should be understood that when the terms "comprises" and/or "comprising" are used in this specification, they specify the presence of stated features, steps, operations, devices, components, and/or combinations thereof;
the first embodiment is as follows:
the present invention will be described in detail with reference to fig. 3 to 10, and specifically, the structure is as follows:
the embodiment provides a method for determining distribution of a transient magnetic field and eddy current loss of a direct current magnetic bias of a transformer, which comprises the following steps: and establishing a transient circuit-magnetic circuit model of the transformer under the DC magnetic biasing. In the circuit model, a direct-current voltage source is added on a primary side winding of a transformer, and a circuit equation containing direct-current voltage is listed; in the magnetic circuit model, different topological structures of a transformer core and the magnetic saturation characteristics of ferromagnetic materials are considered, magnetomotive force generated by eddy current is introduced into the model, and a magnetic circuit equation is listed; and then coupling the magnetic circuit equation with the circuit equation to obtain a nonlinear differential-algebraic equation set, solving the equation set by using a Newton-Raphson method, and determining the transient magnetic field and eddy current loss distribution of the transformer with different topological structures under the direct current magnetic biasing. Specifically, the method comprises the following steps:
firstly, establishing a circuit model:
the no-load equivalent circuit model of the three-phase transformer under DC magnetic bias is shown in FIG. 3, U0Representing a direct voltage u1、u2And u3Is a three-phase alternating voltage, e1a、e2b、e3cInduced electromotive force of the primary side windings, respectively, e1A、e2B、e2CInduced electromotive force, i, of the secondary winding, respectively1a、i2bAnd i3Are respectively primary side winding current, r1a、r2b、r3cAre primary side winding resistances, r1A、r2B、r3CRespectively, secondary side winding resistance, L1a、L2bAnd L3cRespectively, equivalent leakage inductance, L, of only the interlinked primary windings1A、L2BAnd L3CThe equivalent leakage inductance of only the cross-linked secondary side winding is respectively, and a loop formed by three phases in the primary side winding and a neutral wire meets the following circuit equation:
Figure BDA0002324952810000061
wherein p represents a primary side three-phase winding (a, b, c), ΦjRepresents the three legs flux of the transformer, j is 1,2, 3.
If the transformer is a single-phase transformer, j is 1.
Secondly, establishing a magnetic circuit model:
the equivalent magnetic circuit model of the three-phase five-column transformer is shown in FIG. 4, Fa、FbAnd FcMagnetomotive force generated by primary side three-phase winding current (only F if single-phase three-pole transformer)a)，RjAnd phij(j 1-7) is the magnetic resistance and flux of the core limb, the magnetic yoke and the side limb, R8And phi8For leakage magnetic resistance and flux, Fej(j 1-7) is magnetomotive force generated by eddy current (if the transformer is a three-phase three-pole or single-phase three-pole transformer, j 1-5, R6And phi6Leakage magnetic resistance and leakage magnetic flux) and the magnetic path equation is:
Figure BDA0002324952810000071
the magnetic circuit model takes the influence of eddy current into consideration. On a section of magnetic circuit with the cross section area of A and the length of l, a winding with the number of turns of N is connected with current i, the magnetic flux in the magnetic circuit is phi, and the magnetomotive force generated by the eddy current can be expressed as:
Figure BDA0002324952810000072
in the formula, the coefficient keComprises the following steps:
Figure BDA0002324952810000073
where σ denotes an electric conductivity, d denotes a silicon steel sheet thickness, τ denotes a silicon steel sheet width, G-0.1356, H0Representing the internal magnetic potential generated by the magnetic domain.
Defining eddy currents i in branches of the magnetic circuiteThe expression is as follows:
Figure BDA0002324952810000074
the eddy current loss P in this magnetic circuit brancheComprises the following steps:
Figure BDA0002324952810000075
the total core eddy current loss P is determined for single-phase three-column and three-phase three-column transformerseCan be expressed as:
Figure BDA0002324952810000076
for three-phase five-column transformer, the total core eddy current loss PeCan be expressed as:
Figure BDA0002324952810000081
in addition, considering the magnetic saturation characteristics of the ferromagnetic material of the transformer under direct-current magnetic biasing, a B-H nonlinear equation adopts a single-value curve:
B＝f(H) (9)
the nonlinear reluctance in the magnetic circuit is then:
Figure BDA0002324952810000082
respectively substituting formula (7) or formula (8) and formula (10) into magnetic circuit equation (2) of transformers with different topological structures, then combining formula (1) and formula (2), and transforming and sorting the magnetic circuit equation to obtain a matrix form:
RΦj＝F(j＝1,2,3) (11)
and (4) iteratively solving the equation set by using a Newton-Raphson method. By selecting proper step length, the stability and the calculation precision of the solution can be ensured.
Example two:
in the present embodiment, a three-phase five-pole transformer is taken as an example, and the transient magnetic field and the eddy current loss under the presence and absence of dc bias are determined respectively, as shown in fig. 5 to 10.
In FIGS. 5 to 10, H represents the magnetic field intensity, B represents the magnetic induction intensity, and P representseRepresents the eddy current loss; subscripts 1,2,3 denote three core legs of a phase, b phase, c phase, 4, 5 denote yokes, 6, 7 denote side legs.
Fig. 5 and 6 show the transient magnetic field distribution without dc and with dc, respectively. As can be seen from fig. 5, when there is no dc intrusion, the magnetic field intensity in the core is a peak wave (same as the excitation current waveform), the magnetic flux density is a sine wave, and the positive and negative half axes are symmetrical. The magnetic field intensity and the magnetic flux density amplitude are the largest in the three iron core columns, the magnetic yoke is arranged next to the magnetic field intensity and the magnetic flux density amplitude is the smallest in the side columns. As can be seen from fig. 6, when there is dc intrusion, the magnetic field intensity and the magnetic flux density waveform are shifted, the positive and negative half-axes are no longer symmetrical, the amplitude of the positive half-axis is greater than that of the negative half-axis, and the positive amplitude of the waveform of the b-phase pillar is the largest.
Fig. 7 shows the variation of the magnetic field amplitude in the iron core of the three-phase five-limb transformer with the dc current. It can be seen that, as the dc current increases, the magnetic field strength and the magnetic flux density in the iron core increase in amplitude in a non-linear manner. The amplitude of the magnetic field intensity and the magnetic flux density of the phase b core column increases at the fastest speed, then the phase a and the phase c core columns, then the two side columns, and the yoke increases at the slowest speed. The magnetic field of the b-phase core column is influenced the most by the direct current, and the magnetic field of the magnetic yoke is influenced the least by the direct current.
Fig. 8 and 9 are distributions of transient eddy current losses in a three-phase five-limb transformer core without dc and with dc, respectively. As can be seen from FIG. 8, the eddy current loss (P) of the three legs is found in the absence of DCe1、Pe2、Pe3) Eddy current loss (P) from yokee4、Pe5) Equivalent, eddy current losses (P) of the two side legse6、Pe7) And, at a minimum, the amplitude in each half cycle is respectively equal. As shown in FIG. 9, when there is a DC intrusion, the slave Pe1～Pe7The amplitude of each half period is not equal any more, namely the amplitude of one half period is reduced, and the amplitude of the other half period is increased, which indicates that the DC current causes the offset of the eddy current loss.
Fig. 10 shows the eddy current loss amplitude in the iron core of the three-phase five-limb transformer as a function of the direct current. As can be seen, the eddy current loss amplitude is greatest in the yoke, followed by three legs, and the side legs are smallest. Because the length of the yoke is less than the length of the stem, the magnitude of the eddy current loss in the yoke is greater than the magnitude of the eddy current loss in the stem, although the magnitude of the magnetic field in the yoke is smaller. The eddy current loss in the iron core increases in amplitude in a nonlinear manner with the increase of the direct current. Wherein the eddy current loss amplitude in the side leg increases at the fastest rate, followed by the yoke, and the eddy current loss amplitude in the three legs increases the slowest. The direct current has the largest influence on the eddy current loss in the side column and the smallest influence on the eddy current loss in the center column.
The above description is only a preferred embodiment of the present application and is not intended to limit the present application, and various modifications and changes may be made by those skilled in the art. Any modification, equivalent replacement, improvement and the like made within the spirit and principle of the present application shall be included in the protection scope of the present application.

Claims (7)

1.一种确定变压器直流偏磁瞬态磁场和涡流损耗分布的方法，其特征在于，包括以下步骤：1. a method for determining transformer DC bias transient magnetic field and eddy current loss distribution, is characterized in that, comprises the following steps: 建立直流偏磁下变压器的瞬态电路模型和磁路模型，电路模型中，在变压器的一次侧绕组上加直流电压源，并列出含有直流电压的电路方程；磁路模型中，考虑变压器铁心不同的拓扑结构及铁磁材料的磁饱和特性，并将涡流产生的磁动势引入模型，列出磁路方程；The transient circuit model and magnetic circuit model of the transformer under DC bias are established. In the circuit model, a DC voltage source is added to the primary winding of the transformer, and the circuit equation containing the DC voltage is listed; in the magnetic circuit model, the transformer core is considered Different topological structures and magnetic saturation characteristics of ferromagnetic materials, and introduce the magnetomotive force generated by eddy current into the model, and list the magnetic circuit equation; 将磁路方程与电路方程进行耦合，得到非线性微分—代数方程组，利用牛顿—拉夫逊方法求解该方程组，确定直流偏磁下不同拓扑结构变压器的瞬态磁场及涡流损耗分布。The magnetic circuit equation and the circuit equation are coupled to obtain a nonlinear differential-algebraic equation system, which is solved by the Newton-Raphson method to determine the transient magnetic field and eddy current loss distribution of transformers with different topological structures under DC bias. 2.根据权利要求1所述的一种确定变压器直流偏磁瞬态磁场和涡流损耗分布的方法，其特征在于，建立直流偏磁下变压器的空载等效电路模型，并得出相应的电路方程：2. a kind of method of determining transformer DC bias transient magnetic field and eddy current loss distribution according to claim 1, is characterized in that, establishes the no-load equivalent circuit model of transformer under DC bias, and draws the corresponding circuit equation:
Figure FDA0002324952800000011
Figure FDA0002324952800000011
 其中，uj表示交流电压，U0表示直流电压，rjp表示绕组电阻，Ljp表示绕组的等效漏磁感，ijp表示电流，ejp表示感应电动势，Φj表示磁通，p表示一次侧三相绕组。Among them, u j represents the AC voltage, U 0 represents the DC voltage, r jp represents the winding resistance, L jp represents the equivalent leakage inductance of the winding, i jp represents the current, e jp represents the induced electromotive force, Φ j represents the magnetic flux, and p represents the Primary side three-phase winding. 3.根据权利要求1所述的一种确定变压器直流偏磁瞬态磁场和涡流损耗分布的方法，其特征在于，磁路模型中，基于三相五柱变压器得到的磁路方程为：3. a kind of method of determining transformer DC bias transient magnetic field and eddy current loss distribution according to claim 1, is characterized in that, in magnetic circuit model, the magnetic circuit equation that obtains based on three-phase five-column transformer is:
Figure FDA0002324952800000021
Figure FDA0002324952800000021
 其中，j＝1～7，Fa、Fb和Fc为一次侧三相绕组电流产生的磁动势，Rj为铁心柱、磁轭及旁柱磁阻，Φj为铁心柱、磁轭及旁柱磁通，R8为漏磁阻，Φ8为漏磁通，Fej为由涡流产生的磁动势；Among them, j=1～7, F a , F b and F c are the magnetomotive force generated by the three-phase winding current on the primary side, R j is the magnetic resistance of the core column, the yoke and the side column, Φ j is the core column, magnetic The magnetic flux of the yoke and the side column, R 8 is the leakage magnetic resistance, Φ 8 is the leakage magnetic flux, and F ej is the magnetomotive force generated by the eddy current; 对于三相三柱或单相三柱变压器，则j＝1～5，R6为漏磁阻，Φ6为漏磁通。For a three-phase three-column or single-phase three-column transformer, j = 1 to 5, R 6 is the leakage magnetic resistance, and Φ 6 is the leakage magnetic flux. 4.根据权利要求3所述的一种确定变压器直流偏磁瞬态磁场和涡流损耗分布的方法，其特征在于，用B-H非线性方程表示直流偏磁下变压器铁心铁磁材料的磁饱和特性：4. a kind of method for determining transformer DC bias transient magnetic field and eddy current loss distribution according to claim 3, it is characterized in that, represent the magnetic saturation characteristic of transformer core ferromagnetic material under DC bias with B-H nonlinear equation: B＝f(H)B=f(H) 其中，B是磁通密度，H是磁场强度。where B is the magnetic flux density and H is the magnetic field strength. 5.根据权利要求4所述的一种确定变压器直流偏磁瞬态磁场和涡流损耗分布的方法，其特征在于，磁路中的非线性磁阻为：5. a kind of method of determining transformer DC bias transient magnetic field and eddy current loss distribution according to claim 4 is characterized in that, the nonlinear reluctance in the magnetic circuit is:
Figure FDA0002324952800000022
Figure FDA0002324952800000022
 其中，A是铁心截面积，l是磁路长度。Among them, A is the cross-sectional area of the iron core, and l is the length of the magnetic path. 6.根据权利要求5所述的一种确定变压器直流偏磁瞬态磁场和涡流损耗分布的方法，其特征在于，将Rj代入磁路方程，然后与电路方程联立，得到非线性微分-代数方程组，利用牛顿—拉夫逊方法对其进行迭代求解，得到B、H、Φ的瞬时值。6. a kind of method for determining transformer DC bias transient magnetic field and eddy current loss distribution according to claim 5, it is characterized in that, R j is substituted into magnetic circuit equation, then simultaneously with circuit equation, obtains nonlinear differential - The algebraic equations are solved iteratively using the Newton-Raphson method, and the instantaneous values of B, H, and Φ are obtained. 7.根据权利要求6所述的一种确定变压器直流偏磁瞬态磁场和涡流损耗分布的方法，其特征在于，确定变压器铁心瞬态涡流损耗Pe和Φ之间的关系表达式为：7. a kind of method for determining transformer DC bias transient magnetic field and eddy current loss distribution according to claim 6, is characterized in that, it is determined that the relational expression between transformer core transient eddy current loss P e and Φ is: pe＝f(Φ)；p e = f(Φ); 求解得到铁心瞬态涡流损耗Pe的分布。The distribution of the core transient eddy current loss P e is obtained by solving.
CN201911312635.4A 2019-12-18 2019-12-18 Method for determining distribution of direct-current magnetic bias transient magnetic field and eddy current loss of transformer Pending CN111104743A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911312635.4A CN111104743A (en) 2019-12-18 2019-12-18 Method for determining distribution of direct-current magnetic bias transient magnetic field and eddy current loss of transformer

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911312635.4A CN111104743A (en) 2019-12-18 2019-12-18 Method for determining distribution of direct-current magnetic bias transient magnetic field and eddy current loss of transformer

Publications (1)

Publication Number Publication Date
CN111104743A true CN111104743A (en) 2020-05-05

Family

ID=70422255

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911312635.4A Pending CN111104743A (en) 2019-12-18 2019-12-18 Method for determining distribution of direct-current magnetic bias transient magnetic field and eddy current loss of transformer

Country Status (1)

Country Link
CN (1) CN111104743A (en)

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111751765A (en) * 2020-07-06 2020-10-09 华北电力大学 A method for determining AC resistivity of medium-voltage windings of high-frequency step-down transformers
CN111913135A (en) * 2020-08-19 2020-11-10 中国海洋石油集团有限公司 Method for identifying transformer winding deformation fault by Newton iteration method
CN111931310A (en) * 2020-08-28 2020-11-13 西南交通大学 Method for evaluating eddy current loss of wound core interlayer short circuit in consideration of boundary conditions of different magnetic fields
CN112560301A (en) * 2020-11-26 2021-03-26 东南大学 Magnetic material eddy current loss calculation method based on magnetic induction principle
CN112632458A (en) * 2020-12-23 2021-04-09 西南交通大学 Traction transformer eddy current loss solving method considering iron core magnetic saturation
CN112966371A (en) * 2021-02-08 2021-06-15 华北电力大学（保定） Abnormal loss calculation method of ferromagnetic material under alternating current-direct current hybrid excitation
CN113255172A (en) * 2021-07-12 2021-08-13 国网江西省电力有限公司电力科学研究院 Winding real-time temperature rise calculation method under repeated short-time short-circuit working condition
CN114966480A (en) * 2022-06-07 2022-08-30 国网江西省电力有限公司电力科学研究院 Remote direct-current magnetic bias monitoring platform for transformer

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2010072773A (en) * 2008-09-17 2010-04-02 Meidensha Corp Method for analyzing in-magnet eddy current loss for permanent magnet type motor
CN104239652A (en) * 2014-10-13 2014-12-24 国家电网公司 Three-phase five-column transformer modeling analysis method based on EIC principle

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2010072773A (en) * 2008-09-17 2010-04-02 Meidensha Corp Method for analyzing in-magnet eddy current loss for permanent magnet type motor
CN104239652A (en) * 2014-10-13 2014-12-24 国家电网公司 Three-phase five-column transformer modeling analysis method based on EIC principle

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
DONG XIA: "The Distribution of Transient Magnetic Field and Eddy Current Losses of Three-Phase Five-Legged Transformer Under DC Bias", 《2018 IEEE INTERNATIONAL CONFERENCE ON INTERNET OF THINGS (ITHINGS) AND IEEE GREEN COMPUTING AND COMMUNICATIONS (GREENCOM) AND IEEE CYBER, PHYSICAL AND SOCIAL COMPUTING (CPSCOM) AND IEEE SMART DATA (SMARTDATA)》, pages 728 - 732 *
董霞;刘志珍;: "三相变压器直流偏磁仿真分析", 电力自动化设备, no. 07 *

Cited By (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111751765B (en) * 2020-07-06 2021-04-30 华北电力大学 A method for determining AC resistivity of medium-voltage windings of high-frequency step-down transformers
CN111751765A (en) * 2020-07-06 2020-10-09 华北电力大学 A method for determining AC resistivity of medium-voltage windings of high-frequency step-down transformers
CN111913135A (en) * 2020-08-19 2020-11-10 中国海洋石油集团有限公司 Method for identifying transformer winding deformation fault by Newton iteration method
CN111931310A (en) * 2020-08-28 2020-11-13 西南交通大学 Method for evaluating eddy current loss of wound core interlayer short circuit in consideration of boundary conditions of different magnetic fields
WO2022110529A1 (en) * 2020-11-26 2022-06-02 东南大学 Magnetic material eddy-current loss calculation method based on magnetic induction principle
CN112560301A (en) * 2020-11-26 2021-03-26 东南大学 Magnetic material eddy current loss calculation method based on magnetic induction principle
US11790132B2 (en) 2020-11-26 2023-10-17 Southeast University Calculation method of eddy current loss in magnetic materials based on magnetic-inductance
CN112560301B (en) * 2020-11-26 2021-10-08 东南大学 An eddy current loss calculation method for magnetic materials based on the principle of magnetic induction
CN112632458A (en) * 2020-12-23 2021-04-09 西南交通大学 Traction transformer eddy current loss solving method considering iron core magnetic saturation
CN112966371A (en) * 2021-02-08 2021-06-15 华北电力大学（保定） Abnormal loss calculation method of ferromagnetic material under alternating current-direct current hybrid excitation
CN113255172B (en) * 2021-07-12 2021-11-19 国网江西省电力有限公司电力科学研究院 A real-time temperature rise calculation method for windings under repeated short-term short-circuit conditions
CN113255172A (en) * 2021-07-12 2021-08-13 国网江西省电力有限公司电力科学研究院 Winding real-time temperature rise calculation method under repeated short-time short-circuit working condition
CN114966480A (en) * 2022-06-07 2022-08-30 国网江西省电力有限公司电力科学研究院 Remote direct-current magnetic bias monitoring platform for transformer

Similar Documents

Publication Publication Date Title
CN111104743A (en) Method for determining distribution of direct-current magnetic bias transient magnetic field and eddy current loss of transformer
CN100589108C (en) A Digital Simulation Modeling Method of Magnetically Controlled Shunt Reactor
Moses et al. Dynamic modeling of three-phase asymmetric power transformers with magnetic hysteresis: No-load and inrush conditions
Chen et al. Study of analysis and experiment for ability to withstand DC bias in power transformers
Wang et al. Estimation model of core loss under DC bias
CN104849587A (en) Method for analyzing excitation characteristic change of transformer under influence of direct-current magnetic bias
Matsumori et al. Iron Loss Calculation of AC Filter Inductor for Three‐Phase PWM Inverters
Albert et al. Measurement based transformer modelling approach
Shindo et al. High frequency nonlinear modeling of magnetic sheets using polynomial expansions for eddy-current field
Fukumoto et al. Calculating leakage reactance of 9-winding transformer using time-dependent 3D FEM analysis
Brudny et al. Characterization and modeling of a virtual air gap by means of a reluctance network
Del Vecchio et al. Analytic nonlinear correction to the impedance boundary condition
Özüpak Performing Structural Design and Modeling of Transformers Using ANSYSMaxwell
Albert et al. Comparing two topology transformer hysteresis models with power transformer measurements
Sato et al. Study on an accurate iron loss calculation method considering the non-uniformity of the magnetic flux density
Mousavi et al. Comprehensive study on magnetization current harmonics of power transformers due to GICs
CN106532723A (en) Reactive power adjustment scheme of series transformer in unified power flow controller under DC magnetic bias
Zhao et al. The DC bias risk assessment of transformer based on nonlinear magnetic circuit equations
Lin et al. The effects of steel lamination core losses on transient magnetic fields using T-Ω method
CN112467726A (en) Transient characteristic analysis method and system of saturated superconducting current limiter
Ye et al. Analysis of transformer short circuit characteristics based on 3-D finite element method
Sakashita et al. Postcorrection of current/voltage and electromagnetic force for efficient hysteretic magnetic field analysis
Xia The Research of Transient Magnetic Field and Eddy Current Loss of Single Phase Three-Legged Transformer under AC/DC Condition
Susnjic et al. Electromagnetic analysis applied to the prediction of stray losses in power transformer
Zhao et al. Computation and analysis of DC-biased eddy current problems by an efficient fixed-point technique in the harmonic domain

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
CB02 Change of applicant information

Country or region after: China

Address after: 250353 University Road, Changqing District, Ji'nan, Shandong Province, No. 3501

Applicant after: Qilu University of Technology (Shandong Academy of Sciences)

Address before: 250353 University Road, Changqing District, Ji'nan, Shandong Province, No. 3501

Applicant before: Qilu University of Technology

Country or region before: China

CB02 Change of applicant information