Linearizing and Control of a Threephase Photovoltaic System with Feedback Method and Intelligent Control in StateSpace
 Author: Louzazni Mohamed, Aroudam Elhassan
 Publish: Transactions on Electrical and Electronic Materials Volume 15, Issue6, p297~304, 25 Dec 2014

ABSTRACT
Due to the nonlinearity and complexity of the threephase photovoltaic inverter, we propose an intelligent control based on fuzzy logic and the classical proportionalintegralderivative. The feedback linearization method is applied to cancel the nonlinearities, and transform the dynamic system into a simple and linear subsystem. The system is transformed from abc frame to dq0 synchronous frame, to simplify the state feedback linearization law, and make the closeloop dynamics in the equivalent linear model. The controls improve the dynamic response, efficiency and stability of the threephase photovoltaic grid system, under variable temperature, solar intensity, and load. The intelligent control of the nonlinear characteristic of the photovoltaic automatically varies the coefficients K_{p}, K_{i} , and K_{d} under variable temperature and irradiation, and eliminates the oscillation. The simulation results show the advantages of the proposed intelligent control in terms of the correctness, stability, and maintenance of its response, which from many aspects is better than that of the PID controller.

KEYWORD
Threephase photovoltaic , Intelligent control , Fuzzy logic , Classical proportionalintegratorderivate , Feedback linearization

1. INTRODUCTION
The use of renewable energy sources is rapidly developing, and the application of solar energy focusing photovoltaic systems is becoming increasingly popular. The major challenge in photovoltaics is posed by the instability, nonlinearity, and complexity of the currentvoltage characteristic equation. In this paper, we propose an intelligent control of a threephase photovoltaic gridconnected inverter system, which is essential to consider. The behavior of photovoltaic’s in currentvoltage and voltagepower is described by complex and nonlinear analytical equations, and depends on various levels of solar intensity and various cell temperatures. The characteristic of photovoltaic’s is a complex and nonlinear function, and it is difficult to identify a dynamic model. On the other hand, the Lambert Wfunction can be used to find a mathematical equation that is capable of describing the behavior of photovoltaic model, including all related parameters. The lambert Wfunction has been widely used in many publications and fields to reduce the complexity, and describe the behavior of a photovoltaic system, to convert the nonlinear and complex equation into explicit form. Ref. [1] has given many available applications of the Lambert Wfunction, while Ref. [2] used it to develop an analytical compact model for the asymmetric lightly doped MOSFET. Furthermore, the feedback linearization method and state feedback law to transform and eliminate the nonlinearity model of photovoltaic system into a simple equivalent linear model, were used to find a direct relation between the output and the control input, by using inverse dynamics [3,4]. This nonlinear state model transforms the dq0 synchronous frame reference into an equivalent linear system, where the pole placement control loop technique is applied to separate the control, and place the closedloop system pole in the desired location. The system is composed of a photovoltaic array, capacitive DClink, and threephase inverter connected to the grid, which is assumed to be in phase with the inductor current.
The classical PID controller is widely used in 95% of industrial technology. In fact, to assure the stability of the threephase photovoltaic gridconnected with L filters, intelligent adaptive fuzzy logic control and classical PID control using the feedback linearization method, is used to ensure a high quality system. This intelligent control is based on the combination of classical PID, and fuzzy logic controller, and the principle is to find out the fuzzy relationship between the three parameters of PID, error, and error change, by automatically varying the coefficients, K_{p}, K_{i} and K_{d} with variations of radiation, temperature, and load change. This controller combines the advantage of fuzzy control and PID control, and has a high performance under rapidly changing atmospheric conditions. Its effectiveness depends on the experience or knowledge of the right rule, input and output variables, and the membership functions.
Finally, the results obtained from the classical PID controller and intelligent PIDFuzzy logic with feedback linearization method show the advantages of the proposed intelligent control system in terms of its correctness, and the feasibility, stability, and maintenance of its response, which is better than the classical PID controller form many aspects, and which has the advantages of being fast, robust and showing good performance under varying atmospheric conditions.
2. LAMBERT WFUNCTION OF PHOTOVOLTAIC
The Lambert Wfunction is used to find a mathematical equation capable of describing the behavior of a photovoltaic model, including all related parameters. The Lambert wfunction has been widely used in many publications and fields to reduce the complexity and describe the behavior of a photovoltaic system, by converting the nonlinear and complex equation into explicit form [5]. Figure 1 represents a photovoltaic model with a single diode; under illumination the relation between the current and voltage of the photovoltaic single diode is given by [6]:
where, V_{th}=K.T/q is the thermal voltage of the photovoltaic module; I_{ph} is the irradiance current (photocurrent); I_{0} is the cell reverse saturation current (diode saturation current); q is the electron charge; n is the cell ideality factor; K is the Boltzmann constant; T is the cell temperature; and R_{s} and R_{p} represent the cell series and shunt resistance, respectively.
The Lambert Wfunction is the inverse function of
f(w)=wew
Equation (1) is transcendental in nature, and the explicit form
Similarly, the output power of the photovoltaic array can be given in explicit terms, as follows:
Figure 2 presents the characteristic currentvoltage and powervoltage of the photovoltaic cell for different values of temperature and solar incident irradiance.
In the electrical characteristics, two points exist for a specific curve that can generate more power than other points. The maximum output power varies with temperature and irradiation change, so it necessary to extract the maximum available power for any changes.
3. STATE FEEDBACK MODEL OF THE THREEPHASE PHOTOVOLTAIC
Figure 3 shows the general circuit topology configuration of a threephase photovoltaic grid connected inverter connected to the grid connected to the inductance. The system is composed of a photovoltaic array, threephase inverter and capacitive DClink. The photovoltaic array converts solar irradiation into DC current, and the DClink capacitor reduces the harmonic of the DC voltage on the input side of the inverter.
The threephase inverter model of Fig. 3 given in states space coordinates by:
Park’s transformation is applied to equation (5)
We apply the dq transformation to equation (6), and obtain the statespace dq in equation (7)
where, e_{d}, e_{q}, e_{0}, i_{d}, i_{0}, v_{d}, v_{q} and v_{0} are the components of the grid voltage, grid current, and inverter output voltage, respectively, and ω is the angular frequency. The power losses are neglected in inverter switches, and the power of the DCinput and ACoutput is given by:
Where, v_{dc} and i_{dc} are the input voltage and current of the inverter, respectively.
Kirchhoff’s law gives
From equations (7) and (9), the equation of the state model can be grouped as:
3.1 State feedback linearization controller
The main objective of feedback linearization control is to transform and eliminate the nonlinearity model of the photovoltaic inverter system into a simple equivalent linear model, and find a direct relation between the output and the control input, by using inverse dynamics [7,8].
Let us consider a nonlinear control system of the phase state variable vector and phase input vector. The inputoutput state equation of the threephase system is
where,
The outputs are differentiable until the inputs appear for searching for the exact input output. Let r_{i} be a small integer such that at least one of the inputs appears in the r_{i} th derivation of y_{i} . The method is detailed in [9,10].
When
For the first output:
We have u_{2} input from the first derivation and the relative degree is r_{1}=1.
For the second output:
The inputs u_{1} and u_{2} appear after the second derivation and the relative degree is r_{1}=1.
The sum of r_{1}+r_{2}=3 and this is the order of the system. Then we can give an exact linearization of the system.
From equations (14) and (16), the feedback linearization of the system is
where:
Where, β(x) is a non singular matrix, whose determinant is:
The corposant e_{d} of the gridconnected voltage is always different from zero. Then the determinant is not null, and β(x) is nonsingular. The linearization of the control part is
To eliminate the nonlinearity of the system, we substitute equation (21) into (18). So the simple linear relation between outputs yi and the new inputs v_{i} is:
The nonlinearity of the system is cancelled, and Fig. 4 presents the linear equivalent control block diagram of the proposed method.
The controller is designed to obtain the stability of an integral controller, which is added to eliminate the steadystate error due to parameter variations.
where, e_{1}=y_{1}−y_{1ref}, e_{2}=y_{2}−y_{2ref} and y_{ref} is the tracking error, y_{1ref} is the tracking reference of the grid current reference, y_{2ref} the DClink voltage reference, and k_{ij} the gain. Then, the output error is:
The controller gains k_{ij} are determined by Rollscriterion and to assure the characteristic polynomial of the equation (24), these are Hurwitz polynomials, so the error trackings of e_{1} and e_{2} converge to zero.
4. DSIGN OF CONTROLLERS
4.1 Conventional PID controller
The conventional fundamental classical PID controller is often described by the following equation in ideal form for the output voltage of a grid photovoltaic system [1113]:
where, e_{1}(t)=y_{1}−y_{1ref}, e_{2}(t)=y_{2}−y_{2ref} and e_{3}=(t)y_{ref} is the tracking error, in which V_{out}(t) denotes the controller output at time t of the photovoltaic array. K_{p}, K_{i} and K_{d} are known as the proportional constant, integral constant, and derivative constant respectively.
4.2 Structure of the PIDfuzzy selftuning controller
Figure 5 shows the structure of the control system with the proposed PIDFuzzy logic controller. The principle of fuzzy selftuning PID is first to find the fuzzy relationship between the three parameters of PID and error e_{i} (t)(i=1, 2 and 3) and the variation of the error (i=1, 2 and 3), and adjust the PID parameters with fuzzy rules. Fuzzy inference engines modify three parameters, to be content with the online demands of the control system. The inverse reference model of the response of the controlled object is used in the approach to tuning the PIDfuzzy controller.
The input parameters are the deviation e_{i}(t) and the variation of deviation Δe_{i}(t), which are shown in equation (26). The output is the variation of MPPT.
The output parameters are K_{p}, K_{i} and K_{d}, and continuously detect e_{i}(t) and Δe_{i}(t) according to the fuzzy logic rules and regulations in order to satisfy the desire of the two inputs to the parameters of the controller at any time [14].
4.3 Fuzzy gain scheduling
According to the methodology proposed by Zhao [15], the three current PID parameters are determined as follows. We assume the range of K_{p} is [K_{p,min} K_{p,max}], the range of K_{i} is [K_{i,min} K_{i,max}], and the range of K_{d} is [K_{d,min} K_{d,max}]. The range of K_{p} and K_{d} can be normalized into the range, through the linear transformation in equations (27) and (28).
where, K'_{p}, K'_{d,} and α are constants that are determined by means of the fuzzy mechanism. The K_{p,min}, K_{p,max}, K_{d,min}, K_{d,max} are constants adopted to normalize the values of K_{p} and K_{d} which are given by :
The parameters K_{p}, K_{i} , and α are determined by a set of fuzzy rules.
If e_{i}(t) is A_{i} and Δe_{i}(t) is B_{i}, then K'_{p} is C_{i} , K'_{d} is D_{i}, and α=α_{i} (i=1, 2…m). Here, A_{i}, B_{i} , C_{i}, and D_{i} are fuzzy sets on the corresponding supporting sets, while α_{i} is a constant.
The fuzzy rules variables values of the input e_{i}(t) and Δe_{i}(t) are configured for seven membership functions with NB= Big Negative, NM=Medium Negative, NS=Small Negative, ZO=Zero, PS= Small Positive, PM=Medium Positive, and PB=Big Positive, in the format of fuzzy control rules:
{NB, NM, NS, ZO, PS, PM, PB}
with the truth value of {0,6}
Figure 6 shows the memberships of the input variables e_{i}(t) and Δe_{i}(t).
The output K'_{p} and K'_{d} use two exponential membership functions; to calculate K'_{p} and K'_{d}, we need to correct the factors K'_{p} and K'_{d} in the range. The values of the linguistic variable of K'_{p} and K'_{d} are shown in Figs. 9 and 10.
The fuzzy sets C_{i} and D_{i} may be either Big or Small and are characterized by the grade of the membership functions μ, and the variables μ=(K'_{p} or K'_{d}) have the following equations:
For Small:
For Big:
To calculate the integral coefficient K_{i} , correction of the factor is needed, and considered as a fuzzy number, it has a singleton membership function represented by points in Fig. 9, and has 7 linguistic variables, whose values are {α.._{NB}, α_{NM}, α_{NS}, α_{ZO}, α_{PS}, α_{PM}, α_{PB}}, their value being 2, 3, 4, 5.
Through the analysis of the characteristic curve of a photovoltaic array, we can track the fuzzy control rules and correction factor of the classical PID controller parameter, for which e_{i}(t) and Δe_{i}(t) are in different conditions. The fuzzy control rules are shown in Tables 1, 2, and 3.
5. SIMULATIONS AND RESULTS
Theoretical analysis of the proposed hybrid control of the output current and voltage of the threephase photovoltaic system is validated and done by simulation using the Simulink platform. The simulation started when the solar irradiance was at G=1,000 W/m² and T=25℃. Firstly, the threephase gridconnected photovoltaic inverter using the proposed hybrid controller is studied, in comparison with the separated classical PID controller. The inverter is connected to the DC source through a DClink capacitor, with constant voltage reference value used to emulate the output of the maximum power point controller. The DClink capacitor voltage is regulated and controlled by a PIDFuzzy logic control loop, to reduce the oscillation in P_{r} at twice the line frequency due to grid imbalances. The output is the reference power, which must be injected into the grid by the inverter. The simulation control of the sinusoidal output current and voltage obtained with the DClink capacitor replaced by a constant voltage source V_{link}=400 V. The objective is to analyze the behavior of the current and voltage controller of Figs.10, 11, and 12, without influence of the DClink and of the rest the system.
Figures 10, and 11 show the current and voltage, respectively, that are synthesized by the PID controller, and injected into the grid. Figure 14 shows the plot voltage at the load.
Due to the structure of the threephase photovoltaic system converter, the output current and voltage are disturbance and unbalanced operating condition as shown in Figs.10,11, and 12.
The objective now is to present the results of the simulation with all system parts of the photovoltaic array and DCAC converter working together, relative to the PIDfuzzy controller.
From t = 0 s to t = 0.16 s the system drains energy from the AC grid. From t = 0.16 s, the converter injects active power into the grid. At 0.16 seconds the reference current is stepped, and after this time, the converter operates in a steady state condition. Overall, the controller provides an excellent dynamic response.
Figure 16 shows the behavior of the DClink voltage, which initially is not charged. The photovoltaic capacitor is charged from 0 V to 700 V from the photovoltaic array. During the initial charge, the photovoltaic array supplies its maximum current, approximately near 0.16 s; and when the capacitor is charged, the DCAC starts to supply current to the grid. The output sinusoidal current draining power to charge the capacitor suffers a phase inversion near t=0.16 s, and from this time, the DCAC converter begins to deliver active power to the grid, as in Figs. 14 and 15.
6. CONCLUSIONS
This paper presented an intelligent hybrid based on PID control, and a fuzzy logic controller was proposed to control and stabilize the threephases of the output current and voltage of the photovoltaic gridconnected converter. The system is complex and nonlinear, and a Lambert Wfunction and feedback linearization method was proposed. In the analysis, a fully mathematical approach was used to gain insight into the behavior of the photovoltaic model related to the onediode equivalent circuit. The feedback linearization method and state feedback law transform and eliminate the nonlinearity model of the threephase photovoltaic into a simple linear model. Nowadays, it has become evident that complex and nonlinear problems need intelligent systems to be solved. Recently, more than 95% control loops still use the PID controllers. They are used with fuzzy logic control to overcome the problems in the unknown mathematical model of the system. A PIDFuzzy logic hybrid controller was applied to the control power system. Furthermore, fuzzy logic controller is implemented as a gain scheduler to automatically tune the K_{p}, K_{i} and K_{d} parameters for the PID controller. The system was regulated and controlled, and the results show the active power reduced in oscillation, and the current and voltage stabilized. As the results confirmed, initially the DClink voltage capacitor is not charged, and during the charge, the threephase photovoltaic system supplies a maximum current; and when the capacitor charges, the DCAC starts to supply the grid. Compared with conventional controllers, this ensures minimum peak values in the gridinjected currents during unbalanced voltage sags. The proposed intelligent hybrid controller gave fast response and stability without delay, eliminating the oscillation around the operating points. From the simulation results, we concluded that the performance of the intelligent controller is better than the ones obtained with the classical PID controller, since the response time in the transitional state was shortened, and the fluctuations in the steady state were considerably reduced.

1. Corless R. M., Gonnet G. H., Hare D. E. G., Jeffrey D. J., Knuth D. E. 1996 [J. Adv. Comput. Math.] Vol.5 P.329

2. Abebe H., Cumberbatch E., Morris H., Tyree V., Numata T., Uno S. 2009 [J. Semiconduct. Technol. Sci.] Vol.9 P.225

3. Arifa J., Rayb S., Chaudhuria B. 2013 [International Journal of Electrical Power & Energy Systems] Vol.45 P.87

4. Johansen T. A., Hunt K. J. 2000 [IEEE Publication in Decision and Control Conference] Vol.5 P.4467

5. Hwang C., Cheng Y. 2005 [American Control Conference] Vol.6 P.4283

6. Louzazni M., Aroudam E., Yatimi H. 2013 [International Journal of Electrical and Computer Engineering] Vol.3 P.568

7. Charlet B., Levine J. 1989 [Systems & Control Letters] Vol.13 P.143

8. Van Der Schaft A. J. 1984 [Systems & Control Letters] Vol.5 P.27

9. Tall I. A. 2010 [Systems & Control Letters] Vol.59 P.429

10. Monaco S., NormandCyrot D. 2006 [Systems & Control Letters] Vol.55 P.71

11. Chaouachi A., Kamel R. M, Nagasaka K. 2010 [Sol. Energy] Vol.84 P.2219

12. Louzazni M., Aroudam E. 2014 [IEEE Publication, Intelligent Systems: Theories and Applications (SITA14), 2014 9th International Conference on] P.16

13. Louzazni M., Aroudam E. 2014 [IEEE publication, Intelligent Energy and Power Systems (IEPS), 2014 IEEE International Conference on] P.279284

14. Wang Z., Gong Z., Wang Y., Wei G. 2009 [IEEE 9th International Conference on Electronic Measurement & Instruments] P.913

15. Zhao Z., Tomizuka M., Isaka S. 1993 [IEEE Transactions on Systems Man and Cybernetics] Vol.23 P.1392

[Fig. 1.] Photovoltaic model.

[]

[]

[]

[]

[Fig. 2.] Photovoltaic cell electrical curves with constant temperature and different values of irradiance.

[Fig. 3.] Diagram of threephase gridconnected photovoltaic inverter.

[]

[]

[]

[]

[]

[]

[]

[]

[]

[]

[]

[]

[]

[]

[]

[]

[]

[]

[Fig. 4.] Feedback linearization control block diagram.

[]

[]

[]

[Fig. 5.] The structure of the PIDFuzzy logic selftuning control.

[]

[]

[]

[]

[]

[Fig. 6.] Memberships of the input variables ei(t) and Δei(t).

[Fig. 7.] Membership of K'p.

[Fig. 8.] Membership of K'd.

[Fig. 9.] Membership of α.

[Fig. 10.] Threephase output voltage va,b,c of gridconnected using PID.

[]

[]

[Table 1.] The fuzzy control rules for K'p.

[Table 2.] The fuzzy control rules for K'd.

[Table 3.] The fuzzy control rules for α.

[Fig. 11.] Threephase output current ia,b,c of gridconnected using PID.

[Fig. 12.] Plot of the DClink voltage vDC.

[Fig. 13.] Threephase output voltage va,b,c of the photovoltaic gridconnected converter.

[Fig. 14.] Plot of the DClink voltage vDC using the proposed control

[Fig. 15.] Threephase output currents ia,b,c of the photovoltaic gridconnected converter.

[Fig. 16.] Plot of the DClink voltage.

[Fig. 17.] Output power of the photovoltaic array delivered to the AC grid.