Summary of the invention
The technical problem to be solved in the present invention is to provide a kind of method and device that obtains unequal section iron core excitation characteristic, can accurately obtain the excitation property of unequal section iron core.
The embodiment of the present invention provides a kind of method that obtains unequal section iron core excitation characteristic, comprising:
The fitting function of the BH curve of matching unequal section core material; B is magnetic induction density, and H is magnetic field intensity;
The magnetic circuit of unequal section iron core is decomposed, obtain the equivalent circuit of unequal section iron core, the equivalent field excitation branch line in described equivalent circuit is in parallel by a plurality of field excitation branch lines; The corresponding one section of magnetic circuit of field excitation branch line described in each;
Utilize the fitting function of described BH curve to obtain respectively ψ-i curve of field excitation branch line described in each; ψ is magnetic linkage, and i is electric current;
By ψ-i curve of all described field excitation branch lines, obtained ψ-i curve of unequal section iron core.
Preferably, the fitting function of described BH curve comprises linear term and saturated item;
The slope of described linear term equals air permeability μ
0, μ
0=4 π * 10
-7h/m;
Described saturated adopts asinh function, arctan function or ambipolar S compression function.
Preferably, when described unequal section iron core is magnet controlled fault current limiter, and magnet controlled fault current limiter adopts two independently magnetic circuits; The corresponding equivalent field excitation branch line of magnetic circuit described in each, described in each, equivalent field excitation branch line is decomposed into three field excitation branch lines and carries out parallel connection, described three field excitation branch lines respectively corresponding direct current post field excitation branch line, exchange post field excitation branch line and yoke field excitation branch line.
Preferably, when described unequal section iron core is magnet controlled fault current limiter, and magnet controlled fault current limiter adopts three pillar type magnetic circuit; Corresponding three the equivalent field excitation branch lines of described three pillar type magnetic circuit, wherein two described equivalent field excitation branch lines carry out parallel connection for exchanging post field excitation branch line and yoke field excitation branch line, and described in another, equivalent field excitation branch line is direct current post field excitation branch line.
Preferably, the described fitting function that utilizes described BH curve obtains respectively ψ-i curve of field excitation branch line described in each, is specially:
By sectional area and the length of segmentation magnetic circuit, obtain reduction to ψ-i curve of AC winding or DC side winding.
The embodiment of the present invention also provides a kind of device that obtains unequal section iron core excitation characteristic, comprising: fitting function acquiring unit, magnetic circuit resolving cell, segmentation magnetic circuit ψ-i curve acquisition unit and unequal section ψ-i curve acquisition unshakable in one's determination unit;
Described fitting function acquiring unit, for the fitting function of the BH curve of matching unequal section core material; B is magnetic induction density, and H is magnetic field intensity;
Described magnetic circuit resolving cell, for obtaining the equivalent circuit of unequal section iron core, the equivalent field excitation branch line in described equivalent circuit is in parallel by a plurality of field excitation branch lines; The corresponding one section of magnetic circuit of field excitation branch line described in each;
Described segmentation magnetic circuit ψ-i curve acquisition unit, for utilizing the fitting function of described BH curve to obtain respectively ψ-i curve of field excitation branch line described in each; ψ is magnetic linkage, and i is electric current;
Described unequal section ψ-i curve acquisition unshakable in one's determination unit, for ψ-i curve of the ψ-i curve acquisition unequal section iron core by all described field excitation branch lines.
Preferably, the fitting function of described BH curve comprises linear term and saturated item;
The slope of described linear term equals air permeability μ
0, μ
0=4 π * 10
-7h/m;
Described saturated adopts asinh function, arctan function or ambipolar S compression function.
Preferably, when described unequal section iron core is magnet controlled fault current limiter, and magnet controlled fault current limiter adopts two independently magnetic circuits, the corresponding equivalent field excitation branch line of magnetic circuit described in each; Described magnetic circuit resolving cell, is decomposed into three field excitation branch lines for equivalent field excitation branch line described in each and carries out parallel connection, and described three field excitation branch lines are corresponding direct current post field excitation branch line, interchange post field excitation branch line and yoke field excitation branch line respectively.
Preferably,
When described unequal section iron core is magnet controlled fault current limiter, and magnet controlled fault current limiter adopts three pillar type magnetic circuit; Described magnetic circuit resolving cell, be used for corresponding three the equivalent field excitation branch lines of described three pillar type magnetic circuit, wherein two described equivalent field excitation branch lines carry out parallel connection for exchanging post field excitation branch line and yoke field excitation branch line, and described in another, equivalent field excitation branch line is direct current post field excitation branch line.
Compared with prior art, the present invention has the following advantages:
Method provided by the invention and device, carry out magnetic circuit decomposition to the iron core of unequal section, obtains the equivalent circuit of unequal section iron core, and the equivalent field excitation branch line in described equivalent circuit is in parallel by a plurality of field excitation branch lines; The corresponding one section of magnetic circuit of field excitation branch line described in each; Utilize like this fitting function of BH curve of the core material of matching, obtain ψ-i curve of each field excitation branch line; By ψ-i curve of each field excitation branch line, obtained ψ-i curve of whole unequal section iron core, ψ-i curve can represent excitation property unshakable in one's determination again.Method provided by the invention is carried out magnetic circuit decomposition by the iron core of unequal section, is decomposed into a plurality of field excitation branch lines and carries out parallel connection, thereby obtain the excitation property of every section, and then by the excitation property equivalence of every section, gone out the excitation property of whole unequal section iron core.
Embodiment
Below in conjunction with the accompanying drawing in the embodiment of the present invention, the technical scheme in the embodiment of the present invention is clearly and completely described, obviously, described embodiment is only the present invention's part embodiment, rather than whole embodiment.Embodiment based in the present invention, those of ordinary skills, not making the every other embodiment obtaining under creative work prerequisite, belong to the scope of protection of the invention.
For above-mentioned purpose of the present invention, feature and advantage can be become apparent more, below in conjunction with accompanying drawing, the specific embodiment of the present invention is described in detail.
Referring to Fig. 1, this figure is embodiment of the method one process flow diagram of acquisition unequal section iron core excitation characteristic provided by the invention.
The method of the acquisition unequal section iron core excitation characteristic that the present embodiment provides, comprises the following steps:
S101: the fitting function of the BH curve of matching unequal section core material; B is magnetic induction density, and H is magnetic field intensity;
When iron core reach the degree of depth saturated after, owing to there is no rotating magnetic domain, its unit permeance is close to air.Therefore, the limit slope of the core material BH curve degree of depth after saturated is knownly (to equal air permeability μ
0=4 π * 10
-7h/m).
Fitting function can be divided into two parts: linear term and saturated item.The slope of linear term equals air permeability μ
0.The fitting function of saturated can have multiple, as asinh function, arctan function, ambipolar S compression function etc.Wherein ambipolar S compression function is usually used in error back propagation (BP, Error Back Propagation) algorithm.
S102: the magnetic circuit of unequal section iron core is decomposed, obtain the equivalent circuit of unequal section iron core, the equivalent field excitation branch line in described equivalent circuit is in parallel by a plurality of field excitation branch lines; The corresponding one section of magnetic circuit of field excitation branch line described in each;
For prismatic iron core, by BH curve, can obtain ψ-i curve.The number of turn of supposing coil is N, and average cross-section is A, and the length of center line is l, has following relation:
By formula (1), can obtain formula (2):
By formula (1), BH curve can be converted into Ψ-i curve.
But for the iron core of unequal section, first need iron core to be divided into some sections.
By formula (1) and (2), can only obtain Ψ-i curve of every iron leg heart.
By the excitation property of each iron leg heart, how to obtain the excitation property of whole unequal section iron core, be key point of the present invention.
S103: utilize the fitting function of described BH curve to obtain respectively ψ-i curve of field excitation branch line described in each; ψ is magnetic linkage, and i is electric current;
Utilize the fitting function of BH curve and the relation of formula (1) and (2) can obtain ψ-i curve of segmentation magnetic circuit.
S104: the ψ-i curve that is obtained unequal section iron core by ψ-i curve of all described field excitation branch lines.
The iron core of different structure, can adopt the nonlinear field excitation branch line of several parallel connections to represent conventionally.Yet, when employing electromagnetic transient state procedure (EMTP) etc. carries out simulation calculation, if directly a plurality of nonlinear field excitation branch lines are carried out to parallel connection, may produce some problems.For this reason, suggestion is by a plurality of non-linear exciter branch roads of parallel connection equivalent merging in addition.
The method that the present embodiment provides, carries out magnetic circuit decomposition to the iron core of unequal section, obtains the equivalent circuit of unequal section iron core, and the equivalent field excitation branch line in described equivalent circuit is in parallel by a plurality of field excitation branch lines; The corresponding one section of magnetic circuit of field excitation branch line described in each; Utilize like this fitting function of BH curve of the core material of matching, obtain ψ-i curve of each field excitation branch line; By ψ-i curve of each field excitation branch line, obtained ψ-i curve of whole unequal section iron core, ψ-i curve can represent excitation property unshakable in one's determination again.Method provided by the invention is carried out magnetic circuit decomposition by the iron core of unequal section, is decomposed into a plurality of field excitation branch lines and carries out parallel connection, thereby obtain the excitation property of every section, and then by the excitation property equivalence of every section, gone out the excitation property of whole unequal section iron core.
Embodiment bis-:
Introduce the fitting function of the BH curve of the matching unequal section core material that the embodiment of the present invention provides below.
Saturated in the fitting function that the present embodiment provides preferably adopts asinh function.
Consider the problem of fitting precision, saturated in the fitting function in the present invention adopts a plurality of asinh functions.The BH curve of certain core material of take is example, and its fitting function is suc as formula (3), and its raw data and fitting result are as shown in table 1.Fitting effect curve map as shown in Figure 2.
B=μ
0H-0.6299arsh(0.0280H)+0.6925arsh(0.280H) (3)
The fitting result of certain material B-H curve of table 1
Embodiment tri-:
How the unequal section iron core of take below carries out magnetic circuit decomposition as magnet controlled fault current limiter as example introduction.
Magnet controlled fault current limiter has two kinds of typical structures, is respectively to adopt two independently magnetic circuit and employing three pillar type structures.
Referring to Fig. 3, the iron core that this figure is magnet controlled fault current limiter provided by the invention adopts two independently schematic diagram of magnetic circuit.
According to the positive dirction of each physical quantity of stipulating in Fig. 3, there is following winding voltage equation:
Wherein, u
1the voltage at the first interchange winding two ends, u
2the voltage at the second interchange winding two ends, u
3the voltage at DC excitation winding two ends;
Because two windings are series relationship, so the first current i exchanging in winding
1exchange the current i of winding with second
2equate:
i
1=i
2 (5)
According to the Kirchhoff's second law of magnetic circuit, for two iron circuits, can obtain:
H
11l
1+H
12l
2+H
13l
3=N
1i
1+N
3i
3
(6)
H
21l
1+H
22l
2+H
23l=N
2i
2-N
3i
3
Wherein, H
11and H
21be respectively the magnetic field intensity in the first magnetic circuit and AC core limb in the second magnetic circuit (be called for short and exchange post), H
13and H
23be respectively the magnetic field intensity in DC side core limb in the first magnetic circuit and the second magnetic circuit (being called for short direct current post), and H
12and H
22be respectively the magnetic field intensity of yoke part in the first magnetic circuit and the second magnetic circuit; l
1, l
2, l
3be respectively the centerline length that exchanges post, yoke, direct current post.N
1the number of turn of the first interchange winding, N
2the number of turn of the second interchange winding, N
3it is the number of turn of direct current winding; Φ
1the main flux of the first magnetic circuit, Φ
2it is the main flux of the second magnetic circuit; ; i
3it is the electric current in direct current winding.
If arranged:
And meet following constraint condition
N
1i′
1+N
3i
3=0
(8)
N
2i′
2-N
3i
3=0
Can obtain the magnetic potential drop expression formula of two iron circuits:
According to above-mentioned winding voltage equation, can obtain the equivalent circuit of the magnet controlled fault current limiter of two independent magnetic circuit formulas, the equivalent field excitation branch line in described equivalent circuit is in parallel by a plurality of field excitation branch lines, the corresponding one section of magnetic circuit of field excitation branch line described in each.
As shown in Figure 4.Wherein, r
1and L
k1for exchanging resistance and the leakage inductance of winding 1, r
2and L
k2for exchanging resistance and the leakage inductance of winding 2, r
3and L
k3resistance and leakage inductance for DC excitation winding; " equivalent field excitation branch line 1 " and " equivalent field excitation branch line 2 " represents respectively the excitation property of two iron circuits.
Under the unequal condition of each section of sectional area of magnetic circuit, magnetic circuit can be decomposed.Order
And meet following condition:
Like this, each iron core is broken down into three parts: exchange post, direct current post and yoke.The equivalent circuit of equivalent field excitation branch line 1 and equivalent field excitation branch line 2 is respectively as shown in Fig. 5 a and Fig. 5 b.From circuit, equivalent field excitation branch line 1 and equivalent field excitation branch line 2 form by three field excitation branch lines in parallel.Described three field excitation branch lines in parallel are corresponding direct current post field excitation branch line, interchange post field excitation branch line and yoke field excitation branch line respectively.
Be more than to take the magnetic circuit that magnet controlled fault current limiter is two independent magnetic circuit structures to decompose, the magnetic circuit of introducing magnet controlled fault current limiter below and be three pillar type structure decomposes.
Magnet controlled fault current limiter is three pillar type structure, and each physical quantity and positive dirction as shown in Figure 6, according to the positive dirction shown in Fig. 6, have following winding voltage equation:
And have:
i
1=i
2
(13)
Φ
3=Φ
1-Φ
2
Kirchhoff's second law according to magnetic circuit, has
H
11l
1+H
12l
2+H
3l
3=N
1i
1+N
3i
3
(14)
H
21l
1+H
22l
2-H
3l
3=N
2i
2-N
3i
3
Wherein, H
11and H
21be respectively the magnetic field intensity in AC core limb in the first magnetic circuit and the second magnetic circuit, H
3for the magnetic field intensity in DC side core limb, and H
12and H
22be respectively the magnetic field intensity of yoke part in the first magnetic circuit and the second magnetic circuit; l
1, l
2, l
3be respectively the centerline length that exchanges post, yoke, direct current post.N
1the number of turn of the first interchange winding, N
2the number of turn of the second interchange winding, N
3it is the number of turn of direct current winding; Φ
1that the first magnetic circuit exchanges the main flux in post and yoke, Φ
2that the second magnetic circuit exchanges the main flux in post and yoke; Φ
2it is the main flux in direct current post; Φ
k1the leakage flux of the first interchange winding, Φ
k2it is the leakage flux of the second interchange winding; i
3it is the electric current in DC excitation winding.
By formula (14), can be obtained:
H
11l
1+H
12l
2=N
1i
1+(N
3i
3-H
3l
3)
(15)
H
21l
1+H
22l
2=N
2i
2-(N
3i
3-H
3l
3)
Magnetic circuit is decomposed:
And meet following constraint condition:
Can obtain the magnetic potential drop expression formula of all the other magnetic circuits:
According to above-mentioned winding voltage equation, can obtain the equivalent circuit of this restrictor, as shown in Figure 7.
Wherein, " equivalent field excitation branch line 1' " refers to excitation property corresponding to magnetic circuit that left side core limb (having the first interchange winding on it) and upper magnet yoke, lower yoke form, and " equivalent field excitation branch line 2' " represents excitation property corresponding to magnetic circuit that right side core limb (having the second interchange winding on it) and upper magnet yoke, lower yoke form; " equivalent field excitation branch line 3 " represents the excitation property corresponding to magnetic circuit of middle core limb (having direct current winding on it).
Under the inhomogeneous condition of each section of sectional area of magnetic circuit, also can, on the basis of formula (18), magnetic circuit be decomposed.Order
And meet following condition:
Like this, equivalent field excitation branch line 1 ' with equivalent field excitation branch line 2 ' be all broken down into two parts: exchange post field excitation branch line and yoke field excitation branch line.From circuit, be two field excitation branch lines in parallel.Equivalent field excitation branch line 1 ', 2 ' and 3 equivalent magnetic circuit respectively as shown in Fig. 8 a-8c.
From circuit, equivalent field excitation branch line 1 ' and 2 ' by two field excitation branch lines in parallel, forming, correspondence exchanges post field excitation branch line and yoke field excitation branch line respectively.
Certain core material take below as example, introduce ψ-i curve of segmentation magnetic circuit;
ψ-i curvilinear function that can be obtained each section of magnetic circuit by formula (2) and formula (3) is as follows:
According to the length of the sectional area of each section of magnetic circuit and center line, can obtain reduction to ψ-i curve of AC winding or DC side winding.
The iron core of magnet controlled fault current limiter of take below adopt two-way independently structure be introduced as example.
The reduction of take is introduced to direct current winding side as example:
Its number of turn N=868; The sectional area that exchanges post, yoke, direct current post is respectively 0.43m
2, 0.77m
2, 0.86m
2, its centerline length is respectively 2.94m, 5.82m, 2.94m.
According to formula (21), can obtain:
ψ-i the curve that exchanges post can be expressed as:
ψ=0.1385i-235.1arsh(8.267i)+258.5arsh(82.67i) (22)
In like manner, ψ-i curve of yoke can be expressed as:
ψ=0.1253i-421.0arsh(4.176i)+462.8arsh(41.76i) (23)
ψ-i curve of direct current post can be expressed as:
ψ=0.2769i-470.2arsh(8.267i)+516.9arsh(82.67i) (24)
Field excitation branch line is carried out to equivalence: by add certain voltage on three shunt excitation branch roads, solve respectively the branch current on each field excitation branch line, three current summations can be obtained to total current; By branch current, by ψ-i curve, can obtain corresponding magnetic linkage value.
The iron core of magnet controlled fault current limiter of take adopt two-way independently structure be example, the ψ-i characteristic after equivalence is as shown in table 2, corresponding curve as shown in Figure 9.
For comparison curves characteristic, in the embodiment of the present invention, provided the ψ-i characteristic that exchanges post, yoke and direct current post.Can find out, the ψ-i curve after equivalence is more approaching with the ψ-i curve that exchanges post.Because it is minimum to exchange the sectional area of post.When magnetic flux one timing, the degree of saturation that exchanges post is higher, therefore exciting current Main Current is crossed the corresponding field excitation branch line of interchange post.Obviously, the magnetic circuit situation of interchange post winding plays conclusive effect to iron core excitation characteristic.
The above embodiment of the present invention carries out segmentation by unequal section iron core, has proposed the matching of multistage iron core excitation characteristic and equivalent method, thereby obtains the excitation property of whole unequal section iron core.
The method of the excitation property of the acquisition unequal section iron core providing based on above embodiment, the present invention also provides a kind of device that obtains the excitation property of unequal section iron core, below in conjunction with accompanying drawing, introduces in detail its principle of work,
Device embodiment mono-:
Referring to Figure 10, this figure is device embodiment mono-schematic diagram of acquisition unequal section iron core excitation characteristic provided by the invention.
The embodiment of the present invention provides a kind of device that obtains unequal section iron core excitation characteristic, comprising: fitting function acquiring unit 1001, magnetic circuit resolving cell 1002, segmentation magnetic circuit ψ-i curve acquisition unit 1003 and unequal section ψ-i curve acquisition unshakable in one's determination unit 1004;
Described fitting function acquiring unit 1001, for the fitting function of the BH curve of matching unequal section core material; B is magnetic induction density, and H is magnetic field intensity;
When iron core reach the degree of depth saturated after, owing to there is no rotating magnetic domain, its unit permeance is close to air.Therefore, the limit slope of the core material BH curve degree of depth after saturated is knownly (to equal air permeability μ
0=4 π * 10
-7h/m).
Fitting function can be divided into two parts: linear term and saturated item.The slope of linear term equals air permeability μ
0.The fitting function of saturated can have multiple, as asinh function, arctan function, ambipolar S compression function etc.Wherein ambipolar S compression function is usually used in error back propagation (BP, Error Back Propagation) algorithm.
Described magnetic circuit resolving cell 1002, for obtaining the equivalent circuit of unequal section iron core, the equivalent field excitation branch line in described equivalent circuit is in parallel by a plurality of field excitation branch lines; The corresponding one section of magnetic circuit of field excitation branch line described in each;
For prismatic iron core, by BH curve, can obtain ψ-i curve.The number of turn of supposing coil is N, and average cross-section is A, and the length of center line is l, has following relation:
By formula (1), can obtain formula (2):
By formula (1), BH curve can be converted into Ψ-i curve.
But for the iron core of unequal section, first need iron core to be divided into some sections.
By formula (1) and (2), can only obtain Ψ-i curve of every section of magnetic circuit.
By the excitation property of each section of magnetic circuit, how to obtain the excitation property of whole unequal section iron core, be key point.
Described segmentation magnetic circuit ψ-i curve acquisition unit 1003, for utilizing the fitting function of described BH curve to obtain respectively ψ-i curve of field excitation branch line described in each; ψ is magnetic linkage, and i is electric current;
Utilize the fitting function of BH curve and the relation of formula (1) and (2) can obtain ψ-i curve of segmentation magnetic circuit.
Described unequal section ψ-i curve acquisition unshakable in one's determination unit 1004, for ψ-i curve of the ψ-i curve acquisition unequal section iron core by all described field excitation branch lines.
The iron core of different structure, can adopt the nonlinear field excitation branch line of several parallel connections to represent conventionally.Yet, when employing electromagnetic transient state procedure (EMTP) etc. carries out simulation calculation, if directly a plurality of nonlinear field excitation branch lines are carried out to parallel connection, may produce some problems.For this reason, suggestion is by a plurality of non-linear exciter branch roads of parallel connection equivalent merging in addition.
The device that the present embodiment provides, carries out magnetic circuit decomposition to the iron core of unequal section, obtains the equivalent circuit of unequal section iron core, and the equivalent field excitation branch line in described equivalent circuit is in parallel by a plurality of field excitation branch lines; The corresponding one section of magnetic circuit of field excitation branch line described in each; Utilize like this fitting function of BH curve of the unequal section core material of matching, obtain ψ-i curve of each field excitation branch line; By ψ-i curve of each field excitation branch line, obtained ψ-i curve of whole unequal section iron core, ψ-i curve can represent excitation property unshakable in one's determination again.Method provided by the invention is carried out magnetic circuit decomposition by the iron core of unequal section, is decomposed into a plurality of field excitation branch lines and carries out parallel connection, thereby obtain the excitation property of every section, and then by the excitation property equivalence of every section, gone out the excitation property of whole unequal section iron core.
Saturated in the fitting function that the present embodiment provides preferably adopts asinh function.
Consider the problem of fitting precision, saturated in the fitting function in the present invention adopts a plurality of asinh functions.
B=μ
0H-0.6299arsh(0.0280H)+0.6925arsh(0.280H) (3)
Introduce ψ-i curve of segmentation equivalent magnetic circuit below;
ψ-i curvilinear function that can be obtained each section of magnetic circuit by formula (2) and formula (3) is as follows:
According to the length of the sectional area of each section of magnetic circuit and center line, can obtain reduction to ψ-i curve of AC winding or DC side winding.
When described unequal section iron core is magnet controlled fault current limiter, and magnet controlled fault current limiter adopts two independently magnetic circuits; The corresponding equivalent field excitation branch line of each magnetic circuit, this field excitation branch line is decomposed into three field excitation branch lines and carries out parallel connection, and described three field excitation branch lines are corresponding direct current post magnetic circuit, interchange post magnetic circuit and yoke magnetic circuit respectively.
When described unequal section iron core is magnet controlled fault current limiter, and magnet controlled fault current limiter adopts three pillar type magnetic circuit; Three pillar type magnetic circuit is decomposed into three equivalent field excitation branch lines, and wherein two equivalent field excitation branch lines are for exchanging field excitation branch line and field excitation branch line parallel connection corresponding to yoke that post is corresponding, and another equivalent field excitation branch line is the field excitation branch line that direct current post is corresponding.
Magnet controlled fault current limiter has structure and two kinds of typical types of three pillar type structure of two independent magnetic circuits, in embodiment of the method, specifically introduced and how for these two types, to have carried out magnetic circuit decomposition, in device embodiment, no longer specifically introduce, can reference method embodiment content partly.
The above, be only preferred embodiment of the present invention, not the present invention done to any pro forma restriction.Although the present invention discloses as above with preferred embodiment, yet not in order to limit the present invention.Any those of ordinary skill in the art, do not departing from technical solution of the present invention scope situation, all can utilize method and the technology contents of above-mentioned announcement to make many possible changes and modification to technical solution of the present invention, or be revised as the equivalent embodiment of equivalent variations.Therefore, every content that does not depart from technical solution of the present invention,, all still belongs in the scope of technical solution of the present invention protection any simple modification made for any of the above embodiments, equivalent variations and modification according to technical spirit of the present invention.