MPPT (PandO-INC-FUZZY-ANFIS) Simulation and Hardware Implementation.
MPPT (PandO-INC-FUZZY-ANFIS) Simulation and Hardware Implementation.
Cascade
PV dc- dc Multilevel Micro Grid
Array Converter Inverter
MPPT/CVT
Algorithm
NR-SHE
Optimized Algorithm
Firing Angles α
Inferences Drawn Out of Literature Review
i. The multilevel converter technology is a promising technology for high power
electric devices because of its high voltage operation, high efficiency and low
electromagnetic interference (EMI). The desired output of a multilevel converter
may be synthesized by several dc voltages sources. With an increasing number
of dc voltage sources, the converter output voltage waveform approaches a
nearly sinusoidal waveform.
ii. Very little work has been reported in regard of application of a multilevel
inverter with equal and unequal dc voltage source for H-bridge from renewable
zero-pollution energy resources, such as wind, solar, bio and geothermal.
iii. The cascaded H-bridge multilevel inverter topology requires an isolated and
constant dc source for each H-bridges, thus allowing the integration of
renewable energy sources such as wind generators, PV arrays, fuel cells, etc. to
the smart grid.
iv. In literature the switching angles for multilevel inverter have been calculated by
time-consuming equations offline. Alternate approaches were developed to
calculate the switching angle in real time but these approaches were not
extended for unequal dc sources.
v. Switching angles can be calculated offline and stored in lookup table. But for
some operating points the solutions might be missing. This also involves large
memory space for program code.
vi. The performance of cascaded multilevel inverter is compared based on
computation of switching angle using Newton Raphson approach. The genetic
algorithm method has not succeeded to find switching angles for some
modulation indexes which have solutions.
vii. Soft computing research is concerned with the integration of artificial
intelligent tools such as neural networks, fuzzy logic technology,
evolutionary algorithms, etc. in a complementary hybrid framework for
solving real world problems.
viii. A maximum power point tracking (MPPT) strategy is needed for
maximizing the energy harvested from each string of the PV module.
ix. Various MPPT methods, with their merits and demerits are available in
the literature. P and O method has been widely used due to its
simplicity.
x. The ANN based controller reflects better maximum power point tracking
performance as compared to the conventional PID controller and avoids
the tuning of controller parameters.
xi. The harmful effects of non-uniform irradiance (due to partial shading)
on the performance of a PV module are one of the main causes of
overheating of shaded cells which reduces the energy yield of the
module.
xii. Hybrid combination of the fixed and two motor based Sun tracking type
PV system with fractional order controller for maximum power point
tracking also improves the efficiency of PV generation systems.
xiii. MATLABTM/SIMULINKTM have emerged as popular tools and are widely
used for modelling, control and simulation purposes of various electrical
systems, including the component of solar PV generator and under a
variety of operating conditions.
Proposed Work
where,
Ipv : photovoltaic current
I0 : saturation current RsH
V
Rs : equivalent series resistance
Vt =NskT/q : thermal voltage
Rp : equivalent parallel resistance
a : diode ideality constant _
𝑅(1 –𝐷)
𝐿𝐶 = 𝐿 = (10)
2𝑓
1 –𝐷
𝐶𝐶 = 𝐶 = (11)
1𝐿𝑓2
For a switching frequency of 80 KHz and inductance current ripple (∆𝐼) of 10%
the 𝐿𝑐 and 𝐶𝑐 are approximated as 1mH and 100µF respectively
Parameters of Buck Converter
Sr. No. Parameter Value
8 Input voltage 20 V
Radiation Sensor
Buck
PV Current Converter Load
Sensor
PV Array
Output
PV Voltage Voltage
Sensor Dspace PC with Sensor
ds1104 MATLAB/
SIMULINK
(MPPT
Algorithm)
MPPT Algorithms
Membership Function of ANFIS Input (Current) Membership Function of ANFIS Input (Voltage)
Fuzzy inference rules
If (current is in1mf1) and (voltage is in 2mf1) then (duty cycle is out1mf1)
(c ) Fill Factor FF
𝑅𝑠
𝐹𝐹 = 𝐹𝐹0 1 − 𝑉𝑜𝑐 (5)
𝐼𝑠𝑐
𝑣𝑜𝑐−ln(𝑣𝑜𝑐 +0.72)
𝐹𝐹0 = (6)
1+𝑣𝑜𝑐
V
voc = nKToc (7)
q
Response of Different MPPT Algorithms for Variation in Irradiation According to Test Pattern
Proposed ANFIS / Constant Voltage
Tracker (CVT) Algorithm
Output Voltage Corresponding to Irradiance Level
Back
Output Power Corresponding to Irradiance Level
Experimental Setup for Performance Evaluation of PV
Module at Outdoor Condition
Performance Characteristics (a) Efficiency 9.95% (b) Efficiency 14.63 %
Solar Module Analyzer 12A V110906 Software
Experimental Setup for Performance Evaluation of PV
Module at Laboratory Condition
Performance of 37 W PV Module at Laboratory and
Outdoor Conditions
Condition
Voc Im
Angle of PV Irradiation Temperature Isc Vm (V) Pm 𝜂
(V) (mA)
panel Tilt W/m2 oC (mA) (W) (%)
Pmax (W) 37
Vmp (V) 17
Imp (A) 2.25
Voc (V) 21.77
Isc (A) 2.4
Efficiency (%) 12.01
Block Diagram of the Proposed MPPT/CVT System
Subsystem Maximum Power
Prediction Model
in MATLABTM/
SIMULINKTM
Environment
Real Time Data Acquisition using P&O
Real Time Data Acquisition using INC Algorithm
Real Time Data Acquisition using INC
Real Time Data Acquisition using Neural Network
Real Time Data Acquisition using ANFIS
Theoretical)
e (W)
Delay in
dynamic
Response
Max Steady State
Error
Initial
Offset
Theoretical) (W)
Initial
Offset
Initial
Offset
Theoretical) (W)
Initial
Offset
e
Initial
Offset
Initial
Offset
0
5
0.0001
0.0891
0.1781
0.2671
0.3561
0.4451
0.5341
0.6231
0.7121
0.8011
0.8901
0.9791
1.0681
1.1571
1.2461
1.3351
1.4241
1.5131
1.6021
1.6911
1.7801
1.8691
1.9581
2.0471
2.1361
2.2251
2.3141
2.4031
2.4921
2.5811
2.6701
2.7591
2.8481
2.9371
3.0261
3.1151
3.2041
3.2931
3.3821
3.4711
3.5601
12V
3.6491
3.7381
Vout(V)
3.8271
VpvIn(V)
3.9161
MPPT Algorithm Comparison
Delay in Max.
MPPT method Efficiency Over-shoot Settling time dynamic Steady Sensors
( %) ( %) ( sec) response state error used
( sec) ( %)
PO 77.60 to 79.39 No 0.48 0.06 15.14 Voltage,
(ΔD=0.1) Current
PO 81.00 to 81.60 No 0.41 0.039 12.77 Voltage,
(ΔD=0.01) Current
PO 81.23 to 84.37 No 0.40 0.04 12.03 Voltage,
(ΔD=0.001) Current
INC PI 86.32 to 87.25 3.35 1.78 0.001 7.35 Voltage,
Current
NN 87.35 to 2.185 0.6439 0.038 3.88 Voltage,
90.10 Current
ANFIS 87.15 to 6.56 5.35 0 3.55 Voltage,
93.31 Current
≥12V NA 7.28 0.18 0.1 9 Input and
ANFIS & Output
CVT <12V 87.15 to 6.56 5.35 0 3.55 Voltage
93.31 , Current
Back
Conclusion
The MPPT algorithms Perturb and observe , incremental
conductance, neural network, adaptive neuro fuzzy
inference system (ANFIS) and ANFIS & CVT were
discussed, implemented and compared . The modeling
of the PV array was performed in MATLAB/SIMULINK. It
was conclude that ANFIS model gave fast response and
less oscillations compared to Perturb and Observe ,
Incremental Conductance and Neural Network models .
The proposed ANFIS & CVT model tracks 12V for voltage
more than 12V and tracks the maximum power point
according to ANFIS algorithm for voltage less than 12V.
This algorithm is suitable for a PV powered multilevel
inverter which requires an isolated constant dc supply at
its input.
Cascaded H-bridge Inverter
S1 S2
Vdcm
Va
Va[(m-1)/2]
S3 S4
S1 S2
Vdc1
Va
n
S3 S4
(a) (b)
(10)
(11)
(12)
(13)
(14)
(16)
(17)
Block Diagram of the Proposed Harmonic
Elimination System
Ig1(1-4)
Negative Ig2(1-4)
Logic
……
IgL(1-4)
Step 1 Assume any random initial guess for switching angles (say 𝛼0 )
The switching angle matrix is :
𝛼 𝑗 = [𝛼1 𝑗 + 𝛼2 𝑗 + 𝛼3 𝑗 + 𝛼4 𝑗 + 𝛼5 𝑗 ]𝑇 (18)
Step 2 Set modulation index to zero.
𝜕𝑓 𝑗
Step 3 Evaluate the non-linear system matrix 𝐹𝑗
, the Jacobian matrix and
𝜕𝛼
the harmonics amplitude matrix 𝑇 represented below:
The non-linear system matrix,
cos 𝛼1 𝑗 + cos 𝛼2 𝑗 + cos 𝛼3 𝑗 + cos 𝛼4 𝑗 + cos 𝛼5 𝑗
cos 5𝛼1 𝑗 + cos 5𝛼2 𝑗 + cos 5𝛼3 𝑗 + cos 5𝛼4 𝑗 + cos 5𝛼5 𝑗 (19)
𝐹𝑗 = cos 7𝛼1 𝑗 + cos 7𝛼2 𝑗 + cos 7𝛼3 𝑗 + cos 7𝛼4 𝑗 + cos 7𝛼5 𝑗
cos 9𝛼1 𝑗 + cos 9𝛼2 𝑗 + cos 9𝛼3 𝑗 + cos 9𝛼4 𝑗 + cos 9𝛼5 𝑗
cos 11𝛼1 𝑗 + cos 11𝛼2 𝑗 + cos 11𝛼3 𝑗 + cos 11𝛼4 𝑗 + cos 11𝛼5 𝑗
the Jacobian matrix,
NR Algorithms
11 Level Cascaded H-bridge Inverter Applied with NR-SHE Algorithm for 0.8 Value of MI
Current Waveform
Harmonic Spectrum at 0.8 MI for NR-SHE Algorithm for a 11 level Cascaded H-bridge Inverter
Current Spectrum
Magnitude of Harmonic Contents (%) up to 19th Order for 11 Level
Cascaded H-bridge Inverter Applied with NR-SHE Technique
11 Level Cascaded H-bridge Inverter
Magnitude of Harmonic Contents (%) up to 19th Order
Technique Used
Harmonic
Harmonic
Harmonic
Harmonic
Harmonic
Harmonic
Harmonic
Harmonic
Harmonic
Harmonic
Harmonic
Order
Order
Order
Order
Order
Order
Even
Even
Even
Odd
Odd
Odd
0th 0.00 1th 100 0th 0.00 1th 100 0th 0.01 1th 100
2nd 0.00 3rd 0.02 2nd 0.00 3rd 0.60 2nd 0.00 3rd 0.02
4th 0.00 5th 0.09 4th 0.00 5th 0.04 4th 0.00 5th 0.07
NR
6th 0.00 7th 0.08 6th 0.00 7th 0.06 6th 0.00 7th 0.09
8th 0.00 9th 0.06 8th 0.00 9th 3.26 8th 0.00 9th 0.06
10th 0.00 11th 0.10 10th 0.00 11th 0.10 10th 0.00 11th 0.11
12th 0.00 13th 0.02 12th 0.00 13th 0.02 12th 0.00 13th 0.03
14th 0.00 15th 0.09 14th 0.00 15th 1.04 14th 0.00 15th 0.08
16th 0.00 17th 2.65 16th 0.00 17th 2.58 16th 0.00 17th 2.62
Comparison of Magnitude (Peak Value) and THD values of Line Voltage, Phase Voltage
and Current of 7, 9 and 11 Level CHMLI
𝑁𝑑𝑐 ∗𝑉𝑑𝑐
𝑉1 = 𝑀𝑖(4 ∗ )
𝜋
Method
THD
Harmonic
Method
Computational
Level MI Value
𝛂𝟏 𝛂𝟐 𝛂𝟑 𝛂𝟒 𝛂𝟓 Time (s)
(%)
S2 S1 S0 R6 R5 R4 R3 R2 R1 R0 V8 V7 V6 V5 V4 V3 V2 V1 V0
0 1 0 0 0 0 0 0 0 1 0 0 0 0 0 0 1 0 0
CHMLI
Spartan®-3A
DSP FPGA
CHMLI Power Quality Analyzer
3Φ Induction Motor
Experimental Results for 11 Level
Inverter (a) Output Line Voltage (b)
Phase Voltage and (c) Current at
M=0.8 (NR-SHE)
(a)
(b)
(c)
Experimental Results for 11 Level
Inverter (a) Line Voltage FFT
Analysis (b) Phase Voltage FFT
Analysis and (c) Current FFT
Analysis at M=0.8 (NR-SHE)
(a)
(b)
(c)
(b)
Optimum Switching Angles and Minimum THD using NR-SHE, GA-SHE
and PSO-SHE
[email protected]
9417588987
124