Power Electronic Control Device Alzatari - Diss
Power Electronic Control Device Alzatari - Diss
Power Electronic Control Device Alzatari - Diss
genehmigte Dissertation
Von
Mazen Alzatari
aus
Hebron, Palästina
In AC high voltage tests such as partial discharge (PD) measurements and di-
electric test, the test’s voltage source quality can influence the test results accuracy.
Therefore, the international standards for high voltage testing techniques [1,2] deter-
mine a limit of distortion value of a sinusoidal voltage wave shape. √ The ratio of the
voltage peak to it’s root mean square value (RMS) must be within 2 ± 5%, and a
total harmonic distortion (THD) less than 5% as well. Practically due to the pres-
ence of harmonics in power network and using a step up transformer for AC tests, the
distortion of an AC sinusoidal signal could exceed the standard limit.
i
Abstract
Electrical insulators reliability and age degradation estimation can be evaluated
accurately by test. Beside other test procedures high voltage dielectric test(HV-DI)
and partial discharge (PD) measurements are used for this purpose. Since the accuracy
in test is an essential issue, HV-DI test and PD measurements have to be accomplished
in nearly ideal signal quality of the test’s voltage source, as well as controllable dis-
tortions(in presence of harmonics)in order to investigate all possible circumstances.
Performing the test using other frequencies(16,7 Hz, 50 Hz, 60 Hz) are also on market
demand.
For the dimensioning of the power electronic system, the essential components
of a HV test system need to be characterised. The characterisation focus on the
properties relevant for the use cases High Voltage Dielectric Tests and Partial Discharge
Measurements. Based on this characterisation the required parameters for a single
phase DC-AC-Inverter and the corresponding PWM parameters are derived.
The goal of spectral purity disallows the use of standard PWM schemes. Based
on the Naturally Sampled PWM scheme as described in the literature [10], a new
PWM scheme (Enhanced Unipolar PWM) is developed which allows its use for HV
testing. The new driven algorithm is implemented in Digital Signal Processor (DSP)
unit which is controlling an H-Bridge and an attached 100 kV, 5 kVA test transformer.
ii
For all those who have respect and love to the prophets and messengers, especially for
Noah, Abraham ,Moses, Issa(Jesus) and Mohammed Peace be upon all of them
iii
Acknowledgements
“ Be grateful to Me and to both your parents; to Me is the eventual coming”(
An approximate meaning,the Qur’an-Luqman 14). I am so grateful to the GOD Allah,
for giving me the strength to do this work until the end.
This thesis was written while I was in the Institute of Electrical Power Transmis-
sion (ETS), in the Faculty of Engineering Sciences at the University of Duisburg Essen
in Germany. I would like to express my sincere gratitude to Prof. Dr.-Ing. Holger
Hirsch, the head of the institute for his wise guidance,assistance and golden advices
which led to this successful work. I would like to thank Prof. Dr.-Ing. Frank Jenau
(second supervisor) for his efforts and interest in this work research.
Many thanks to all of the ETS team for creating a friendly research atmosphere,
the scientific discussions and helping in lab experiments, especially for Dr.-Ing. Joerg
Honerla, Ir. Budi Sudiarto, Ir. Aji Nur Widyanto and M.Sc. Christoph Schwing.I feel
also grateful to my colleague M.Sc. Rasha Almazedi for the help in text editing.
Mazen Alzatari
Duisburg Sep. 2015
iv
Contents
Dedication iii
Acknowledgements iv
List of Tables x
Abbreviations xi
Symbols xii
1 Introduction 1
1.1 Thesis principles and aims . . . . . . . . . . . . . . . . . . . . . . . . . 1
1.1.1 Signal quality in high voltage test . . . . . . . . . . . . . . . . . 1
1.1.2 Safety and protection in dielectric test . . . . . . . . . . . . . . 2
1.1.3 Frequency conversion . . . . . . . . . . . . . . . . . . . . . . . . 2
1.1.4 Generation of harmonics . . . . . . . . . . . . . . . . . . . . . . 4
1.2 Requirements of a power electronic control device for high voltage test
transformers . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4
Bibliography 83
List of Figures
3.1 Comparison between the simulated results and the theoretical calculated
from[11 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 45
3.2 Compare the harmonics of V ohP W M spectrum in three cases: (1) simu-
lation with dt=0.(2)simulation with dt=0.3 us (3) measured data from
Pec-HV . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 48
3.3 THD evaluation by system simulation at dt=0.3µs and zero intended
harmonics . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 52
3.4 Measured THD for the three fundamental frequencies at zero intended
harmonics . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 52
3.5 Simulated THD of the three fundamental frequencies with intended har-
monics . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 54
3.6 Measured THD for the three fundamental frequencies with intended
harmonics . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 56
3.7 System dynamic behaviour at various ma . . . . . . . . . . . . . . . . 57
3.8 Highest voltage harmonics peaks of the measured V ohP W M in frequency
range 99.8 kHz to 500 kHz . . . . . . . . . . . . . . . . . . . . . . . . 60
3.9 Highest voltage harmonics peaks of the simulated (50 Hz) V ohP W M in
frequency range 99.8 kHz to 500 kHz at various ma . . . . . . . . . . . 61
x
Abbreviations
xi
Symbols
xii
Symbols xiii
Dielectric tests (DT) are used for the validation of the selection of a suitable
insulation materials, and Partial Discharge (PD) measurements allows to estimate the
ageing behaviour of the materials. In classical way in order to perform a dielectric or
PD test a step up transformer is used to generate high voltage, where the Equipment
Under Test (EUT) is connected to the transformer’s secondary winding. The step up
transformer’s primary voltage is variable, that can be done with a motor driven or
manually operated variable transformer.
Unfortunately the voltage from the grid nowadays is not a pure sinusoidal signal
due to the high number of nonlinear devices ( e.g switched mode power supplies)
connected to the grid. Nonlinear loads are the reason behind the appearance of a
low voltage signals appear as multiple frequencies of the fundamental(main sinusoidal
frequency) , these are known as Harmonics. Harmonics will adversely effect the test
accuracy, that is why an international standards[1,2] has put a limit for the Total
Harmonic Distortion (THD) for AC high voltage test. It has to√be less than 5% and
the ratio of sinusoidal peak voltage to its RMS must be within 2 ± 5%. Harmonics
are not stable in the power grid, the amount of THD will change according to the load.
The Union of the Electricity Industry - EURELECTRIC has made a study on Europe
power quality[5,6] which shows that the THD is increasing every year, more over it
may exceed the standard limit( 5%). Moreover, the step up transformers in HV tests
will amplify the harmonics present in the voltage source and they produce additional
harmonics due to the non-linearity of the iron core[7,13,14,29,46]. In Lab experiments
the distortion factor could be more than 5% especially when the applied voltage is
over 23 from the transformer’s primary rated value. Therefore, the requirements of the
standards are hard to fulfil in such cases.
In [29,30,8] devices based on power electronics are presented. These devices are
intended to do PD measurements and dielectric tests, they and many other devices
described in companies websites e.g[53] or in other literature are focusing on the output
signal quality of the device itself which is still not pure sine(THD roughly 3 − 5%) .
Furthermore those devices do not allow the compensation of the generated harmonics
by the test transformer itself, which in some cases is a real challenge to fulfil the
HV-tests standards requirements.
1
Chapter 1. Introduction 2
In dielectric tests which are used for the insulation investigations, the voltage
applied to the dielectric material is gradually increased. At a certain point a breakdown
occurs. The breakdown voltage is the amplitude of the voltage immediately before the
breakdown. After the breakdown event an electric arc is built within the breakdown
channel, which is fed by the power supplied through the test transformer. In case of
solid state materials the arc leads to larger areas of damage so that further investigation
on the break down development is nearly impossible. Therefore, the power feeding into
the test transformer needs to be interrupted within a short time.
The concept of the system protection in dielectric tests is based on a fast detec-
tion of voltage fall time and/or current rise time in the breakdown channel. When
a breakdown is detected a Solid State Circuit Breaker(SSCB) is trigged in order to
switch off the voltage source, that will limit the injected energy in the breakdown
channel.
A previous development work of an external fast switch-off units was done at ETS
institute at University of Duisburg Essen[24]. It consists of two IGBTs as switching
off elements and the break down detector which is based on the voltage fall time. The
switching off time for this device is roughly 500 µs (relatively high) and required extra
semiconductor elements for construction, i.e. using the same semiconductors (IGBTS)
of the H-bridge inverter will be an advantage.
The power frequency used worldwide is either 50 Hz or 60 Hz. For some railway
systems the power frequency of 16.7 Hz is used. Therefore according to the market
demand, performing HV-DT and PD tests with the power frequencies 50 Hz, 60 Hz
and 16.7 Hz are requested.
Using a static frequency converter as a voltage source for HV tests and PD mea-
surements is a great step forward in the field of HV tests. The test can be accomplished
with variable frequencies(16,7Hz, 50Hz, 60Hz) as well as controllable amplitudes. The
test source output signal quality still have two challenges; THD and background noise
which are adversely influence on the test accuracy [17,8,18,22,16,29,22].
In order to analyse the reasons behind these challenges, the source of harmonics
distortions and the background noise has to be comprehensively understood (from the
frequency converter side).
• DSP which includes the PWM software algorithm and controls the switching
elements through a specific driver.
To generate a unipolar PWM signal, two sinusoidal signals (desired output) have
to be compared with a triangle carrier frequency (switching frequency), the result
is a voltage pulse with variable width proportional to the amplitude of sinusoidal
signal. In the ideal case the switching has to occur simultaneously with the intersection
points between the sinusoidal and triangle signals, this is known as Naturally Sampled
PWM (NS-PWM). On the other hand, practically the sinusoidal signal is stored as
samples in the DSP memory, sampled sine signal instead of analogue continuous one
has an influence on the signal output quality, especially the THD and the filtration
of the background noise, thus it has been considered as a different switching scheme
(Regularly Sampled PWM /RS-PWM)[10], RS-PWM is the conventional method of
implementing the naturally sampled unipolar PWM [11,31,34,35], more details are
given in clause 2.5.
Different methods and researches were done to eliminate the harmonics and back-
ground noise of PWM converters such as:
• controlling the pulses timing in order to eliminate harmonics especially with low
switching frequency (<500 Hz)[15,19,20,21]
• adding shielding layers on the test transformer’s windings in order to reduce the
common mode electromagnetic interference and thus the background noise in PD
measurements.[17]
In addition, these methods are applied for conventional frequency converters, but
in this work the new Pec-HV is able to perform signal preconditioning and to generate
a controllable harmonics.
PD measurements are made in the labs with nearly pure signals. Nevertheless,
harmonics are present in real installations and may influence the PD behaviour of the
high voltage insulation. In order to investigate this influence a sources is needed which
is able to generate intended controllable(amplitudes/phase shift) harmonics.
As shown above the solutions available in the market or described in the literature
are individually designed for specific application fields. Typically they can not be used
for another application. There is so far no holistic solution. Therefore this dissertation
aims to the realisation of a holistic approach. The contributions of this dissertation
and the goals for the device are:
improved to be faster and reliable based on voltage fall time and current raise
time.
• High quality output: In order to have a pure sinusoidal signal for HV-DI tests or
PD measurements, the signal has to be free from low order harmonics(frequencies
below 2 kHz) and acceptable filter complexity/cost for high frequency region(100
kHz to 500 kHz). That is achieved by bridging the gap between the regularly
sampled and naturally sampled PWM (Ideal case), thus eliminating the low order
harmonics and reducing the filtration complexity of the background noise. The
new Enhanced Unipolar PWM algorithm based on calculating the intersection
points between the desired and the carrier frequencies using programmed Matlab
model. According to the mentioned intersection points the PWM signal duty-
cycle values are calculated. The outcome, which is actually a duty-cycle values,
are stored in the DSP memory, multiplying the duty-cycle values with a factor
will control the output voltage amplitude, while the time of processing each
duty-cycle value will determine the desired output frequency (16,7 Hz, 50 Hz, 60
Hz).
• Generation of intended harmonics: The power electronic control device for high
voltage test transformers (Pec-HV) is able to generate intended harmonics for
lab investigations, these harmonics (3rd and 5th) are fully controllable; phase-
shift and amplitude. Intended harmonics are generated based on the Enhanced
Unipolar PWM algorithm (explained in previous paragraph), where the duty-
cycle values of the fundamental frequency and the intended harmonics are added
together(simplified concept) to generate the wanted harmonic. The challenge
is to keep the real measured THD close to the desired one, and to keep the
background noise around 1 pC. A Matlab model is used for optimizing the En-
Unipolar-PWM parameter such as the dead time and switching frequency to
achieve the desired goals. The difference between desired THD and the measured
in case of intended harmonics generation is ± 1 % and the background noise
remains below 1 pC. Furthermore the generated harmonics can be employed
eliminate the dead time problem in the unipolar single phase PWM inverter.
2. Conceptual Framework of
controlling HV test transformers
2.1 General
AC high voltage test quality is the major interesting point in this research. Beside
other issues, the test quality depends on the test’s voltage source. The purity of the
test source signal and its ability to compensate the generated distortion by the test
transformer’s iron core are the two main challenges. Since the used voltage test source
in this research is a DC-AC inverter (the developed Pec-HV), the implemented PWM
switching scheme will influence the quality of Pec-HV.
In order to achieve the goals of this research(clause 1.1) the effort is focused on the
PWM switching schemes starting from choosing the best PWM method from quality
point of view (less THD), and then filling the gaps in the selected PWM scheme to get
the optimal quality and to enhance it to generate intended controllable harmonics to
support the required Pec-HV features and functionality. Figure 2.1 shows the area of
this research which is focused on the test signal quality.
6
Chapter 2. Conceptual Framework 7
Nowadays the dielectric stress to materials can be calculated with the appropriate
computation tools. From the understanding of the dielectric behaviour of the material
and the mechanism of breakdowns, the maximum permissible electric field strength can
be estimated. However, the experience shows that breakdowns are strongly influenced
by impurities, local inhomogeneities and surface effects (roughness, dust, humidity).
Since these items are hard to implement in calculation models, each high voltage device
used in the power grid will be tested (type tests or routine tests) .
Insulators are able to withstands voltages higher than rated until a specified
limit, then due to increasing of electric field stress the insulator will collapse and a
current will flow though it. That is the basic description of insulator voltage break
down. The breakdown test can be done by applying a power frequency voltage (50
Hz, 60 Hz, 16.7 Hz) through a step up transformer(up to 100 kV in this work ). Two
electrodes are connected to the secondary winding of the test transformer, where the
tested insulator specimen is fixed, then gradually increase the electric filed(Voltage)
until a breakdown is detected, and then the applied voltage has to be switched off(see
figure 2.2)[42,43,44,1,2].
Figure 2.3 shows a simulated breakdown process of the above shown schematic.
Before the breakdown occurs, the transformer’s primary current is relatively low (few
amps), but after the breakdown where the voltage falls to nearly zero, the current
raise up immediately to a destructive values (few hundreds Amps). Detecting the
breakdown and switching off the applied voltage very fast is a significant issue in such
test.
Chapter 2. Conceptual Framework 8
Figure 2.3: Breakdown, voltage and current characteristics, the diagram shows
the voltage at the output of HV divider and the primary current
In this study, The experimental investigations were done using two spheres as
electrode, and the air as an insulating material. The new device (Pec-HV) has the
ability of interrupting the applied voltage within roughly 50 µs in case of breakdown.
The device also allows a lab dielectric breakdown tests with a pure voltage sinusoidal
signal as well as with different power frequencies and an optional generated intended
harmonics in order to investigate a real conditions where the harmonics are present in
the power grid.
The well known theories of Townsend mechanism and Streamer mechanism de-
scribe the behaviour of charged gaseous particles under the influence of electric field
and how that leads to gaseous breakdown. For simplicity, gaseous breakdown can be
explained as shown in figurer 2.4, where the three main steps are[42,43,44]:
1. Ohmic behaviour
3. The current will start to flow between the electrodes (avalanche) which leads to
the insulator breakdown.
The best definition for Partial discharge is from the International standard (IEC
60270)[45] “Partial discharge (PD) is a localized electrical discharge that only partially
bridges the insulation between conductors and which may or may not occur adjacent
to a conductor ”
PD could occur :
• Between sharp edges and ground, where high voltage is applied toward the ground
in gaseous media , which is known as corona discharge(figure 2.5).
So PD occurs due to high stress electric field on insulators where high voltage
is applied. These discharges appears as small voltage dips(pulses) in a range of few
millivolts as amplitude and a few nanoseconds as width. Thus a conventional voltage
or current measurement devices are not able to measure PD pulses. That is why a
special device is needed.[45,44,17,57,56]
Figure 2.7 shows one of the common circuits which is used to measure PDs. A
coupling capacitor is connected in series with a coupling device, where both are in
parallel to the tested object. The coupling device contains a passive element (R,L,C)
which has relative low impedance (roughly 100Ω or less). The coupling device and
the capacitor also work as high pass filter. The connecting cable between the coupling
device and the Measurement Instrument (MI) could be an optical fibre or a coax
cable.[45,44,55]
1. Integration along the pulse boundary in time domain: The PD is a narrow current
pulse (few nano seconds), where the integration of the current pulse is the charge
q which is expressed in pC. Figure 2.8(a).
Chapter 2. Conceptual Framework 11
Figure 2.8: (a) PD pulse in time domain (b) PD pulse in frequency domain
Band pass filter is needed to limit the measurements bandwidth. The interna-
tional standard IEC60270 [45] has determined the recommended values for wide band
PD measurements as the following:
• 30kHz ≤ f 1 ≤ 100kHz
• f 2 ≤ 500 kHz
The PD-MI works at a high frequency range as mentioned before(f 1 and f 2) and
the PD patterns are evaluated in parallel with time domain (phase angle) to identify
the type of PDs[44,57,55]. Taking in account that the typical frequency converters
synchronise the carrier with the sine signal(desired), therefore the PWM switching time
points appear at constant phase angles, thus the PD-MI has no chance to discriminate
the signal from PWM switching, neither in frequency nor in time domain.
In order to eliminate this interference the PWM output voltage has to be filtered
with a proper attenuation ratio, especially for the bandwidth (100 kHz to 500 kHz).The
PWM harmonics magnitude in the interference area as shown in figure 2.9 are measured
in milli volts. The relation between the two units mV and pC can’t have a general
formula [56], while it depends on the impedance of coupling between the tested object
and MI, but as it is clear from [56] and [61] in addition to lab work experience that the
range of PD patterns can be expressed in volt starting from roughly 1 millivolt and
less. According to experimental investigations, it is found that 1 pC would correspond
to roughly 35 to 5 µV (magnitude of the PWM harmonics in the interference area).
That means, the key of a proper filter design is to attenuate the PWM spectrum
in the range (100 kHz to 500 kHz) to be below roughly 5 µV. Figure 2.10 shows a
real PD measurements performed in the ETS lab in case of improper filtering of PWM
converter, where the PWM switching side-band frequency interference with PD pattern
is obviously clear.
2.4.1 Harmonics
The total harmonic distortion factor is measured by calculating the square root
of summation squares of harmonic amplitudes divided by fundamental frequency am-
plitude. as shown in equation 2.1 below.
p
V32 + V52 + V72 + .....V492
T otal harmoinc distortion T HD = (2.1)
V1
change the load impedance( increasing in case of inductive load, or decreasing the in
case of capacitive load). Thus harmonics have different effects depending on the device
which is exposed to them, but in general harmonics (if not eliminated) can over heat
the devices[58]. They must be mitigated, otherwise they may cause damage to the
electrical system.
In AC high voltage tests such as dielectric Test and PD measurements, the inter-
national standard requirements (IEC 60060-1 2010) for the test-source voltage signal
is to be sinusoidal, √
a small deviation from sinusoidal is accepted if ratio of peak to
r.m.s. values equals 2 within ± 5 % , and the THD has to be taken in account too.
The IEEE standard for High-Voltage Testing Techniques[2] put a limit for the THD
to be less than 5 %.
In practice these requirements are difficult to fulfil in some cases and the stan-
dards limits values can be exceeded. That is due to the following :
• Harmonics are presence in the voltage test source whether it is from mains 50
Hz or conventional frequency converters.
• Test transformers are nonlinear devices, they will increase the THD while they
produces extra harmonics in addition to amplifying the voltage source harmon-
ics.(clause 2.4.2 presents the details)
The flux φ will flow through the iron core. According to the Faraday’s law the
induced voltage will depend on the flux changing rate dφ
dt
and number of the coil turns
(N).
dφ dφ
up = .N p and f or the secondary usec = .N sec (2.2)
dt dt
Figure 2.14 explains the nonlinear behaviour of an iron core. Assuming a sinu-
soidal voltage u is applied to a coil on an iron core transformer, the flux φ is 90 degree
shifted. The flux φ will start to increase in approximately linear behaviour in time
intervals (t0 to t1) and (t1 to t2) as well as the current i. When the core is nearly
saturated (t2 to t3), the iron permeability µ is no longer linear as well as the current
and the magmatic field H. The induced i is non sinusoidal, it contains harmonics.
These current harmonics will goes back to the voltage source, the amount of voltage
distortion will depend on the source impedance. By applying Fourier series analysis,
only odd harmonics will be seen in voltage, since the distorted current is an symmetry
odd function.
In high voltage tests where a step up transformer is used to generate a high volt-
age, the tested objects or measurement instruments are almost have a high impedance,
thus the secondary current is normally less than 100 mA. The main current loop is the
primary current(few Amps or few 10s Amps) as shown in figure 2.15 . This current Ip
is a nonlinar due to the nonlinarity of the iron core and contains odd current harmonics
which will appears as a voltage drop throw the source impedance, thus Vp=Vs+Vz
will be induced to the secondary windings and will be amplified with the contained
Chapter 2. Conceptual Framework 16
harmonics. This interprets the reasons for a distorted sinusoidal voltage signal in figure
2.12 (clause 2.4.1).
Figure 2.15: Nonlinearity-core’s effect on the current loop in the Test Transformer
A measurement in lab was setup in a similar way of figure 2.15, were the voltage
source is the Pec-HV and a harmonics network analyser is connected to measure the
harmonics at Vp. This experiment shows that the 3rd and the 5th harmonics are
the dominants. Figure 2.16 show the result of this experiments where the core non-
linearity is clear since the THD is 7.1 %, which is more than 5 % in this case, thus a
solution has to be found for this problem.
Figure 2.17 shows an examples of high voltage test transformers, where insulation
issue is more significant and they should have less flux density within the core to avoid
high order harmonics .[47,44].
Chapter 2. Conceptual Framework 17
Figure 2.16: Measured 3rd and 5th harmonics which are generated due to core
non-linearity
The interesting point for this thesis is the inverter’s switching algorithms, espe-
cially how it can be modified and implemented in order to get a pure clean signal or
combined with intended harmonics output which can be efficiently used as a voltage
source for HV tests.
The well known inverter topologies are half bridge and full bridge inverter as
shown in figure 2.19. Since the full bridge topology is more efficient for real inverters
design, it will be adopted in this thesis[10,11,12].
As seen in figure 2.18, the main switching element is an IGBT, it is efficient reli-
able semiconductor and can be driven with a simple driver circuit. Full bridge IGBTs
are available now in advanced complete design including the drivers and a protection
circuits, they are well known as power modules (PM).A digital signal processor (DSP)
is used to control the IGBTs switching process through an implemented switching
algorithm, where a modulation of a desired signal (Vsine ) and carrier signal (Vtri ) is
executed.
The PWM switching algorithms specifies the inverter’s output signal quality and
shape. There are many switching algorithms/bridge topologies that can be applied for
a full bridge inverters design, some of them are too complex to be implemented for a
high level of power(few 10s kW) such as multi stage inverters. Practically bipolar and
unipolar PWM switching strategies are efficient, controllable and not complex to be
implemented. In this thesis a Unipolar PWM strategy is developed and enhanced to
achieve the goals. In the next parts of this clause one can find an explanation about the
unipolar concept; starting from the general PWM concept and approaches, and then
the half bridge topology with bipolar PWM switching strategy, then the full bridge
bipolar PWM, and then the full bridge unipolar PWM .
• If Vsine > Vtri then VP W M is high and when Vsine < Vtri then VP W M is low.
Chapter 2. Conceptual Framework 20
• The peak voltage of V sine ( V psine ) should be less or equal to the peak of the
carrier signal V ptri
Vp sine
ma = ≤1 (2.3)
Vp tri
where ma is the modulation index.
• mf is the frequency ratio of the carrier frequency ftri and desired frequency fsine
.
ftri
mf = (2.4)
fsine
• mf should be an even integer(in case of using Unipolar PWM scheme) in order
to get half wave symmetry of the periodic desired output voltage, therefore less
harmonics in the total output voltage (Vo) [11,12].
• For optimal system behaviour, in order to produce quality output with less har-
monic distortions, Vsine & Vtri should be synchronised(mf is an integer) and the
modulation index ma (figure 2.21) should be in the linear region[11,12].
Assuming that the VP W M signal (figure 2.20) is fed to the leg A of the full bridge
in figure 2.19(b) where VP W M is directly fed to SA and through an inverter to SA` .
Va(t) will have the same characteristics of VP W M but having the VDC value for the
high value and 0 for the low value.
Chapter 2. Conceptual Framework 21
In such a system the output voltage Va(t) is a result of the modulation process
of the two signals Vsine and Vtri , where :
• Vsine (t) = V psine sinwsine t . where, wsine is the angular frequency of the desired
output voltage.
• and Vtri (t) = 2∗Vπptri arcsin [sin(wtri t− π2 )] . where , wtri is the angular frequency
of the carrier triangle signal as shown in figure 2.20.
The well known method to analyse any signal to its frequency components is the
Fourier Series, in such system where two periodic functions ( Vsine and Vtri ) determine
the characteristics of Va(t), since there is no rotational relation between them, a single
Fourier Series is not applicable, therefore a Double Fourier Series method will be
applicable to calculate the system spectrum[10,50,51]. This method is used first by
Bennett and Black in communication modulation theory[51,10].
Equation 2.5 [52,51,10] shows the general form of Double Fourier Series which
has four main parts:
A00
• The first part 2
is a constant which represents the DC offset value.
• The second part (summation of n’s) represents the fundamental frequency and
it’s baseband harmonics (see figure 2.23) where y is the desired signal frequency
which is wsine t.
• The 3rd part (summation of m’s) represents the carrier harmonics where x is the
carrier frequency wtri t
∞ ∞
A00 X X
f (x, y) = + [A0n cosny + B0n sinny] + [Am0 cos mx + Bm0 sin mx]
2 n=1 m=1
∞
X ∞
X
+ [Amn cos(mx + ny) + Bmn sin(mx + ny)]
m=1 n=−∞,n6=0
(2.5)
where Z π Z π
1
Amn = 2 f (x, y)cos(mx + ny) dx dy (2.6)
2π −π −π
Z π Z π
1
Bmn = 2 f (x, y)sin(mx + ny) dx dy (2.7)
2π −π −π
Chapter 2. Conceptual Framework 22
or
Z π Z π
1
Amn + jBmn = 2 f (x, y) ej(mx+ny) dx dy (2.8)
2π −π −π
π
Z π Z (1+ ma sinwsine t)
1 2
Amn + jBmn = 2 V DC ej(mx+ny) dx dy (2.9)
2π −π − π2 (1+ ma sinwsine t)
As a final result of applying the previous explained method for a full bridge
inverter, the PWM spectrum of a sine-triangle function will be as in figure 2.23 [
52,51,10]. A generalised spectrum harmonics amplitude values of both Unipolar and
bipolar switching schemes are clearly presented in [11], table 2.1 contains a samples of
these values.
Chapter 2. Conceptual Framework 23
From table 2.1 and figure 2.23 the harmonics always appear as the fundamental
frequency h1 , carrier frequency and its sideband frequencies which can be described
as the following[52,51,10,11,12] :
h = imf ± j (2.10)
where:
j is even for odd i and vice versa. (i.e j=2,4,6... for i= 1,3,5... and j=1,3,5.... for i=
2,4,6...)
For example if mf =250 ,i=2 and j=1 so h = 2*250 ± 1 = 499 and 501 , if the
fundamental frequency is 50 Hz then the harmonic number 499 is at the frequency of
499*50 =24.950 kHz and has an amplitude of 0.314 V if ma=0.8 and VDC=1 V .
The task of the low pass filter in figure (2.19 b) is to filter out these high fre-
quencies and to let the fundamental and the baseband harmonicas ( called intended
in this thesis) to pass.
Chapter 2. Conceptual Framework 24
In this thesis the desired sinusoidal output voltage will contain an intended 3rd
and 5th harmonic. In such case the double integral of equation 2.9 cannot be simplified
using the conventional method (Bessel function) [52]. Furthermore in case of exact
integer carrier ratio mf , simulation studies of PWM harmonic investigations using the
well known Fast Fourier Transform(FFT) is effective[10]. From this point of view FFT
will be used in this thesis to investigate the harmonic content of the Matlab model
system which will be explained in clause 3.3 .
PWM approaches
• Naturally Sampled approach: This approach represent the ideal case(or the so-
called naturally sampled PWM), where both signals the desired sine signal and
the carrier signal are analogue. In such a case the VP W M signal switching will
occur simultaneously with the intersection points between the two signals (Vsine
and Vtri ) , see figure 2.24 (a). The analysed VP W M signal in frequency domain
will not contain any baseband harmonics[10,40], if a low pass filter is applied to
the inverter’s output in order to filter out the carrier and its sideband frequencies
, then a zero voltage THD at the filter output is expected. Unfortunately that
is not the case in real inverters design due to implementation challenges such
regular sampled PWM and the dead time effect.
Dead Time
Dead Time (dt ) means a delay time between switching the two IGBTs which are
working inversely as seen in figure 2.25, where a delay time is generated at the raising
edge of SA and SA` . This delay can be generated optionally by the DSP, the main
point is to generate a waiting time until the switching process of an IGBT is completed,
then to switch the complementary IGBT. That is important to avoid any high current
through both IGBTs during switching, the delay time must be equal or larger to the
IGBT switching on or off time (the typical value from data sheet is roughly 2 µ s).
The voltage error due to the dead time distortions is proportional to the switching
frequency ftri , the dead time dt and VDC [49]. In [49] the modulation index ma whether
it is ignored or accepted as 1, in addition to that it has been mentioned in [60,48] that
the voltage error is inversely proportional to the modulation index ma . Thus the
complete description for the dead time voltage error is the following :
ftri ∗ d t ∗ 2VDC
Dead time voltage error = (2.11)
ma
The Matlab model which will be explained in chapter three is used to illustrate
equation 2.11 by the following example :
Assuming a unipolar PWM switching scheme has the following parameters; ma=0.4
, dt=0, ftri = 800 Hz ,VDC = 300 ,Vsine & Vsine`= 50 Hz, in this case the fundamental
output voltage V op is 119.7 V . If the same system has a 2µs dead time, then the
fundamental output voltage V op is 117.4 V , that means the voltage error is 2.3V,
which is in good agreement with the results (2.4 V) if equation 2.16 is used for the
calculation.
Chapter 2. Conceptual Framework 26
Many solutions such as using extra hardware elements or applying the control
theory to design a control loop are applied to eliminate/compensate the dead time
distortions [32,48,49]. In this thesis it will be eliminated by the concept of intended
harmonic generation, more details will be shown in the next chapter.
As mentioned before half bridge with bipolar strategy inverter(refer to figure 2.19
a) is the simplest one. The aim is to get a sinusoidal output voltage with predetermined
frequency. Figure 2.26 shows how to generate the VP W M signals (SA and SA`) by the
DSP, to simplify the explanation the naturally sampled PWM approach is adopted.
The desired sine signal Vsine is compared with the carrier Vtri , as a result the signal
SA is generated, and then by inverting SA the signal SA` is generated. Feeding the
signal SA through the IGBT driver to IGBT 1 and the signal SA` to IGBT 2 will let
the output voltage V oP W M to switches between VDC/2 and -VDC/2 . Filtering the
signal V oP W M by a Low Pass Filter(LPF) will get out the desired signal Vo as seen in
figure 2.27.
In half bridge case the maximum peak voltage of the filtered Vpwm ( V op ) is half
of the VDC at ma =1. This is the main disadvantage of half bridge inverter.
V DC
V op = ma (2.12)
2
Chapter 2. Conceptual Framework 27
In this case the amplitude of the low pass filtered output voltage V op is:
V op = ma ∗ V DC (2.13)
Unipolar PWM strategy (only applicable for full bridge topology) can be achieved
by comparing the carrier frequency Vtri with two sinusoidal signals which have 180˚
phase shift between them, these two signals represent the desired output frequency.
The results from a simple LTspice simulation visualise this modulation scheme withma =0.9
and mf =8 (Vsine =50 Hz) as an example. The naturally sampled PWM approach is
adopted for simplicity (Figure 2.29).
1. Two signals of the desired output frequency are 180˚out of phase ( Vsine and
Vsine`) are compared with the triangle carrier frequency Vtri .
2. As a result of comparing Vsine with Vtri the SA signal is generated, and by
comparing Vsine`with Vtri the SB signal is generated. SA`is generated by inverting
SA and as well SB ` by inverting SB. The four signals controls the switching
process of the four IGBTs , where IGBT 1 fed by SA , IGBT 2 fed by SA` ,
TGBT 3 fed by SB and IGBT 4 by SB`.
3. The outcome of the previous switching process is V oP W M , which is actually VA-
VB, this subtraction will eliminate some of the harmonics content and doubling
the switching frequency[11,12] . The amplitude of the low pass filtered output
signal Vo is V op where:
V op = ma ∗ V DC. (2.14)
4. Here one can see the frequency domain of the output voltage VP W M using FFT .
In order to be in a good agreement with naturally sampled PWM, the simulation
time step is chosen to be 10ns and FFT is applied for 10 sine cycles.
As a Comparison to the previous explained full bridge Bipolar PWM ,the full
bridge Unipolar PWM has the following advantages over Bipolar PWM :
Chapter 2. Conceptual Framework 29
• Less harmonics content: Unipolar PWM has harmonics only at even multiplier
of mf .[11,12,32,33]
• More power efficiency : IGBT produce the most power losses during the transi-
tion on-off or vice versa. Therefore, the higher the switching frequencies the more
loss power averaged over time. On the other hand higher switching frequency
required to reduce the complexity and cost of the low pass filter . The unipolar
PWM strategy doubling the switching frequency at the outputVpwm while keep-
ing the switching frequency the same for each IGBT . That is a great advantage
for the Unipolar PWM.[11,12,32,33]
High quality inverter output voltage signal, free from baseband harmonics or
with the generation of intended harmonics is an important issue for this thesis. In
the Unipolar PWM, the timing of switching from high to low or vice versa is crucial
process. It is the main key which links between the time domain and the frequency
domain (harmonics), any deviation in the time of switching(from the ideal case) will
yield to harmonics. On the other hand, justifying the timing process of the pulses (time
of switching) according to a required frequency domain characteristics (with/without
intended harmonics) of Unipolar PWM is the main concept of Alfa( α ) Approach.
The following mathematical analysis is come out through an effort which has
been put in order to employ the theoretical mathematics for practical solutions for the
Unipolar PWM switching scheme[10,11,12,19,20]. To simplify the problem, a Unipolar
PWM switching scheme is assumed which has mf =8 and ma=0.8 , in such a case (figure
2.30) the signal of Unipolar PWM ( V oP W M ) is quarter-wave symmetry odd function.
• n= 1,2,3... mf
∞
X
f (t) = av + ak cos kw0 t + bk sin kw0 t (2.15)
k=1
where :
1 RT
• av is the average DC value f (t)dt
T 0
2π
• wo =
T
where T again is the period of the function (figure 2.29) and k=1,3,7,9...
Now applying equation 2.16 for the V oP W M (t) where αn = tnw0 , and then
substitute the result in the general form of Fourier series(equation 2.20) , the result
can be written like the following form :
4 ∗ VDC
hv(k) = [cos kα1 −cos kα2+cos kα3 −cos kα4+cos kα5 −cos kα6+cos kα7 −cos kα8]
k∗π
(2.17)
where
• k=1,3,5,7 ......2*mf -1 .
Chapter 2. Conceptual Framework 31
• hv represent the amplitude of the kth harmonic where hv1 is the fundamental
frequency.
2∗mf −1
4 ∗ VDC X
hvk = − (−1)n cos kαn (2.18)
k∗π n=1
where n=1,2,3,4....mf
This method is used to eliminate the harmonics from the PWM signal [10,11,12,19,20],
For this thesis the main idea is to employ it to determine the amplitudes of the intended
harmonics.
3.1.1 Entrance
Assuming that a general sinusoidal signal compared with a carrier triangle signal
to generate PWM as in figure 3.1 , where t1 t2 t3 .....t2mf are the intersection points
of the triangle and sinusoidal functions and Ttri is the period of one triangle cycle.
32
Chapter 3. Pec-HV 33
where Dcy(n) is the PWM duty cycles matrix [Dcy(1) Dcy(2) Dcy(3)...Dcy(mf )
] and n=1,2,3,4....mf .
At the implementation level, if the Dcy(n) is stored in the DSP instead of a sam-
pled sinusoidal, and if a proper DSP with high resolution counter for the triangle signal
is chosen( > 1.9 ns in this thesis) , then such a switching scheme can be considered as
Naturally Sampled PWM. In this switching scheme, the DSP counter starts counting
from zero until its value equals to the corresponding stored Dcy(n) value and simulta-
neously switching on/off the PWM signal as seen in figure 3.1. Actually that will lead
to a small deviation in the switching time ∆t compared with the original t1 and t2,
but since both times(t1 and t2) are shifted by the same amount (because the carrier
is a equal two sides triangle), the duty cycle in the total triangle cycle (Ttri ) remains
exactly the same. If mf is large enough (roughly mf > 30) the ∆t is very small thus
this timing deviation can be neglected, therefore the spectrum of this method is in
agreement with Naturally sampled PWM. (System validation in clause 3.4.2 )
As a backward solution, if the Dcy(n) values are known and for simplification
the values of ∆t are neglected, then the corresponding switching time t1, t2, t3...2mf
can be calculated using the following equations:
Ttri
t(2n−1) = (Dcy(n) ∗ ) + (Ttri ∗ (n − 1)) (3.2)
2
Ttri
t(2n) = (−Dcy(n) ∗ ) + (n ∗ Ttri ) (3.3)
2
where n=1,2,3,...mf
The desired signal Vsine in this switching algorithm (as seen in figure 3.2) consist
of the fundamental frequency Vsineh1 and a 3rd harmonic Vsineh3 (as an example, adding
more harmonics are possible ), since the switching strategy is based on the Unipolar
Chapter 3. Pec-HV 34
PWM the signal Vsine`is required which is also consist of Vsineh1`and Vsineh3`, where all
of the signals are compared to one carrier frequency Vtri .
The concept based on the calculation of the intersection points between each
of the fundamental frequency and its harmonic signals( Vsineh1 , Vsineh1` , Vsineh3 ,
Vsineh3` ) and the carrier frequency, then to calculate the PWM duty cycles values
according to the calculated intersection points, and then adding the duty cycle values
of the fundamental frequency and the 3rd harmonic (again, adding more harmonics
is possible) in order to find out α values for the desired signal one has to perform
calculation steps as the following :
Chapter 3. Pec-HV 35
• Finding the intersection points(t1h1 , t2h1 .... t2mfh1 ) of the signals Vsineh1 and
Vtri , and then calculating the duty cycle values Dcyh1(n) using equation 3.1
substituting t(n)h1 instead of t(n).
• Finding the intersection points(t1h3 , t2h3 .... t2mfh3 ) of the signals Vsineh3 and
Vtri , and then calculating the duty cycle values Dcyh3(n) using equation 3.1
substituting t(n)h3 instead of t(n).
• Adding the matrices Dcyh1 and Dcyh3, considering that it is not allowed to have
duty cycle more than 1 (100%), thus a subtraction of 0.5 (50%) from each duty
cycle value of the harmonics is required except for the fundamental frequency.
So the total duty cycle is Dcyh1h3 =Dcyh1 +(Dcyh3-0.5) . As a general Form :
• Substitute the elements of Dcyh1h3(n) in the equations 3.2 and 3.3 in order to
find the switching times corresponding to the duty cycles of Dcyh1h3 , the result
matrix is t(n)=[t1 , t2 ,t3 ....t2mf ]. Amusing that tk is the odd values of tn and
wo is the angular frequency of the desired signal Vsine then the corresponding α
values of the signal SAh are
mf
where n=1,2,3,4... 2
and k=1,3,5,...mf − 1
and
• Applying the explained above bullets for the signals Vsineh1`Vsineh3`in order to find
out Dcyh1h3 ` and then the signal SBh. The result from substituting Dcyh1h3 `
in the equations 3.2 and 3.3 is t0 (n) = [t01 , t02 ,t03 ....t02mf ]. Amusing that t0k is the
odd values of t0n and wo is the angular frequency of the desired signal Vsine then
the corresponding α values of the signal SBh are
where n=1,2,3,4... mf
2
and k=1,3,5,...mf − 1
and
α(4n − 2) = tk w0 and α(4n − 1) = tk+1 w0 (3.8)
Chapter 3. Pec-HV 36
mf mf
where n= 2
, 2
+ 1 , mf
2
+ 2 ....mf and k=mf + 1 , mf+3, mf+5,....2*mf-1.
For example if mf =8 , then t0 (n)=[t01 , t02 ,t03 .....t016 ] and the corresponding α(n)
values for SBh are [α1 , α4 , α5 , α8 , ........, α19 , α22 , α23 , α26 , α27 ,α30 , α31 ]
• Using the signals SAh and SBh instead of SA and SB(figure 2.29) respectively
in order to generate V oP W M as explained in clause 2.5.4 ,where V oP W M in this
case contains the intended harmonics amplitudes, thus it is renamed to V ohP W M .
The α values of the signal( V ohP W M =SAh-SBh ) are in agreement with equation
2.18 , as the En-Unipolar-PWM is validate in clause 3.4 .
The inputs for the Matlab model are the signal’s parameters of figure 3.2(Vtri ,
Vsineh1 ,Vsineh1 `, Vsineh3 , Vsineh3 `) and in addition to the 5th harmonic parameter. The
system signals equations are:
•
2 ∗ V ptri π
Vtri (t) = arcsin [sin(wtri t − )] (3.9)
π 2
where , wtri is the angular frequency ( 2πftri ) of the carrier triangle signal and
V ptri is the peak amplitude(in Volt) of the triangle signal.
•
Vsineh1 (t) = V psineh1 sin w0 t (3.10)
•
Vsineh1`(t) = −V psineh1 sin w0 t (3.11)
•
Vsineh3 (t) = V psineh3 sin( 3w0 t + φh3) (3.12)
where φh3 is the phase shift between the fundamental frequency Vsineh1 and the
3rd harmonic Vsineh3
•
Vsineh3`(t) = −V psineh3 sin( 3w0 t + φh3) (3.13)
Chapter 3. Pec-HV 37
•
Vsineh5 (t) = V psineh5 sin( 5w0 t + φh5) (3.14)
where φh5 is the phase shift between the fundamental frequency Vsineh1 and the
5th harmonic Vsineh5
•
Vsineh5`(t) = −V psineh5 sin( 5w0 t + φh5) (3.15)
As seen in figure 3.3, the first step is to determine (V ptri and V psineh1 ) where these
values are used to calculate ma, and then the desired signal frequency (fsineh1 ) and
the carrier frequency (ftri ) where they used to calculate mf . The intended harmonic
amplitudes (V psineh3 and V psineh5 ) in addition to V psineh1 must not exceed V ptri . Phase
shift (φh3 and φh5) can be added as well as the dead time dt which can be determined
according to the switching element (IGBT) manufacturer datasheet.
The next step is to use the input parameters to find out the intersection points(see
figures 3.2 and 3.3 ) and then duty cycle values: Dcyh1 , Dcyh3-0.5 ,Dcyh5-0.5 and
then add them together in order to calculate Dcyh1h3h5 which is used to find SAh
α values. Applying the same concept to find out Dcyh1h3h5 ` and the corresponding
SBh α values.
Instead of using the equations 3.5, 3.6, 3.7 and 3.8 to find out V ohP W M α values
one can simply combine and sort the two matrices α’s of SAh and α’s of SBh. The
result up to this step is a matrix which contains the V ohP W M α values.
Dead time as explained before is a waiting time during the switching on/off
between each pair of the H-bridge IGBTs , for the Unipolar PWM strategy it is
considered as a positive time shift for the odd α values of V ohP W M . In this model dt
is denote to the dead time.
This PWM switching scheme was developed to allow the use of PWM for PD
measurement and dielectric tests, it can be implemented in any 16-bit DSP which
includes PWM modules.
The time of timer1 will determine the fundamental signal frequency where timer1*mf
is the time of a complete desired sine cycle. Timer2 is 3 times faster than timer1, there-
fore it is able to process 3 complete sine cycles compared with timer1 to generate the
3rd harmonic. Timer3 is 5 times faster than timer1 where it is able to process 5
complete sine cycles compared with timer1, thus it generates the 5th harmonic.
Two counters are programmed inside each ISR, one to count forward for the
generation of SAh and the another in backward for the generation of SBh. Multiplying
the Duty cycle value by a factor from 0.1 to 0.99 will determine the amplitude where
selecting the starting point (n, n3 ,n5 and n`, n3`, n5`) in each counter will determine
the phase shift. As an example if mf =250 and 90˚phase shift between the desired sine
and the 3rd harmonic is required, then the algorithm has to wait until zero crossing
(normally at one complete cycle of the fundamental desired frequency) and then to
update the values of n3 and n3`according to the following :
360˚
• = 1.44˚which the resolution of phase shift.
250
90˚
• = 62.5 will be rounded to 63 , so n3=63 and n3`= 250-63 = 187
1.44
The voltage amplitude factor (Vf in Timer1) can have a value from 0.0001 to
0.9999 which will determine the modulation index ma, as well as Vf3 and Vf5 in
Timer2 and Timer3 respectively can have a value from 0.001 to 0.999 will determine
the amplitude of the 3rd and 5th using the shown equations in each timer where D5
denote to 50% duty cycle.
4
While Timer1 just pass its period, the duty cycle registers PDC1 and PDC2
5
of the DSP PWM module will update their values according to shown equations in
Timer1 (figure 3.4), where more details regarding a programming options in the DSP
PWM module are taken in account in order to avoid any timing errors, so each time
PDC1 and PDC2 were updated their values at the exact required values of Dcyh5,
Dcyh5`as well as Dcyh3 and Dcyh3`.
As seen in figure 3.5, the duty cycles table stored in the DSP which controls the
switching process of the En-Unipolar-PWM. Since the main target of En-Unipolar-
PWM is to allow the usage of PD measurements as well as HV-DI tests with a specific
amplitudes of intended harmonics, the analysis of the output signal V ohP W M in fre-
quency domain(up to 500 kHz) is required for the validation test and for the low pass
filter investigation (clause 3.5).
Time domain EMI measurements are alternatives for the methods working in
frequency domain which can be performed with a measurement receiver. The usually
possible bandwidth with a measurement receiver is the B=200 Hz, in order to have a
frequency resolution other than the 200 Hz (e.g. 34 to 120 Hz), the time domain EMI
measurements method is applied in this research work.
In order to meet the previous conditions, an oscilloscope with 12 bit ADC and
high memory capacity is used as acquisition unit for the test, the chosen sampling
frequency is 10 MHz, thus according to equation 3.16 the dynamic range is 74 dB
which is fulfilled the previous explained conditions. The measured data is analysed
then using Matlab FFT function.
• 12.47 kHz (1/ftri = 80.1875 us) for the fundamental of about 50 Hz.
Chapter 3. Pec-HV 43
• 12.47 kHz for the 16.7 Hz , that can be achieved by repeating each duty cycle
value for 3 times . i.e the duty cycle value will be updated once each 3 triangle
cycles (3* 80.1875 µ s) .
The dead time is according to the IGBTs data sheet is chosen to be 1.57 µs, but
due to a deference propagation delay time between the raising edges and falling edges
of the used optocouplers in the IGBTS drivers ( according to data sheet information)in
addition to another issues, the real effective dead time is reduced to much lower values
(explained in the discussion at the end of this clause).
A- Matlab Model validation test at zero dead time and zero amplitudes of
intended harmonics
• The signal V ohP W M has a steps time resolution of 100 ns and was simulated for
one complete cycle of the fundamental frequency (20.047 ms).
Figure 3.6 shows the spectrum of simulated V ohP W M . It is obvious that the re-
sults are in agreement with the explained Naturally Sampled(NS) PWM especially the
baseband harmonics are not present, further more the carrier sidebanded frequency
has a similar peak amplitude compared with the presented data in table 2.1 as well as
reference [11](see table 3.1), since the baseband area is harmonics free , it is obviously
clear that the calculated α values by the EN-Unipolar PWM are in agreements with
previous explained α approach, in other words if a set of equations are founded ac-
cording to equation 2.18 and α values of this simulated exampled are substituted , the
result which is the amplitudes of hvk will equal approximately zeros for the ks from
k=3 until k=499.
Table 3.1 shows the compared sidebanded harmonics peak values of the simulated
En-Unipolar PWM with the previously mentioned parameters with the calculated val-
ues from reference [11]. The first column is the harmonic order which is always a mul-
tiple of mf (equation 2.15) and by multiplying its values by the fundamental frequency
(49.88 Hz) will give the exact frequency of the harmonic order(the second column).
The numbers in the 3rd and 4th column are almost the same, the small differences
is due to calculation round off error. That is clearly proofs that En-Unipolar-PWM
algorithm passes the first validation test.
Chapter 3. Pec-HV 45
Table 3.1: Comparison between the simulated results and the theoretical calculated
from[11]
• ma = 0.8 , mf =250 , dt= 1.57 µs for the measured data and dt =0.3 us for
the simulated data , V psineh3 =V psineh5 = 12.5% of V psineh1 which is 0.1 V at
VDC=1V (for the simulation) and 10 V for the measurements of V ohP W M which
is produced by the device Pec-HV, so the expected harmonics amplitudes are 1
V in this case for the 3rd and 5th.
• The simulated signal V ohP W M has a steps time resolution of 100 ns and it’s
spectrum was calculated for one time domain cycle of the fundamental frequency.
• The measured signal V ohP W M has a sampling time of 100 ns and recorded for
two complete cycles of the fundamental frequency.
The model MEn-Unipolar-PWM is used to simulate this test results for the
following two cases; (1) with zero dead time. and (2) with 0.3 us dead time. The
Chapter 3. Pec-HV 46
spectrum of the simulated signal for both cases are compared to the measured V ohP W M
(3rd case) from the device Pec-HV. The measurements are performed more than once
to assure that the signal is stationary.
Figures 3.7, 3.8 and 3.9 are present the results of this test for the three taken
cases, where table 3.2 presents a selected harmonics peaks values in volt of the 3 cases
up to roughly 500 kHz.
Figure 3.7: 1st case, Simulated spectrum of V ohP W M at dt=0 with generated 3rd
and 5th harmonics at 50 Hz
As seen in the figure 3.7, the 1st case results shows that the spectrum contains
the expected 3rd and 5th harmonics with the amplitudes of 0.0999 V and 0.0998 V
respectively, that is so closed from the desired amplitude (0.1 V), it is also clear that
the baseband area is free from any other harmonics rather than the intended ones.
In the second case (figure 3.8) where the dead time is taken into account, the
baseband area contains some undesired harmonics, that is due to the dead time effect.
This result is matched with the briefly explained in [11,12] about the dead time effects.
In the 3rd case(figure 3.9), which is the experimentally measured, the signal
V ohP W M spectrum shows a very good desired results, the shape is close to the desired
one. The amplitudes of the intended harmonics almost match the expected (0.9116 and
0.9234 V) of the 3rd and 5th harmonics respectively, as well as the amplitude of the
fundamental frequency(7.349 V) which is not far from the desired one(8 V). Further
more the baseband area is free from undesired harmonics(one of the targets). In spite
of this satisfying results, a clear interpretation is needed to explain the differences
between the selected dead time values(1.57 us in 3rd case and 0.3 us in the second
case). Before answering this question a clear reading of the three cases is presented as
numerical values in table 3.2 .
Chapter 3. Pec-HV 47
Figure 3.8: 2nd case, Simulated spectrum of V ohP W M at dt=0.3 us with generated
3rd and 5th harmonics at 50 Hz
Figure 3.9: 3rd case, Measured spectrum of V ohP W M at dt=1.57 us with generated
3rd and 5th harmonics at 50 Hz
Chapter 3. Pec-HV 48
Table 3.2: Compare the harmonics of V ohP W M spectrum in three cases: (1)
simulation with dt=0.(2)simulation with dt=0.3 us (3) measured data from Pec-HV
Chapter 3. Pec-HV 49
Before starting the expatiations of table 3.2 one has to take in account that: in
order to compare the simulation results of (case 1 and 2) of the harmonics amplitudes
with the measured one, the simulated values have to be multiplied by 10, since the
relation between VDC and the amplitudes of the harmonics is linear (clause 2.5.21).
The table (3.2) is divided vertically to five parts according to the frequency range,
the following is briefly reading and analysis of the table’s data:
1. Baseband area : the harmonics order from the first (fundamental) to the ninth
are selected because they have the highest amplitudes among the baseband har-
monics. The third part of the system validation test(clause 3.4.2-C ) is focused
only on this range of frequency for the three fundamentals (16.7 Hz 50 Hz and
60 Hz)
2. Carrier sideband frequencies around 2*mf : here all the harmonics peaks are
considered, the number of peaks and their shapes are almost identical in the
three cases, the differences are in the amplitudes. The measured harmonics
amplitudes mostly are less than the simulated, that is due to the difference of
dv/dt of the signal V ohP W M in time domain, where the measured pulse (will
be explained later) shows larger dv/dt and thus less harmonics amplitudes in
frequency domain.
3. Carrier sideband frequencies around 4*mf : here also all the peaks are considered,
the same as the previous point the signal shape in the three cases are almost
identical but with roughly +20 Hz shifted(case 3 shifted from case 1 and 2) , the
measured harmonics amplitudes are slightly less than the both simulated case
except for the harmonics 4*mf ± 1 for is between the fist case and the second.
4. Carrier sideband frequencies around 6*mf : here just some of the sideband fre-
quencies are considered, the measured harmonics amplitudes(case 3) are more
closed to the first case, the measured amplitudes are also shifted by roughly +26
Hz from the simulated harmonics. This shift is due to the slightly difference
between the simulation SWF and the measured one.
5. 100 - 500 Khz: here only the highest peaks are selected,while this frequency range
is important for the low pass filter investigations. Also the measured peaks shows
less amplitudes than the case 1 and 2 .
According to the shown results above, one can notes the following :
1. The dead time in case 2 is 0.3 us where it is 1.57 us in the 3rd case.
2. The baseband area of the measured data are more closed to case 1 , where no
dead time is considered in the simulation.
Chapter 3. Pec-HV 50
3. Most of the measured carrier sideband harmonics amplitudes are slightly less
than the simulated ones ( when the simulated are multiplied by 10).
In order to find out the reason behind that, the time domain signal of V ohP W M
is examined for case 1 and 3 of table 3.2. Figure 3.10 shows the 7th cycle (as an
example) pulse for the 1st and 3rd cases. Actually most of the pulses are examined,
as result it is found that the measured pulses are having a real dead time value ranges
from 0.1 us to 0.87 us where most of them has a real dead time value of roughly 0.3
us(the selected for case 2), while the real programmed dead time value is 1.57 us. That
can be understood by analysing the hardware prorogation delay behaviour, starting
from the DSP PWM signal outputs pins, then going through the optocoupler of the
IGBT’s driver and ending by the IGBT switching on/off process. The used optocoupler
(according to datasheet) has more propagation delay by roughly 600 ns for the falling
edges of pulses than for the raising edge, where the programmed dead time is creating
a delay for the raising edge, in other words it reduces the dead time effect by roughly
600 ns by means of delaying the falling edges.
Another cause which reduces the effect of dead time is the dv/dt of the IGBTs
switching on/off characteristics, where larger raising/falling time reduced the pulse
width (measured at 80% of the pulse amplitudes) since the switching on/off ranges
from 0.3 to 1.8 us. With respect to the explained above, the dead time of 0.3 us is
adopted in MEn-Unipolar-PWM for the simulations.
Since the measured pulses of V ohP W M has a larger dv/dt than the simulated
,thus it will have less harmonics amplitudes for the carrier sideband frequencies . On
the other hand, the simulation has to be performed at high time resolution (100 ns) in
order to keep the accuracy of switching at α values, as well as to use to same number
of data points for the FFT function.
One of the main goals of this thesis is to have a high quality output voltage signal
from the Pec-HV. In this test the THD will be evaluated for the three fundamental
frequencies (16.7 Hz, 50 Hz and 60 Hz) for both; the modelled system (MEn-Unipolar-
PWM) and for the device Pec-HV. The measurements are performed according to
the standard IEC61000-4-7 (Testing and measurement techniques General guide on
harmonics and inter-harmonics measurements and instrumentation, for power supply
systems and equipment connected thereto). The signal V ohP W M from Pec-HV is
filtered by a low pass filter, where the filter is needed here to keep the measurement
device save from high frequency disturbances.
• The VDC value is chosen to be 125 V for the simulation and 180 V for the
experimental measurement of Pec-HV(the VDC values are selected to simplify
result’s reading).
• The modulation indexma =0.8 for both. The dead time is modified to 0.3 µs for
simulations.
This test is performed in two parts, the first with zero amplitudes of the intended
harmonics and the second with 12.5 % of the 3rd and 5th harmonics. The ideal THD
for the fist test part is 0% , where it 17.68% for the second test part.
First Part:
Table 3.3 shows the results of using the Matlab model MEn-Unipolar-PWM to
find out the spectrum of V ohP W M for each fundamental frequency(16.7 Hz 50 Hz and
60 Hz), and then to calculate the THD up to the harmonic order 63 using equation 2.1
. It is clear that all the results are so closed from the ideal case (THD= 0%) as well
as the amplitudes of the fundamental frequencies (ideally 0.8*125= 100 V). The three
fundamental frequencies has slightly different from the ideal values (16.63 Hz 49.883
Hz and 60.01 Hz), because the model MEn-Unipolar-PWM is built in an accurate way
to simulate the implemented En-Unipolar-PWM where an integer numbers has to be
filled in registers to perform the required PWM switching process. Having zero dead
time means zero THD up to harmonic order 63 as proofed in clause 3.4.21, according
to various performed simulations, the THD amount depend on the dead time value (
direct proportional), but the THD in case of fixed dead time is inversely proportion to
the SWF (the carrier triangle period) , that is why the 16.63 Hz has the lowest THD
where 60.01 has the highest one.
Chapter 3. Pec-HV 52
Table 3.3: THD evaluation by system simulation at dt=0.3µs and zero intended
harmonics
Table 3.4 and the figures 3.11, 3.12, 3.13 shows the measured THD of the three
fundamental frequencies. The harmonics are measured by the power quality analyser
(janitza UMG511) which is directly connected at the output terminals of the device
Pec-HV, where the measured signal here is Vo . The results are so closed from the
simulated data in table 3.2, In other words the system Matlab model as well as the de-
vice Pec-HV are achieved what is required until this point(high quality output voltage
signal).
Table 3.4: Measured THD for the three fundamental frequencies at zero intended
harmonics
Figure 3.11: Measured THD for the 16.7 Hz at zero intended harmonics
Chapter 3. Pec-HV 53
Second Part :
This test is performed with the same parameters of the first part but a 12.5%
of 3rd and 5th harmonics were added. High accuracy (theoretical THD must be so
closed from measured/simulated) is required for generating an intended harmonic in
order to achieve one of the main goals of this thesis.
Table 3.5 shows the simulated THD with the generation of 12.5% 3rd and 5th
harmonics, the expected THD in this case is 17.68%. The simulated data for the
three fundamental are so closed from the expected THD, the maximum difference is
acceptable as well (0.3% at 60 Hz).
Table 3.5: Simulated THD of the three fundamental frequencies with intended
harmonics
The measured data which presented in the figures 3.14, 3.15, 3.16 and table 3.6
are in agreement with the simulated as well as the expected THD, also the maximum
deviation between simulated and measured (0.79% at 60 Hz) is within an acceptable
range.
Figure 3.14: Measured THD for the 16.7 Hz with intended harmonics
Chapter 3. Pec-HV 55
Table 3.6: Measured THD for the three fundamental frequencies with intended
harmonics
The results of the three validation tests(clause 3.4.2- A,B and C) are leads to the
conclusion that the system Pec-HV and MEn-Unipolar-PWM are fulfils their require-
ments and both are valid.
Chapter 3. Pec-HV 57
The device Pec-HV has a various dynamic output signal behaviour where the
voltage amplitude, frequency and the amplitudes/phase-shift of the generated intended
harmonics are controllable within a wide range of digitised steps. (for example ma of
the fundamental 50 Hz can have a value from 0.0790 to 0.9999 within a step resolution
of 0.0005).
Table 3.7 shows a simulated THD up to the harmonic order 63 at different mod-
ulation index ( ma ) without a generation of intended harmonics. Through performing
this simulation, it has been noticed that the amplitudes of the 3rd, 5th, 7th and
9th harmonics are roughly fixed around (0.4%, 0.3%, 0.2% and 0.1 %)respectively of
the maximum fundamental frequency(atma =1, mf =250 and dt=0.3us) regardless of
changing ma value later, that means if ma is decreased then the fundamental ampli-
tude will decrease too, thus the THD will increase and vice versa. These undesired
harmonics are produced due to the dead time effect and they are roughly fixed since
the dead time has a fixed value. i.e if ma=1 at VDC=1 V then the amplitudes of the
3rd, 5th, 7th and 9th are 0.004 V, 0.003 V, 0.002 V and 0.001 V respectively, even if
ma value is changed, these harmonics amplitudes remains almost the same.
ma THD[%]
0.2 2.2
0.3 1.51
0.4 1.41
0.5 0.91
0.6 0.77
0.7 0.66
0.8 0.57
0.9 0.51
0.99 0.47
Since the system is able to to generate intended 3rd and 5th harmonics(various
amplitudes and phase-shifts), therefore it is possible to compensate the undesired 3rd
and 5th (generated due dead time) by generating intended ones with an adequate
amplitude/phase-shift to reduce the THD, thus having a high quality output with an
THD less than 0.4%. Figure 3.17 shows a measured 50 Hz at (ma =0.8) where the 3rd
and 5th are eliminated by the generation of intended 3rd and 5th having the same
amplitudes of the undesired one but with 180˚as phase shift .
Chapter 3. Pec-HV 58
Figure 3.17: Eliminate the effect of dead time at the output voltage of Pec-HV
From the previous presented simulated and measured results, one can expect the
system behaviour for all of the possible output voltage signal shape of the fundamental
frequencies(16.7 Hz , 50 Hz and 60 Hz) with the exact desired THD, where difference
between the desired THD and the produced by Pec-HV is roughly 1% at the worst
case.
As explained in clause 2.3, a proper attenuation ratio of the low pass filter has
to be determined in order to attenuate the carrier sideband harmonics of the sig-
nal V ohP W M in the frequency range of (100 kHz to 500 kHz), thus eliminating the
background noise interference in the PD measurements. According to an experimen-
tal investigations, the attenuation has to be enough to keep any harmonic amplitude
roughly in the range of 35 uV to 5 uV at the high voltage terminal of the tested object
for the PD (back to figure 2.7).
1. Adjust the Pec-HV parameters and then measuring the signal V ohP W M in the
frequency domain using a measurement receiver according to CISPR-16-1-1. (fig-
ure 3.18)
2. Keeping the same setting of the Pec-HV and connect the device to a coupling
capacitor and PD MI as seen in figure 3.18.
3. Measuring the PD in the frequency domain at different selected bandwidth.
Chapter 3. Pec-HV 59
The result of this experiments shows that the sideband harmonics have the value
of roughly 1 pC if its amplitude is ranges from 5 to 35 uV within a selected measurement
bandwidth from 9 kHz to 300 kHz.
From the presented results in clause 3.4.2-B, the highest harmonic amplitude
appears around the harmonic order of 2*mf , and then start decreasing at higher fre-
quencies. For the filter investigation, it is important to know the highest harmonic
amplitude in the range (100 kHz to 500 kHz), which is expected around 100 kHz. The
highest peak has to be determined at a various values of ma or with/without intended
harmonics for the 3 fundamental frequencies.
Table 3.8: Highest voltage harmonics of the measured V ohP W M in frequency range
99.8 kHz to 500 kHz
Table 3.8 shows the measured highest harmonic amplitude for the three funde-
mantal frequencies using the same measurements parameters(ma=0.8 at VDC=10)
and method which are mentioned in clause 3.4.2-B As seen in the table the highest
peak is 0.3805 V, since the relation is linear between VDC and the harmonics am-
plitudes (clause 2.5.2.1) so the maximum is (320/10)*0.3805=12.176 V , where 320 is
the maximum VDC. The complete measured spectrum of the presented results shows
that adding an intended harmonics will slightly influence on the amplitude of carrier
sidebanded frequencies.
ma highest peak
V
0.2 0.030
0.3 0.0646
0.4 0.039
0.5 0.0594
0.6 0.0478
0.7 0.0543
0.8 0.0462
0.9 0.0485
0.99 0.0457
Table 3.9: Highest voltage harmonics peaks of the simulated (50 Hz) V ohP W M in
frequency range 99.8 kHz to 500 kHz at various ma
In order to design a proper Low Pass Filter, the test transformer’s transfer func-
tion has to be known, especially for the frequency range from 100 kHz to 500 kHz. That
can be achieved by building up a complete model for the transformer such presented
in the references [67, 68] or by lab measurements investigations. Since the modelling
methods are complex, the measurements investigations is chosen for this thesis.
• Input impedance : The input in such usage for the transformer is the LPF and
a low impedance voltage source (less than 50 Ω) whether it is the Pec-HV or the
mains.
• The output impedance: Normally in the PD measurements the load at the trans-
former secondary is a capacitive voltage divider in addition to the PD coupling
capacitor and the PD coupling device which is actually a high pass filter with 50
Ω.
• A variable transformer is used to control the mains voltage amplitude at the pri-
mary and the LPF to separate the high frequency component from the mains in
Chapter 3. Pec-HV 62
• A signal generator, a power amplifier and a current clamp are used to inject a
sinusoidal signal in the frequency range of 100 kHz to 500 KHz.
To investigate the effect of the core non-linearity, the first measurements are
performed as seen in figure 3.19 but without the PD coupling device and coupling
capacitors, where the secondary voltage was measured through the capacitor voltage
divider. The measurements are performed at selected frequency values from 10 kHz
to 1.4 MHz where the mains voltage is 10 Vrms and then raised up to 190 Vrms.
The results of this measurements are as seen in figure 3.20 which shows that the core
non-literary is increasing the attenuations.
where:
• The red points are the selected frequency values among the interested range.
The measurement are performed again at the worst attenuation case (10 Vrms
of the mains voltage) with connected the PD coupling device and capacitor as shown
in figure 3.19. The results are shown in figurer 3.21.
The curve in the figure 3.21 shows an attenuation of roughly -65 dB in frequency
range of 100 kHz to 180 kHz then the attenuation is reduced to - 43 dB at roughly 300
kHz. The shown curve is measured at 10 Vrms of the mains 50 Hz voltage. Increasing
the 50 Hz voltage will increase the attentions as measured before. The curve is different
from figure 3.20 due to a different output impedance at the secondary of the test
transformer (50 Ω of the PD coupling device). This result would be a positive point
for the filter design complexity and cost.
Chapter 3. Pec-HV 64
According to the presented information in the clauses 2.3 and in this clause, a
LPF as seen in figure 3.22 is selected to allow the usage of En-Unipolar-PWM for PD
measurements. The filter consists of a sine-filter followed by multiple stages of low pass
filters in order to achieve the required attenuations, that is in addition to transformer’s
attenuation in order to reduce or eliminate the V ohP W M signal’s harmonics interference
with the PD measurements.
Figure 3.23 shows the measured voltage-frequency response of the selected LPF,
the frequency response curve were measured at the mains (low impedance) as input
impedance and 1 MΩ as output impedance. High output impedance is selected because
it is more close to the real usage of the filter where it is normally connected to the test
Chapter 3. Pec-HV 65
Adding the two curves (figure 3.21 and 3.23) shows a satisfying voltage attenua-
tion for the frequency range of 100 kHz to 500 kHz as seen in figure 3.24.
Figure 3.24: The total transfer function Vout/ Vin of the LPF and the transformer
The harmonics amplitudes of the signal V ohP W M are getting smaller with the
respect to increasing the frequency, i.e the highest harmonic amplitude in the interest-
ing range (100 kHz - 500 kHz) is at roughly 100 kHz , then higher frequency harmonics
amplitudes are getting decrease (figure 3.9) , while the attenuation due to both (LPF
and transformer) are increasing with respect to the frequency as seen in the figure
Chapter 3. Pec-HV 66
The filter is tested for all the expected dynamic singles output which the Pec-HV
can produce. Figure 3.25 shows one of this results where the background noise is less
than 1 pC.
As mentioned in clause 1.1.2 and 2.2, the voltage breakdown is an expected result
of the DI-HV test, where the current could raise up to a destructive values especially
if the test voltage source has not been interrupted very fast (within few 10’s of us). To
avoid any kind of damage in the test system, a protection unit based of the current raise
up time (di/dt) at the primary side of the test transformer and on the voltage fall time
(dv/dt) at the secondary. Figure 3.26 shows how the protection unit is integrated with
the Pec-HV. A special algorithm is programmed and implemented in the DSP of the
protection unit(PU-DSP) ,the PU-DSP will interrupt the IGBTs switching process in a
case of voltage breakdown and/or if the current exceeds a certain limit. An operational
amplifier circuit is designed to match between the current(I) measuring point at the
primary part and the ADC of the PU-DSP circuit.Also Another operational amplifier
circuit is designed to match between the divided voltage (V) at the secondary and the
ADC of the PU-DSP.
Figure 3.26: Overview schematic of the Pec-HV and the Protection Unit
taking more than three samples consume more calculation time and less than three
will increase the detection error. If the calculated step is less than a dynamic pre-
calculated threshold, the process will be repeated again, and if the step is higher
than the threshold, a breakdown is detected and a trigger signal will be sent to the
main DSP(figure 3.26) in order to stop the IGBT’s switching process. If the main
DSP receives the trigger signal from the protection unit, then IGBTs 1 and 3 will be
opened and 2 and 4 will be closed ,thus reducing the injected energy in the breakdown
channel.
The time of the detection process and thus the total switching off/interrupting the
test voltage source is depending on the sampling time of the ADC and a propagation
delay of the trigger signal as well as the IGBTs gate signals to close/open the desired
IGBTs as explained before. The sampling time is experimentally adjusted to be 10us
where faster time gives more detection error, as well as the consumed time of calculating
Chapter 3. Pec-HV 68
the step value must be less than the sampling rate, in other words no samples will be
missed. The overall measured switching off time of the H-bridge in case of a detected
breakdown is roughly 50 us.
The protection against any high raise up current is required, not only at the
case of breakdown, but as a general protection from any unexpected cases. The device
Pec-HV can be used to accomplish the HV-DI test and PD measurements where a step
up transformer is installed to generate the high AC voltage. The maximum current
of the primary’s transformer coil could be less than the capable one of the Pec-HV,
that is in addition to the destructive current value which can raise up rapidly in case
of voltage Breakdown. Therefore a programmable limit (max Ilmit and min Ilimit)
as seen in figure 3.28 has to be determined, where max I is the maximum current
positive peak value of the H-bridge module which is (50 A) and min I is the negative
peak. Interpreting the switching process will executed in case of exceeding the value
of max Ilmit or beneath the value of min Ilimit .The interruption is executed in the
same way as explained for the voltage-breakdown case.
4. Experimental Procedures and
results
The previous presented THD measurements in the figures 3.11 to 3.17 were per-
formed at the output terminals of the Pec-HV (no load was connected). In this clause
as well as the following clauses, all measurements/test were accomplished while a high
voltage step up transformer is connected to the Pec-HV.
Figure 4.1 shows the Pec-HV test’s setup in the high voltage lab at the University
of Duisburg Essen. Pec-HV is connected to the 100 kV, 5 kVA transformer, the
Network Harmonic Analyser (NWHA) is connected in parallel to Pec-HV to measure
the voltage THD of the test’s signal, the same THD will transferred to the secondary
windings, since the secondary can be approximately considered as an open circuit. For
the PD measurements the Omicron device MPD600 was connected through a fiber-
optic, quadripole as a coupling device and 100 pF as a coupling capacitor to the HV
transformer. The system is calibrated for the proper PD measurements.
69
Chapter 4. Experimental Procedures 70
Figure 4.2: Undesired harmonics produced by the non-linear core of the test trans-
former
Since the Pec-HV is able to generate intended 3rd and 5th harmonics, it is possible
to compensate the undesired generated ones and thus preconditioning the test’s signal
quality by reducing the THD to roughly 0.5%(Figure 4.3). The compensation can be
done by generating the same undesired harmonics amplitudes, but with 180˚phase
shift
Figure 4.4 shows how the signal preconditioning improves the test’s signal quality
as seen by the PD measurements device, where 4.3 (a) shows the THD due the iron
core and (b) shows the test signal after the compensation of the undesired harmonics.
The second aim of this thesis is a reliable fast interrupting of the test voltage
source in case of breakdown during the DI-HV tests. Figure 4.5 shows the experiment
test setup in HV-lab, where the primary current I is observed through a current probe
by CH2 of an oscilloscope and the secondary divided voltage V by CH1. The point V is
connected as well to the dv/dt input of the Pec-HV’s protection unit. Two spheres was
connected to the secondary to simulate a breakdown in the air as dielectric material.
The experiment was performed at various distances between the spheres as range of
(5 mm to 40 mm), the machine Pec-HV was started up and voltage was increased
gradually, depending on the adjusted distance the breakdown was detected at nearly
24 kV peak up to nearly 100 kV peak. This range of breakdown is tested at the
fundamental frequencies of 50 Hz and 60 Hz, while at 16.7 Hz the maximum possible
voltage to apply is 13 of the nominal transformer’s primary voltage (220/3) V , thus
the maximum HV at the secondary is roughly 34 kV RMS.
Figure 4.6 presents one of the performed tests at 50 Hz where the distance be-
tween spheres was 5 mm. Both figures are for the same experiment but at different
time scale. The voltage breakdown occurs at 24 KV peak, due to that the current
raised up rapidly. Simultaneously the protection unit detected the breakdown volt-
age and current. At the moment of detecting the breakdown, the protection unit will
send a trigger in order to interrupt the H-Bridge switching process. The detection
and interrupting process as seen in figure 4.6 takes roughly 50µs. The voltage signal
oscillated rapidly at the moment of breakdown, that is due to the capacitive voltage
divider (figure 4.5). since the voltage has to sink to zero value within nano seconds,
the capacitors are not able to be discharged in this fast time, thus the voltage will
oscillate within few micro seconds unit the capacitors are completely discharged.
Chapter 4. Experimental Procedures 73
Regarding the time 50 µs :as explained in 3.7 , the dv/dt algorithm has a sampling
time of 10µs , that means the programmed algorithm takes 30 µs after the breakdown
to read another 3 samples, another nearly 10 µs to calculate the difference absolute
value, where the rest 10 µs is trigger propagation delay and switching off/on the IGBTs
of the H-bride.
The third aim is to reduce the PD measurements background noise when the
measurements are fed by the Pec-HV as a power source. The test was setup in the
same way as in clause 4.1 (figure 4.1). The background noise has to be evaluated
for the three selectable frequencies (50 Hz, 60 Hz, 16.7 Hz). In this test the signal
preconditioning was not performed while the main point from the test is to observe
the background noise level.
The figures 4.7 (a, b and c) presents the measured (in pC) level of the back-
ground noise according to the standard IEC 60270 for the three selected fundamental
frequencies.
The background level is measured within the bandwidth (100 to 200 kHz) and
ranges from (0.4 to 1.2 pC) for the three cases(50 Hz, 60 Hz, 16.7 Hz) . It is obvious
as seen in figure 4.7 that the noise is proportional with the voltage level (if measured
within the same bandwidth). As explained in clause 3.6, the highest carrier sidebanded
frequencies are more close to 100 kHz, that means if the measurements band width is
shifted (e.g : 200 kHz to 500) the background noise will be reduced to be less than 1
pC. In this case the measurements bandwidth is still within the recommended range
according to (IEC 60270), but f 1 (figure 2.8-b) is out of the recommended range. Since
the PD pulses has a wide spectrum range and through many lab experiment, it has
been noticed that this shift of f 1 will not change the value of PD in pC in case of
having real PD(adding a pin at the high voltage terminal).
Chapter 4. Experimental Procedures 75
Pec-HV is able to generate 3rd and 5th harmonics. The harmonics amplitude
and phase shift (with respect to the fundamental frequency) are controllable.
Again the same test setup of figure 4.1 was used for this test experiment. The
test was accomplished for the three fundamental frequencies at various amount of 3rd
and 5th harmonics amplitudes. The voltage peak of the generated signal (fundamental
+ intended harmonics) must be within the linear region of the system (figure 2.20). In
other words the peak voltage of the signal Vsine in figure 3.2 must be less than the peak
of Vtria . Through these experiments it was noticed that generating the harmonics dose
not effect on the background level.
Figure 4.8 shows one the performed PD measurements where a pin is fixed at the
high voltage terminal to simulate a real PD pulses. The machine Pec-HV is started up
without the generation of intended harmonics in as seen in figure 4.8-(a), where a 5th
harmonic of 10% of the fundamental frequency’s amplitudes was generated as seen in
Figure 4.8 (b) .
The interested point for this thesis, is the quality of the generated signal, in
other words the deviation between the required THD and the generated one. Figure
4.9 shows the quality of the generated 10% of the 5th harmonics where the measured
THD is 10.02 % .
Chapter 4. Experimental Procedures 77
Modifying the Pec-HV switching algorithm in order to produce the required high
frequency signal is a challenge, since the impeded instructions in the DSP need to be
processed much faster in this case (i.e 500 µ s at 2 kHz) in order to complete the sine
cycle.
Figure 4.10: Modifying the switching algorithm to produce high frequency sine
80µs. The duty cycle value remains the same during the 80 µs. Therefore less duty
cycle values will be processed when increasing the sine frequency. This will not sharply
reduce the quality of the required sine signal. The signal quality was tested using the
a computer sound card scope program since the THD can be measured up to 20 kHz.
The signal worst THD is roughly 3.5% (through the frequency range 50 Hz to 2 kHz)
1
ISR T ime = (4.1)
F ∗ mf
where :
• ISR Time is the required time to update the duty cycle register in second.
Figure 4.12 shows as example a 2 kHz sinusoidal signal produced by the HF-
Unipolar-PWM measured at the output terminals of the Pec-HV, where the same
characterised LPF is employed to filter the carrier frequency and its side-band har-
monics.
The measurements are performed in the Lab of the University of Duisburg Essen.
Figure 4.13 shows the measurements setup, where a current transformer is used to
produce the required high current for the measurements, since Pec-HV has a maximum
current capability of 6 A the maximum current at the transformer secondary is roughly
Chapter 4. Experimental Procedures 80
200 A. The shunt resistor is connected at the secondary of the current transformer,
while a voltage drop at the shunt is measured using a low voltage measurement device
and the current is measured by a current probe and a scope. The shunt impedance
is then measured at a selected frequency values to cover the range from 50 Hz to 2
kHz. Figure 4.14 shows the result of this measurement , where the shunt impedance
is slightly varied with the frequency.
5.1 Conclusions
In this dissertation, the use of pulsed power electronic device for high voltage
testing is investigated and demonstrated experimentally.
The En-Unipolar-PWM algorithm was embedded in a suitable DSP and all the
required hardware elements were designed/combined in order to accomplish the Pec-
HV. Electromagnetic compatibility issues such, EMI and EMC filters, grounding and
shielding was taken in account.
An integrated protection unit was devolved for fast switching off the device in
case of breakdown as well as current surge.
• Produce High quality (T HD < 0.5%) one phase controllable amplitude sinu-
soidal signal with different selectable frequencies (16,7 Hz, 50 Hz, 60 Hz) as well
as the ability of signal preconditioning for HV tests (where a step up transformer
is used).
81
Chapter 5. Conclusions and Future work directions 82
• Fast switching off the applied voltage in case of voltage Breakdown (within
roughly 50 µs).
• Generate harmonics (3rd and 5th) of fundamental selected frequency with ability
of controlling the harmonics amplitude and phase shift. Signal distortion in this
case (< 1%) from the intended one.
• The device power efficiency is relatively high, power losses in switching elements
(IGBTs) are 1.6% from the input power.
• The Pec-HV is able to produce high frequency and high power sinusoidal signal
up to 2 kHz and 1.32 kW receptively.
[1] IEC 60060-3:2006: High-voltage test techniques - Part 3: Definitions and require-
ments for on-site testing, .
[2] IEEE STD98158 4-2013. IEEE Standard for High-Voltage Testing Techniques, .
[3] V.M. Catterson, S. Bahadoorsingh, S. Rudd, S.D.J. McArthur, and S.M. Row-
land. Identifying harmonic attributes from online partial discharge data. In IEEE
Transactions on Power Delivery, volume 26, pages 1811–1819, July 2011.
[8] F. Martin. A. Thiede. Power frequency inverters for high voltage tests.
In HighVolt Prueftechnik Dresden GmbH. Universitaet Karlsruhe, Institut
fuer Elektroenergiesysteme und Hochspannungstechnik, 2007. Downloaded on
Mar. 2012. URL https://www.ieh.kit.edu/rd_download/Martin_HighVolt_
Kolloquium_2007.pdf.
83
Bibliography 84
[10] Lipo Holmes. Pulse width modulation for power converters:principles and practice.
Wiley-IEEE Press. 2003.
[11] Ned Mohan, Tore M. Undeland, and William P. Robbins. Power electronics con-
verters, applications,and design. THIRD EDITION, JOHN WILEY and SONS,
INC., 2003.
[12] Ph.D. Muhammad H. Rashid. Power electronics handbook devices, circuits, and
applications. Third Edition,Elsevier Inc, 2011.
[13] W. L A Neves and H.W. Dommel. On modelling iron core nonlinearities. IEEE
Transactions on Power Systems,, volume 8(2):p 417–425, May 1993. ISSN 0885-
8950.
[14] M. Vakilian and R.C. Degeneff. A method for modeling nonlinear core character-
istics of transformers during transients. IEEE Transactions on Power Delivery,,
volume 9(4):pages 1916–1925, Oct 1994. ISSN 0885-8977.
[15] Niklas Erik Rueger and Nuernberg. Zur anwendung eines bandseparierenden mod-
ulationsverfahrens mit niedrigen taktzahlen in der leistungselektronik. Leibniz
Universitaet Hannover. Institut fuer Antriebssysteme und Leistungselektronik,
Verein Deutscher Ingenieure, Verlag VDI Verlag, Dissertation, 2012.
[16] Florian Martin and Thomas Leibfried. An universal high-voltage source based on
a static frequency converter. university of karlsruhe, institute of electric energy
systems and high voltage technologies (ieh). brochure. 2006. Downloaded on Mar.
2012. URL https://www.ieh.kit.edu/rd_download/ISEI_2006_Martin.pdf.
[17] K. Moessner, T. Leibfried, J. Stolz, and M. Gamlin. Reduction of the back-
ground noise for partial discharge measurements in transformer test circuits fed
by static frequency converters. In Proceedings of 2011 46th International Univer-
sities Power Engineering Conference (UPEC),, pages 1–6, Sept 2011.
[18] F. Martin, T. Leibfried1, and M. Steger. Comparison between a voltage sourced
and a current sourced static frequency converter. XV th International Symposium
on High Voltage Engineering, Aug 2007.
[19] J.N. Chiasson, L.M. Tolbert, K.J. McKenzie, and Zhong Du. A complete solution
to the harmonic elimination problem. IEEE Transactions on Power Electronics,,
volume 19(2):pages 491–499, March 2004. ISSN 0885-8993.
[20] S. Ramkumar, V. Kamaraj, and S. Thamizharasan. Ga based optimization and
critical evaluation she methods for three-level inverter. In 1st International Con-
ference on Electrical Energy Systems (ICEES), 2011, pages 115–121, Jan 2011.
[21] Ayong Hiendro and JalanJend. A. Yani. Multiple switching patterns for shepwm
inverters using differential evolution algorithms. International Journal of Power
Electronics and Drive System (IJPEDS), Vol No.2(2088-8694), Dec. 2011.
[22] T. Taufik, M. McCarty, M. Anwari, and A.S. Prabuwono. Optimization of op-
erating parameters in a unipolar pwm inverter. In Applied Power Electronics
Colloquium (IAPEC), 2011 IEEE, pages 57–62, April 2011.
Bibliography 85
[23] Kai Mueller. Entwicklung und Anwendung eines Messsystems zur Erfassung von
Teilentladungen bei an Frequenzumrichtern betriebenen elektrischen Maschinen.
Universitaet Duisburg Essen. Dissertation 2003.
[24] Mazen Alzatari. Realization and testing of a fast switch-off unit for dielectric high
voltage tests. Universitaet Duisburg Essen, ETS, Master Thesis 2010.
[25] Fast switch-off circuit with control. no. 1.59/1., March HIGHVOLT Prftechnik
Dresden GmbH, Brochure 2008. Downloaded on Mai.2013. URL http://www.
highvolt.de/portaldata/1/Resources/HV/Downloads/1-59-2.pdf.
[28] A.K. Murugan and R.R. Prabu. Modeling and simulation of thirty bus system
employing solid state circuit breaker. In International Conference on Emerging
Trends in Electrical Engineering and Energy Management (ICETEEEM), 2012,
pages 98–103, Dec 2012.
[29] A. Winter and A. Thiede. A new generation of on-site test systems for power
transformers. In Conference Record of the 2008 IEEE International Symposium
on Electrical Insulation. ISEI 2008., pages 478–482, June 2008.
[30] A. Thiede Y. Huang and S. Schierig. Transformer test system based on static
frequency converter. Proceedings of the 16th International Symposium on High
Voltage Engineering, ISBN 978-0-620-44584-9, 2009.
[31] Chih-Yuan Chen, Yu-Hsien Lin, Jiann-Fuh Chen, and Ray-Lee Lin. Design and
implementation of dsp-based voltage frequency conversion system. In Interna-
tional Symposium on Computer Communication Control and Automation (3CA),
volume 2, pages 435–438, May 2010.
[32] Zhongwei Guo and F. Kurokawa. A novel pwm modulation and hybrid control
scheme for grid-connected unipolar inverters. In Applied Power Electronics Con-
ference and Exposition (APEC), 2011 Twenty-Sixth Annual IEEE, pages 1634–
1641, March 2011.
[33] S. Halasz, G. Csonka, A. A M Hassan, and B.T. Huu. Analysis of the unipolar pwm
techniques. In 8th Mediterranean Electrotechnical Conference, 1996 MELECON
96, volume 1, pages 353–356 vol.1, May 1996.
[34] O. Pop, G. Chindris, and A. Dulf. Using dsp technology for true sine pwm
generators for power inverters. In The 27th International Spring Seminar on
Meeting the Challenges of Electronics Technology Progress, 2004., volume 1, pages
141–146 vol.1, May 2004.
Bibliography 86
[35] S.M. Mohaiminul Islam and G.M. Sharif. Microcontroller based sinusoidal pwm
inverter for photovoltaic application. In 1st International Conference on the De-
velopments in Renewable Energy Technology ( ICDRET), pages 1–4, Dec 2009.
[36] M.A. Rojas-Gonzalez and E. Sanchez-Sinencio. Two class-d audio amplifiers with
89/90thd+n consuming less than 1mw of quiescent power. In International Solid-
State Circuits Conference Digest of Technical Papers,IEEE . ISSCC 2009., pages
450–451,451a, Feb 2009.
[37] Myoung-Il Choi Jeong-Chay Jeon, Hyun-Jae Jeon and Chee-Hyun Park. The
design and implementation of a 5 kw programmable three-phase harmonic gen-
erator. Journal of Electrical Engineering and Technology, Vol. 3, No. 2,:162–166,
2008.
[38] Jorge Cerezo. Class d audio amplifier performance relationship to mosfet pa-
rameters. International Rectifier IOR. Application Note AN-1070, 2005. Down-
loaded on Mar. 2012. URL http://www.irf.com/technical-info/appnotes/
an-1070.pdf.
[40] J. Hamman and F.S. Van der Merwe. Voltage harmonics generated by voltage-fed
inverters using pwm natural sampling. IEEE Transactions on Power Electronics,,
volume 3(3):297–302, July 1988. ISSN 0885-8993.
[41] M. Balda. Nonlinear least squares. April 2013. Downloaded on May.2013. URL
http://www.mathworks.com/matlabcentral/fileexchange/17534.
[42] Prof. Holger Hirsch. High Voltage Engineering. Universitaet Duisburg Essen ,ETS.
Lecture script 2010.
[45] IEC 60270: 2000 : High-voltage test techniques - Partial discharge measurements,
.
[47] Prof. Dr.-Ing. Dirk Peier. Erzeugung und Messung hoher Wechsel- und Stoss
spannungen, Praktikumsversuch, Hochspannungstechnik und elektrische Anlagen
,Universitaet Dortmund. Nov. 1999.
Bibliography 87
[48] Yen-Shin Lai and Fu-San Shyu. Optimal common-mode voltage reduction pwm
technique for inverter control with consideration of the dead-time effects-part i:
basic development. Transactions on Industry Applications, IEEE, volume 40(6):
pages 1605–1612, Nov 2004. ISSN 0093-9994.
[49] Lihua Chen and F.Z. Peng. Elimination of dead-time in pwm controlled inverters.
In Applied Power Electronics Conference, APEC 2007 - Twenty Second Annual
IEEE, pages 306–309, Feb 2007.
[51] S.R. Bowes and B.M. Bird. Novel approach to the analysis and synthesis of mod-
ulation processes in power convertors. Proceedings of the Institution of Electrical
Engineers,, volume 122(5):507–513, May 1975. ISSN 0020-3270.
[52] Iannelli Luigi (Eds.) Vasca, Francesco. Dynamics and control of switched elec-
tronic systems. Springer, 2012, XIV.
[53] H+H High Voltage Technology GmbH. product/Compact Test Systems, Com-
pany Brochure, Downloaded on Feb. 2012. URL http://www.hundh-mk.com/
download/pruefsysteme.pdf.
[54] S. Bahadoorsingh and S.M. Rowland. Modeling of partial discharges in the pres-
ence of harmonics. In Electrical Insulation and Dielectric Phenomena, CEIDP
09. IEEE Conference, pages 384–387, Oct 2009.
[55] user manual for the device MPD600, Brochure, Downloaded on Mai.2015. URL
https://www.omicronenergy.com/fileadmin/user_upload/pdf/literature/
MPD-600-Brochure-ENU.pdf.
[56] Oldrich Sekula Thomas Strehl Hendrik Elze Sacha Markalous Werner Weis-
senberg, Toni Wunderlin. Uhf-pd-monitoring and on-site-commissioning-test of
400 kv xlpe-insulated cable circuits. Jicable 2007.
[58] J.W. Wood, H.G. Sedding, W. K. Hogg, IJ. Kemp, and H. Zhu. Partial discharges
in hv machines; initial considerations for a pd specification. Science Measurement
and Technology, IEE Proceedings A, 140(5):409–416, Sep 1993. ISSN 0960-7641.
[59] Nikunj Shah. Harmonics in power systems causes, effects and con-
trol. Siemens Industry, Inc., May 2013 downloaded in Sep. 2013. URL
http://www.industry.usa.siemens.com/drives/us/en/electric-drives/
ac-drives/Documents/DRV-WP-drive_harmonics_in_power_systems.pdf.
[60] F. Koeslag, H.T. du Mouton, H.J. Beukes, and P. Midya. A detailed analysis of
the effect of dead time on harmonic distortion in a class d audio amplifier. In
AFRICON 2007, pages 1–7, Sept 2007.
Bibliography 88
[61] Heun-Jin Lee Kwang-Hwa Kim, Sang-Hwa Yi and Dong-Sik Kang. Setup of
standard pd calibrator and its uncertainties. Journal of Electrical Engineering
and Technology, Vol. 6, No. 5, pp. 677-683, 2011.
[62] Florian Krug and P. Russer. The time-domain electromagnetic interference mea-
surement system. Transactions on Electromagnetic Compatibility, IEEE, volume
45(2):330–338, May 2003. ISSN 0018-9375.
[63] Martin Aidam, Stephan Braun, Lutz Dunker, Jean-Claude Nickel, and Man-
fred Stecher. Proposed amendments of cispr 16 parts 1-1, 2-x and 3. Sta-
tus of FFT methods for emission measurements, May 2006 . downloaded
on Aug. 2013. URL http://cispra.iec.ch/FFT_JTF_Documents/CISPR_A_
WG1(DE-adhocFFT)06-01.pdf.
[65] Multiple Access Communications Ltd. Time domain emc emissions measure-
ment system, final report. Company Brochure. May 2004. Downloaded on
Sep. 2013. URL http://stakeholders.ofcom.org.uk/binaries/research/
technology-research/timedomainemc.pdf.
[66] Peter Russer. Emc measurements in the time-domain. Institute for Nanoelec-
tronics, Technische Universitaet Muenchen, Brochure, 2011 . downloaded on Sep
2013. URL http://www.ursi.org/proceedings/procga11/ursi/et-1.pdf.
[67] W. H. Tang and Q. H. Wu. Condition monitoring and assessment of power trans-
formers using computational intelligence, power systems,. Springer-Verlag London
Limited, 2011.
[69] Maximilian Heindl, Stefan Tenbohlen, Juan Velsquez, Alexander Kraetge, and
Ren Wimmer. Transformer modelling based on frequency response measurements
for winding failure detection. International Conference on Condition Monitoring
and Diagnosis,, September 6-11, 2010, Tokyo, Japan.
[71] Walt Kester. Analog-digital conversion. Analog Devices, 2004, ISBN 0-916550-
27-3, Chapter 2.