Comparison Between Classical P&O' Algorithm and FLC of MPPT For GPV Under Partial Shading
Comparison Between Classical P&O' Algorithm and FLC of MPPT For GPV Under Partial Shading
Corresponding Author:
Youcef Abdelaziz
L.G.E, Laboratory, Departement of Electrical Engineering
Dr Moulay Tahar University of Saida, Algeria
Email: [email protected]
1. INTRODUCTION
The use of renawble energy has become a major necessity for the future of our planet. The MPPT
becomes complex when the photovoltaic generator is in the presence partial shading. The research was
directed towards optimizing the performance of a solar photovoltaic system by extracting its maximum
power for all atmospheric conditions. The GPV is formed by a series/parallel combination of photovoltaic
modules to converge on the desired GMP. Several MPPT algorithms have been proposed by researchers [1]-
[9], the most used are: Perturb and Observe, Incremental Conductance, Constant Voltage. These methods are
based on finding the maximum of the carateristic P-V of the PV solar system and are effective only when the
curve has only one maximum. When passing cloudy or obstacle on a GPV, multiple points appear on the
power-tension characteristic with an overall GMP maximum and LMP points. Conventional methods do not
converge on the desired GMP [10]-[13].
Figure 1. We found that this approach to the 'P&O' algorithm leads to an error in tracking the global power
point when it changes position from Pm3 to Pm2 for values under 200 W/m² in figure 2. The special
categorization and terminology given by [26] is used for this purpose. Figure 1 shows a GPV composed of
100 modules in a mixed group (10x10).
2
Pm1
Local maximum
points
0
0 100 200
Voltage across PV array [V]
Figure 1. GPV (10x10) under parial shading divided Figure 2. P-V curve of GPV under partial shading of
into three groups Figure 1
V nV
+ N s Rs I
V+N RI
I = I sc − I 0 e t N s − 1 − s s
(1)
N R
s p
I sc _ std V nV
+ N s Rs I
V +N R I
f (V , I , G) = I −
+ I 0 e t N s − 1 + s s
(2)
G N s Rp
3. RESULTS AND DISCUSSION
This simulation was made for the configuration of Figure 1. The non-shaded modules receive
irradiation level of G=1kW/m². The shaded modules are under G=0.1kW/m². The characteristic P-V
illustrates the presence of three peaks with a single GMP and two local maximums in Figure 2. Figure 3
brings together all the P-V curves for different irradiations changes in steps of 0.1kW/m² and a constant
temperature T=25°C of the Figure 1. From the simulation of Figure 3, the power increases as the irradiation
increases. We also note the change in the position of the GMP from Pm3 to Pm2 for low irradiation values.
Table 1 shows the effect of partial shading on the ‘P&O’ MPPT algorithm for a grid of values of G. It gives
the ‘P&O’ tracked power and the maximum available power.
9
G= 1.0
Power output from PV array [KW]
8
G=0.9
G= 0.8
G= 0.7
6 G= 0.6
G= 0.5
G= 0.4
4
2
G= 0.1 G= 0.3 G= 0.2
0
0 100 200 250
Voltage across PV array [V]
Figure 3. GPV under partial shading given by Figure 1 for different shading irradiation levels
Int J Pow Elec & Dri Syst, Vol. 12, No. 2, June 2021 : 1131 – 1138
Int J Pow Elec & Dri Syst ISSN: 2088-8694 1133
Figure 3 shows that for a sudden change in irradiation for shaded modules of 1000 to 50 and 100
W/m², the classical "P&O" algorithm does not have the intelligence to converge the GMP. In this case the
power losses are 08,24% and 08.99%. Figure 4 shows that for a sudden change in the shading irradiation
level (from 1000 to 100 W/m²), the MPP tracked value using the P&O algorithm converges to 2886 W and
not towards the global MPP who is 3208 W. The idea is to optimize the maximum power at the exit of the
GPV, in this case GMP by inserting the FLC control whose gains are regulated by PSO in order to optimize
the maximum output power of the GPV.
(W/m²) (W) Power (W) (%) P-V curve of the partially Pm3=2886 W
shaded PV array
50 3099 2587 08,24 4
100 3208 2886 08,99
200 3488 3488 0
300 4093 4093 0 2
Pm2=3208 W
400 4698 4698 0
500 5300 5300 0
600 5896 5896 0 0
0 40 80 120 160 200
700 6480 6480 0 Voltage across PV array [V]
800 7045 7045 0
900 7570 7570 0
1000 7997 7997 0 Figure 4. Classical ‘P&O’ tracking process from a
uniform standard irradiation to the partial shading
configuration of figure 1 with a shading levelof
0.1kW/m2
4. SEARCH FOR THE GMP OF A GPV UNDER PARTIAL SHADING BY THE FLC
Explaining in this section we are going to propose a FLC command to give it the ability to track the
global MPP of a partially shaded PV array. In the presence of uniform irradiation, a single GMP appears on
the P-V chacarteristic and it is easy to track it down by the classical 'P-O'. In the event of partial shading, the
power output from GPV will decrease significantly and several multiple points may appear in the P-V
characteristic. In this case, the P-O algorithm has not the potential to find the GMP unless the operating point
does not change position. Our model is set with a limited number of rules [27]-[29]. The FLC command used
in the search for the GMP performs better, namely its robustness and simplicity. The main parts of FLC,
fuzzification, rule-base, inference and defuzzification, are show in Figure 5 [30]. The strategy to simulate our
configuration of Figure 1 is to insert between the GPV and the load a DC-DC converter (boost) controlled by
the FLC algorithm.
The principle of a FLC command is based on two input variables which are the error E and the error
change dE and an output variable dD (variation of the duty cycle). The value of the output variable, which
drives the static converter to find the GMP, is determined using a truth table and changing input parameters.
error E is defined as the difference between dP/dV(k) and the desired value dP/dV (k)=0. This last value
corresponds to the unique extremum of the curve P=f(V). This extremum is a maximum. The more positive E
is, the higher the value of P. Conversely, the more negative E is, the lower the value of P. Finally, when E
tends to 0, the value of P tends to its maximum, the MPP. It can be likened to the slope of P=f (V).
𝑑𝑃 𝑑𝑃
𝐸(𝑘) = (𝑘) − 0 = (3)
𝑑𝑉 𝑑𝑉
The variation in the deviation dE indicates in which direction and in what proportion the deviation
changes as the algorithm proceeds. So, when dE tends to 0, the system stabilizes. The first, second input
variable and the only one output for the FLC are divided into five fuzzy sets: NB (negative big), NM
(negative medium), Z (zero), PM (positive medium) and PB (positive big). The rules are formed as shown in
Table 2 and Figure 7 [31].
Figure 7. Input (E), (dE) and output (dD) membership Figure 8. Characteristic surface of the FLC
functions for proposed FLC command
Int J Pow Elec & Dri Syst, Vol. 12, No. 2, June 2021 : 1131 – 1138
Int J Pow Elec & Dri Syst ISSN: 2088-8694 1135
The FLC command depends heavily and necessarily on the correct selection of certain design
parameters namely input gains and output gain to improve the performances of the structure of the FLC
method. To guide us in searching for the appropriate gains for the FLC controller, we used the particle swarm
optimization method to optimize the parameters of the belonging functions associated with the input and exit
variables. The idea is to optimize our maximum power delivered by the GPV in this case GMP inserting the
command FLC by using a PSO program for the research of the appropriate gains of the FLC applied to the
pursuit of the MPP. PSO as one of meta-heuristic optimization techniques depends on searching the optimal
solution within the search area based on the exchange of experiences among particles in the population. The
particles in the swarm modify their positions in the next iteration based on individual local best position and
global best position of swarm. Each particle represents a solution for the control variables of the entire
optimization problem [32]. PSO is used to optimize the gains of FLC algorithm which optimize the power
output from PV array (MPPT) [33]-[35]. The proposed system simulated in MATLAB environment. To this
end a 'fuzzy _optmize_PSOMATLAB' program has been designed for the search for the gains of the FLC.
The appropriate gains for this optimization approach by running this PSO program as presented in Table 3.
function Fuzzy_Optimize_PSOMATLAB
warning off
tic
rng default
options =
optimoptions('particleswarm','PlotFcns',@pswplotbestf,'MaxIter',20,'SwarmSize',20);
lb = 0.00005*ones(1,2);
ub = 0.009*ones(1,2);
[x,fval,exitflag,output] = particleswarm(@func,2,lb,ub,options);
toc
end
function fitness=func(x)
assignin('base', 'Ta', [x(1) x(2)])
sim('Panneau_Shading_MPPT',0.1);
fitness=max(Sys_Out1.Data)
end
Figure 10. Bloc simulink of the GPV under partial shading with MPPT methods
8000
8
Power output from PV array (KW)
7000
X: 0.2035
X: 0.5729 Y: 7995
Power output from PV array (W)
Y: 7.997 6000
6
5000
X: 1.522
X: 1.612
4 Y: 3.208 4000
Y: 2880
3000
2 2000
1000
0 0
0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2
Time (S) Time (S)
(a) (b)
Figure 11. Curve of P(t) with variation of G (1000W/m² for100W/m²) of Figure 1 with, (a) ‘P&O’ and
(b) FLC
Power output from PV array(W)
8
Global MPP f or 1000W/m²
Figure 12. FLC tracking process from a uniform standard irradiation to the partial
shading configuration of Figure 1 with a shading level of 0.1kW/m2
Int J Pow Elec & Dri Syst, Vol. 12, No. 2, June 2021 : 1131 – 1138
Int J Pow Elec & Dri Syst ISSN: 2088-8694 1137
6. CONCLUSION
In this article, we show the inability of the classical 'P&O' algorithm to track the overall power point
for <200W/m2 values causing losses for a GVP in the presence of partial shading. We then proposed an
optimization approach to a fuzzy logic MPPT command that moves the operating point imposed by the
partial shading configuration to a point near the global. The results showed a perfect ability of the proposed
algorithm to track the overall MPP for the shading configuration of Figure1 for a low level of irradiation.
Finally, these results lead us to say that the proposed FLC control method can solve the overall problem of
tracking MPs with the classical 'P&O' algorithm.
APPENDIX
REFERENCES
[1] Dezso Sera, Laszlo Mathe, Tamas Kerekes, Sergiu Spataru and Remus Teodorescu, "On the parturb-and-observe
and Inc.Cond MPPT methods for PV systems," IEEE Journal of Photovoltaics, vol. 3, no. 3, pp. 1070-1078,
2013, DOI: 10.1109/Jphotov.2013.2261118
[2] T.-Y. Kim, H.-G. Ahn, S. K. Park, and Y.-K. Lee, "A novel MPPT control for photovoltaic power system under
rapidly changing solar radiation," IEEE Int. Symp. Ind. Electron., 2001, pp. 1011-014, DOI:
10.1109/ISIE.2001.931613
[3] Liu X., and Lopes L.A.C., "An improved ‘P&O’ MPPT algorithm for PV arrays," Conf. Rec. PESC’04, Aachen,
Germany, June 2004, pp. 2005-2010, DOI: 10.1109/PESC.2004.1355425
[4] Liu C., Wu B., and Cheung R., "Advanced algorithm for MPPT control of photovoltaic systems," Conf. Rec. of
the first Canadian Solar Buildings Conference, Montreal, August 2004.
[5] Femia, G. Petrone, G. Spagnuolo, and M. Vitelli, "Optimization of ‘P&O’ MPPT method," IEEE Trans. Power
Electron., vol. 20, no. 4, pp. 963-973, Jul. 2005, DOI: 10.1109/TPEL.2005.850975
[6] W. Xiao, J. Lind, W. Dunford, and A Capel, "Real-time identification of optimal operating points in photovoltaic
power systems," IEEE Transactions on Industrial Electronics, vol. 53, no. 4, pp. 1017-1026, August 2006, DOI:
10.1109/TIE.2006.878355
[7] J. C. Hernández, O. G. Garcia and Francisco Jurado, "Photovoltaic devices under partial shading conditions,"
International Review on Modelling and Simulations, vol. 5, no. 1, pp. 414-425, 2012.
[8] R. Ramaprabha and Dr.B.L. Mathur, "Comparative study of series and parallel configurations of solar PV array
under partial shaded conditions," International Review on Modelling and Simulations (IREMOS), vol. 3, no. 6,
pp. 1363-1371.
[9] M. Mostefai, A. Miloudi and Y. Miloud, "An intelligent MPPT for photovoltaic systems based on neural
network," International Review on Modelling and Simulations (IREMOS), vol. 6, no. 5, October 2013.
[10] K. H. Hussein and I. Muta, "MPPT: An algorithm for rapidly changing atmospheric conditions," Proc. Inst.
Electr. Eng. Gener., Transmiss. Distrib. vol. 142, no. 1, pp. 59-64, Jan. 1995.
[11] S. Jain and V. Agarwal, "A new algorithm for rapid tracking of approxi-mate MPP in photovoltaics systems,"
IEEE Power Electron. Lett., vol. 2, no. 1, pp. 16-19, Mar. 2004, DOI: 10.1109/LPEL.2004.828444
[12] L.A. Hecktheuer, A. Krenzinger, C.W.M. Prieb, "Methodology for photovoltaic modules characterization and
shading effects analysis," Journal of the Brazilian Society of Mechanical Sciences, vol. 24, no. 1, pp. 26-32,
2002, DOI: 10.1590/S0100-73862002000100004
[13] K.S. Parlak and H. Can, "A new MPPT method for PV array system under partially shaded conditions," 2012 3rd
IEEE International Symposium on Power Electronics for Distributed Generation Systems (PEDG), June 2012,
pp. 437-441, DOI: 10.1109/PEDG.2012.6254039
[14] Mohamed Amine Abdourraziq and Mohamed Maaroufi, "Experimental verification of the main MPPT techniques
for photovoltaic system," International Journal of Power Electronics and Drive Systems (IJPEDS), vol. 8, no. 1,
pp. 384-391, 2017, DOI: 10.11591/ijpeds.v8.i1.pp384-391
[15] Yaichi, M., Fellah, M. K, and Mammeri, A, "A neural network based MPPT technique controller for photovoltaic
pumping system," International Journal of Power Electronics and Drive Systems (IJPEDS), vol. 4, no. 2, pp.
241-255, 2014, DOI: 10.11591/ijpeds.v4i2.5875
[16] Ishaque, K and Salam, Z, "A review of MPPT techniques of PV system for uniform insolation and partial shading
condition," Renew. Sustain. Energy Rev, vol. 19, pp. 475-488, 2013, DOI: /10.1016/j.rser.2012.11.032
[17] H. Zheng, S. Li, R. Challoo and J. Proano, "Shading and bypass diode impacts to energy extraction of PV arrays
under different converter configurations," Renewable Energy, vol. 68, pp. 58-66, 2014, DOI:
10.1016/j.renene.2014.01.025
[18] M. A. M. Ramli, S. Twaha, K. Ishaque and Y. A. Al-Turki, "A review on MPPT for photovoltaic systems with
and without shading conditions," Renewable and Sustainable Energy Reviews, vol. 67, pp. 144-159, 2017, DOI:
10.1016/j.rser.2016.09.013
[19] K. Latha Shenoy, C. Gurudas Nayak and Rajashekar P Mandi, "Effect of partial shading in grid connected solar
PV system with FL Controller," International Journal of Power Electronics and Drive Systems (IJPEDS),
vol. 12, no. 1, pp. 431-440, 2021, DOI: 10.11591/ijpeds.v12.i1.pp431-440
[20] Liqun Liu, Xiaoli Meng and Chunxia Liu, "A review of MPPT methods of PV power system at uniform an partial
shading," Renew. Sustain. Energy Rev., vol. 53, pp. 1500-1507, 2016, DOI: 10.1016/j.rser.2015.09.065
[21] T. Sudhakar Babu, N. Rajasekar and K. Sangeetha, "Modified particle swarm optimization technique based
MPPT for uniform and under partial shading condition," Appl. Soft Comput. J., vol. 34, pp. 613-624, 2015, DOI:
[22] Kadhim Hamzah Chalok and Mohammad Faridun Naim Tajuddin, "Optimal extraction of photovoltaic energy
using FLC for MPPT technique," International Journal of Power Electronics and Drive Systems (IJPEDS),
vol. 11, no. 3, pp. 1628-1639, September 2020, DOI: 10.11591/ijpeds.v11.i3.1628-1639
[23] Dadjé, Abdouramani, et al. "MPPT methods for photovoltaic systems operating under partially shaded or rapidly
variable insolation conditions: a review paper," International Journal of Sustainable Engineering, vol. 9, no. 4,
pp. 224-239, 2016.
[24] Mohammad Syahir Bin Ishak, Rahmatul Hidayah Salimin, Ismail Musirin and Zulkiffli Abdul Hamid,
"Development of PV array configuration under different partial shading condition," International Journal of
Power Electronics and Drive Systems (IJPEDS), vol. 10, no. 3, pp. 1263-1269, Sep 2019, DOI:
10.11591/ijpeds.v10.i3.pp1263-1269
[25] Youcef and A. Miloudi, "Optimization of partially shaded PV array using a modified P&O MPPT algorithm,"
Leonardo Electronic Journal of Practices and Technologies, vol. 15, no. 28, pp. 179-196.
[26] H. Patel and V. Agarwal, "MATLAB-based modeling to study the effects of partial shading on PV Array
Characteristics," IEEE Trans. on Energy Conversion, vol. 23, no. 1, pp. 302-310, Mar 2008, DOI:
10.1109/TEC.2007.914308
[27] Adel Haddouche, Mohammed Kara and Lotfi Farah, "MPPT using FLC with reduced rules," International
Journal of Power Electronics and Drive Systems (IJPEDS), vol. 9, no. 3, pp. 1381-1390, Sep 2018
[28] Ahmad Saudi Samosir, Herri Gusmedi, Sri Purwiyanti and Endah Komalasari, "Modeling and Simulation of FLC
based MPPT for PV Application," International Journal of Electrical and Computer Engineering (IJECE), vol. 8,
no. 3, p. 1315-1323, 2018, DOI: 10.11591/ijece.v8i3.pp1315-1323
[29] Boualem Boukezata, Abdelmadjid Chaoui, Jean-Paul Gaubert and Mabrouk Hachemi, "An improved FLC MPPT
based P&O method to solve fast irradiation change problem," J. Renew. Sustain. Energy, vol. 8, no. 4, 2016,
DOI: 10.1063/1.4960409
[30] Masafumi Miyatake, Mummadi Veerachary, Fuhito Toriumi, Nobuhiko Fujii and Hideyoshi Ko, "MPPT of
Multiple Photovoltaic Arrays: A PSO Approach," IEEE Transactions on Aerospace and Electronic Systems,
vol. 47, no. 1, pp. 367-380, 2011, DOI: 10.1109/TAES.2011.5705681
[31] Yasser Ahmed and Ayman Hoballah, "Adaptive filter-FLC integration for torque ripples minimization in PMSM
using PSO," International Journal of Power Electronics and Drive Systems (IJPEDS), vol. 10, no. 1, pp. 48-57,
2019, DOI: 10.11591/ijpeds.v10.i1.pp48-57
[32] Ayman Hoballah and Istvan Erlich, "PSO-ANN approach for transient stability constrained economic power
generation," IEEE Bucharest Power Tech, Bucharest, 2009, pp. 1-6, DOI: 10.1109/PTC.2009.5281926
[33] Yen-Shin Lai and Juo-Chiun Lin, "New hybrid fuzzy controller for direct torque control induction motor drives,"
IEEE Transactions on Power Electronics, vol. 18, no. 5, pp. 1211-1219, 2003, DOI: 10.1109/TPEL.2003.816193
[34] Kashif Ishaque, Zainal Salam, Muhammad Amjad, and Saad Mekhilef, "An Improved PSO-Based MPPT for PV
with Reduced Steady-State Oscillation," IEEE Transactions on Power Electronics, vol. 27, no. 8, pp. 3627-
33638, 2012, DOI: 10.1109/TPEL.2012.2185713
[35] K. Sundareswarana, V. Vignesh Kumara and S. Palanib, "Application of a combined PSO and perturb and
observe method for MPPT in PV systems under partial shading conditions," Renewable Energy, vol.75, March
2015, DOI: 10.1016/j.renene.2014.09.044
Int J Pow Elec & Dri Syst, Vol. 12, No. 2, June 2021 : 1131 – 1138