Modelling of Grid Connected Single Phase
Modelling of Grid Connected Single Phase
Modelling of Grid Connected Single Phase
Author . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
Certified by . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
Jorge Garcia Garcia
Associate Professor
Thesis Supervisor
Certified by . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
Cristina Gonzalez Moran
Associate Professor
Thesis Supervisor
2
Modelling of Grid connected Single Phase Two Stage Full
Bridge Converter for Photovoltaic Application
by
Manish Patil
Abstract
This Master thesis presents modelling of Single phase two stage grid connected full
bridge inverter. The whole system includes the PV modules, a boost converter and a
full bridge Inverter. Reference datasheet of Monocristalline Silicon SLK72M6L have
been used for the photovoltaic array. PV modules were made up of a number of PV
cells connected in series and parallel.
One array consists of six PV panel. Then two arrays are connected in paral-
lel. These PV arrays are used for generating electricity in PV generation systems.
Thereafter, PV array have been simulated under varying environmental condition
like varying irradiance, temperature and partial shading conditions. PV array evalu-
ated by varying number of series, parallel array and by physical parameters like Rsh
(shunt) and Rse (series) resistance. I-V and P-V operating curves studied closely
under open loop of PV array. Observed open circuit current and voltage are around
9.02 A and 44.5 V respectively.
In the second part, circuit have been designed such that it delivers constant and
step up DC voltage to load. Open loop boost converter has been studied by using
constant DC voltage source and observed output voltage and current for boost con-
verter. Thereafter, our study have been progressed by coupling PV array with boost
converter. The characteristics of I-V and P-V curves of photovoltaic module coupled
with boost converter were tuned by means of simulations to achieve Maximum Power
Point (MPP) with optimal voltage and current Vmpp and Impp respectively. To
achieve MPP, MPPT algorithm has been used.
In this system, output voltage of boost converter have been observed around 400
V to charge the DC link. In third model, full bridge inverter have been simulated in
open loop to achieve required power at load. LCL filter has been designed at load
side and achieved maximum power transfer. When power at input and output were
made almost equal, load has been replaced with grid. Simulated open loop model
has been replaced with closed loop model. Grid angle for phase lock loop have been
used for alpha beta and D Q conversion. Finally it has been verified. By means of
3
step by step simulation, power loss for grid connected two stage full bridge inverter
have been calculated. Finally full bridge inverter has been designed with parameter
almost similar to Ingecon@Sunlite 5TL model (transformer less).
4
Acknowledgments
I wish to express my sincere appreciation to all the people who helped me in many
ways throughout my Master candidature at University of Oviedo. I would like to pay
my greatest gratitude and appreciation to my supervisors, Dr Jorge Garcia Garcia
and Dr Cristina Gonzalez Moran of University of Oviedo, for their generous technical
support provided while troubleshooting the thesis project.
The personal and administrative support provided by LEMUR group, Dr Islam, Bas-
sam at various stages of my Master thesis degree are acknowledged with gratitude.
Timely personal and technical assistance received from Major Umesh & Sujay, co-
operate to commission the experimental setup of my master thesis within short period
of time. I am very much thankful to them. Thanks also go to my colleagues Umer,
Saikot for their valuable input to my research work.
Last but not least, I would like to thank my parents, my sister Swati and rest of
family for their unconditional love and continuous support. I would not have been
able to complete Thesis without them.
5
6
Contents
1 Introduction 17
1.1 Statement of the Problem . . . . . . . . . . . . . . . . . . . . . . . . 17
1.2 Research Objectives and Methodologies . . . . . . . . . . . . . . . . 17
1.3 Outline of the Thesis . . . . . . . . . . . . . . . . . . . . . . . . . . . 18
2 Literature Review 19
2.1 Introduction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 19
2.2 Solar Photovoltaic Energy . . . . . . . . . . . . . . . . . . . . . . . . 20
2.2.1 Types of Photovoltaic cell and Equivalent circuit of a PV cell 20
2.3 Power Converters for Integrating PV Systems to the Grid . . . . . . . 22
2.4 The Modulation Techniques . . . . . . . . . . . . . . . . . . . . . . . 23
2.5 Operation and control of VSC . . . . . . . . . . . . . . . . . . . . . . 26
2.5.1 Basic Operation of a VSC . . . . . . . . . . . . . . . . . . . . 26
2.5.2 Control Systems in a VSC . . . . . . . . . . . . . . . . . . . . 26
2.6 Grid-Connected PV System . . . . . . . . . . . . . . . . . . . . . . . 26
7
3.3 DC-DC Boost Converter Operation . . . . . . . . . . . . . . . . . . . 36
3.3.1 Mode 1:Operation of the Boost Converter . . . . . . . . . . . 36
3.3.2 Mode 2:Operation of the Boost Converter . . . . . . . . . . . 36
3.4 Modulation Technique of Single Phase Full Bridge VSC . . . . . . . . 37
3.5 LCL Filter . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 38
3.6 DC-link Capacitor . . . . . . . . . . . . . . . . . . . . . . . . . . . . 39
3.7 Control system in VSC . . . . . . . . . . . . . . . . . . . . . . . . . . 40
5 Simulation Result 55
5.0.1 Simulation of PV with variation in Temperature & Irradiance 55
5.0.2 E↵ect of Change in Irradiance for Photovoltaic output Power,
Voltage, Current, DC-link voltage & Active-Reactive Power . 56
A 63
A.1 Boost Converter . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 63
A.1.1 Selection of Inductor . . . . . . . . . . . . . . . . . . . . . . . 63
A.1.2 Selection of Input Capacitor . . . . . . . . . . . . . . . . . . . 64
A.1.3 Selection of output Capacitor . . . . . . . . . . . . . . . . . . 64
B 67
B.1 LCL Filter . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 67
B.1.1 Basic Equation . . . . . . . . . . . . . . . . . . . . . . . . . . 67
8
B.1.2 Designing Procedure for 5 kVA Single-Phase VSC . . . . . . . 69
C 71
C.0.1 Selection DC-link Capacitor . . . . . . . . . . . . . . . . . . . 71
D 73
D.1 Data Sheet SLK72M6L . . . . . . . . . . . . . . . . . . . . . . . . . . 73
9
10
List of Figures
3-1 I-V characteristic of a solar array for a fixed temperature & irradiance
[31] [38] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 32
3-2 P-V characteristic of a solar array for a fixed temperature & irradiance
[31] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 33
3-3 Flow chart for P& O [8] . . . . . . . . . . . . . . . . . . . . . . . . . 34
3-4 Principle of P& O [8] . . . . . . . . . . . . . . . . . . . . . . . . . . . 35
3-5 Circuit diagram of boost converter (mode 1) [32]. . . . . . . . . . . . 36
3-6 Circuit diagram of boost converter (mode 2)[32]. . . . . . . . . . . . . 37
3-7 Wave form of boost converters operation[32]. . . . . . . . . . . . . . . 38
3-8 Output current for S1, S2 ON; S3, S4 OFF for t1 < t < t2 [32]. . . . . 39
3-9 Unipolar PWM generator [12]. . . . . . . . . . . . . . . . . . . . . . . 39
3-10 Output Voltage Waveform [32]. . . . . . . . . . . . . . . . . . . . . . 40
3-11 shows Waveform for SPWM with Unipolar voltage switching [25]. . . 40
11
3-12 shows Carrier and Reference Waveform for SPWM with Unipolar volt-
age switching [25]. . . . . . . . . . . . . . . . . . . . . . . . . . . . . 41
3-13 shows LCL filter [30]. . . . . . . . . . . . . . . . . . . . . . . . . . . . 41
3-14 shows Schematic VSI including voltage and current control loop [43]. 42
3-15 shows Current control loop [43]. . . . . . . . . . . . . . . . . . . . . . 43
3-16 shows DC Voltage control loop [43]. . . . . . . . . . . . . . . . . . . . 43
3-17 shows MATLAB/SIMULINK results for real and imaginary phase wave-
forms [36]. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 44
3-18 shows Grid angle Calculation [36]. . . . . . . . . . . . . . . . . . . . . 44
3-19 shows Construction of the imaginary Orthogonal Circuit Variables [41]. 45
3-20 shows Real Circuit and its Imaginary Orthogonal Circuit [41]. . . . . 46
3-21 shows DQ transformation using constructed Imaginary Orthogonal Cir-
cuit Variable for fixed fundamental frequency application [41]. . . . . 47
5-1 I-V characteristic of a solar array for a fixed temperature but varying
irradiance [31] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 55
5-2 P-V characteristic of a solar array for a fixed temperature but varying
irradiance [31] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 56
5-3 I-V Characteristic of a PV array under a fixed irradiance but varying
Temperatures [31] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 56
5-4 P-V Characteristic of a PV array under a fixed irradiance but varying
Temperatures [31] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 57
5-5 shows Step change in the Photovoltaic Power Output [31] . . . . . . . 57
12
5-6 shows Step change in the Photovoltaic PV array Output Voltage [31] 58
5-7 shows Step change in the Photovoltaic PV array Output Current [31] 58
5-8 shows PV array current Ipv reaches steady state after 0.002 s [31] . . 58
5-9 shows Step change in the Single Phase Full Bridge Inverter Grid Cur-
rent [31] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 59
5-10 shows Step change in the Single Phase Full Bridge Inverter Grid Volt-
age [31] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 59
5-11 shows change in the Single Phase Full Bridge Inverter DC link Voltage
[31] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 59
5-12 shows change in the Single Phase Full Bridge Inverter for Active-
Reactive Power [31] . . . . . . . . . . . . . . . . . . . . . . . . . . . . 60
13
14
List of Tables
15
Glossary
16
Chapter 1
Introduction
17
simulation model of a PV system representing the dynamic behaviour of an actual
PV system connected to an LV power grid. Therefore, the dynamic behaviour of all
the physical components and all control loops should be accurately modelled. Data
sheet for Siliken PV panel given in appendix D.1 and IngeconSun 5TL model have
been refered to design simulation model [39].
18
Chapter 2
Literature Review
2.1 Introduction
In most of power generating plants, fossils fuels are main source of energy to generate
electricity. Till 19th century energy has been created largely by burning hydrocarbon,
oil, coal and natural gas. It produces huge carbon emissions. Emission of greenhouse
gases, limited coal resources, environmental distortion, rising prices of fossils fuels and
scarcity of foreign exchanges reserves responsible for developing new energy economy.
However, energy generation from renewable sources like wind, solar, hydro, biomass,
tidal, geothermal, wave, bio-fuels & bagasse does not a↵ect the environment that
adversely. Government and electric utility around world has introduced subsidies
and significant change in their traditional policies to allow investors to integrate
their own renewable power generating systems to electric power grid to promote
sustainable energy technologies. Further, commercially attractive rebates, grants and
credit schemes will continue to make significant changes to the traditional electrical
power grid.
Electrical power generated from a PV system solely depends on factors such as
solar irradiance and ambient temperature. These factors are subjected to sudden
or slow variations, resulting in changes in the output power. As PV systems are
interfaced to LV power grids through power electronic converters (inverters) and as no
mechanical inertia is involved, changes in power flow can take place in relatively short
19
periods. Such changes can cause rapid voltage fluctuations at the point of common
coupling (PCC) of a PV system and can lead to overvoltage, under-voltage and/or
flicker problems. Further, switching frequency and low-order harmonic levels may
increase with the addition of PV systems to LV power grids, possibly deteriorating
the quality of power supply.
Simulation studies can be performed to investigate the terminal characteristics
of a PV system under di↵erent operating conditions and to analyse how LV power
grids will behave when multiple PV units are integrated. In order, to perform such
simulation studies, the development of accurate simulation models of PV systems is
necessary. A detailed simulation model of a PV system that represents the dynamic
behaviour of all the physical components of a grid-connected PV system and the
associated control loops are essential to investigate the dynamic behaviour of a single
or multiple grid-connected PV system under di↵erent operating conditions.
All photovoltaic (PV) cells consist of two or more thin layers of semi-conducting
material. It is mostly silicon. Electrical charges are generated, when semiconductor
exposed to sunlight. It produces direct current when multiple cells are connected.
20
silicon panels silicon is deposited in continuous process onto base material giving
fine grained, sparkling appearance. Then it is encapsulated in transparent insulating
polymer with tempered glass cover and then bound into metal framed module. In
this case, monocrystalline panel SLK72M6L of SILIKEN has been used D.1 [38] [39].
A panel is 72 series cells. The peak power is between 285 Wp to 310 Wp.
The equivalent circuit of PV cell is shown in figure 2-1 b. Internal losses of current
in model has been ignored. A diode is connected in anti-parallel with current source.
When PV cell exposed to sunlight DC current is generated is denoted by I. By using
Kirchho↵ law output current, I can be obtained by equation 2.1,
V +IRs
I = Iph ID = Iph Isat [e Vt 1] (2.1)
In case of short circuit situation, the equation presented in 2.1 can be rewritten
as:
Isc = Iph ID (2.2)
where
Isc : Short circuit current (A)
Iph : Photo current (A)
ID : Diode current (A)
Isat : Saturation current (A)
Vt : Thermal voltage (v)
Rs : Series resistance (Ohm)
V : Terminal voltage (V)
I : Output current (A)
21
Figure 2-1: shows Equivalent circuit of a solar cell and PV device [3].
22
direction of average power flow.
To control the power flow through a CSC, it should be designed from fully con-
trollable bipolar power electronic switches. The power flow of a VSC can be fully
controlled if the converter is made of reverse conducting power electronic switches
IGBT or MOSFET. These switches basically used for low power converter applica-
tion. That is reason VSC converter become most widely used type power electronics
converter to interfaces low power PV system to grid.
In this dissertation, single phase two stage DC-DC-AC converter has been used.
DC-AC inverter in single phase full bridge transformer less inverter is usually preferred
for relatively small domestic rooftop PV installations.
Single phase voltage source converter has been used for H bridge bipolar or unipolar
modulation method. By controlling the magnitude and the angle of voltage at the
output of the H-bridge, v, of the VSC with respect to the grid side voltage, vg , amount
of power injected or absorbed has been controlled. Basic structure of Single-phase
VSC has been shown in figure 2-2.
23
signal(m) is proportional to magnitude, phase angle and frequency of voltage that
should be generated at the output of the H-bridge. Sinusoidal signal is compare
with high frequency triangular signal using relational operator to get pulse-width
modulated switch-ON and switch-OFF pattern. It is shown in figure 2-3 a and b.
The switching pattern is then applied to the power electronic switches of the H-bridge
shown in figure 2-4 and figure 2-5. The current injected by the VSC to the grid is
highly distorted. Filtering mechanism has been used in order minimise the distortion
in current.
Figure 2-4: shows Output voltage of a H-Bridge of a VSC that uses unipolar-SPWM
technique [30].
24
Figure 2-5: shows Output voltage of a H-Bridge of a VSC that uses bipolar-SPWM
technique [30].
due to reduced losses during zero voltage states and low losses in the output filter.
The advantage of bipolar modulation is that there is no switching harmonics, very
low leakage current and EMI
Other technique which is largely used in power electronics is space vector Pulse
width modulation (SVPWM). Space vector PWM control scheme for three-phase
PWM inverter in PV generation system discussed in Jiyong Li et al. [34]. From
Vdc input voltage of the inverter, we can calculate each state voltage. Once voltage
in reference frame is determined, duration period T1 and T2 time for each vector in
each PWM cycle has been calculated. Low total harmonic distortion (THD), constant
switching frequency, well-defined output harmonic spectrum, optimum switching pat-
tern and excellent dc-link voltage utilization are advantages with SVPWM. In Q. Zeng
et al. [6] explained drawbacks for SVPWM [31].
25
2.5 Operation and control of VSC
To reduce ripple in output current, VSC is integrated to grid via. coupling inductor.
It also minimises the e↵ect of the PWM switching voltage developed at the output
of a VSC. To improve limited capability of attenuating switching current ripple of
VSC, LCL filters often used in grid-connected VSCs. LCL filter which shown in figure
2-6, attenuates high frequency switching ripple current of VSC. LCL filter design also
a↵ect dynamic performance of control system of grid connected VSC. While designing
LCL filter accuracy should be high. In detail procedure of LCL filter also explain in
appendix B.
Based on control method, VSC can be classified as either voltage or current controlled.
By adjusting the phase angle and the magnitude of the output voltage of the VSC
relative to the PCC voltage, real and reactive power output are controlled. Due to lack
of close loop to control its output current, there is no over current protection in VSC.
But in case of current-controlled VSC the output current of the VSC is controlled. The
control of the active and reactive power outputs is achieved by controlling the output
current. By adjusting the magnitude and phase angle of the output voltage of the VSC
with respect to the grid voltage, current has been controlled. In a current controlled
VSC, the closed-loop current controller and grid synchronisation mechanism are the
two main control systems being used to integrate PV systems to the grid.
26
Figure 2-6: shows Grid connected VSC with LCL filter and PV System [30].
27
28
Chapter 3
Equivalent circuit of solar cell with one diode and some equation has been discussed
in chapter 2. Equation which responsible for output current I and diode current Id ,
discussed in chapter 2. In this chapter, equation and mathematical equivalent circuit
has been derived for photovoltaic array.
The topologies are characterized to calculate several unknowns. It classified with
respect to total number of unknowns. Normally the single diode model is known
as five parameter model. And other one double diode is known as seven-parameter
29
model. Both model based three characteristic points of the current-voltage (I-V)
curve i.e. open-circuit condition, short-circuit condition, and MPP [26] [37]. Five
characteristics has been considered midway between MPP and open-circuit condition,
at one-half of the open circuit condition [18]. For these models, mathematical-based
equation has been used and cannot be used directly in a circuit simulator. Considering
the environmental conditions (ECs), mathematical models of a PV cell has been
derived. While designing model, parameters has been used usually available in Siliken
SLK72M6L manufacturers PV module datasheet [38].
Correct value of the parameters Rs , Rp , Isat , Iph under the actual ECs are not
available at all in the datasheets of the PV modules. This parameter can be used in
equation 3.2, after determining correct value of above parameter. Under standard test
condition short-circuit current (Isc ), open-circuit voltage (Voc ), voltage and current
in MPP (Vmpp & Impp ), maximum power (PM P P ), coefficients ( Isc ) (temperature
coefficient for the short circuit current) and ( Voc ) (temperature coefficient for the
open circuit voltage) these parameters are defined. Diode current can be neglected
and can be very small in comparison with Isc . Photo current Iph is equal to Isc current
can be expressed as each irradiation as cell temperature as,
V +IRs
I = Iph ID = Iph Isat [e Vt 1] (3.1)
30
Ga
T = (Top ) + (N OCT (Tn )) ⇤ ( ) (3.3)
Gan
where,
Voc : Open Circuit Voltage (V)
Vocs : Open Circuit Voltage in STC (V)
Voc : Open Circuit Voltage Coefficient with Temperature (V/K)
Voc (T )
ID (Ga , T ) = Iph (Ga , T ) = Isat (Ga , T ) ⇤ [e Vt (T ) 1] (3.5)
T
Vt (T ) = A ⇤ K ⇤ (3.6)
q
Where equation 3.6 is the thermal temperature (V ), A the diode ideality factor,
q the charge of an electron and k the Bolztman constant. With equations 3.1 to 3.6
the operating point for each pair of irradiance Ga and room temperature Top can be
determined, once A and Rs are known.
Manufacturer of the solar module gives another parameter needed to model the
solar cells. The datasheet which gives the electrical characteristics is calculated under
standard test condition STC when the temperature T is 25 and the irradiance G is
1000 W/m2 .
Matlab Simulink I-V model has been constructed using all equation. Technical
characteristics sheet provided by manufacturer for monocrystalline SLK72M6L SI-
LIKEN panel given in table 3.1. Panel has series of 72 cells. The peak power per
31
Table 3.1: Technical Characteristics for SILIKEN SLK72M6L
Parameter Value
Short circuit current Isc (A) 9.02
Open circuit voltage Voc 44.5
Current at maximum power point Imp 8.51
Voltage at maximum power point Vmp 35.8
Number of cells in series ns 6
05
Temperature coefficient of Isc 1.9780*10
05
Temperature coefficient of Voc -5.8608*10
Pmax (Wp) 305
panel 305 (W p) per panel, short circuit current coefficient, open circuit voltage co-
efficient with temperature and other parameter of panel measured in Standard Test
Condition (STC) given in table 3.1. For this experiment, 12 panel i.e. six set of
two array connected in parallel. In this case, array voltage and current have been
calculated from equation 3.1. This array voltage has been used as input to boost
converter.
Figure 3-1: I-V characteristic of a solar array for a fixed temperature & irradiance
[31] [38]
The characteristic I-V curve 3-1 shows that there are two regions in the curve: one
is the current source region and another is the voltage source region. In the voltage
32
Figure 3-2: P-V characteristic of a solar array for a fixed temperature & irradiance
[31]
source region (in the right side of the curve), the internal impedance is low and in the
current source region (in the left side of the curve), the impedance is high. Irradiance
& temperature plays an important role in predicting the I-V characteristic and e↵ects
of both factors should be considered while designing the PV system. Whereas the
irradiance a↵ects the output current, temperature mainly a↵ects the terminal voltage
[5] [30].
Nonlinear solar cell I-V characteristic changes with irradiation and temperature. It is
difficult to track unique point for the maximum power (PM P P ). With help of optimal
voltage Vmpp and the optimal current Impp , maximum power (PM P P ) has been located.
MPPT draws maximum power from the solar array regardless of weather or variations
in load conditions. There are diverse MPPT algorithms that have been manufacturing
for PV systems [9] [27].
Because of the simplicity in the algorithm, perturb and observe (P&O) method
has been identified as a simple MPPT algorithm. It needs only one sensor. The
P&O method has been identified as a slow tracking algorithm compared to other
available MPPT. In the occurrence of rapid variations of the environmental conditions,
it fails sometimes [8]. P&O method is not working properly in case of continuous
oscillations around the MPP. Other MPPT algorithm is the incremental conductance
33
(InC) method. As compared to P&O method, the InC method tracks the MPP quickly
and accurately [11]. InC method is working properly in case of continuous oscillations
around the MPP as compared to P&O method. It requires output current and voltage
feedback of the PV array needed to implement. The InC method is accurate, less
complex and easily implementable MPPT algorithm.
In P&O or InC method, firstly output voltage and current of PV module has
been measured. Then control parameter needs to be increased or decreased, once
PV power has been observed. Reference signal (voltage and current) or duty ratio
as controller for switching signal has been used for DC/DC converter. P&O MPPT
algorithm used in this thesis, does not require previous knowledge of the PV module
characteristics. DC-link capacitor has been calculated correctly to work MPPT al-
gorithm correctly. Switching frequency and step size has also been selected properly
as it a↵ecting algorithm. In section 3.2.1, searching algorithm, their flow chart and
theory has been discussed.
34
Perturb & Observe algorithm flow chart has given in the figure 3-3. It states that
when the operating voltage of the PV panel is perturbed by a small increment and if
the resulting change in power P is positive, then direction of MPP followed. Or else
keep on perturbing in the same direction. If P is negative, then direction of MPP not
followed. The sign of perturbation supplied must be changed in this case.
In this algorithm, voltage for PV array has been used as sensor. So, the cost
of implementation is less and hence easy to implement. The time complexity of this
algorithm is very less. By reaching very close to the MPP it does not stop at the MPP
and keeps on perturbing in both the directions. When the algorithm has reached very
close to the MPP and appropriate error limit has been set or can used a wait function.
It ends up increasing the time complexity of the algorithm.
In figure 3-4, on the left of the MPP, while incrementing the PV voltage, power of
the PV increases. By decrementing the PV voltage, the power of the PV decreases.
On the right of MPP, incrementing the voltage decreases the power and decrementing
the voltage increases the power. To extract the maximum power from the PV module,
MPPT controller has been implemented. Detail explaination regarding Incremental
& Conductance method explained in reference [22] [8].
35
3.3 DC-DC Boost Converter Operation
Output voltage of boost converter is always greater than input voltage so it is also
called step up converter. Boost converter has high efficiency only for reasonable gains
and for low power values.
Inductor gets charged through the DC source i.e. PV panel and stores the energy,
when switch is closed as shown in figure 3-5. Inductor charges linearly because output
DC source is constant due to large C. The diode blocks the flowing current and so
the load current remains constant. Current is being supplied to load due to the
discharging of the capacitor. Selection of input capacitor, inductor, output capacitor
and load resistor explain in appendix A.1.
In mode 2, figure 3-6 the switch is open. So, the diode becomes short circuited. The
energy stored in the inductor gets discharged through opposite polarities. Discharging
of inductor charge the capacitor. Throughout the operation, the load current remains
constant. Boost converter works properly because output C is very large. Due to large
C, variations in current ripples do not a↵ect output voltage. Figure 3-7 shows the
waveforms for a boost converter.
36
Figure 3-6: Circuit diagram of boost converter (mode 2)[32].
Figure 3-8 & 3-10 shows a simplified circuit diagram and output Voltage Waveform
of a Single-Phase H Bridge inverter. The logic behind the switching of the devices is
given below,
37
Figure 3-7: Wave form of boost converters operation[32].
The values of the inductors and the capacitor of the LCL filter of the PV system
3 3
shown in figure 3-13 were determined as Li = 3 ⇤ 10 H, Lg = 1.5 ⇤ 10 H and
6
Cf = 2.6 ⇤ 10 F, by following the design guidelines established in [1]. Resistances
of inductances Lf c and Lf g are assumed as Rf c = 0.5 Ohm. The total inductance of
the designed LCL filter is less than 0.1 pu (using a base - Sb = 5e4 kVA, Vb = 230 V)
and the resonant frequency, fres , is 10.7 kHz; less than 0.5fsw . A damping resistor
Rd = 1.5 Ohm, calculated as one third of the impedance of Cf at fres , is added in
series with Cf to improve the stability of the current controller [21]. Detail design for
LCL filter explain in appendix B.
38
Figure 3-8: Output current for S1, S2 ON; S3, S4 OFF for t1 < t < t2 [32].
The output capacitance of a power converter is a vigorous part of the overall system.
It maintains the output voltage. It a↵ects stability of the regulator. The component
parameters are part of the feedback control loop. Using DC link capacitor value,
reference d and q axis current has been designed. These references has been used for
further PWM signal. The output capacitance determines the position the dominant
pole. It therefore a↵ects the stability boost converter. Output voltage ripple of the
converter has been a↵ected by output capacitor. The output voltage peak-to-peak
ripple Vdc given by equation 3.7.
P
Vdc = (3.7)
(Cdc ⇤ Vdc ⇤ !)
39
Figure 3-10: Output Voltage Waveform [32].
Figure 3-11: shows Waveform for SPWM with Unipolar voltage switching [25].
voltage ripple in DC link voltage of single phase VSC. It calculated to limit peak-peak
ripple of DC-link voltage. V dc is average 5 or 2 % of V dcavg . An excessively large
capacitor at the DC-link can reduce the DC-link voltage ripple almost to zero as per
3.7. But that is practically impossible and there is a limitation to the size of Cdc .
Hence there is a voltage ripple in the DC-link voltage of a single-phase VSC. The
capacitance of the DC-link capacitor of a single-phase VSC has been determined to
limit the peak-peak ripple of the DC-link voltage. The detailed design procedure of
the DC link capacitor described in appendix A.1.2. Vdcavg = 400 V, Vdc =0.05 * 400
=20 V. Selected Cdc =2.5 mF.
40
Figure 3-12: shows Carrier and Reference Waveform for SPWM with Unipolar voltage
switching [25].
power. This method is used to synchronize the inverter output current with the grid
voltage, to obtain a unity power factor. It has been developed to create a current
reference that consists of active and reactive components.
The design of the control system for the inverter has been divided into current
controller, DC voltage controller and grid synchronization. A schematic block diagram
of the controller is shown in figure 3-14. Current controller control AC current injected
into the grid and voltage controller control DC voltage. The single-phase feedback
current loop is used to regulate the grid current. The current controller loop and the
model as shown in figure 3-15.
Stationary reference frame component alpha has been created by delaying grid
voltage/current component using integer delay block in Simulink. And orthogonal
component beta has been created by further delaying alpha component. Grid angle
is calculated as shown in figure 3-18. Grid angle has been used in park transformation
41
Figure 3-14: shows Schematic VSI including voltage and current control loop [43].
for calculation of DQ rotating reference frame component [43]. Figure 3-17 shows the
two orthogonal stationary component waveforms formed in the single-phase system.
It has been observed that V is e↵ective after 0.025 sec which represents a quarter-
cycle delay from V↵ .
DC-link film capacitor producing significant double line frequency ripple. Low
pass filter with 100 Hz cut o↵ frequency has been used to filter Vdcmsr . LPF min-
imises ripple. Error between measured and reference DC linked voltage has been pass
through PI to get Id reference current. Error between Idref and Idact has been again
passed through PI to get Vdref d-axis voltage. Vdqref dq-axis voltage has been further
used as feed forwarder. After final addition, Vdq(total) passed through Inverse Park
42
Figure 3-15: shows Current control loop [43].
transformation block i.e. DQ-↵ block. V ↵ is real part of this voltage has been used
for Unipolar SPWM. The feed-forward signal employed while generating SPWM to
minimize drawback of slow dynamic response of cascade control.
In thesis PI with pole zero cancellation has been designed in a DQ synchronous
frame. It can achieve zero steady state error at the fundamental frequency. It also im-
proves its dynamic response. But this method is not readily applicable to single-phase
power converters because there is only one phase variable. AB to DQ transformation
needs at least two orthogonal variables which do not have in case of single phase
power converter. Imaginary orthogonal circuit concept has been introduced to get
other variable. The imaginary Orthogonal Circuit has the same circuit components
and parameters e.g. power switches, inductors and capacitors. Imaginary circuit vari-
able like inductor current and capacitor voltage maintain 90 phase shift with respect
to their counterparts in Real circuit.
1
So, for simulation, 4
cycle delay method has been used. Stationary reference frame
have been constructed using delay method explain in figure 3-19. And then this alpha
and beta has been used as real and imaginary signal respectively to convert it into
43
Figure 3-17: shows MATLAB/SIMULINK results for real and imaginary phase wave-
forms [36].
complex variable. This complex signal has been used to get grid angle and magnitude.
This grid angle has been used further for alpha-beta stationary reference signal to
rotating DQ reference signal. Grid angle has also been used in Park transformation
and Inverse Park tranaformation. At the fundamental frequency, controller can be
designed and easily implemented in the rotating DQ frame to achieve a theoretically
infinite control loop gain.
44
Figure 3-19: shows Construction of the imaginary Orthogonal Circuit Variables [41].
Imaginary Orthogonal Circuit does not physically exist; orthogonal Circuit variables
can be constructed by using the Real Circuit variables. Imaginary Orthogonal circuit
variable and the real Circuit variables has fixed 90 phase shift. Imaginary signal,
delay by quarter cycle using integer delay block in Matlab Simulink as shown in
figure 3-19. Real & Imaginary circuit variable is denoted by alpha & beta component
respectively.
Figure 3-21 shows that Imaginary Orthogonal Circuit is estimated from the Real
Circuit and it introduces the dynamics of a quarter cycle delay in this construction
process. In result shows that delay does not a↵ect the control accuracy in steady
state. It has guaranteed zero steady state error. As compared to conventional ap-
proach, dynamic response of the control is still much faster. In case of unknown or
variable fundamental frequency application, phase lock loop can be used to track the
line frequency. In case of Imaginary Orthogonal Circuit, PLL dynamics need to be
included.
DQ rotating frame for controller design is like that of DC-DC converters. Imag-
inary Circuit variables are build using the Real Circuit variables. Complex real and
imaginary block has been used in Matlab Simulink. Then complex to magnitude-
angle block, grid angle has been calculated which used in DQ transformation block.
Then, the variables in the RI stationary frame are transformed into the DQ rotating
frame. Then voltage and/or current signal designed in the rotating DQ frame has
45
Figure 3-20: shows Real Circuit and its Imaginary Orthogonal Circuit [41].
been converted to real and imaginary stationary frame. For this conversion, inverse
transformation matrix has been used. Imaginary frame signal has been discarded.
Real signal as reference signal has been used for Unipolar modulation scheme. Trian-
gular signal has been used as carrier signal. Frequency of triangular signal is around
10 Khz.
46
Figure 3-21: shows DQ transformation using constructed Imaginary Orthogonal Cir-
cuit Variable for fixed fundamental frequency application [41].
47
48
Chapter 4
This initial set up is used to obtain I-V, P-V curve. Then same model has been
used to obtain I-V and P-V curve by keeping irradiance or temperature constant. Ef-
fect on curve using both conditions has been observed. P-V model has been attached
to boost as well as Single phase full bridge inverter. Best value of Rs and Rp has
49
Figure 4-2: shows I-V characteristic model [9].
been calculated to get maximum power point of PV. Model of photovoltaic array has
been implemented given in figure 4-2. Matlab Simulink model uses control current
source. Initially specified value of irradiance and temperature has been used. Six
panel in series has been used and set of six panel connected in parallel. So total 12
panel has been used. This has internal series cell around 72 per panel. In this number
of series and parallel set are denoted by ns and nm. Current I used for controlled
current source. Controlled current source using best value of Rs nad Rp for boost
converter has been used. Rs and Rp is the intrinsic series and shunt resistance re-
spectively. Value for Rs is very small and Rp is of very high value [3]. Irradiance and
temperature has been used as input for solar array to get output voltage and current.
50
operation is limited due to core loss. The PWM control signal for the dc converter is
used to get required output. Initially for open loop, pulse generator has been used to
control duty cycle. And in case of close loop boost converter, MPPT algorithm has
been used. Figure 4-3 shows PV array boost converter with MPPT algorithm.
Figure 4-3: shows Simulink model of the photovoltaic system with MPPT controller
& boost converter [31].
By tracking maximum power, efficiency for PV system has been increases at dif-
ferent environmental condition. MPPT has been used to extract the maximum power
from the solar PV module and transfer that to the load. That will achieve by con-
trolling duty cycle. In this case load impedance seen by source is varied which match
point of peak power with source. MPPT is not a mechanical tracking system that
physically moves panel. It is a fully electronic system that varies the electrical oper-
ating point of the modules to deliver maximum available power.
51
Figure 4-4: shows Simulink model of the Full Bridge VSC [31].
to the amplitude of a sine wave and is evaluated at the centre of the same pulse.
The frequency of reference signal determines the inverter output frequency and its
peak amplitude reference signal. Then it controls the modulation index m and Vrms
output voltage Vo . Carrier frequency control number of pulses per half cycle. In
PWM switching technique, DC input voltage is usually constant in magnitude. And
produced output AC voltage where magnitude and frequency can be controlled. Fre-
quency for trianguler wave is same as switching frequency and is kept constant.
Switching frequency of triangular wave controls speed at which the inverter switches
are turned o↵ and on. Modulation ratio ma is used to control duty cycle of one of
inverter switches. Duty cycle control based on comparison of reference (Vcontrol ) and
triangular carrier signal (Vtri ). The various frequency triangular carriers with di↵er-
ent amplitude modulation ratio SPWM signal had been programmed and tested in
single phase inverter. The two switches are never o↵ at the same time which results
to avoid output voltage fluctuation between ± V2dc .
Because of Unipolar switching, output voltage level changes between either 0 to -
Vdc or from 0 to +Vdc , this scheme has the e↵ect of doubling the switching frequency as
far as the output harmonics are concerned, compared to the bipolar switching scheme.
To generate an AC waveform in single-phase inverter, the switches S1, S2 ON or S3,
S4 on for respective period to get alternate positive as well as negative voltage. Two
52
switches in same leg should not be turn on at same time. Final sinusoidal fundamental
component as shown in figure 4-5. The load voltage is calculated by
cos(n↵)
Vn = 4 ⇤ Vdc (4.1)
n⇤⇡
n=1,2,3....
Where
↵: The blanking time along with the scale of each harmonic rests on it.
The angle ↵ decides magnitude of the harmonic. Harmonic has been reduced by
PWM technique. High frequency signal has been used to reduce number of filter.
The resulting output voltage include fundamental alternating component. Voltage
includes high frequency harmonics. The harmonics can be filtered with LCL filter
or by implementing appropriate modulation techniques. The instantaneous inductive
load current is given in 4.2 [35].
Figure 4-6: shows Simulink model Single Phase Two stage Full Bridge VSC [31].
54
Chapter 5
Simulation Result
Figure 5-1: I-V characteristic of a solar array for a fixed temperature but varying
irradiance [31]
The I-V & P-V curves shown in figure 5-1 & 5-2 respectively for various irradiance
but a fixed temperature (25 ). It has been observed that the open circuit voltage
decreases with increase in irradiance at a fixed temperature. Moreover, from the I-V
& P-V curves shown in figure 5-3 & 5-4 respectively at a fixed irradiance at 1000
W/m2 , it has been observed that the open circuit voltage decreases with increase in
temperature. From the I-V curve shown in figure 5-1, short circuit current has been
increases with increase in irradiance at a fixed temperature.
55
Figure 5-2: P-V characteristic of a solar array for a fixed temperature but varying
irradiance [31]
Figure 5-3: I-V Characteristic of a PV array under a fixed irradiance but varying
Temperatures [31]
The simulation is run at 0.6s. At the beginning, the irradiation is set at G=1000
W/m2 and at t=0.3 s a step change of irradiation to 600 W/m2 is performed. The
output power of the PV varies from 3.947 kW to 2.383 kW as shown in figure 5-5.
Though there is a variation of the irradiance, PV array operates at maximum power.
Figure 5-6 shows step change in output voltage of PV array with respect to change
in irradiance. The output voltage oscillates around the maximum power due to P&O
algorithm. Figure 5-7 shows step change in output current of PV array with respect
to change in irradiance. The PV array current Ipv reaches steady state after 0.002 s
shown in figure 5-8. Moreover, the current oscillates around the MPP. The current
56
Figure 5-4: P-V Characteristic of a PV array under a fixed irradiance but varying
Temperatures [31]
Figure 5-5: shows Step change in the Photovoltaic Power Output [31]
drops from 17.02 A to 10.41 A with respect to irradiance at 0.3s. At 600 W/m2
irradiances, current 10.41 A which is the maximum current of the PV array. Maximum
power almost track instantly by MPPT.
Figure 5-10 & 5-9 shows step change in grid voltage & current respectively of
Single Phase Full Bridge Inverter with respect to change in irradiance. Figure 5-11 &
5-12 shows change in DC link voltage & Active-Reactive power respectively of Single
Phase Full Bridge Inverter with respect to change in irradiance.
57
Figure 5-6: shows Step change in the Photovoltaic PV array Output Voltage [31]
Figure 5-7: shows Step change in the Photovoltaic PV array Output Current [31]
Figure 5-8: shows PV array current Ipv reaches steady state after 0.002 s [31]
58
Figure 5-9: shows Step change in the Single Phase Full Bridge Inverter Grid Current
[31]
Figure 5-10: shows Step change in the Single Phase Full Bridge Inverter Grid Voltage
[31]
Figure 5-11: shows change in the Single Phase Full Bridge Inverter DC link Voltage
[31]
59
Figure 5-12: shows change in the Single Phase Full Bridge Inverter for Active-Reactive
Power [31]
60
Chapter 6
6.1 Conclusion
With help of PV & IV curve shown in simulation, characteristics for voltage, current,
power has been verified and is the same as the characteristics given from the data
sheet. PV models output voltage & current vary with respect to change in irradiance
or temperature. Maximum power has been achieved irrespective of operation condi-
tion with help of Perturb & observe algorithm which can track maximum power point
of the PV array.
Whole single phase two stage full bridge grid connected DC-DC-AC inverter with
maximum point has been analysed and simulated. Finally, DC voltage generated by
PV array could produce AC sinusoidal current at output of inverter. Finally, PV
power decides, amplitude of current. SPWM (sinusoidal pulse width modulation)
unipolar technique and designed LCL filter has been used to get grid voltage and
current with less harmonics and distortion. To control DC link voltage and current,
current and voltage control loop has been used.
61
system should be done. To keep the boost converter output voltage constant, Perturb
and observation MPPT controller has been used. The fast-inner current loop with a
slower outer voltage loop are designed with help of LCL filter to eliminate harmonic.
For further improvement, this method can be used.
⇤ Single phase grid connected system can be used as a basic guideline when de-
signing the grid connected PV inverters and their controls 50Hz and 240V point of
common coupling (PCC) voltage.
⇤ With better control on active and reactive power, research on power management
can also be conducted for future work.
⇤ More advance voltage control techniques and operating and coordination strate-
gies developed with help of Thesis can be evaluated by carrying out simulation stud-
ies on PV systems integrated to three-phase four wire systems and to Distributed
resources (DR) for further development.
⇤ When PV system is integrated with a battery storage system, power curtailing
may not be necessary in certain situations. In future, it is beneficial to integrate
developed PV system with a battery storage systems.
62
Appendix A
The inductor determines the stability of the current sensor loop and the inductor
current ripple. The peak-to-peak inductor ripple current is inversely proportional
to the inductor value. There are two ways to calculate current ripple. To calculate
inductor, IL has been assumed to 5 % of the rated input current of the DC-DC
boost converter.
63
A.1.2 Selection of Input Capacitor
Vpv D
Cpv = 2
(A.2)
4 ⇤ Vpv fdc Ldc
215
D=1 400
= 0.4625 Duty cycle of the DC-DC boost converter,
Vpv = 215 V Average voltage across the PV array,
Vpv =20 V Assume Switching ripple voltage of the PV array,
Ldc =19 ⇤ 10 3 H
fdc =10 kHz.
Calculated Cpv =0.621 µF implemented in PV system in order to limit Vpv .
Selected Best Value Cpv =0.584 µF
Equation A.3 has been used for selection of boost converter output capacitor. It de-
pends on minimum ripple voltage. This output capacitor acts as DC link capacitor.
Since,
Pin = 215 ⇤ 17.02 = 3659.3W = Pout
Pout 3659.3
Iout = Vout
= 400
= 9.14825 A
Vload 400
Rload = Iout
= 9.14825
= 44 Ohm
64
where Vload : output voltage of the boost converter
Vload : output ripple voltage
Due to variation in temperature and irradiance, output voltage of PV array also
varies. DC link capacitor is installed between the PV and the inverter to compensate
the variation of the output voltage of the PV. Proper selection of DC link capacitor
reduces voltage ripple and o↵er energy storage for short period of time.
65
66
Appendix B
In LCL filter includes two inductors and capacitor connected in parallel i.e. per phase.
In above LCL filter Rf and Lg are grid side resistor and Inductor respectively. And
Ri and Li is inverter side inductor. In this Cf is parallel capacitor and Rd is resistor
for passive damping [15].
The voltage at the output of the VSC, the voltage across Cf and Rd and the grid
voltage are v, vc and vg respectively. The output current of the VSC and the current
injected to the grid by the VSC are i and ig respectively. The grid voltage appears as
a short circuit for all the harmonic voltages created by the VSC as grid voltage. It is
assumed as a pure sinusoid at fundamental frequency, 0 f 0 . The harmonic number is
denoted by “\nh00 .
(ig (s)) Rd ⇤ C f s + 1
GLCL (s) = = (B.1)
v(s) As + Bs2 + Cs + D
3
Where A = Lf c ⇤ Lf g ⇤ Cf
B= (Lf c Rf g + Lf g Rf c + (Lf c + Lf g ) ⇤ Rd ) ⇤ Cf
C = (Rf g Rf c + Rf c Rd + Rf g Rd )Cf + Lf c + Lf g
D = Rf c + Rf g
67
After neglecting resistance of each inductor,
(ig (s))
GLCL (s) = v(s)
Rd ⇤ C f s + 1
GLCL (s) = (B.2)
s(Lf c ⇤ Lf g ⇤ Cf s2 + (Lf c + Lf g ) ⇤ Rd ⇤ Cf s + Lf c + Lf g )
ig (s) 1
GLCL (s) = = (B.3)
v(s) s(Lf c ⇤ Lf g ⇤ Cf s2 + Lf c + Lf g )
ig (s) 1
= (B.4)
v(s) s(Lf c ⇤ Lf g ⇤ Cf s2 + Lf c + Lf g )
Where
2 1
ZLC = Lf g ⇤Cf
and
2 Lf c +Lf g
!res = Lf c ⇤Lf g ⇤Cf
LCL filter has been designed to reduce switching ripple in grid current. In below
procedure, relationship between the converter side current and voltage has been de-
rived. There are two ways to design LCL filter any frequency. In first, assume filter
capacitance should be with very small impedance. Capacitor acts as short circuit as
frequency is higher than switching frequency of VSC, fsw .
i(s) 1
= (B.5)
v(s) sLf c
ig (s) Z2
= 2 LC 2 (B.6)
i(s) S + !res
In other way, it has been derived in equation B.6. Grid voltage considered as short
68
Table B.1: Specifications of the VSC
Parameter Value
Pr 5.4 KW
Vr 230
Vdc 400
f 50 Hz
fsw 25 KHz
circuit for harmonics. So, frequency is higher than switching frequency of VSC, fsw .
ig (s) Z2
= 2 LC 2 (B.7)
i(s) S + ZLC
It has been proved while analysing the performance of a designed LCL filter.
Defining base quantities for the VSC system.
Vr2
Zb = (B.8)
Pr2
1
Cb = (B.9)
! ⇤ Zb
Detail design procedure & constraints of LCL Filter explained in reference [24]
[23] [13] [2] [30].
Base values can be calculated as Zb = 9.8 and Cb = 325 F by using B.8 and B.9
respectively. The peak switching ripple current, ip through Lf c to a given value can
be obtained from B.10 [14] [10] to limit Lf c inductance. VSC has been modulated
with sinusoidal unipolar PWM and with the modulation index, ‘m’ is 0.5 < m < 1.
Then only equation B.10 will be valid.
Vdc
Lf c = (B.10)
16 ⇤ fsw ip
Lf c can be calculated as 300 µH using B.10. Switching ripple current attenuated
to 10 % of the rated peak current of the VSC i.e. 10%*33=3.3 A.
69
Table B.2: LCL designed Filter Parameters
Parameter Value
Lf c 3 mH
Lf g 1.5 mH
Cf 2.6µF
Rd 1.5 ohm
f 50 Hz
fsw 25 KHz
70
Appendix C
P
e(t) = ⇤ Sin(2!t) (C.1)
2⇤!
The di↵erence in the maximum and the minimum energy stored in Cdc , emax that
is the energy di↵erence when (2!t)= ⇡2 and (2!t)= - ⇡2 is given by
P
emax = (C.2)
!
The maximum and the minimum energy stored in Cdc causes the maximum and the
minimum voltages at Cdc since the voltage across a capacitor is proportional to the
stored energy in the capacitor. Therefore emax can be written in terms of the
maximum and the minimum DC-link voltages that are Vdcmin and Vdcmax respectively
as given in
1 2 2
emax = ⇤ Cdc ⇤ (Vdcmax Vdcmin ) (C.3)
2
An expression for the peak-peak 100 Hz DC-link voltage ripple, Vdc , can be derived
as given in (C.6) from (C.2) and (C.3). Equation (C.5), is the average DC-link voltage
of the VSC.
Vdc = Vdcmax Vdcmin (C.4)
(Vdcmax + Vdcmin )
Vdcavg = (C.5)
2
71
P
Vdc = (C.6)
(Cdc ⇤ Vdc ⇤ !)
The capacitor at the DC-link should be large to reduce the DC-link voltage ripple
almost to zero as per (C.6). Practically it is not possible and there is a limitation
to the size of Cdc . But voltage ripple at DClink voltage of a single-phase VSC has
been controlled. The Vdc is 5 % of Vdcavg i.e. 20 V. With reference to value given
in table B.1, Cdc calculated 2.2 mF. Cdc selected best value 2.5 mF. Detailed design
procedure for DC-link capacitor has been explained in [30] [42] [16] [20] [33].
72
Appendix D
!"#$%&'"!%()*+,*-*./0*+,
"#$%&'(!('&)*+(!,'-'.*/(0)&/-'(
1+*/2/.)$'!+-!+&!,%-$'!*+)&!3!4.*+$/0)$'!)!+(.)&)!,%-$/)&
!""#$%!&%#' !"#$%&'(%)*+*,-./
!"#$%&'()*+%'(,"-.+/0(1'-$%+%"2#("-/(3/$#($40+./(5/&/(/%#-.#&(0/$
-#1#$+./.#$(/(0/&6'(50/7'(.#(0/(+-."$%&+/($'0/&8
9'-(,:$(.#(;<<(=>(+-$%/0/.'$(#-(%'.'(#0(,"-.'?(@#,'$(1&#/.'
"-/(,/&1/(&#$5#%/./(/0('A#&(.#(,/-#&/(1'-%+-"/./("-(5&'."1%'
.#(1/0+./.(1'-(&#-.+,+#-%'(1'-$%/%/.'8
!"#$%&/(.#1+.+./(/5"#$%/(5'&(0/(BCDC+(-'$(5#&,+%#("-/(&#."11+4-
.#(1'$%#$(2(.#(,#E'&/(.#(0/(#A+1+#-1+/(5#&,/-#-%#$8(F$+,+$,'?(-'$
!"# $
&'(
%%
%%%
!"#$%&,(*-.&
5#&,+%#(,#E'&/&(0'$(5&'."1%'$(2/(#*+$%#-%#$(2(.#$/&&'00/&(-"#G/$
%#1-'0'6H/$8
!"#$%&'()*+&
!("#$%&''()#*%&+(),('-
B BB B BB B BB B
!(./-0(1('-'()#2*"3.*45657
!"#$%&'()*(+
!(8&$9*,&*$&:/+(,-,2*"3.*45;<=
BB B B
BBBB
0 4
1 22
222334
5
/
!"#$%&'
()#*$+$,-.$
&,#?(2,->*1
"*-,1"
01$1.231%)*+*,-.
^'$( ,4."0'$( _+0+\#-( $'-( 1#&%+A+1/.'$( .#( /1"#&.'( 1'-( 0/$( -'&,/$
+-%#&-/1+'-/0#$(`^(#(Ba98(^/(/0%/(1/0+./.(2(A+/3+0+./.(/$#6"&/-("-(&:5+.'
4%567-+-.2-%2#+-$1.7*1%8-%9#2-.7*1%:;<=>
&#%'&-'(.#(0/(+-G#&$+4-(#-(1"/0b"+#&(+-$%/0/1+4-(A'%'G'0%/+1/?(2/($#/
4%?=%1@#A%8-%B1$1.231%8-%9$#8'72#
&#$+.#-1+/0?(+-."$%&+/0('(5/&b"#($'0/&8
4%CD%1@#A%8-%B1$1.231%8-%9#2-.7*1%+*.-1+
9/./(,4."0'($#($',#%#(/($"$(A/$#$(.#(#-$/2'?(+-10"+./("-/(&+6"&'$/
4%567-97*#.1+%7#E9#$21E*-.2#%7#.%F1G1
G#&+A+1/1+4-(.#(0/(5&'."11+4-(.#(#-#&6H/(2(0'$(#-$/2'$(.#(&#$+$%#-1+/
%%%+'E*.#A*818
#-($+,"0/.'&#$(.#(&/.+/1+4-8(F0(1'-%/&(1'-(5&'."1%'$(2($'0"1+'-#$
4%H-$2*I*717*J.%KLM%NOP%-%Q.2-$2-,%91$1%19+*717*#.-A
+--'G/.'&#$(+-$%/0/.'$($/%+$A/1%'&+/,#-%#(#-(0'$(a$%/.'$(`-+.'$(2(#-
%%-.%2#8#%-+%E'.8#
%'.'(#0(,"-.'?(_+0+\#-(#$("-('-'1+.'(0H.#&(,"-.+/0(#-(#-#&6H/($'0/&8
4%RJ8'+#A%7#.%'.1%-I*7*-.7*1%8-%S1A21%-+%?T>
4%U-A9-2'#A#%7#.%+1A%-E*A*#.-A%8-%71$F#.#
!"#$%&'"(#)*"*+,-#,.#,"#-/$,0(#'-(#&'01-2,#,"#34546
7,-,01# '-# 89:;# $<.# &,# ,-,07=1# >(-# 0,.?,>2(# 1"
?0($,&*(#&,"#>(-@'-2(#&,#"(.#$%&'"(.#,.2'&*1&(.#A
'-#539B;#>(-#0,.?,>2(#1"#$=-*$(C
!"#"$%&'!()((I(J'-./(B$//1(K#&/0(2(9/3/00#&'?(LM(I(MNOP<(K/%#&-/(Q(R/0#-1+/(Q(*+,-.-'I(S#08T(UC;MV(O<W(ML(WW(;;(I(X/*T(UC;MV(ON(<Y<(OW(NZ(I(+-A'[$+0+\#-81',(I(]]]8$+0+\#-81',
73
!"#$%&'"!%()*+,*-*./0*+,
"#$%&'(!('&)*+(!,'-'.*/(0)&/-'(
%"&#
'#(#)$
0"/")#(/1%#,)"%&)$+%#/2)#,3"%
=2I3*1*PI3
B)(%"7$)"*&7%#&*$&7%#8(*5&*"[:N(*5&*"7$<"N&
B)0&7()87&(*5&*<8(
Q*%"<"5#8(*5&*"7$<"N&
/(&0123 /
%(456153 %
*
""#
"$#
*(7(,((840
=2S? 9
9(2:";"0!
,
!"N"*5&*$87&1)87&(*+/*S3
'(<3"53(=10!813
'
E8<&#"7$)"*f2*00
b7)5"5&(*&7*00
_%)p:&%"*)5&7%)6)$"$)H7 E+_GGJ
K58I (":&/Z/[T/WX(DY
<EFG48/H8I8J84IK378/L2IK6F=M
@I=6F=MF/BP8I8J84I=37/<8FG4,
!"#$%&'#(%)*+%,-./&01%1210(-#0#(+%3"1'
14,.$1/%(.%2#5"(6%7.81-%*0(#81%$&-9*01
L,8I *:'*/]/WX(DY
PU>/Z.79./[/Z.79:4/P30@.7
!""#$%#&'!
!""#$%#&'!
M57/M9.42/85;;.8-95;,)/3,./09;9030/P5>/%#/QRS/8566.7/T97.,/9;,34:-.2/G57
:19.-1%"*'/2#'5%;&'0(#.'%:.46%<"1'
)./&21$=%"*>*-/.&$%8.2(*51%)*+:1
09;9030/5G/"#U?/V%"'UMW>/?4:,,/?/G97./7:-9;D)/:,/2.-.709;.2/:G-.7/XE/G4:0.
,9,-I !$$$/Z/[T
"*'/2#'5%("1%,".(.8.2(*#0%)./&216
'!!$&($!)*
+4.8-798:4/B:43.,/3;2.7/.,-A;2:7/85;29-95;,)/:885729;D/L+?/&#"#'
,-./&01/6%?1*/%0.),21(12+%*'/
0*-19&22+%("1%&$1-@$%5&#/1%:19.-1
K:457.,/.4N8-7985,/@:C5/85;29895;.,/.,-A;2:7)/,.DO;/L+?/&#"#'
MMI \C:'*D
0.''10(1/%#'%$1-#1$%3#("%.("1-
!"#$%$&'
!*
K0:<
&'1&#/+Q>I':/\\aH/=6F/]<!1C
!$$$/+Q9':/'CaT/b/]<!1C
$(20>0(68(23?8@53?8!
+,-./000123435,61789:/,;<=34./36>8?23435,61789
-.,->/Y:<9030/,.79.,/G3,./7:-.2/%#/H?/Q06,>
J.,5/[/R.9D\-I !&:$$/)c2/Q/#!:*&/4-2
<8F1. @A)"$B"A/'($+/CD
'!!$&($!)*
$
2.-.;92:0.;-./:;-.,/2./0:;9634:7/.4
A7.:/:8-9B:/6:7:/0:;9634:7/4:/8:C:/2.
,.79./85;/5-75,/012345,/,./:48:;=:;
B:457.,/2./-.;,91;/6.49D75,5,>/E.:/.4
83:;25/,./.<65;./:/4:/43=>/?3@797/.4
85;.<95;.,>/?3:;25/,./85;.8-:/.;
H=-M37=F8/V8M/Q/<=6G>=7IGM,F/-b
NBM,7=G73E6O
+,-./012345/675238./.4.8-79892:2
0:;3:4/2./9;,-738895;.,/.;-.7:/F
K066I '&:C/Z/[T/WX(DY
J066I '($:$/+V/WXCDY
\:\&/]/WX(DY
H90>I !"#$%&&$%#$
ST/U)
@3435,6/<8FG4,2/@1A1^1
^_VJ;WI ;$1#(DQaT
P]?^I #"X'aT
012345/G5-5B54-A985>
BQR. $!'(*!
T_H./`;&\\##*""
L066I
$!'(*!
!8##)&7%&**VJY
/8%&7$)"*V@Y
/099A/099*(%$
-8<%"N&*V-Y -8<%"N&*V-Y
+##"5)"7$)"*V@A02Y
!"#$%&'()*+,)$% 0(/#,9,)"5$%
B)0&7()87&(*VW1J1XY =US?*1*UU?*1*Q?*00 G&[)(%#"58*&7*bW*& bW*dGB,!=R?T,?=*A*bW=R?T
/&(8 2T*Z[ +7%&#%&Z
W87[)%:5&(*5&*$"O<&*()0;%#)$"(*5&*=I2S0I*eQ*002I*58O<&*$"9"*")(<"7%&I*<)O#&*5&*]"<H[&7:#8(I !"%&[8#'"*)[7'6:[" !<"(&*!
!"O<&(*5&*("<)5" #&()(%&7%&*"*<"*#"5)"$)H7*b- !&#%)6)$"$)H7*El- +_!*S=2=3*A*+_!*S=RT?*A*S=R?=*&7(".8*7)&O<"*("<)7"
!"N"*5&*$87&1)87&( +/,S3*$87*5)858(*O.9"((**5&*9#8%&$$)H7 +Dd*U??=F2??? cL*_D?PA3=R?
K"#$8 J<&"$)H7*5&*"<:0)7)8*"785)\"58*5&*=3*0)$#"(*5&*&(9&(8#*%)98*S?ST*ES +Dd*=Q??= cL*_D?UAS32?
-)5#)8*5&<"7% -)5#)8*%&09<"58*5&*TI2*00*$87*O"N8*$87%&7)58*5&*])&##8*$87*&<&^"5"*$"9"$)5"5*5&*%#"7(0)()H7 B&$<"#"$)H7*5&*$8768#0)5"5*!_*V0"#$"*!_Y
!;<:<"(*(8<"#&( R2*$;<:<"(*0878$#)(%"<)7"(*=3S*1*=3S*00
0$+5,),$+(%&5(&$6(/"),7+&)$'6/$8"5"%
!"#$%&(-.)#/,)$% E&09&#"%:#" ,Q?*L!*"*`P3*L!
/8%&7$)"*0C1)0"*"*DE!*V`TA?*aY /09*V@9Y 2P3 2U? 2U3 T?? T?3M T=?M !"#["*&(%C%)$" 2Q??*/"
_6)$)&7$)"*"*DE! *VaY =QIR =QIU =3I2 =3I3 =3IR =SI? !"#["*0C1)0" 3Q??*/"
X"$%8#*5&*<<&7"58 XX ?IRS= ?IR3U ?IR3U ?IR3P ?IRS ?IRST G&()(%&7$)"*"<*)09"$%8 +09"$%8*98#*[#"7)\8*e23*00*"*2T*0A(
-8<%"N&*"*98%&7$)"*0C1)0" -09*V-Y T3IR T3IR T3IR T3IP T3IP T3IU
!8##)&7%&*"*98%&7$)"*0C1)0" +09*VJY RIUP PI=2 PI23 PITP PI3= PIST 4"/"+#1"&5(&6/$52)#$
-8<%"N&*5&*$)#$:)%8*"O)&#%8 -8$*V-Y QQI2 QQIQ QQIQ QQIQ QQI3 QQI3 =?*"k8(*5&*["#"7%'"*&7*0"%&#)"<&(*.*0"78*5&*8O#"
-8<%"N&*5&*$8#%8$)#$:)%8 +($*VJY PIQR PIS= PIR3 PIU= UI?2 UI=T
-8<%"N&*0C1)08*bW*A*+_! -0"1*V-Y*bWA+_! S??*A*=??? :;&"<$%&5(&="/"+#1"&5(&6$#(+),"&-,+("-
!8&6)$)&7%&*%&09&#"%:#"*5&*/09 EZ/09*VaAL!Y ,?IQ= /#)0&#*"k8F*URa*5&*<"*98%&7$)"*&%)p:&%"5"
!8&6)$)&7%&*%&09&#"%:#"*5&*-8$ EZ-8$*VaAL!Y ,?ITS? B&<*"k8*2*"<*23F*#&5:$$)H7*0C1)0"*5&*98%&7$)"*5&*?IRa*9>">
!8&6)$)&7%&*%&09&#"%:#"*5&*+($ EZ+($*VaAL!Y `?I?S
E&09&#"%:#"*5&*$;<:<"*78#0"<*5&*89&#"$)H7 cd!E*VL!Y QRf2
X:()O<&(*5&*9#8%&$$)H7 J == !A*K"(("0"[#&<<I*=T*m*/8<>*+75*Wno8#%&%">
B)858(*j.,9"(( JA- =2AQ? QS=TP*G"6&<O:7.8<*,*-"<&7$)"*,*_(9"k"
1%-/'!2344
!GCGNEA(!O0O*g*G875"*+(""$*/&#"<*.*!"O"<<I*=Q*g*QSUP?*/"%"*,*-"<&7$)"*,*8P.HQH(g*E&<>F*V`TQY*U?2*Q=*22*TT*g*X"1F*V`TQY*US*?R?*U2*S3*g*)768h()<)Z&7>$80*g*iii>()<)Z&7>$80
74
Bibliography
[1] L[eslie] A. Aamport. The gnats and gnus document preparation system.
G-Animal’s Journal, 41(7):73+, July 1986. This is a full ARTICLE entry.
[2] C. Bao. Step-by-step controller design for lcl-type grid connected inverter
with capacitor current-feedback active-damping. IEEE Trans. Power Electron,
29(3):1239, March 2014.
[3] Habbati Bellia, Ramdani Youcef, and Moulay Fatima. A detailed modeling of
photovoltaic module using matlab. NRIAG Journal of Astronomy and Geo-
physics, 3(1):53–61, 2014.
[4] Sergio Augusto Oliveira da Silva, Leonardo Bruno Garcia Campanhol, Alessan-
dro Goedtel, Claudionor F Nascimento, and Denilson Paiao. A comparative
analysis of p-pll algorithms for single-phase utility connected systems. In Power
Electronics and Applications, 2009. EPE’09. 13th European Conference on, pages
1–10. IEEE, 2009.
[5] Debashis Das. Modeling and simulation of PV array with boost converter: An
open loop study. PhD thesis, NATIONAL INSTITUTE OF TECHNOLOGY
ROURKELA, 2011.
[6] Sairaj V Dhople, Ali Davoudi, and Patrick L Chapman. Dual-stage converter to
improve transfer efficiency and maximum power point tracking feasibility in pho-
tovoltaic energy-conversion systems. In Applied Power Electronics Conference
and Exposition (APEC), 2010 Twenty-Fifth Annual IEEE, pages 2138–2142.
IEEE, 2010.
[7] JHR Enslin and PJM Heskes. Harmonic interaction between a large number of
distributed power inverters and the distribution network. In Power Electronics
Specialist Conference, 2003. PESC’03. 2003 IEEE 34th Annual, volume 4, pages
1742–1747. IEEE, 2003.
[8] Trishan Esram and Patrick L Chapman. Comparison of photovoltaic array maxi-
mum power point tracking techniques. IEEE Transactions on energy conversion,
22(2):439–449, 2007.
75
[9] Cristina Gonzalez-Moran, P Arboleya, G Diaz, and J Gomez-Aleixandre. Mod-
eling photovoltaic dc primary sources as grid connected inverter supplies consid-
ering non linear e↵ects. In Electrical Power Conference, 2007. EPC 2007. IEEE
Canada, pages 50–55. IEEE, 2007.
[10] D Grahame Holmes and Thomas A Lipo. Pulse width modulation for power
converters: principles and practice, volume 18. John Wiley & Sons, 2003.
[15] Hea-Gwang Jeong, Kyo-Beum Lee, Sewan Choi, and Woojin Choi. Performance
improvement of lcl-filter-based grid-connected inverters using pqr power trans-
formation. IEEE Transactions on Power Electronics, 25(5):1320–1330, 2010.
[16] CENGİZ KAHRAMAN, BAŞAR ÖZTAYŞI, and SEZI CEVIK ONAR. Photo-
voltaics type selection using a projection model-based approach to intuitionistic
fuzzy multicriteria decision making. In Uncertainty Modelling in Knowledge En-
gineering and Decision Making: Proceedings of the 12th International FLINS
Conference, pages 924–929. World Scientific, 2016.
[17] Tamás Kerekes, Remus Teodorescu, Pedro Rodrı́guez, Gerardo Vázquez, and
Emiliano Aldabas. A new high-efficiency single-phase transformerless pv inverter
topology. IEEE Transactions on Industrial Electronics, 58(1):184–191, 2011.
[18] David L King, Jay A Kratochvil, and William Earl Boyson. Photovoltaic array
performance model. United States. Department of Energy, 2004.
[19] Soeren Baekhoej Kjaer, John K Pedersen, and Frede Blaabjerg. A review of
single-phase grid-connected inverters for photovoltaic modules. IEEE transac-
tions on industry applications, 41(5):1292–1306, 2005.
76
[21] Marco Liserre, Antonio Dell’Aquila, and Frede Blaabjerg. Stability improve-
ments of an lcl-filter based three-phase active rectifier. In Power Electronics
Specialists Conference, 2002. pesc 02. 2002 IEEE 33rd Annual, volume 3, pages
1195–1201. IEEE, 2002.
[22] Fangrui Liu, Yong Kang, Yu Zhang, and Shanxu Duan. Comparison of p&o
and hill climbing mppt methods for grid-connected pv converter. In Industrial
Electronics and Applications, 2008. ICIEA 2008. 3rd IEEE Conference on, pages
804–807. IEEE, 2008.
[23] W. Xiao M. Hanif, V. Khadkikar and J. Kirtley. Two degrees of freedom active
damping technique for an lcl filter based grid connected pv systems. IEEE Trans.
Ind. Electron., 16(6):2795, July 2013.
[24] F. Blaabjerg M. Liserre and Ste↵an Hansen. Design and control of an lcl-filter
based three-phase active rectifier. Int. J. Electron.
[25] Sachin Maheshri and Prabodh Khampariya. Simulation of single phase spwm
(unipolar) inverter. International Journal of Innovative Research in Advanced
Engineering (IJIRAE), 1(3), 2014.
[26] Yousef A Mahmoud, Weidong Xiao, and Hatem H Zeineldin. A parameterization
approach for enhancing pv model accuracy. IEEE Transactions on Industrial
Electronics, 60(12):5708–5716, 2013.
[27] K Manohar and P Sobha Rani. Mppt and simulation for a grid-connected pho-
tovoltaic system and fault analysis. The International Journal of Engineering
And Science (IJES), 1(2):158–166, 2012.
[28] Ned Mohan and Tore M Undeland. Power electronics: converters, applications,
and design. John Wiley & Sons, 2007.
[29] Anuja Namboodiri and Harshal S Wani. Unipolar and bipolar pwm inverter.
International Journal for Innovative Research in Science & Technology, 1(7):237–
243, 2014.
[30] Brian K Perera. Modelling of inverter interfaced renewable energy resources to
investigate grid interactions. 2015.
[31] Falinirina F Rakotomananandro. Study of photovoltaic system. PhD thesis, The
Ohio State University, 2011.
[32] Dr Zainal Salam. Power electronics and drives (version 2). Pdfheart. com, 2002.
[33] Fritz Schimpf and Lars Norum. Parallel power decoupling for single-phase pv-
inverters–making film capacitors applicable.
[34] Dezso Sera, Remus Teodorescu, and Pedro Rodriguez. Pv panel model based on
datasheet values. In Industrial Electronics, 2007. ISIE 2007. IEEE International
Symposium on, pages 2392–2396. IEEE, 2007.
77
[35] R Sha↵er. Power electronic harmonic nalysis and fundamentals of power elec-
tronics with matlab, power quality in electrical systems. Da vinci engineering
series, Thomson learnin Inc, USA, 2007.
[38] http://www.coloradosolarpower.net/images/brochures/p100094/SLK60P6L
[39] https://autosolar.es/pdf/Ingeteam-Ingecon-Sun-Lite.pdf.
[41] Richard Zhang, Mark Cardinal, Paul Szczesny, and Mark Dame. A grid simula-
tor with control of single-phase power converters in dq rotating frame. In Power
Electronics Specialists Conference, 2002. pesc 02. 2002 IEEE 33rd Annual, vol-
ume 3, pages 1431–1436. IEEE, 2002.
[42] Xiangdong Zong. A single phase grid connected DC/AC inverter with reactive
power control for residential PV application. PhD thesis, University of Toronto,
2011.
[43] Xiangdong Zong and Peter W Lehn. Reactive power control of single phase grid
tied voltage sourced inverters for residential pv application. In IECON 2012-
38th Annual Conference on IEEE Industrial Electronics Society, pages 696–701.
IEEE, 2012.
78