Harmonics Study
Harmonics Study
Harmonics Study
POWER SYSTEM
HARMONICS
CIRILO C. CALIBJO, Ph.D., PEE, ACPE, ASEAN Engr.
2019 – IIEE National President
2018 – IIEE Vice President for Internal Affairs
2017 – IIEE National Secretary
2016 – IIEE National Treasurer
2012 - IIEE-WV Regional Governor
Dean, CEA, Angeles University Foundation
OUTLINE
• Introduction
• Theory and Computation Techniques
• Quantifying Harmonics
• Sources
• Effects and Symptoms
• Standards
• Solutions
Introduction
Utility Power to a Facility
of a fundamental 80
-0.5
50
Actual
frequency -1
-1
60
Desired
0 0.5 1 1.5
Volts
Current
• Harmonic sine and 0
Transient
Non-periodicDCSteady
Input State:aResponse
from Interharmonic
Diode RectifierComponent
cosine functions sum up
40
%%
11
to a periodic (non-
1
20
0.5
-50
sinusoidal) function 0.5
0
0 Specified
function;
00 0.5
0.5 11 1.5
1.5
0 0.5 Time 1 1.5
Time
Time
Definition: Harmonic Analysis
What is it?
– Principles, properties and methods for expressing
periodic functions as sum of (harmonic) sine and cosine
terms:
• Fourier Series
• Fourier Transform
• Discrete Fourier Transform
Definition: Harmonic Analysis
Where is it used?
– Obtain the response of a system to arbitrary periodic
inputs; quantify/assess harmonic effects at each
frequency
– Framework for describing the quality of the system input
and output signals (spectrum)
Definition: Generalization-Fourier Series
• Rectangular Form
1 ∞ 2𝜋
𝑓 𝑡 = 𝐴
2 0
+ 𝑛=1 𝐴𝑛 cos 𝑛𝜔𝑡 + 𝐵𝑛 sin(𝑛𝜔𝑡) , 𝜔 = 𝑇
𝑡𝑜 +𝑇 𝑡𝑜 +𝑇
2 2
where 𝐴0 = 𝑓 𝑡 𝑑𝑡 𝐴𝑛 = 𝑓 𝑡 cos(𝑛𝜔𝑡)𝑑𝑡
𝑇 𝑡𝑜 𝑇 𝑡𝑜
𝑡𝑜 +𝑇
is the zero frequency 2
𝐵𝑛 = 𝑓 𝑡 sin(𝑛𝜔𝑡)𝑑𝑡
or average value of 𝑇 𝑡𝑜
f(t).
Definition: Generalization-Fourier Series
• Polar Form
1 ∞
𝑓 𝑡 = 𝐴 + 𝑛=1 𝐶𝑛 sin(𝑛𝜔𝑡 + 𝜃𝑛 )
2 0
where
𝐶𝑛 = 𝐴2𝑛 + 𝐵𝑛2
𝐴𝑛
tan 𝜃𝑛 =
𝐵𝑛
Example: Square Wave Pulse
0.8
0.6
0.4
0.2
−1 −0.50 ≤ 𝑡 ≤ −0.25
0
𝑓 𝑡 = 1 −0.25 ≤ 𝑡 ≤ 0.25
-0.2
−1 0.25 ≤ 𝑡 ≤ 0.50
-0.4
-0.6
𝑃𝑒𝑟𝑖𝑜𝑑, 𝑇 = 1
-0.8
-1 2𝜋
-0.5 -0.4 -0.3 -0.2 -0.1 0 0.1 0.2 0.3 0.4 0.5 ω= 2𝜋𝑓 = = 2𝜋
𝑇
Example: Square Wave Pulse
−0.25 0.25
−1 cos 2𝑛𝜋𝑡 𝑑𝑡 + 1 cos 2𝑛𝜋𝑡 𝑑𝑡
−0.5 −0.25
𝐴𝑛 = 2 0.5
+ −1 cos 2𝑛𝜋𝑡 𝑑𝑡
0.25
Example: Square Wave Pulse
−1 1
sin −0.5𝑛𝜋 − sin −𝑛𝜋 + sin 0.5𝑛𝜋 − sin −0.5𝑛𝜋
𝐴𝑛 = 2 𝑛𝜋 2𝑛𝜋
−1
+ sin 𝑛𝜋 − sin 0.5𝑛𝜋
2𝑛𝜋
−1 1
sin −0.5𝑛𝜋 − sin −𝑛𝜋 + sin 0.5𝑛𝜋 − sin −0.5𝑛𝜋
𝐴𝑛 = 2 2𝑛𝜋 2𝑛𝜋
−1
+ sin 𝑛𝜋 − sin 0.5𝑛𝜋
2𝑛𝜋
Example: Square Wave Pulse
−1 1
−sin 0.5𝑛𝜋 + sin 𝑛𝜋 + sin 0.5𝑛𝜋 + sin 0.5𝑛𝜋
𝐴𝑛 = 2 2𝑛𝜋 2𝑛𝜋
−1
+ sin 𝑛𝜋 − sin 0.5𝑛𝜋
2𝑛𝜋
1 sin 0.5𝑛𝜋 − sin 𝑛𝜋 + sin 0.5𝑛𝜋 + sin 0.5𝑛𝜋
𝐴𝑛 =
𝑛𝜋 − sin 𝑛𝜋 + sin 0.5𝑛𝜋
1
𝐴𝑛 = 4 sin 0.5𝑛𝜋 − 2sin 𝑛𝜋
𝑛𝜋
𝟒
𝑨𝒏 = 𝒔𝒊𝒏 𝟎. 𝟓𝒏𝝅
𝒏𝝅
Example: Square Wave Pulse
4
𝐴𝑛 = sin 0.5𝑛𝜋 n A B
𝑛𝜋
4 4 1 4/p 0
𝑛 = 1: 𝐴𝑛 = sin 0.5𝜋 =
𝜋 𝜋 2 0 0
4 3 -4/3p 0
𝑛 = 2: 𝐴𝑛 = sin 𝜋 = 0
2𝜋
4 4
4 0 0
𝑛 = 3: 𝐴𝑛 = sin 1.5𝜋 = −
3𝜋 3𝜋 5 4/5p 0
4 6 0 0
𝑛 = 4: 𝐴𝑛 = sin 2𝜋 = 0
4𝜋
n ±4/np 0
𝑡𝑜 +𝑇
2
𝐵𝑛 = 𝑓 𝑡 sin 𝑛𝜔𝑡 𝑑𝑡 = 0
𝑇 𝑡𝑜
Example: Waveform Recovery
1
0.5 n=1-3
n=1-9
n=1-7
n=1-5
n=1
0
-0.5
-1
-0.5 -0.4 -0.3 -0.2 -0.1 0 0.1 0.2 0.3 0.4 0.5
1
0.5
E 0 as N
0
-0.5
-1
-0.5 -0.4 -0.3 -0.2 -0.1 0 0.1 0.2 0.3 0.4 0.5
Example: Approximation of the Rectified sine
v(t)
Vm
v (t ) Vm sin(t )
2p
T
vo(t)
Vm
T/2 T t
d=T/2
Example: Average Value
vo(t)
Vm
T/2 T t
d=T/2
d T /2
1 2
A0 Vm sin( t ) dt V m sin( t ) dt
d 0 T 0
p
Vm Vm 2
sin( t ) dt p cos(t )
t p
Vm
p 0
t 0
p
Example: Harmonic Terms
vo(t)
Vm
T/2 T t
d=T/2
2p
d
2
Ak Vm sin( t ) cos(k t ) dt
d 0 d
T /2
4Vm
T sin( t ) cos(2kt ) dt
0
p
Vm 4Vm
p sin( 2k 1)t sin( 2k 1)t dt
0
p (2k 1)( 2k 1)
Example: Summary
v(t) 2 4Vm 1
Vm
v (t ) Vm sin(t )
vo (t )
p
Vm
p k 1, 2, (2k 1)( 2k 1)
cos( 2kt )
2p
T
t ( n 2k )
2 4Vm 1
vo(t)
vo (t )
p
Vm
p n 2, 4, (n 1)( n 1)
cos( nt )
Vm
T/2 T t
Example: Numerical Problem-DC Drive
+ E=150 V
is Ra=1 W
La=5 mH
vo
s 2 240 sin(2p 60 t ), V
-
Example: Procedure to obtain response
+ E=150 V
is Ra=1 W
La=5 mH
vo
Input Voltage
200
v,V
0
s
s 2 240 sin(2p 60 t ), V -200
-
0 0.005 0.01 0.015 0.02 0.025
Output Voltage
Input Voltage
300
200
200
v,V
v,V
0
s
100
-200
0
0 0.005 0.01 0.015 0.02 0.025
0 0.005 0.01 0.015 0.02 0.025 Time, s
Output Voltage
300
200
v,V
o
Example: Input Harmonic Approximation
Input Voltage
200
s 2 240 sin( 2p 60 t ) 339.4 sin( 377 t ) Average or dc component
v,V
2
s
-200
Vo ,dc
2
Vm 216.1 V Vo ,dc Vm 216.1 V
0 0.005 0.01
p 0.015 0.02 0.025 p
Output Voltage
300
Harmonic Expansion
200
v,V
o
4Vm 1
100
o Vo,dc cos(2p 60n t )
0
0 0.005 0.01 0.015 0.02 0.025
p n2,4, (n 1) (n 1)
Time, s
o 216.1 (144.1 cos(2p 120 t ) 28.8 cos(2p 240 t ) 12.3 cos(2p 360 t ))
Example: Equivalent Circuit
io
+ E=150 V
is Ra=1 W
La=5 mH
vo
s 2 240 sin(2p 60 t ), V
-
o 216.1 (144.1 cos(2p 120 t ) 28.8 cos(2p 240 t ) 12.3 cos(2p 360 t ))
Example: Equivalent Circuit
io
Vo ,dc 216.1 V
+ La= 0.005 H
o , 2 144.1 cos( 2p 120 t ) V
-
+
o, 4 28.8 cos( 2p 240 t ) V Ra=1 W
-
+ E=150 V
o,6 12.3 cos( 2p 360 t ) V
-
o 216.1 (144.1 cos(2p 120 t ) 28.8 cos(2p 240 t ) 12.3 cos(2p 360 t ))
Example: Superimpose Sources: DC Source
io
Vo ,dc E
Vo ,dc 216.1 V
I o ,dc 66.1 A
+
Ra
La= 0.005 H
o , 2 144.1 cos( 2p 120 t ) V
-
+
o, 4 28.8 cos( 2p 240 t ) V Ra=1 W Ra=1 W
-
Vo ,dc 216.1 V
+ E=150 V
o,6 12.3 cos( 2p 360 t ) V
- E=150 V
Example: Superposition: n=2, f=120 Hz
Vo ,dc 216.1 V
Vo , 2 144.1
I o, 2 p 75 36.9(p 75) A
+ Za 3.9 La= 0.005 H
o , 2 144.1 cos( 2p 120 t ) V
-
+
iX a i (2p 120) La
+ o, 4 28.8 cos( 2p 240 t ) V i3.77 W Ra=1 W Z a (2) Ra iX a
Vo , 2 144.1p V
-
1 i3.77 3.9 75 W
-+ Ra=1 W E=150 V
o,6 12.3 cos( 2p 360 t ) V
-
Example: Superposition: n=4, f=240 Hz
io
Vo ,dc 216.1 V
Vo, 4 28.8
+ I o, 4 p 82L.a4=0.005
3.78
H (p 82.4) A
o , 2 144.1 cos( 2p 120 t ) V Za 7.61
-
+ iX a i (2p 240) La
o, 4 28.8 cos( 2p 240 t ) V
+ i 7.54 W Ra=1 W Z a (4) Ra iX a
- V 28.8p V
o,4 1 i 7.54 7.6182.4 W
+- Ra=1 W E=150 V
o,6 12.3 cos( 2p 360 t ) V
-
Example: Superposition: n=6, f=360 Hz
io
Vo ,dc 216.1 V
Vo , 6 12.3
I o,6 p 85 1.08(p 85) A
+ Za 11.35
La= 0.005 H
o , 2 144.1 cos( 2p 120 t ) V
- iX a i (2p 360) La
++ i11.31 W Z a (4) Ra iX a
o, 4 28.8 cos( 2p 240 t ) V Ra=1 W 1 i11.31 11.3585 W
- Vo , 6 12.3p V
- Ra=1 W
+ E=150 V
o,6 12.3 cos( 2p 360 t ) V
-
Example: Summary
200 vo
o
200
100
V , V and I , A
100 0 150
-
v , V and i , A
oh
0 Output Current
100
oh
100
s
-100
i ,A
-20050 50
o
-300
0
0
00 0.005
0.005 0
0.01
0.01 0.015
2 0.015 0.02
4
0.02 0.025
0.025 6
Time,
Time, s
sHarmonic Number
Quantifying
Harmonics
Root Mean Square
𝑖 𝑡 = 𝐼𝑛 sin(𝑛𝜔𝑡 + 𝜃𝑛 )
𝑛=1
then
2 2 2 2 1 ∞ 2
𝐼𝑟𝑚𝑠 = 𝐼1,𝑟𝑚𝑠 + 𝐼2,𝑟𝑚𝑠 +𝐼3,𝑟𝑚𝑠 +⋯ =2 𝑛=1 𝐼𝑛
Individual Harmonic Distortion (IHD) Level
𝐼𝑛
𝐼𝐻𝐷𝑛 = 𝑥100%
𝐼1
Total Harmonic Distortion
• Since
∞
2 2
𝐼𝑟𝑚𝑠 = 𝐼𝑛,𝑟𝑚𝑠
𝑛=1
2
∞
𝐼
𝑛=2 𝑛 𝐼2 2 + 𝐼3 2 + 𝐼4 2 + ⋯
𝑇𝐷𝐷 = =
𝐼𝐿 𝐼𝐿
Average Power
∞ ∞
𝑉𝑛 𝐼𝑛
𝑃𝑎𝑣𝑔 = cos(∅𝑛 − 𝜃𝑛 ) = 𝑉𝑛,𝑟𝑚𝑠 ∗ 𝐼𝑛,𝑟𝑚𝑠 ∗ 𝑑𝑝𝑓𝑛
2
𝑛=1 𝑛=1
Current Maximum
THD pftrue
20% 0.98
50% 0.89
100% 0.71
K Factor
∞ 2 2
𝑛=1 𝑛 𝐼𝑛
𝐾= ∞ 2
𝑛=1 𝐼𝑛
Example:
If
∞
𝑖 𝑡 = 𝐼𝑛 sin(𝑛𝜔𝑡 + 𝜃𝑛 )
𝑛=1
Phase Shift
When delayed by ΔT
∞
𝑖 𝑡 − ∆𝑇 = 𝐼𝑛 sin(𝑛𝜔(𝑡 − ∆𝑇) + 𝜃𝑛 )
𝑛=1
𝑖 𝑡 − ∆𝑇 = 𝐼𝑛 sin(𝑛𝜔𝑡 + 𝜃𝑛 − 𝑛𝜔∆𝑇)
𝑛=1
𝑖 𝑡 − ∆𝑇 = 𝐼𝑛 sin(𝑛𝜔𝑡 + 𝜃𝑛 − 𝑛𝜃1 )
𝑛=1
Phase Shift
𝑖𝑎 𝑡 = 𝐼𝑛 sin(𝑛𝜔𝑡 + 𝜃𝑛 )
𝑛=1
∞
2𝜋
𝑖𝑏 𝑡 = 𝐼𝑛 sin(𝑛𝜔𝑡 + 𝜃𝑛 − 𝑛 )
3
𝑛=1
∞
2𝜋
𝑖𝑐 𝑡 = 𝐼𝑛 sin(𝑛𝜔𝑡 + 𝜃𝑛 + 𝑛 )
3
𝑛=1
Phase Sequence
2𝜋 4𝜋
𝑖𝑏 𝑡 = 𝐼1 sin 𝜔𝑡 + 𝜃1 − + 𝐼2 sin 2𝜔𝑡 + 𝜃2 −
3 3
6𝜋
+𝐼3 sin(3𝜔𝑡 + 𝜃3 − )
3
2𝜋 2𝜋
𝑖𝑏 𝑡 = 𝐼1 sin 𝜔𝑡 + 𝜃1 − + 𝐼2 sin 2𝜔𝑡 + 𝜃2 +
3 3
+𝐼3 sin(3𝜔𝑡 + 𝜃3 )
Phase Sequence
2𝜋 2𝜋
𝑖𝑐 𝑡 = 𝐼1 sin 𝜔𝑡 + 𝜃1 + + 𝐼2 sin 2𝜔𝑡 + 𝜃2 −
3 3
+𝐼3 sin(3𝜔𝑡 + 𝜃3 )
Phase Sequence
Compact Fluorescent
White LED
Classical Nonlinear Loads: Arc Furnace
• 6 Pulse VFD
Power Electronic Loads: Variable Frequency Drives
• 6 Pulse VFD
Power Electronic Loads: Voltage-Source
Converters
• For applications
below 500HP,
voltage source
converters
employing pulse-
width modulators
with turn-on/turn-
off switches on the
motor side are often
the choice for ASDs.
Power Electronic Loads: Switched-Mode Power
Supplies.
• These power supplies
are the "front-end" of
single phase 120V loads
such as PCs and home
entertainment
equipment.
Power Electronic Loads: Switched-Mode Power
Supplies.
• Typically, they have a full-wave diode rectifier
connected between the AC supply system and a
capacitor, and the capacitor serves as a low-ripple
“battery” for the DC load.
• Unfortunately, low ripple means that the AC system
charges the capacitor for only a fraction of each
half-cycle, yielding an AC waveform that is highly
peaked
Other Nonlinear Loads
• Measure of heating in
conductor is product of
squared current I and
impedance Z (I2Z)
• Harmonics increase
heating in conductors
due to
– additional load current
flowing through conductor
– current redistribution inside
conductor (skin effect)
Overheating Distribution Transformers
𝐼𝑟𝑚𝑠 = 1 + 𝑇𝐻𝐷𝐼2
• Increased maintenance
– Excessive heat burdens electrical infrastructure, from
transformers, cables, bussing, to across the line motors.
2 R sys
kVLL X
R sys cos tan 1 , R sys (33)
MVA sc R 2
α
2
kVLL 1 X Xsys
Xsys sin tan , Xsys (34)
MVA sc R α 2
α = The turns ratio of the transformer at PCC
2
1000 kVLL (35)
R tr R pu
kVA tr
2
1000 kVLL (36)
X tr X pu
kVA tr
Example
The impedance Zin looking into the system from the load, consists
of the parallel combination of source impedance R tot jX totand the
capacitor impedance
Zin
R tot jωL tot j / ωC
1 (43)
R tot jωL tot j
ωC
1 1 (44)
ωo L tot , fo
ωo C 2πωo
The equation for Zincan be used to determine the equivalent system
impedance for different frequencies. The harmonic producing loads
can resonate (parallel resonance), the above equivalent circuit.
Designating the parallel resonant frequency by ωo(rad/sec) or f o(HZ)
and equating the inductive and capacitive reactances.
Example
- Harmonic current components that are close to the parallel
resonant frequency are amplified.
- Higher order harmonic currents at the PCC are reduced because
the capacitors are low impedance at these frequencies.
- The figure below shows the effect of adding capacitors on the 400
Volts bus for power factor correction.
Example
• This figure shows that by adding some typical sizes of power
factor correction capacitors will result in the magnification of the
5th and 7th harmonic components, which in turns makes it even
more difficult to meet the IEEE 519-1992 harmonic current
standards .
• Power factor correction capacitors should not be used without
turning reactors in case the adjustable speed drives are >10% of
the plant load.
Example
R sys
202
250
cos tan 12.4 0.6154Ω
Xsys
202
250
sin tan 12.4 1.4769Ω
20
α 50
0.4
Example
1000 0.42
R tr 0.01 0.00160Ω
1000
1000 0.42
X tr 0.07 0.0112Ω
1000
R tot = 0.000246 + 0.0016 = 0.001846 Ω
X tot 0.000591 0.0112 0.011791Ω
0.011791
L tot 37.55 10 6 H
2 π 50
Example
Case a:
1000 0.42
Xc 0.8 Ω
200
1
C 3.98 103 F
2 π 50 0.8
1
fo 412.18HZ
2 π 37.50 10 6 3.98 103
For 200 KVAR, the harmonic order at which parallel resonance occurs is:
h 412.18 50 8.24
Example
Case b: Case c:
1000 0.42 1000 0.42
Xc 0.4 Ω Xc 0.267 Ω
400 600
C 7.96 103 F C 11.94 103 F
f o 291.45HZ f o 237.97HZ
h 5.83 h 4.76
Example
Case d:
1000 0.42
Xc 0.2 Ω
800
C 15.92 103 F
f o 206.08HZ
h 4.12
It is clear for the above system that in the 600 KVAR case,
there exists a parallel resonant frequency f o close to the 5th
harmonic.
Solutions
Harmonic Solutions