This paper presents the design and realization of a digital PV simulator with a Push-Pull Forward (PPF) circuit based on the principle of modular hardware and configurable software. A PPF circuit is chosen as the main circuit to restrain the magnetic biasing of the core for a DC-DC converter and to reduce the spike of the turn-off voltage across every switch. Control and I/O interface based on a personal computer (PC) and multifunction data acquisition card, can conveniently achieve the data acquisition and configuration of the control algorithm and interface due to the abundant software resources of computers. In addition, the control program developed in Matlab/Simulink can conveniently construct and adjust both the models and parameters. It can also run in real-time under the external mode of Simulink by loading the modules of the Real-Time Windows Target. The mathematic models of the Push-Pull Forward circuit and the digital PV simulator are established in this paper by the state-space averaging method. The pole-zero cancellation technique is employed and then its controller parameters are systematically designed based on the performance analysis of the root loci of the closed current loop with
ki
and
RL
as variables. A fuzzy PI controller based on the Takagi-Sugeno fuzzy model is applied to regulate the controller parameters self-adaptively according to the change of
RL
and the operating point of the PV simulator to match the controller parameters with
RL
. The stationary and dynamic performances of the PV simulator are tested by experiments, and the experimental results show that the PV simulator has the merits of a wide effective working range, high steady-state accuracy and good dynamic performances.
I. INTRODUCTION
In recent years, due to the serious problems of energy shortages, environmental pollution and energy security throughout the world, the development and utilization of clean renewable energy is becoming an important way to improve the energy structure, reduce environment pollution and persistently ensure the energy supply all over the world
[1]
,
[2]
. The worldwide interest in renewable energy has greatly promoted the development of PV generation technology and industry
[3]
-
[5]
.
In the research and development of PV generation systems, if a real PV array is used to test and verify the correctness and performance of the power electronic equipment and control systems, it would be both costly and time-consuming. Therefore, the output characteristics of a PV array under various operating conditions are always simulated by a PV simulator, which can be made by an analog circuit, a digital circuit or an analog and digital hybrid circuit. There are several methods to implement a PV simulator, including:
-
1) amplifying the output of the current and voltage of a sample PV cell or photosensor diode with a controlled light source imitating sunlight or an LED light emission circuit to simulate variations in sunshine intensity[6],[7]
-
2) constructing a PV array by taking a PV cell equivalent circuit with analog circuits as a basic unit[8]
-
3) using a DC power supply controlled by a personal computer (PC) with a data acquisition card or a special real time digital simulator (such as RT-LAB)[9],[10]
-
4) adopting a unidirectional DC-DC converter with a basic topology controlled by a DSP, a PC, a microprocessor, a microcontroller, or some other digital control technology as the main circuit of the digital PV simulator to simulate the output characteristics of the PV array, where the DC input voltage is supplied by a DC power supply or an uncontrolled diode rectifier[11]-[14]
-
5) adopting the topology of a PWM rectifier with a bidirectional compound DC/DC converter controlled by a DSP or a PC[15]
-
6) using multiple multiphase DC/DC converters controlled by a PC or a DSP as the main circuit of the PV simulator[16],[17]
-
7) selecting a high frequency PWM rectifier for a three-phase voltage source controlled by the host computer as the main circuit of the PV simulator[18].
Despite achieving adequate simulative effects, most of the PV simulators still exhibit some problems, such as a limited effective working range, bigger ripples of the output voltage and current, lack of systematic modeling and design methods for the controller parameters, lower development efficiency, and so on.
Considering the above drawbacks, the digital PV simulator presented in this paper adopts modular hardware, configurable software, systematic modeling and design methods for the controller parameters so as to obtain the merits of a wide effective working range, high steady-state accuracy, good control performances, high reliability and development efficiency. A comparison of the characteristics of the different methods for PV simulators is listed in
Table I
.
CHARACTERISTICS COMPARISON OF DIFFERENT METHODS FOR PV SIMULATOR
CHARACTERISTICS COMPARISON OF DIFFERENT METHODS FOR PV SIMULATOR
II. CONFIGURATION AND WORKING PRINCIPLE OF THE DIGITAL PV SIMULATOR
The overall configuration of the digital PV simulator is shown in
Fig. 1
. It mainly consists of a DC power supply, a DC-DC converter, a digital controller, a PWM signal generation circuit, a measuring and signal conditioning circuit, a driving and isolating circuit, and so on. The main operation principle is as follows: the DC power supply provides a steady DC input voltage for the DC-DC converter. According to the output voltage measured in real-time, the digital controller calculates the reference current signal by the I-V characteristic model and then calculates the duty ratio of the converters’ PWM signals by a current closed-loop regulator. After signal conversion, the duty ratio is converted into a signal-level control voltage which is sent to the PWM signal generation circuit to generate PWM signals. The PWM signals control the states of the power switches so as to adjust the converter’s output voltage and current after signal isolation and power driving. Finally, the PV simulator can always work at a certain operational point at which the output voltage and output current of the converter can satisfy both the load characteristics and the I-V curve of the PV array. Thus, the output characteristics of the PV array can be easily simulated by the digital PV simulator.
Overall configuration of the digital PV simulator.
III. WORKING MECHANISM OF THE PUSH-PULL FORWARD CIRCUIT
- A. Circuit Structure
The main circuit of the DC-DC converter adopts a Push-Pull Forward (PPF) circuit in the PV simulator. A full-wave rectifier is applied in the secondary windings of the high frequency transformer so as to reduce the conduction losses of the rectifying diodes with a low-output voltage and a high-output current. An LC filter is used to provide a continuous load current and to reduce harmonics.
Fig. 2
shows the topological structure of the DC-DC converter. It can be seen from the figure that the PPF circuit is different from a conventional Push-Pull converter in that a clamping capacitor is connected in series between the two primary windings of the high frequency transformer and between the two power switches. As a result, it has the advantages of both the Push-Pull converter and Forward converter, and it can restrain the magnetic biasing of the core for the DC-DC converter, magnetize the magnetic core bi-directionally and reduce the spike in the turn-off voltage across each switch. Therefore, the PPF circuit is becoming a kind of preponderant circuit topology for low voltage and high current output applications.
Diagram of Push-Pull Forward circuit.
- B. Analysis of the Operating Modes
In steady state, the PPF circuit has 8 operating modes in a switch cycle. Each power switch matches with the 4 operating modes of one cycle. In the following analysis, it is assumed that all of the power switches and diodes are ideal components and that the on-state voltage drops for all of the power electronic components are negligible
[19]
. The equivalent circuits of the former four operating modes are shown in
Fig.3
.
Equivalent circuits of the former four operating modes in a switch cycle.
1) Mode [t0-t1]:
Prior to
t0
, the power switches S
1
and S
2
are off. The primary current freewheels through
Uin(+)
→
Np2
→
Cb
→
Np1
→
Uin(-)
and forms a circular current which is
Ia
=
ip1
=
ip2
. During this period, the rectifying diodes of the secondary windings turn on simultaneously, and the current through each of the rectifying diodes is
iDR1
=
iDR2
=
Io
/2. At
t0
, the filtering inductor current runs down to the minimum value
ILfmin
.
At
t0
, S
1
turns on.
Uin
and
ucb
are applied to the leakage inductances
Lσ
of
Np1
and
Np2
, respectively.
ip1
increases rapidly due to
Uin
being in the same direction as
ip1
, and
ip2
decreases rapidly due to
ucb
being in the opposite direction as
ip2
. Corresponding with
ip1
and
ip2
,
iDR1
and
iDR2
increase and decreases respectively because of the affection of the magnetic circuit. During this period, the filtering inductor current increases gradually from
ILfmin
. This mode ends when
iDR1
is equal to the load current, both
ip2
and
iDR2
are decreased to 0 and the commutation process ends. The voltage across S
1
is always 2
Uin
during this period.
2) Mode [t1-t2]:
In this mode,
Uin
and
ucb
are applied to the magnetizing inductor
Lm
and the primary reduced inductance of the filtering inductances
Lf
of
Np1
and
Np2
, respectively. Each half of the changing rate of the magnetizing current and load current is shared by the primary winding.
ip1
continues to increase in the original direction defined as the reference direction whereas
ip2
increases gradually from 0 in the opposite direction. This mode lasts from
t1
to
t2
and S
1
turns off at
t2
. Because the circuit works like two forward converters in parallel, the circuit is called Push-Pull Forward circuit or converter. At the instant of S
1
turning off,
ip1
runs up to the maximum value, because both the magnetizing current and primary reduced current of the load current flow through S
1
. The filtering inductor current also runs up to the maximum value accordingly. The voltage across S
1
is also 2
Uin
during this period.
3) Mode [t2-t3]:
After S
1
turns off, the body diode
DS2
of S
2
is forced to turn on to continue the leakage inductance current because
ip1
is always larger than
ip2
before
t2
. The energy of the leakage inductance is released to charge the clamping capacitor
Cb
through the low impedance loop
Np1
→
DS2
→
Cb
.
Uin
and
ucb
are applied to the leakage inductances of
Np2
and
Np1
, respectively.
ip1
decreases rapidly, and
ip2
reduces to 0 gradually and then changes direction to increase rapidly in the reference direction.
iDR1
and
iDR2
decrease and increase, respectively. When
iDR1
and
iDR2
reach the same value
Ia
at
t3
, this mode ends. During this period, the load current decreases gradually from
ILfmax
.
4) Mode [t3-t4]:
During this period, both S
1
and S
2
are off, and the leakage inductance current is freewheeling through
Uin(+)
→
Np2
→
Cb
→
Np1
→
Uin(-)
and forms a circular current in the primary windings. Due to the two primary windings being series-opposing connections, the voltages across them are both 0 and the voltages across the power switches are both
Uin
. At the same time, magnetizing current also forms a circular current in the secondary side of the high frequency transformer. This mode ends when S
2
turns on at
t4
.
5) Mode [t4-t5][t5-t6][t6-t7][t7-t8]:
At
t4
, S
2
turns on and the circuit goes into the second half of a switch cycle. The modes [
t4-t5
], [
t5-t6
], [
t6-t7
] and [
t7-t8
] in the second half of a switch cycle correspond to [
t0-t1
], [
t1-t2
], [
t2-t3
] and [
t3-t4
], respectively. The operating principle of the second half of a switch cycle works in the same way as the first half, except that the magnetizing currents of the two half cycles are opposite in direction. During the second half of a switch cycle, demagnetization of the high frequency transformer is completed.
The key operating waveforms are shown in
Fig.4
.
Key waveforms of Push-Pull Forward circuit.
- C. Input-Output Relationship
Suppose the primary winding turn of the transformer is
NP1
=
NP2
=
w1
and its secondary winding turn is
Ns1
=
Ns2
=
w2
, then the turns ratio of the secondary to primary is
n
=
w2
/
w1
. The duty ratio D is defined as 2
Ton
/
Ts
, where
Ton
is the on duration of each power switch (S
1
or S
2
) and
Ts
is the switch time period which is equal to the reciprocal of the switching frequency
fs
.
During steady-state operation, the filtering inductor current is of the triangle waveform, and it varies periodically between
ILfmin
and
ILmaxf
. The inductor current increase Δ
ILf(+)
, during the on duration of S
1
or S
2
, is equal to its decrease Δ
ILf(-)
while S
1
and S
2
are both off. The relationship between them can be expressed as:
Thus:
That is:
where
Uin
and
Uo
are the average input voltage and output voltage, respectively.
If the power losses of the circuit are neglected, the average input current can be expressed as :
where
Io
is the average load current,
Io
=(
ILfmin
+
ILfmax
)/2.
IV. MODELING OF THE DIGITAL PV SIMULATOR
Fig. 5
shows a control block diagram of the digital PV simulator. There are two control loops in it: the outer-reference loop and the inner-current loop. The former measures the output voltage
uo
of the PV simulator and feeds it back to the I-V characteristic mode so as to generate the reference current
iref
; while the latter controls the filtering inductor current
iLf
.
RL
represents the equivalent load resistance of the PV simulator’s load characteristics;
G
(
s
) is the transfer function of the duty ratio to the filtering inductor current;
Gc
(
s
) is the transfer function of the PI controller,
Gc
(
s
)=
kp
+
ki
/
s
;
Gfi
(
s
) is the equivalent transfer function of the current measuring, conditioning and filtering,
Gfi
(
s
)=1/(
Tf1s
+1); and
Gfu
(
s
) is the equivalent transfer function of the voltage measuring, conditioning and filtering,
Gfu
(
s
)=1/(
Tf2s
+1). The I-V characteristics can be described from an engineering analytical model of the PV array. Because the changing rate of the filtering capacitance voltage
uCf
is very small in steady-state operation, the current
iCf
of the filtering capacitor branch is negligibly small and
iLf
is almost equal to
io
. Thus,
iLf
replaces
io
for use as the feed current in the actual control system. The modeling and controller parameter design of the PV simulator are carried out mainly based on the inner-current loop.
Control block diagram of the digital PV simulator.
According to the working principle and the input-output relationship of the PPF circuit, when S
1
(or S
2
) is on, the relationship between the primary and secondary is:
where
iin
,
iLf
and
uF
correspond to the instantaneous value of input current, the filtering inductor current and the pulse voltage rectified by the rectifying diodes, respectively.
When S
1
and S
2
are both off, the relationship is:
By the above relationship and Kirchhoff’s voltage and current laws, considering the effect of the inductance parasitic resistor
Rf
, the state-variable equations of the PPF circuit with S
1
(or S
2
) being on are:
where
Ton
=
t0-2
or
Ton
=
t4-6
,
t0-2
is the duration from
t0
to
t2
for each switch cycle, and
t4-6
is the duration from
t4
to
t6
for each switch cycle. The defining methods for the following similar variable are similar to
t0-2
and
t4-6
.
Similarly, the state-variable equations of the PPF circuit when S
1
and S
2
are both off are:
where
Toff
=
t2-4
or
Toff
=
t6-8
.
Combining the above two switch statuses and according to the working mechanism of the simulator, the state space averaging equations of the PPF circuit in a switch time period are:
Due to
Ton
+
Toff
=
Ts
/2, D=2
Ton
/
Ts
, the above equations can be rewritten in state space expression as:
Assuming that the initial conditions are zero, the transfer function of the duty ratio to the filtering inductor current by Laplace transforms of equation (10) can be obtained as:
V. CONTROLLER PARAMETER DESIGN OF THE DIGITAL PV SIMULATOR
- A. Ideas of the Controller Parameter Design
The open loop transfer function of the current loop is:
There are two parameters (
kp
and
ki
) that need to be determined in equation (12).
kp
and
ki
are usually designed to match with the transfer function
G
(
s
). However,
G
(
s
) varies with
RL
, and
RL
varies with the operating point and has a wide variation range. As a result, together with the high dimension of
Gok
(
s
), it is very difficult to determine the two controller parameters. If the pole of
Gfi
(
s
) is cancelled by the zero of the PI controller, namely,
kp
/
ki
=
Tf1
, the design difficulty will be reduced to a certain extent. However,
ki
remains difficult to determine due to the uncertainty of
RL
.
Under the condition of
kp
/
ki
=
Tf1
, the performances of the root loci with
ki
and
RL
as variables are analyzed for designing the controller parameters. On the basis of the analytical method, several sets of PI controller parameters matched with different
RL
are designed. In order to match the controller parameters with
RL
, a fuzzy PI controller based on the Takagi-Sugeno fuzzy model is applied to self-adaptively regulate the controller parameters so as to satisfy the control demands as the operating point changes.
- B. Performance Analysis of the Root Loci with ki
The open loop transfer function of the current loop with
ki
as a variable is:
The root loci of the current loop with
ki
as a variable for different values of
RL
are a family of curves. It is found by analyzing the root loci that the root loci show three different types of motion curves with
RL
=1.287Ω and
RL
=1.403Ω as cut-off points. The typical root loci are shown in
Fig. 6
. It can be seen from the figure that the system is always stable no matter how
RL
changes. However, the locations of the zeros and poles vary with
RL
. For different locations of the zeros and poles, the performances of the current loop show large differences. At a given
RL
, different closed-loop poles have different values of
ki
and performance indexes, such as the damping ratio
ζ
and the overshot. It is obvious that the value of
ki
determines the location of the closed-loop poles and the performances of the system.
Root loci of the current loop with ki as variable( (a) RL=1Ω, (b) RL=1.3Ω, and (c) RL=50Ω).
- C. Performance Analysis of the Root Loci with RL
According to the characteristic equations of the closed current loop system, its open loop transfer function with
RL
as a variable is rewritten as:
The root loci of the current loop with
RL
as a variable for different values of
ki
are also a family of curves. Although the locations of the zeros and poles change with
ki
, the form of the root loci is basically the same. Typical root loci are shown in
Fig. 7
.
Root loci of the current loop with RL as variable.
It can be seen from the figure that the system is always stable no matter how
ki
changes. When
ki
kept constant, different values of
RL
have different performance indexes. Apparently,
RL
is also a key factor to influence the system’s performances.
- D. Design of Controller Parameters
From the analysis presented above, it can be seen that only suitable values of
ki
matched with
RL
can satisfy the performance requirements of the system because
RL
varies with the operating point.
For this purpose, by combining the time domain performance indexes of the equivalent 2-order system with the output characteristics of the PV array, typical operating points of the segmented operating ranges are selected to determine the closed-loop dominant poles and the corresponding
ki
based on the root loci analysis. The closed-loop dominant poles are selected close to an imaginary axis and the absolute values of the real parts for the other close-loop poles and zeros are three times more than those of the closed-loop dominant poles. Once the closed-loop dominant poles are decided, the corresponding values of
ki
can be calculated. For the given
Tf1
=0.01s, the value of
kp
can be obtained according to the relationship
kp
/
ki
=
Tf1
. The PI controller parameters of the typical operating points in each of the segmented operating ranges are listed in
Table I
of the Appendix.
- E. Self-adaptive Regulation of the Controller Parameters
It can be seen from the above analysis that the controller parameters must be regulated timely and appropriately according to
RL
in order to obtain satisfactory dynamic and stationary performances. For these reasons, a fuzzy PI controller based on the Takagi-Sugeno fuzzy model is applied to regulate the controller parameters according to the change of the operating point of the PV simulator. The structure of the fuzzy PI controller is shown in
Fig. 8
.
Structure of fuzzy PI controller.
In this paper,
RL
is selected as the reference factor of the input variable, the input space is divided into four fuzzy subspaces based on the above design results of the controller parameter and every fuzzy subspace corresponds to the linear PI controller designed above. The models of the linear PI controllers are joined together smoothly by membership functions. With overlapping among the other inference rules, the Takagi-Sugeno fuzzy model achieves the nonlinear global mapping and makes the PI controller parameters of each rule have a different weighting factor. Therefore, with the fuzzy inference of the Takagi-Sugeno fuzzy model, the PI controller parameters can self-adaptively regulate according to the change of the operating point and
RL
. The value of R
i
c
of the Takagi-Sugeno fuzzy model’s fuzzy implication relationship has the following form
[20]
:
-
ifRLisAithenkp=kip,ki=kii,i=1, 2, …,m
where,
i
is the
i
-th rule;
m
is the number of fuzzy subspaces,
m
=4;
Ai
is the fuzzy subsets of the
i
-th rule; and
kip
and
kii
are the outputs of the
i
-th rule.
The total output of the fuzzy implication relation is given by:
where,
μAi
(
RL
) is the membership function of
RL
and the application degree of the
i
-th rule and is decided by the membership degree for all of the input subsets in the rule.
The membership function of
RL
is the Gaussian type and it is shown in
Fig. 9
.
Membership function of RL.
VI. HARDWARE AND SOFTWARE DESIGN OF THE DIGITAL PV SIMULATOR
- A. Push-Pull Forward Circuit
The design parameters of the PV simulator are shown in
Table II
of the Appendix. According to the above design requirements, the duty ratio of the PPF circuit and the turns ratio of the secondary to the primary are confirmed firstly. Considered the extreme case where the maximum output voltage prescribed by the design indexes can be obtained under the minimum input voltage and the effective duty ratio is increased as much as possible in order to decrease the turns ratio of the secondary to the primary, the turns ratio of the secondary to the primary can be calculated by:
where
Dmax
is the maximum duty ratio of the PPF circuit,
Dmax
=0.9.
The other parameters of the high-frequency transformer are calculated by using the Area Product method. According to the working mechanism of the PPF circuit and considered a certain safety margin, 1MBH60D-100 IGBT modules and MUR3040PT diodes are selected as the power switches and output rectifying diodes, respectively.
According to the working mechanism of the PPF circuit, the averaging voltage of the clamping capacitor is equal to
Uin
. Due to the finite value of the clamping capacitor, there is a certain pulsating quantity of voltage Δ
ucb
, and it can be expressed as:
Based on the estimation of the above formula and experimental research, a 10
μF
high-frequency non-inductive capacitor is selected as the clamping capacitor.
The key parameters of the PPF circuit are listed in
Table III
of the Appendix.
- B. PWM Signal Generation, Driving and Isolating Circuit
In the PV simulator, a SG3525 is adopted as the pulse width modulation chip to generate two push-pull PWM signals according to the input control voltage. The two PWM signals generated by the SG3525 are amplified by two EXB840 modules to drive the power switches. The output signal of the EXB840’s pin5 is sent to pin10 of the SG3525 by photocoupler so as to shut down the PWM signals for protecting the power switches when an overcurrent occurs. In order to prevent the power switches from turning on simultaneously, the two PWM signals are interlocked in the driving circuit so as to enhance the reliability of the system. A diagram of the PWM signal generation, driving and isolating circuit is shown in
Fig.10
.
Diagram of PWM signal generation, driving and isolating circuit.
- C. Signal Measuring and Conditioning Circuit
The output voltage and inductor current are measured and converted into 0~2.5V voltage signals by Hall sensors. The signals are converted into 0~10V voltage signals by a signal conditioning circuit and then sent into a data acquisition card for A/D conversion. Taking the voltage measuring as an example, a diagram of the signal measuring and conditioning circuit is shown in
Fig. 11
.
Diagram of signal measuring and conditioning circuit.
- D. Control and I/O Interface
The control and I/O interface of the PV simulator, which consists of a personal computer and a PCI1710 multifunction data acquisition card, can conveniently achieve the data acquisition and configuration of the control algorithm and interface by the abundant software resources of the computer.
- E. Control Software Design
The modules of the Real-Time Windows Target are loaded to become a seamless integration of the Simulink environment and the external equipment under the external mode of Simulink, which can make Simulink become a real-time configuration development environment and a graphical control platform. In order to construct and adjust the models and parameters of the PV array, sunshine intensity, PV array’s temperature and controller conveniently and to improve the software development efficiency and reliability, the control software of the PV simulator is designed and developed based on this configuration environment and control platform.
A diagram of the real-time configuration and control program for the PV simulator is shown in
Fig. 12
. To increase the readability of the program, models of the I-V characteristics and the PI controller are both encapsulated into subsystems. Sampled-data control in equal intervals of time is adopted in the system and its sampling period is 1ms.
Block diagram of real-time configuration and control program for digital PV simulator.
In order to solve the PV model conveniently, the engineering analytical model of the PV array, which can accurately describe the output characteristics of the PV array with only several electrical parameters under the standard test conditions provided by the manufactures and is convenient for engineering applications, is adopted to build the I-V characteristics subsystem in the MATLAB/Simulink environment. The engineering analytical model of the PV array can be expressed as
[21]
:
where:
DT
=
T
-
Tref
,
VOC
is the open-circuit voltage,
ISC
is the short-circuit current,
Vm
is the optimum operating voltage and
Im
is the optimum operating current, which gives the maximum power at that time.
Rs
is the series internal resistor,
Q
and
Qref
are actual sunshine intensity and the reference sunshine intensity of 1000
W/m2
, respectively.
T
and
Tref
are the actual module temperature and the reference temperature of 25
℃
, respectively.
α
and
β
are the current and voltage temperature coefficients under the reference sunshine intensity, respectively.
A diagram of the I-V characteristic subsystem is shown in
Fig.13
. It can be seen from the diagram that the models and parameters of the sunshine intensity, the PV array’s temperature, etc. can be built and adjusted conveniently.
Subsystem of I-V characteristics.
A diagram of the fuzzy PI controller subsystem is shown in
Fig.14
. In order to ensure real time control, the Takagi-Sugeno fuzzy model is realized with the Embedded MATLAB Function block of Simulink in the subsystem. The block can achieve complicated function flexibility with the MATLAB programming language and it can execute simulations and generate code for a Real-time Workshop target.
Subsystem of fuzzy PI controller.
VII. EXPERIMENT RESULTS AND ANALYSIS
- A. System Implementation
Based on the principle of modular hardware and configurable software, the functional modules of the PV simulator’s hardware which consists of the PPF circuit, the PWM signal generation, the driving and isolating circuit, the measuring and signal conditioning circuit, the control and I/O interface, and so on, are all implemented and tested. After the module level test, all of the functional modules are connected with software to form a complete system and then tested in a whole system. In this paper, a DC programmable electronic load of IT8516C by ITECH is adopted to simulate the load characteristics. The whole system is shown in
Fig. 15
, and the parameters of the PV array used to verify the performance of the PV simulator are listed in
Table IV
of the Appendix.
Experimental setup of PV simulator.
- B. Steady-state Experiment
In the experiment, the load mutation test is done to verify the steady-state performances of the different operating points. The typical waveforms of the output voltage and filtering inductor current are shown in
Fig. 16
. In the test,
RL
is changed suddenly from 20Ω to 5Ω and then changed suddenly from 5Ω to 20Ω after remaining a 5Ω for some time. The experimental waveforms indicate that the ripples of the voltage and current are both very small and that the system can work steadily at the operating point corresponding with
RL
. The PV simulator has a fast dynamic response and a high steady-state accuracy.
Waveforms of output voltage and filtering inductor current.
Based on the above results, the static output characteristics of the PV simulator are simulated under 25℃ and different sunshine intensities as well as 1000W/m
2
sunshine intensity and different temperatures. Due to limitations on the length of this paper, only the static output characteristics under 25℃ and different sunshine intensities are given. The experimental results are shown in
Fig.17
.
Static output characteristics of PV simulator under 25℃ and different sunshine intensities ((a) Current-voltage characteristics and (b) Power-voltage characteristics).
The steady-state experimental results show that the output voltage and current of the PV simulator have a wide adjusting range, and that its current-voltage characteristics and power-voltage characteristics comply well with the theoretical output characteristics.
- C. Dynamic Experiment
Due to variations of the operating point of the PV simulator with a load, the test of the power variation of the PV simulator with the sunshine intensity changed under natural conditions is done to verify the dynamic characteristic of the PV simulator. In the test, the sunshine intensity change under natural conditions is simulated by a combination of the positive and negative ramp function, and the other models. A DC programmable electronic load is selected in the constant-voltage mode to maintain a constant output voltage of the PV simulator so as to measure its output power. The curves of the sunshine intensity change and the output power for the PV simulator are shown in
Fig. 18
. The test result shows that the power of the PV simulator varies linearly with the sunshine intensity change and it has a good dynamic response.
Dynamic output characteristics of PV simulator ( (a) curve of sunshine intensity change and (b) curve of output power of PV simulator).
- D. System Efficiency
The PV simulator’s efficiency as a function of different output voltages under 1000W/m
2
sunshine intensity is measured. The efficiency curve is shown in
Fig. 19
. The efficiency curve shows that the PV simulator has a high working efficiency within a wide operating range and that it can satisfy the design requirements.
Efficiency curve of the PV simulator in 1000W/m2 sunshine intensity.
VIII. CONCLUSIONS
In this paper, the hardware of a digital PV simulator is designed and implemented based on the principle of modular hardware. The control structure of the system is based on a PC and a multifunction data acquisition card. A PPF circuit is adopted as the main circuit of the PV simulator to restrain the magnetic biasing of the core for the DC-DC converter and to reduce the spike of the turn-off voltage across every switch. Based on an analysis of the operating mode and the input-output relationship of the PPF circuit, mathematic models of the PPF circuit and PV simulator are built by the state-space averaging method. Aiming at the problem where the load resistance
RL
changing with the operating point of the PV simulator leads to the difficult controller parameter design, the inertial element of the equivalent transfer function of current measuring, conditioning and filtering is cancelled by the pole-zero cancellation technique so as to reduce the design difficulty of the controller parameters. Then the controller parameters are systematically designed based on the performance analysis of the root loci of the closed current loop with
ki
and
RL
as variables. In order to match the controller parameters with
RL
, a fuzzy PI controller based on the Takagi-Sugeno fuzzy model is applied to self-adaptively regulate the controller parameters.
The real-time configuration and the control program for the digital PV simulator are developed in the real-time environment of Matlab/Simulink which loads the modules of the Real-Time Windows Target and works under the external mode of Simulink so that the models and parameters of the system can be adjusted conveniently and the system can run in real-time. Lastly, the stationary and dynamic performances of the PV simulator are tested by experiments. The experimental results show that when the controller parameters designed in this paper according to the state-space averaging model are applied to the PV simulator, satisfactory dynamic and stationary performances are obtained. This makes it have the merits of a wide effective working range, a high steady-state accuracy and good dynamic performances. Therefore, it can be seen that the modeling and controller parameter design of the digital PV simulator based on the PPF circuit proposed in this paper is correct and reasonable. It also has a significant reference value for the modeling and controller parameter design of many other types of PV simulators and converters.
BIO
Jike Zhang was born in Hohhot, China. He received his B.S. and M.S. degrees in Control Science and Engineering from the Inner Mongolia University of Technology, Hohhot, China, in 2000 and 2003, respectively. He is presently working toward his Ph.D. degree at the Inner Mongolia University of Technology. Since 2005, he has been a Lecturer in the Department of Automation, Inner Mongolia University of Technology. His current research interests include power electronics and control in renewable energy.
Shengtie Wang was born in Hohhot, China. He received his B.S. and M.S. degrees in Control Science and Engineering from the Chengdu University of Science and Technology, Chengdu, China, in 1985 and 1988, respectively. He received his Ph.D. degree in Electrical Engineering from the Huazhong University of Science and Technology, Wuhan, China, in 1999. Since 2000, he has been a Professor in the Department of Automation, Inner Mongolia University of Technology, Hohhot, China. From September to December of 2003, he was a Researcher at Mie University, Tsu, Japan, where he was engaged in research on the control of small wind turbines. From 2007 to 2008, he worked as a Visiting Scholar at the Norwegian University of Science and Technology (NTNU), Trondheim, Norway, where he was engaged in research on the control of large wind turbines. His current research interests include power electronics and control in renewable energy, intelligent control and applications.
Zhihe Wang was born in Baotou, China. He received his B.S. and M.S. degrees in Control Science and Engineering from Chongqing University, Chongqing, China, in 1985 and 1988, respectively. From 1991 to 1996, he worked for the Automation Research Institute, Inner Mongolia University of Technology, Hohhot, China, where he was engaged in research on electric drives, motion control, and power quality control. Since 2003, he has been a Professor in the Department of Automation, Inner Mongolia University of Technology. His current research interests include power electronics and converter control.
Lixin Tian was born in Baotou, China. He received his B.S. degree in Electronic Science and Technology from Inner Mongolia University, Hohhot, China, in 1989. From 1992 to 1996, he worked for the Automation Research Institute, Inner Mongolia University of Technology, Hohhot, China, where he was engaged in research on power electronics and power quality control. Since 2001, he has been a Assistant Professor in the Department of Electric Power, Inner Mongolia University of Technology. His current research interests include power electronics and energy conversion.
Ahmed E. M.
,
Shoyama M.
2012
“Scaling Factor design based variable step size incremental resistance maximum power point tracking for PV systems,”
Journal of Power Electronics
12
(1)
164 -
171
DOI : 10.6113/JPE.2012.12.1.164
Thangaprakash S.
2012
“Unified MPPT control strategy for z-source inverter based photovoltaic power conversion systems,”
Journal of Power Electronics
12
(1)
172 -
180
DOI : 10.6113/JPE.2012.12.1.172
Jang D.-H.
,
Han S.-K.
2010
“Low cost high power density photovoltaic power conditioning system with an energy storage system,”
Journal of Power Electronics
12
(3)
487 -
494
DOI : 10.6113/JPE.2012.12.3.487
Ji S.-K.
,
Jang D.-H.
,
Hong S.-S.
2012
“Analog control algorithm for maximum power trackers employed in photovoltaic applications,”
Journal of Power Electronics
12
(3)
503 -
508
DOI : 10.6113/JPE.2012.12.3.503
Lee S.
,
Kim J.-E.
,
Cha H. J.
2010
“Design and implementation of photovoltaic power conditioning system using a current-based maximum power point tracking,”
Journal of Electrical Engineering & Technology
5
(4)
606 -
613
DOI : 10.5370/JEET.2010.5.4.606
Shen Y. L.
1997
“A photovoltaic array simulator,”
Acta Energiae Solaris Sinica
(in Chinese)
18
(4)
448 -
451
Nagayoshi H.
2004
“I-V curve simulation by multi-module simulator using I-V magnifier circuit,”
Solar Energy Materials & Solar Cells
82
(1-2)
159 -
167
DOI : 10.1016/j.solmat.2004.01.014
Li Q.
,
Li Y. D.
,
Wang K.
,
Rao J. Y.
,
Shen Z. Y.
2008
“Design and development of solar array simulator,”
Electric Drive
(in Chinese)
38
(5)
8 -
10
Bun L.
,
Raison B.
,
Rostaing G.
,
Bacha S.
,
Rumeau A.
,
Labonne A.
2011
“Development of a real time photovoltaic simulator in normal and abnormal operations,”
in Proc. IECON 2011
867 -
872
Nanakos A. C.
,
Tatakis E. C.
2008
“Static and dynamic response of a photovoltaic characteristics simulator,”
in Proc. EPE-PEMC 2008
1827 -
1833
Li Y.
,
Lee T. W.
,
Peng F. Z.
,
Liu D. C
2009
“A hybrid control strategy for photovoltaic simulator,”
in Proc. APEC 2009
899 -
903
Koran A.
,
Sano K.
,
Kim R.-Y.
,
Lai J.-S.
2010
“Design of a photovoltaic simulator with a novel reference signal generator and two-stage LC output filter,”
IEEE Trans. Power Electron.
25
(5)
1331 -
1338
DOI : 10.1109/TPEL.2009.2037501
Vijayakumari A.
,
Devarajan A. T.
,
Devarajan N.
2012
“Design and development of a model-based hardware simulator for photovoltaic array,”
International Journal of Electrical Power and Energy Systems
43
(1)
40 -
46
DOI : 10.1016/j.ijepes.2012.04.049
Di Piazza M. C.
,
Pucci M.
,
Ragusa A.
,
Vitale G.
2010
“Analytical versus neural real-time simulation of a photovoltaic generator based on a DC-DC converter,”
IEEE Trans. Ind. Appl.
46
(6)
2501 -
2510
DOI : 10.1109/TIA.2010.2072975
Wang F. S.
,
Wen J. Y.
,
Zhang X.
,
Zhao W.
2012
“An investigation on digital photovoltaic array simulator,”
Acta Energiae Solaris Sinica
(in Chinese)
33
(3)
494 -
499
Lee J.-P.
,
Min B.-D.
,
Kim T.-J.
,
Kim J.-H.
,
Ryu M.-H.
,
Baek J.-W.
,
Yoo D.-W.
,
Yoo J.-Y.
2009
“Development of a photovoltaic simulator with novel simulation method of photovoltaic characteristics,”
in Proc. INTELEC 2009
1 -
5
Martin-Segura G.
,
Lopez-Mestre J.
,
Teixido-Casas M.
,
Sudria-Andreu A.
2007
“Development of a photovoltaic array emulator system based on a full-bridge structure,”
in Proc. EPQU 2007
1 -
6
Zhang C. Z.
,
Chen Q. F
,
Zou J. Z.
2011
“Research on large capacity PV array simulator,”
Acta Energiae Solaris Sinica
(in Chinese)
32
(10)
1461 -
1465
Zhang F. H.
,
Qin H. H.
,
Wang H. Z.
,
Yan Y. G
2003
“Freewheeling current in Push-Pull Forward converter,”
in Conf. Rec. IEEE-PESC Annu. Meeting
Vol. 1
353 -
358
Zhang J. K.
,
Wang S. T.
,
Wang Z. H
2010
“Nonlinear fuzzy PID control strategy for nonlinear plant with large inertia,”
in Proc. WCICA 2010
2537 -
2541
Su J. H.
,
Yu S.-J.
,
Zhao W.
,
Wu M.
,
Shen Y.
,
He H.
2001
“An investigation on engineering analytical model of silicon solar cells,”
Acta Energiae Solaris Sinica
(in Chinese)
22
(4)
409 -
412