www.ietdl.org
Published in IET Power Electronics Received on 31st January 2010 Revised on 14th July 2010 doi: 10.1049/iet-pel.2010.0040
ISSN 1755-4535 Transformer winding loss caused by skin and proximity effects including harmonics in pulse-width modulated DC–DC flyback converters for the continuous conduction mode D. Murthy-Bellur1 N. Kondrath2 M.K. Kazimierczuk3 1School of Engineering, Pennsylvania State University, The Behrend College, Erie, PA 16563, USA 2Department of Electrical and Computer Engineering, University of Minnesota Duluth, Duluth, MN 55812, USA 3Department of Electrical Engineering, Wright State University, Dayton, OH 45435, USA E-mail: [email protected]
Abstract: High-frequency (HF) transformers used in pulse-width-modulated power converters conduct periodic non-sinusoidal HF currents, which give rise to additional winding losses caused by harmonics. This study presents expressions for winding power losses in a two-winding transformer with periodic non-sinusoidal primary and secondary current waveforms for flyback converters operating in the continuous conduction mode (CCM). Dowell’s equation, which takes into account both the skin and proximity effects, is used to determine the primary and secondary winding resistances as functions of frequency. Fourier series of the primary and secondary current waveforms and the winding resistances are used to determine the winding power losses at harmonic frequencies. The harmonic primary and secondary winding loss factors FRph and FRsh are introduced for CCM operation. The theory is illustrated by performance evaluation of the transformer in a flyback converter in CCM over the entire range of converter operation. Power losses in the primary and secondary transformer windings are illustrated as functions of the output power and the DC input voltage.
1 Introduction [16], the AC resistances for a number of arbitrary waveforms were determined using the rms values of the current High-frequency (HF) transformers are an integral part of pulse- waveforms and the rms values of its derivatives. Methods to width-modulated (PWM) converters, which are widely used in design inductors and transformers for PWM converters are many applications such as power-factor correctors, DC–DC available in [30–32]. However, the transformer design in [30] power converters and electromagnetic interference (EMI) filters and other publications do not account for winding losses [1–5]. For low-power applications, flyback PWM DC–DC caused neither by HF skin and proximity effects nor by converter is a popular choice because of its low parts count harmonics. Until now, the Fourier expansion of the flyback when compared with other isolated DC–DC converter transformer primary and secondary current waveforms in topologies. The HF transformer in a flyback converter is used PWM converters operating in the continuous conduction to provide galvanic isolation and also to store the magnetic mode (CCM) as shown in Fig. 1 has not been presented in the energy required for power conversion. Operating the flyback literature. converter at high switching frequency reduces the size of the The objectives of this paper are (i) to develop winding power transformer, but increases the transformer winding losses loss equations for a two-winding transformer conducting caused by skin and proximity effects. Additionally, since periodic non-sinusoidal HF currents containing DC transformers in isolated PWM converters conduct periodic non- components and harmonics such as those in flyback PWM sinusoidal currents in steady state, the spectra of the primary converters [33–35] operating in CCM, using Dowell’s theory and secondary current waveforms of the transformer are very and the methods developed in [36–38], (ii) to give a step-by- rich and the losses caused by the harmonics further increase the step procedure to design a two-winding transformer for winding power losses. The harmonic losses in transformer flyback converter operated in CCM, using the area product windings at HF also increase the operating temperature. method and (iii) to study the effect of harmonics on winding HF effects in windings and cores were studied in several losses in the transformer of a flyback PWM DC–DC papers [6–29]. The winding AC resistance caused by the skin converter in CCM over the entire range of converter operation. and proximity effects for sinusoidal current in a winding Sections 2 and 3 present the derivations of core area was presented in [6–8]. The theoretical analyses presented in product and winding power loss equations for a flyback [6–8] were further extended in [9–15]. In another method transformer in CCM conducting periodic non-sinusoidal
IET Power Electron., 2011, Vol. 4, Iss. 4, pp. 363–373 363 doi: 10.1049/iet-pel.2010.0040 & The Institution of Engineering and Technology 2011 www.ietdl.org
where Awp is the bare wire cross-sectional area of the primary winding in m2. The window area limited by the maximum current density in both primary and secondary windings is given by
N A + N A = p wp s ws Wa (4) Ku
where Ns is the number of turns in the secondary winding, Aws is the secondary winding bare wire cross-sectional area in 2 m and Ku is the window utilisation factor. Assuming the winding allocation is such that NpAwp ¼ NsAws, the window area is
2N A = p wp Wa (5) Ku
where the window utilisation factor is defined as
2N A = ACu = p wp Ku (6) Wa Wa
where ACu is the total copper cross-sectional area in the window in m2. The maximum energy stored in the = 2 magnetic field of the transformer is Wm (1/2)LmILm max, where Lm is the magnetising inductance on the primary side and ILm is the current through the magnetising Fig. 1 Flyback converter circuit and two-winding flyback inductance Lm. Assuming that Lp ≃ Lm and Ipmax ≃ transformer current waveforms for CCM operation ILm(max), and using (2), (3) and (5), the area product of the a Flyback PWM DC–DC converter circuit core is defined as b Primary winding current ip c Secondary winding current is 2L I 2 4W A = W A = p p max = m (7) p a c K J B K J B current waveforms, respectively. The procedure to design a u m pk u m pk two-winding transformer for a flyback PWM DC–DC converter operated in CCM is given in Section 4. In Section In (7), Bpk , Bs, where Bs is the core saturation flux density at 5, the performance evaluation of the designed flyback the highest operating temperature. The area product Ap is a transformer is given. Experimental and simulation results measure of energy handling capability of a core and can be are presented in Section 6. Finally, conclusions and used to select a suitable core using the manufacturer’s discussions are given in Section 7. datasheets.
2 Core area product for a two-winding 3 Copper power loss for flyback transformer flyback transformer conducting periodic non-sinusoidal current in CCM The instantaneous magnetic flux linkage of the primary winding is given by In transformers used in isolated PWM DC–DC power converters in CCM, the current waveforms through the l(t) = N f(t) = N A B(t) = L i (t) (1) primary and secondary windings are pulsating. Therefore p p c p p the amplitudes of the harmonics are high, resulting in high winding losses. This situation is different in non- where f(t) is the instantaneous magnetic flux linking Np turns isolated PWM converters in CCM because the current of the primary winding of the transformer, Ac is the core 2 waveform of the inductor is nearly constant and the cross-sectional area in m , Lp is the inductance of the amplitudes of harmonics are low as compared to the DC primary winding, ip(t) is the instantaneous current in the component [38]. primary winding and B(t) is the instantaneous magnetic flux density in T. Equation (1) for the peak values becomes 3.1 Primary winding power loss l = N f = N A B = L I (2) pk p pk p c pk p p max Assume that the primary winding current for CCM is a rectangular waveform with a duty ratio D as shown in The maximum current density of the primary winding wire is Fig. 1b. The primary winding current waveform is given by I ≤ = p max ≃ Ip max, for 0 , t DT Jm (3) ip (8) Awp 0, for DT , t ≤ T
364 IET Power Electron., 2011, Vol. 4, Iss. 4, pp. 363–373 & The Institution of Engineering and Technology 2011 doi: 10.1049/iet-pel.2010.0040 www.ietdl.org
The DC component of the primary winding current waveform harmonic loss factor is defined as the total primary winding for CCM is loss consisting of DC and AC power losses normalised with respect to the primary winding DC power loss given 1 T I DT by I = i dt = p max dt = DI (9) pdc T p T p max 0 0 P 1 2 = wp = + sin(npD) FRph 1 2 FRpn (15) The Fourier series of the current waveform in the primary Pwpdc n=1 npD winding is The skin depth of the winding conductor at the nth 1 sin(npD) harmonic is i = I 1 + 2 cos(nv t) (10) p pdc npD s n=1 d d = √w1 (16) wn n and the amplitudes of the fundamental component and the harmonics of the primary current waveform for CCM is = where dw1 r/pm0f is the conductor skin depth at the fundamental frequency f, which is equal to the 2Ip max 2Ipdc I = sin(npD) = sin(npD) (11) converter switching frequency fs. The winding conductor pn np npD thickness normalised with respect to the conductor skin depth is The AC resistance factor F obtained from Dowell’s one- R √ dimensional solution in Cartesian coordinates [7] is √ = h = h n = An A n (17) + dwn dw1 = Rw = + = sinh(2A) sin(2A) FR Fs Fp A − Rwdc cosh(2A) cos(2A) where the normalised winding conductor thickness at (12) 2(N 2 − 1) sinh(A) − sin(A) the fundamental frequency for a rectangular conductor is + l [39] 3 cosh(A) + cos(A) h w where the first term represents the skin effect winding AC A = (18) d p resistance factor FS and the second term describes the w1 proximity effect winding AC resistance factor FP. The power loss in the primary winding for CCM is for a square conductor is [39] ⎡ ⎤ 2 1 1 R I = h h P = R I 2 ⎣1 + wpn pn ⎦ A (19) wp wpdc pdc dw1 p 2 n=1 Rwpdc Ipdc ⎡ ⎤ 2 and for a round conductor is [39] 1 1 I = P ⎣1 + F pn ⎦ wpdc 2 Rpn I n=1 pdc p 3/4 d d A = (20) 1 2 4 dw1 p = + sin(npD) Pwpdc 1 2 FRpn n=1 npD In (18), (19) and (20), w/p, h/p and d/p are called the layer = porosity factors, w is the width of the rectangular bare PwpdcRwpn (13) conductor, h is the thickness of the square and rectangular bare conductors, d is the diameter of the round bare where the normalised primary winding resistance at the conductor and p is the distance between the centres of fundamental frequency and harmonic frequencies obtained two adjacent conductors, called the winding pitch. from (12) is 3.2 Secondary winding power loss R √ = wpn = + = FRpn FSn FPn A n Assume that the secondary winding current for CCM is a Rwpdc √ √ rectangular waveform as shown in Fig. 1c. The secondary + winding current is given by × sin(2A√n ) sin(2A √n) − cosh(2A n) cos(2Ap n) 0, for 0 , t ≤ DT √ √ i ≃ (21) 2 s ≤ 2(N − 1) sinh(A n) − sin(A n) Is max, for DT , t T + l √ √ (14) 3 cosh(A n) + cos(A n) The DC component of the secondary winding current waveform for CCM is The factor FRpn takes into account the skin and proximity effects in the winding. Note that FR and FRpn are different T T = 1 = Is max = − as the latter is dependent on the nth harmonic of the Isdc isdt dt (1 D)Is max (22) primary winding current waveform. The primary winding T 0 T DT
IET Power Electron., 2011, Vol. 4, Iss. 4, pp. 363–373 365 doi: 10.1049/iet-pel.2010.0040 & The Institution of Engineering and Technology 2011 www.ietdl.org
Table 1 Two-winding transformer design for flyback PWM DC–DC converter in CCM
Step Parameter Equation Value number
Converter specifications
1 input voltage VI 50 + 10 V
2 output voltage VO 24 V
3 switching frequency fs 100 kHz
4 output power POmax(POmin) 30 W (5 W) Design equations for the transformer in flyback converter 2 5 maximum load resistance = VO 115.2 V RL max PO min 6 minimum DC voltage transfer function = VO 0.4 MVDC min VI max 7 maximum DC voltage transfer function = Vo 0.6 MVDC max VI min ∗ 8 transformer primary-to-secondary turns ratio = nconvDmax 1.5 nT − (1 Dmax)MVDC max 9 minimum duty cycle = nTMVDC min 0.4 Dmin + nTMVDC min nconv 2 − 2 10 minimum magnetising inductance for CCM ≃ = nTRL max(1 Dmin) 466 mH Lp Lm( min ) 2f1 pick 500 mH L 11 secondary winding inductance = p 222 mH Ls 2 nT − 12 maximum peak-to-peak primary winding current = = nTVo(1 Dmin) 0.43 A DiLp( max ) DiLm( max ) fsLp
13 maximum DC input current = MVDC maxIO max 0.83 A II max hconv − 14 maximum peak primary winding current = IO max + nTVO(1 Dmax) 1.84 A Ip max − nT(1 Dmax) 2fsLD Area product
15 window utilisation factor Ku 0.3
16 peak flux density Bpk 0.31 T 2 17 maximum current density of the winding wire Jm 6A/mm 18 maximum energy stored in the transformer magnetic = 1 2 0.85 mJ Wm LpIp max field 2 4 19 core area product (calculated) = = 4Wm 0.61 cm Ap WaAc KuJmBpk Core selection 20 selected core magnetics ferrite p-type round-slab (RS) core 0P-43019 4 21 area product of the selected core Ap 0.63 cm 2 22 core cross-sectional area Ac 1.23 cm 2 23 core window area Wa 0.51 cm
24 mean magnetic path length lc 4.56 cm
25 mean length of single turn lT 6.05 cm 3 26 core volume Vc 5.61 cm
27 core permeability mrc 2500 + 25%
28 core saturation flux density Bs 0.49 T = a b 29 coefficients in core power loss density equation (for p- Pv k(fs in kHz) (10Bm) type ferrite material and 100 kHz ≤ fs , 500 kHz) 29(a) k 0.0434 29(b) a 1.63 29(c) b 2.62 Wire selection r (208C) d = cu 30 skin depth of copper wire (at 208C and fs ¼ 100 kHz) w1 0.209 mm pfsm0
31 selected copper wire AWG No. 26
32 bare wire diameter of the strand wire dis 0.40 mm
33 insulated wire diameter of the strand wire dos 0.45 mm 2 34 bare wire cross-sectional area of the strand Awsi 0.12 mm 2 35 insulated wire cross-sectional area of the strand Awso 0.16 mm
36 DC resistance of the strand per unit length Rwdc/lw 0.1345 V/m Primary winding design I 2 37 cross-sectional area of the primary winding wire = pmax 0.30 mm Awp Jm Continued
366 IET Power Electron., 2011, Vol. 4, Iss. 4, pp. 363–373 & The Institution of Engineering and Technology 2011 doi: 10.1049/iet-pel.2010.0040 www.ietdl.org
Table 1 Continued
Step Parameter Equation Value number A 38 number of strands in the primary winding = wp 2.38 sp Awsi pick 2
39 number of turns of the primary winding = Ku(Wa/2) 29.82 Np SpAwsi pick 30
= 40 total length of the primary winding wire Iwp NplT 181.5 cm pick 187 cm Primary winding losses 41 primary strand DC and low-frequency resistance = Rwdc 251.5 mV Rwpdcs lwp lw R 42 primary winding wire DC and low-frequency resistance = wpdcs 125.8 mV Rwpdc Sp = 2 43 primary winding DC and low-frequency power loss Pwpdc RwpdcIs max 0.087 W 44 primary winding harmonic loss factor FRph (Use (15) @ VImin, POmax) 5.195 = 45 primary winding total power loss Pwp FRphPwpdc 0.453 W Secondary winding design = 46 maximum peak secondary winding current Is max nTIp max 2.81 A 2 47 cross-sectional area of the secondary winding wire = Is max 0.46 mm Aws Jm 48 number of strands in the secondary winding = Aws 3.63 Ss Awsi pick 4
N 49 number of turns of the secondary winding = p 19.88 Ns nT pick 20
= 50 total length of the secondary winding wire Iws NslT 121 cm pick 130 cm Secondary winding losses 51 secondary strand DC and low-frequency resistance = Rwdc 0.174 V Rwsdcs lws lw 52 secondary winding wire DC and low-frequency = Rwsdcs 43.71 mV Rwsdc resistance Ss = 2 53 secondary winding DC and low-frequency power loss Pwsdc RwsdcIO max 0.068 W 54 secondary winding harmonic loss factor FRsh (Use (26) @ VImin, POmax) 5.19 = 55 secondary winding total power loss Pws FRshPwsdc 0.354 W Total winding losses = + 56 total DC and low-frequency power loss in both Pwdc Pwpdc Pwsdc 0.155 W windings = + 57 total winding power loss Pw Pwp Pws 0.808 W Core losses m A N2 58 air-gap length = 0 c p − lc 0.26 mm lg Lp mrc
59 maximum peak magnetic flux density m0NpIp max 0.25 T B = pk l + (l /m ) g c rc m N Di 60 magnetic flux density AC component = 0 p Lp 0.024 T Bm + lg (lc/mrc) 2 = a b 3 61 core power loss density Pv k(fs in kHz) (10Bm) 2.07 mW/cm = 62 core loss Pc VcPv 10.99 mW Total losses, temperature rise and efficiency = 63 total transformer power loss Pcw PcPw 0.819 W 2 64 total core surface area At 31.95 cm 65 surface power loss density P 0.025 W/cm2 c = cw At 66 temperature rise of the transformer 21.798C DT = 450c0.826 (N s + N )A 67 window utilisation factor (recalculated) = p p s wsi 0.35 Ku Wa P 68 transformer efficiency at full power and VImin = O 97.34% ht + PO Pcw
∗ Assume converter efficiency h ¼ 0.9 and Dmax ¼ 0.5
IET Power Electron., 2011, Vol. 4, Iss. 4, pp. 363–373 367 doi: 10.1049/iet-pel.2010.0040 & The Institution of Engineering and Technology 2011 www.ietdl.org
Table 2 List of symbols
A winding conductor thickness normalised with respect to conductor skin depth 2 Ac transformer core cross-sectional area in m
ACu total wire cross-sectional area in the core window in m2 4 Ap area product in m
Awp bare wire cross-sectional area of the primary winding in m2
Aws bare wire cross-sectional area of the secondary winding in m2
Bpk peak magnetic flux density in T
Bs saturation magnetic flux density in T D converter duty cycle
FR winding AC resistance factor
FRph primary winding harmonic loss factor
FRsh secondary winding harmonic loss factor
ILm magnetising inductance current in A
Ip primary winding current in A
Is secondary winding current in A 2 Jm winding wire maximum current density in A/m
Ku window utilisation factor
Lm magnetising inductance of the transformer in H
Lp inductance of the primary winding in H
Ls inductance of the secondary winding in H
Np number of turns of the primary winding of the transformer
Ns number of turns of the secondary winding of the transformer
Pwp primary winding power loss in W
Pws secondary winding power loss in W T time period of one switching cycle of the converter in s 2 Wa core window area in m
Wm maximum energy stored in the magnetic field of the transformer in J dw skin depth of the winding conductor in mm l magnetic flux linkage in Wb-turns w magnetic flux in Wb
The Fourier series of the current waveform in the secondary winding is