Denver PELS 20070410 Hesterman Magnetic Coupling
Denver PELS 20070410 Hesterman Magnetic Coupling
Denver PELS 20070410 Hesterman Magnetic Coupling
SLIDE #
SLIDE #
Presentation Outline
Introduction Modeling magnetic coupling with electric circuit equations Measuring electric circuit model parameters Equivalent circuits for transformers and coupled inductors Magnetic circuit modeling overview Tips for creating magnetic circuit models Deriving electric model parameters from magnetic model parameters Matrix theory requirements for coupling stability Examples
SLIDE #
SLIDE #
SLIDE #
SLIDE #
SLIDE #
Time domain
di v=L dt
Frequency domain
v = jLi
SLIDE #
L11 = Self inductance of winding 1 L22 = Self inductance of winding 2 L12 = L21 = Mutual inductance
SLIDE #
v2 = L12
v1 L11 v L 2 = 21 v N LN 1
L12 L 22 LN 2
L1N i1 L 2 N d i2 dt LN N iN
d [v] = [L] [ i ] dt
SLIDE #
10
v1 = jL11 i1 + jL12 i2
v2 = j L21i1 + j L22 i2
SLIDE #
11
v N = j LN 1i1 + j LN 2 i2 + j LN 3i3 ++ j LN N iN
v1 L11 v L 2 = j 21 v N LN 1
L12 L 22 LN 2
L1N i1 L 2N i2 LN N iN
[v] = j[L][ i ]
SLIDE #
12
Lqr = Lrq
L. O. Chua, Linear and Nonlinear Circuits. New York: McGraw-Hill, 1987, pp. 771-780 This principle can be derived from Maxwells equations: C. G. Montgomery, R. H. Dicke, and E. M. Purcell, Ed., Principles of Microwave Circuits. New York: Dover Publications, 1965 It can also be derived from a stored energy argument: R. R. Lawrence, Principles of Alternating Currents. New York: McGrawHill, 1935, pp. 187-188
SLIDE #
13
k=
L12 = L11L22
L21 = L11L22
M L1L2
1 k 1
SLIDE #
14
k qr =
1 k13 = L 21 k 23 L 22 L11 1 L 31 L L 33 11
SLIDE #
15
L12 =
Laid Lopp 4
L2
k=
L1
series-aiding, Laid
series-opposing, Lopp
SLIDE #
16
v1 = j L11i1 + j L12 i2
v2 = j L12 i1 + j L 22 i2
1 k12 1
Write down what is known for the series-aiding configuration
i2 = i1
i1
L1
i2
L2
vaid = v1 + v2
vaid = j Laid i1
series-aiding, Laid
SLIDE #
17
v1 = j L11i1 + jL12 i1
v2 = j L 21i1 + j L22 i1
Recall Substitute Simplify Recall Therefore
vaid = v1 + v2
SLIDE #
18
i2 = i1
i1
L1
i2
L2
vopp = v1 v2
vopp = j Lopp i1
series-opposing, Lopp
v1 = j L11i1 jL12 i1
v2 = j L21 i1 j L22 i1
Recall Substitute
vopp = v1 v2
SLIDE #
19
Lopp = L11 2 L12 + L 22 Laid = L11+ 2 L12 + L 22 Laid Lopp = 4L12 L12 =
k=
Laid Lopp
L12 L11L22
Therefore
k=
SLIDE #
20
k=
voc 21 voc12 vd 12 vd 21
Measurements can be corrupted by capacitance and loading Particularly suited for measurements between windings on different legs of three-phase ungapped transformers
SLIDE #
21
v1 = j L11 i1 + jL12 i2
v2 = jL12 i1 + j L22 i2
Find i1 and v2 when winding 1 is driven and winding 2 is unloaded
v1 = vd 12 = j L11 i1
v2 = voc 21 = j L12 i1
Substitute i1, simplify and rearrange
i1=
vd 12 j L11
voc 21 = j L12
vd 12 j L11
SLIDE #
22
v2 = vd 21 = j L22 i2
v1 = voc12 = j L12 i2
Substitute i2, simplify and rearrange Recall
i2 =
vd 21 j L22
vd 12 j L22
k=
voc 21 voc12 vd 12 vd 21
Multiply
SLIDE #
23
voc 21 voc12 = vd 12 vd 21
L12 =k L11L22
Therefore
voc 21 voc12 k= vd 12 vd 21
SLIDE #
24
k 21 = 1
Lleak 21 L22
SLIDE #
25
v1 = j L11i1 + jL12 i2
v2 = jL12 i1 + j L 22 i2
Find the leakage inductance at winding 1 when winding 2 is shorted, L12s Write down what is known for that condition
v1 = j Lleak 12 i1
v2 = 0 = j L12 i1 + j L22 i2
Simplify
0 = L12 i1 + L22 i2
SLIDE #
26
Solve for i2
L12 i2 = i1 L 22
Simplify
Substitute value of i2
L2 Lleak 12 = L11 12 L 22
SLIDE #
27
Divide by L11
Rearrange
SLIDE #
28
SLIDE #
29
SLIDE #
30
SLIDE #
31
SLIDE #
32
k = 1
Lleak 12 L11
Generally use at least four significant digits when the coupling coefficient is close to 1 (k is known to more significant digits than the measured inductances it was derived from)
SLIDE #
33
1 k 1
Pin-for-pin equivalent electrical behavior
SLIDE #
34
k12 k 23 k31 0
SLIDE #
35
k
L11 L22
1
La Lb
1: Ne
2
Ideal
La = Lleak12
Measure the inductance at winding 1 with winding 2 shorted
Lb = L11 La
Measure the self inductance of winding 1 and subtract La
L 22 Ne = Lb
Measure the self inductance of winding 2 and calculate above equation
La Lb
1: Ne
2
Ideal
If winding 2 is shorted, then the inductance measured at winding 1 is simply La because the ideal transformer shorts out Lb. Thus La is equal to the leakage inductance measured at winding 1.
La = Lleak12
To convert the coupled inductor model to the cantilever model, recall that Therefore
Lleak12 = L11 1 k 2
La = L11 1 k 2
k
L11 L22
1
La Lb
1: Ne
2
Ideal
Lb = L11 La
To convert the coupled inductor model to the cantilever model, recall that Substitute Therefore
La = L11 1 k 2
k
L11 L22
1
La Lb
1: Ne
2
Ideal
The inductance of winding 2, L22 , is equal to Lb reflected through the ideal transformer
L 22 = Ne2 Lb
L 22 Ne = Lb
Lb = k 2 L11
1 Ne = k
L 22 L11
SLIDE #
39
Ll1
N1 :N 2
Ll 2
LMag
Ideal
N1 N 2
SLIDE #
40
v1 = j L11i1 + jL12 i2
v2 = j L12 i1 + j L 22 i2
Write down the open-circuit primary voltage due to a current in the secondary
N2 v1 = j L12 i2 = j LMag N i2 1
N2 L12 = LMag N 1
SLIDE #
41
N1 L Mag = N 2
L12
2
Recall
N2 L 22 = Ll 2 + LMag N 1
N2 Ll 2 = L 22 LMag N 1
12
Recall the equation for Lleak12 from the derivation of the Self and Leakage Inductance model L2
Lleak 12 = L11
L 22
By symmetry,
SLIDE #
42
Recall
Therefore
Also,
SLIDE #
43
SLIDE #
44
R2 R1 R13 R7 R8
R3 R15
R4
R5 R6
R1
MMF1
MMF3
R17
R18 R16
R14
R18 R13 R16 R7 R8 R9 R10 R11 R12 R14
R9
Winding 1
Winding 2
Winding 3
SLIDE #
45
R1
MMF1
MMF2 R19 MMF3
Rth1
MMF1
N2
R17 R18
N1
R13
N3
R14
N1 short short
N12 L1 = th1
Nx = turns of winding x
SLIDE #
46
Leakage Inductances
Replace the MMF source for each winding not being considered with a short circuit Replace the MMF source for the shorted winding with an open circuit Determine Thvenin equivalent reluctance, th , at the MMF source for the winding where the leakage inductance is to be determined Inductance = turns squared divided by the Thvenin equivalent reluctance
R2 R3 R4 R5 R6 R15 R1
MMF1
Rth1
MMF1
N2
R17 R18
MMF2 R19
MMF3
N1 open
N1
R13
N3
R14
short
Lleak 13
N12 = th1
Nx = turns of winding x
SLIDE #
47
SLIDE #
48
Energy Storage
The magnetic energy stored in one inductor is:
1 2 WM 1 = Li 2
The energy stored in a set of N coupled windings is:
WM N
1 T = [i] [L ][i] 2
SLIDE #
49
1 k 1
1 k12 1
1 k 23 1
1 k31 1
Yilmaz Tokad and Myril B. Reed, Criteria and Tests for Realizability of the Inductance Matrix, Trans. AIEE, Part I, Communications and Electronics, Vol. 78, Jan. 1960, pp. 924-926
SLIDE #
50
k12 1 k 1 21 K= k N 1 k N 2
k1N k2 N 1
k qr = k rq
A set of coupled inductors is passive if and only if all of the eigenvalues of K are non-negative There are N eigenvalues for a set of N windings Eigenvalue calculations are built-in functions of programs like Mathcad and Matlab
SLIDE #
51
Consistency Checks
The eigenvalue test can let you know if there are measurement errors, but it wont help identify the errors The ratios of magnetizing inductances on the same core leg should be approximately equal to the square of the turns ratios Set up test simulations to verify that the models match the test conditions Test leakage inductances by shorting one winding and applying a signal to the other winding (Compute the inductance indicated by the voltage, current and frequency) Check for coupling polarity errors, especially if there are windings on multiple core legs Compare test data with multiple windings shorted to simulations or calculations with multiple windings shorted
SLIDE #
52
= L1
Each diagonal element is equal to the reciprocal of the inductance of the corresponding winding when all of the of the other windings are shorted Example application: You have a four-winding transformer. What is the inductance at winding 1 when windings 3 and 4 are shorted? Set up an inductance matrix, Lx, by extracting all of the elements from the total inductance matrix that apply only to windings 1, 3 and 4
SLIDE #
53
L14 L 24 L 34 L44
x = L x1
SLIDE #
54
SLIDE #
55
Assign node numbers The reference node is node 0 Assign currents for each branch Compute mutual inductances from the self inductances and coupling coefficients
L12 = k L11 L 22
SLIDE #
56
Write a KCL equation using the node-to-datum voltages as variables for each node other than the datum node, unless the node has a fixed voltage with respect to the datum node. In that case, write an equation that assigns the node voltage to the fixed value.
SLIDE #
57
v 1 = vs i1 = j C1(v1 v2 ) = j C1(vs v2 )
1 v3 i2 = + j C 2 v3 = + j C 2 v3 R R1 1
Current variables must be maintained for each coupled inductor Other current variables can be replaced by equivalent expressions
SLIDE #
58
Write an equation for the voltage dropped across each coupled inductor with the node variables on one side and the inductance and current terms on the other side
v2 = j L11i1 jL12 i2
v3 = j L12 i1 j L22 i2
Note that the minus signs are due to the fact that i2 was assigned to be flowing out of the dotted end of the winding instead of into it
SLIDE #
59
i1 = j C1(vs v2 )
1 v2 + i1 = vs j C1
1 i2 = + j C 2 v3 R 1
v2 = j L11 i1 jL12 i2
1 + j C 2 v3 i2 = 0 R 1 v2 j L11 i1 + j L12 i2 = 0
v3 = j L12 i1 j L22 i2
v3 j L12 i1 + j L22 i2 = 0
SLIDE #
60
G X =U
1 v2 + i1 = vs 1 j C1
3 1
1 2 R + j C 2 v3 i2 = 0 1
4
v2 j L11 i1 + jL12 i2 = 0
v3 j L12 i1 + j L22 i2 = 0
2 3 4
1 0 1 0
1 + j C 2 0 R 1 0 j L11 1 jL12
1 j C1
0 v2 v s v3 0 1 = i1 0 jL12 i 2 0 j L22 i2
SLIDE #
61
SLIDE #
62
Suppose we have a three-winding inductor with a resistor terminating each winding as is shown in Fig. 1. If we externally force an initial condition of currents in the windings, and then let the inductor "coast" on its own, the currents should all exponentially decay to zero as the stored energy is dissipated. If the coupling coefficient matrix is not positive definite, however, at least one of the eigenvalues for the system will be not be negative.
SLIDE #
63
As is shown below, this produces an unstable situation where more energy can be extracted from the inductor than was initially stored there. When the initial energy is dissipated, the stored energy becomes negative, and the inductor continues to deliver power. This, of course, is impossible. Therefore, a set of coupling coefficients in which the coupling coefficient matrix is not positive definite describes a coupled inductor that is not physically realizable. Because Spice and other circuit simulators allow such inductors to be specified, it is instructive to examine what happens when nonphysical sets of coupling coefficients are specified.
SLIDE #
64
We first examine a physically-realizable case. After that, we show what happens when the coupling coefficients are improperly chosen. Enter the values of the inductances and coupling coefficients L1 := 10 H L2 := 11 H L3 := 10 H k12 := 0.96 k23 := 0.99 k13 := 0.98
SLIDE #
65
All of the eigenvalues are positive. As is shown below, this leads to negative eigenvalues for the time response of the system.
SLIDE #
66
L1 K L1 L2 K L1 L3 1, 2 1, 3 L2 K L2 L3 L := K1 , 2 L1 L2 2, 3 K L1 L3 K L2 L3 L3 2, 3 1, 3
R1 := 1
R2 := 1
R3 := 1
SLIDE #
67
I1i := 1 A
I2i := 0 A
I3i := 0 A
I1 I I2 I 3 I1i Ii := I2i I 3i 1 Ii = 0 A 0
We can think of this as a situation in which the current in winding 1 is externally forced to be 1 Ampere, and then the circuit is left to dissipate the stored energy starting at time t = 0.
SLIDE #
68
Ii L Ii = 1.000 10
The time response of the circuit of Fig. 1 can be described by the following equation
I1 R1 I2 R2 I R 3 3
d L I dt
(1)
:= L
SLIDE #
69
(2)
(3)
SLIDE #
70
Define G
R1 0 0 G := 0 R2 0 0 0 R 3
d I dt G I
(4)
(5)
:= eigenvals ( G)
SLIDE #
71
0.295 0.77 0.565 -1 := eigenvecs ( G) = 0.501 0.628 0.595 s 0.814 0.107 0.571 e1 t The solution to (5) will have the form c1 0 0 I 0 c2 0 e 2 t (6) 0 0 c 3 3 t e
SLIDE #
72
The values of the constants c1, c2 and c3 can be determined from the initial conditions
I0
c1 c2 c 3
c1 1 c2 := Ii c 3
c1 1 c2 = 0 A c 0 3
1 Ii = 0 A 0
SLIDE #
73
c1 0 0 C := 0 c2 0 0 0 c 3
(Note that the minimum index value for matrices in this document is set at 1, not zero.) Define current functions I1( t) := C I2( t) := C I3( t) := C e e e
1 t
1, 1
+C +C +C
1, 2
e e e
2 t
+C +C +C
1, 3
e e e
3 t
2, 1
1 t
2, 2
2 t
2, 3
3 t
3, 1
1 t
3, 2
2 t
3, 3
3 t
SLIDE #
74
Check the current values at a few points I1( 0 s ) = 1 A I2( 0) = 0 A I3( 0) = 0 A I1( 0.1 s ) = 0.803A I2( 0.1 s ) = 0.021A I3( 0.1 s ) = 0.214A I1( 10 s ) = 0.23A I2( 10 s ) = 0.242A I3( 10 s ) = 0.233A
SLIDE #
75
t := 0 , 10
,
s .. 20 10
0.5 0
5 .10
1 .10 t
1.5 .10
2 .10
I1 drops quickly as the current builds up in the other two windings. All of the currents then decay.
SLIDE #
76
I1( t) I( t) := I2( t) I ( t) 3
W ( t) :=
1 2
( I( t) ) L I( t)
SLIDE #
77
t := 0 , 10
5 .10 4 .10 W( t ) 3 .10 2 .10
7
6
s .. 10
2 .10
4 .10
6 .10 t
8 .10
1 .10
SLIDE #
78
We can now examine a case where the coupling coefficient are improperly specified. Are relatively small change in one of the couple coefficients is all that is required to create an unstable configuration.
L1 := 10 H L2 := 11 H L3 := 10 H
SLIDE #
79
One of the eigenvalues is negative. As is shown below on, this makes one of the eigenvalues for the time response of the system positive.
SLIDE #
80
K L1 L2 K L1 L3 L1 1, 2 1, 3 L2 K L2 L3 L := K1 , 2 L1 L2 2, 3 K L1 L3 K L2 L3 L3 1, 3 2, 3
R1 := 1
R2 := 1
R3 := 1
SLIDE #
81
I1i := 1 A
I2i := 0 A
I3i := 0 A
I1i Ii := I2i I 3i
1 Ii = 0 A 0
Ii L Ii = 1.000 10
SLIDE #
82
:= L
SLIDE #
83
Define G
R1 0 0 G := 0 R2 0 0 0 R 3
(4)
SLIDE #
84
0.414 0.713 0.566 -1 := eigenvecs ( G) = 0.395 0.701 0.594 s 0.82 0.023 0.572 e1 t The solution to (5) will have the form c1 0 0 I 0 c2 0 e 2 t (6) 0 0 c 3 3 t e
SLIDE #
85
The values of the constants c1, c2 and c3 can be determined from the initial conditions
I0
c1 c2 c 3
c1 1 c2 := Ii c 3
c1 1 c2 = 0 A c 0 3
1 Ii = 0 A 0
SLIDE #
86
c1 0 0 C := 0 c2 0 0 0 c 3
1, 1
e e e
+C +C +C
1, 2
e e e
+C +C +C
1, 3
e e e
2, 1
1 t
2, 2
2 t
2, 3
3 t
3, 1
1 t
3, 2
2 t
3, 3
3 t
SLIDE #
87
I1( 0.01 s ) = 4 10 A
5 5
SLIDE #
88
t := 0 , 10
11
s .. 4.2 10
4 0 5 .10
10
1 .10 t
1.5 .10
I1 rises instead of decaying. The other currents start at zero and then build up.
SLIDE #
89
I1( t) I( t) := I2( t) I ( t) 3
1 2
T
W ( t) :=
( I( t) ) L I( t)
SLIDE #
90
5 .10
W( t )
5 .10
1 .10
2 .10
3 .10 t
4 .10
5 .10
Acknowledgement
Thanks to my manager, Dr. David J. Christie, for reviewing the presentation and providing helpful comments and suggestions Also, thanks to Dr. James H. Spreen for his continuing guidance
SLIDE #
92