Complexity

Complexity / 2019 / Article

Research Article | Open Access

Volume 2019 |Article ID 6141607 | 15 pages | https://doi.org/10.1155/2019/6141607

Super-Twisting Sliding Mode Control for Gearless PMSG-Based Wind Turbine

Academic Editor: Hassan Zargarzadeh
Received22 Jan 2019
Revised26 Mar 2019
Accepted11 Apr 2019
Published22 Apr 2019

Abstract

In recent years, the complexities of wind turbine control are raised while implementing grid codes in voltage sag conditions. In fact, wind turbines should stay connected to the grid and inject reactive power according to the new grid codes. Accordingly, this paper presents a new control algorithm based on super-twisting sliding mode for a gearless wind turbine by a permanent magnet synchronous generator (PMSG). The PMSG is connected to the grid via the back-to-back converter. In the proposed method, the machine side converter regulates the DC-link voltage. This strategy improves low-voltage ride through (LVRT) capability. In addition, the grid side inverter provides the maximum power point tracking (MPPT) control. It should be noted that the super-twisting sliding mode (STSM) control is implemented to effectively deal with nonlinear relationship between DC-link voltage and the input control signal. The main features of the designed controller are being chattering-free and its robustness against external disturbances such as grid fault conditions. Simulations are performed on the MATLAB/Simulink platform. This controller is compared with Proportional-Integral (PI) and the first-order sliding mode (FOSM) controllers to illustrate the DC-link voltage regulation capability in the normal and grid fault conditions. Then, to show the MPPT implementation of the proposed controller, wind speed is changed with time. The simulation results show designed STSM controller better performance and robustness under different conditions.

1. Introduction

Recently, wind energy is a rapidly growing source of electricity [1, 2]. Installed wind power generation is expected to exceed 760 GW by 2020 in the world [3]. Therefore, the significant amount of power generation capacity belongs to wind power. In the past, wind farms could be disconnect from the grid in the grid side fault conditions, but nowadays, the modern grid codes, in addition to not allowing the separation of the grid to wind turbines, require them to perform some auxiliary tasks in grid fault conditions. For instance, Danish grid code forces wind farm to follow low-voltage ride through (LVRT) requirements in fault condition according to Figure 1 [4].

Among different types of wind turbines, doubly fed induction generators (DFIGs) and permanent magnet synchronous generators (PMSGs) are of attractive types in wind farms [5]. Due to advantages, developments in semiconductor switching devices and high-reliability and efficiency, the use of PMSG is growing.

In traditional control of back-to-back converter, machine side converter (MSC) performs the maximum power point tracking (MPPT) and the grid side converter (GSC) regulates the DC-link voltage [6, 7]. Hence, when grid faults occurred, the DC-link voltage may be increased because MSC does not sense grid fault and GSC loses its own control on DC-link voltage regulation [8]. To overcome this problem, the GSC and MSC tasks can be changed [9, 10]. Due to the nonlinear relationship between DC-link voltage and input control signal, different nonlinear controls are presented [11, 12].

Among the modern nonlinear control, sliding mode control has good feature because it is robust with respect to system parameter uncertainties and external disturbances [13]. The main drawback of conventional sliding mode control is chattering. To reduce chattering, several methods are introduced. One of the attractive chattering-free sliding modes is high order sliding mode (HOSM) [14, 15]. But the main problem in the implementation of HOSM algorithms is the increasing information demand. The super-twisting sliding mode (STMC) is one of HOSM that does not need additional information [16, 17].

Several papers use the HOSM to enhance LVRT capability in PMSG-based wind turbines [1820]. Nevertheless, all of these works use traditional control of back-to-back converters. In fact, they cannot keep the DC-link voltage in safe range in deep voltage drops. Therefore, they need to use external devices which result in the overall cost increases. As mentioned before, by exchanging the task of converters, the DC-link voltage is kept in safe range, MSC regulates DC-link voltage, and GSC controls MPPT. In present paper, the super-twisting sliding mode control method is employed in the new control structure.

It can be mentioned that the new contribution of the proposed method is implementing a HOSM controller for the gearless high inertia PMSG-based wind turbine, in which no additional instrument is required to implement LVRT without imposing stress to the DC-link capacitor. In fact, the grid voltage drop and the wind variations are considered as an external disturbance. Furthermore, the chattering is eliminated and the controller does not need additional information.

In this paper, wind energy conversion system model, including modeling of the wind turbine, PMSG, DC-link, and the grid is presented in Section 2. In Section 3, first, the brief review of STSMC is presented; then full STSM controllers for back-to-back converters are designed. Section 4 proves the improvement of the new control strategy by comparing with the conventional approach on the MATLAB/Simulink. Finally, in Section 5, conclusions are presented.

2. Description of Wind Energy Conversion System

The schematic of the PMSG wind power system is shown in Figure 2. The modeling of each section has been introduced in the following subsections.

2.1. Modeling of Wind Turbine

The mechanical output power of wind turbine is given by the following equation:The turbine power coefficient is defined by the following equations [21]:The tip-speed ratio (λ) depends on the shaft speed () and the wind speed as given below:From (1) and (4), the turbine torque is given asThe equation of motion of single-mass modeling of the mechanical system is expressed aswhere is the electromagnetic torque of the generator.

2.2. Modeling of PMSG

A dynamic mathematical model of a surface-mounted PMSG is presented in the synchronous d-q equivalent circuits as [22]

The electromagnetic torque from the generator is given as

2.3. Modeling of DC-Link

The DC-link is the interface between the generator and the grid. The dynamic DC-link voltage equation can be expressed bywhere and are the generator and the back-to-back converters losses and the grid power, respectively. In this study, it is supposed that converters are lossless.

2.4. Modeling of Grid

The GSC is connected to the point of common coupling (PCC) by RL grid filter and a coupling transformer. In this study, it is modeled by RL that the state equations in the d-q reference frame can be expressed aswhere ωf is the angular frequency.

The equations of instantaneous injected active and reactive powers to the PCC are expressed as

If the PCC voltage space vector is oriented on d-axis, , then

Now, let us introduce the state variables , , , , , .

The state-space equations obtained up to now are put together to obtain a state-space model of the whole system. It is supposed that the MSC and GSC are ideal and they produce the desired voltages. In fact, the average model is used for GSC and MSC. For convenience, the state-space models of the whole system are rewritten as

3. Control of Back-to-Back Converters

In this section, new controllers for the back-to-back converters are designed using second-order sliding mode. Among several proposed modern techniques, the sliding mode provides more advantages such as robustness and high-accuracy solution especially for nonlinear systems under uncertainty conditions and external disturbances [23]. The main disadvantage of the classic sliding mode is the chattering problem. Though, the second-order sliding mode is appropriate to obtain chattering-free control and guarantees a finite-time reaching phase [24]. In Section 3.1, second-order sliding modes are briefly introduced, then, in Section 3.2, super-twisting controller that is one of attractive second-order sliding mode controllers will be presented.

3.1. A Brief Introduction of Second-Order Sliding Modes

As mentioned before, the main disadvantage of first-order sliding mode control (FOSMC) is a chattering problem. This effect consists of the oscillation of the system variables around the sliding surface, which causes a discontinuous control signal. This effect can disturb or damage the physical system. One of the interesting ways to eliminate or reduce chattering effect is using HOSM. In this technique, the higher order time derivatives of sliding surface keep to zero [24]. As a result, this action mitigates the chattering effect. But, the increasing information demand is the main problem in the implementation of HOSM algorithms. Generally, r-th order sliding controller requires the knowledge of (r-1)-th order of time derivative. Thus, a r-th order sliding mode is determined by

Hence, the second-order sliding mode controller is simple. In fact, it needs the knowledge of the first-order time derivative of the sliding surface. These controllers require real-time measurement of or at least of .

Among several second-order sliding mode controllers, the super-twisting controller has a good feature because it does not need the knowledge of . The super-twisting controller can be used instead of the first-order (conventional) sliding mode using the same available information.

3.2. Super-Twisting Controller

Consider the nonlinear dynamic system given by [24]where , , and the smooth functions , , are unknown. From (24), the time derivative of the sliding surface can be expressed aswhere and are smooth functions that for some positive constants , , , , and the following relations can be held.and the input control signal can be defined aswhere and should be satisfied with the following conditions:

Figure 3 shows the phase portrait in the super-twisting controller.

3.3. Generator Side Converter Controller

As mentioned in the Introduction, to suppress DC-link overvoltage in a grid fault condition and implement LVRT, DC-link voltage is controlled by MSC. In this method, GSC implements MPPT.

The reference value of d-axis current of the generator is zero to reduce the copper loss and to avoid demagnetization of permanent magnets. The q-axis voltage of generator (q-axis input control signal) is produced by DC-link voltage surface.

3.3.1. D-Axis Controller Design

The goal of the d-axis controller is keeping the d-axis current of the generator to zero. Hence, the control object can be expressed by the sliding variable:where is a positive constant. By taking the first derivative of and using (18), we have

Suppose that the control signal, , can be divided into known and unknown terms. Equation (30) can be rewritten aswhere is known term of the controller and is the unknown term of the controller. By selecting as the following form: and substituting (32) into (31), we have where

According to (33), we can use super-twisting sliding mode control. Hence, the unknown term of the d-axis controller can be expressed bywhere the values of control parameters are given in the Appendix.

3.3.2. Q-Axis Controller Design

The q-axis controller is regulating DC-link voltage in suitable constant value in all conditions such as wind variation condition and grid fault condition that the DC-link voltage may be varied. Also, the relationship between DC-link voltage and the q-axis controller signal is nonlinear. Hence, the control object can be expressed by the sliding surface:

Let us introduce sliding variable aswhere is derivative of . and are positive constants. By dividing known and unknown terms and taking the derivative of and substituting (19) and (20) in it, we have: where

According to (39), we can use super-twisting sliding mode control. Hence, the unknown term of q-axis controller can be expressed by

As a result, q-axis control law becomes

3.4. Grid Side Converter Controller

The GSC has two main tasks: first, MPPT control of wind turbine, second, reactive power supporting the grid in different situations such as grid voltage drop. The d-axis controller performs MPPT control and q-axis controller supports reactive power.

3.4.1. D-Axis Controller Design

As mentioned above, MPPT is implemented by the d-axis controller of GSC. There are several approaches to implement MPPT that among them optimal power control (OPC) can be applied to PMSG-based wind turbines [25]. In this design, the OPC method is used.

Considering (4) and (5) and replacing the optimal value of and the wind turbine mechanical torque can be given as a function of , as follows:where .

The optimal mechanical power of the turbine can be obtained by multiplying both sides of (43) by . Thus

The PCC power reference can be obtained by subtracting the losses of optimal mechanical wind turbine power aswhere is the RL filter loss. From (15) and (45), the d-axis reference current is obtained as

It should be noted that, in voltage drop conditions, Danish grid code enforces wind power plants to inject reactive current to PCC according to Figure 1. Hence, the d-axis current should be limited to produce reactive current by the q-axis controller because the current capacity of power electronic converter is limited. The value and upper limit of d-q axis currents in different conditions becomeswhere , , and are the per unit of rated current of the converter, the per unit of the upper limit of d-axis current, and the per unit of the upper limit of q-axis current, respectively. , , and are the per unit of d-axis voltage in PCC, the per unit of d-axis current, and the per unit of q-axis current, respectively. Now, by calculating the d-axis current from (46) and applying (47), the d-axis controller can be designed.

The control objective can be expressed by the sliding variable:where is a positive constant. By taking the first derivative of and using (21), we have

Suppose that can be divided into known and unknown terms. Equation (49) can be rewritten aswhere is known term of the controller and is an unknown term of the controller. By selecting as the following form and substituting (51) into (50), we havewhere ,

According to (52), we can use super-twisting sliding mode control. Hence, the unknown terms of the d-axis controller can be expressed by

As a result, the q-axis control law becomes

3.4.2. Q-Axis Controller Design

The main task of the q-axis controller is reactive power supporting the grid in different conditions. The q-axis current reference can be produced by reactive power demand of PCC by (16) in normal condition and by (47) in the grid voltage drop conditions. The design procedure is similar to the d-axis controller design. Hence, to avoid repeating, it is written only sliding variable and q-axis control low equation.

As a result, q-axis control law becomes

4. Simulation Results

To evaluate the performances of designed super-twisting sliding mode control in the PMSG-based wind turbines, the several simulations have been carried out on the MATLAB/Simulink software. Then this method is compared with the PI controller [8] and FOSMC [13]. The parameters of a 1.5 MW PMSG-based wind turbine and grid characteristics are given in Table 1.


Turbine parameters PMSG and Grid parameters
ParameterValueParameterValue

36.63.174
1.2253.07
48000000.23
2007.017
8.1p80
0.481500

4.1. Operation with Symmetrical Grid Faults

As mentioned above, the relationship between the DC-link voltage and the generator q-axis input control signal is nonlinear. In controlling and regulating the DC-link voltage at a constant value in all conditions such as grid fault situation, three controllers (PI, FOSMC, and STSMC) are compared. To evaluate the performances of these controllers, symmetrical voltage drop similar to the Danish grid code (Figure 1) at t=5s is applied as shown in Figure 4.

Figure 5 shows the DC-link voltage variation with three controllers. As shown in Figure 5, the STSMC has the best performances in DC-link voltage control in both starting time and grid voltage drop. According to Figure 5(a), the FOSMC has higher chattering in all time. In addition, STSMC has the fast starting time (Figure 5(b)). Furthermore, in the PI controller, the peak of the DC-link voltage reaches 1556 V while in STSMC, it is less than 1510 V (Figure 5(c)).

Figure 6 shows the q-axis current of the generator. This component of current controls the DC-link voltage. As shown in Figures 6(a) and 6(b), the STSMC has fast response in the grid voltage drop. Hence, this subject can reduce the DC-link overvoltage. The FOSMC has chattering and it can fatigue the DC-link capacitor.

The d-axis current of the GSC is shown in Figure 7. According to Figure 7, PI controller has worst performance in the grid fault condition.

In Figure 8, q-axis current of the GSC is shown. Although all of the controllers satisfy the injecting reactive current according to the grid code, similar to Figure 7, the PI controller does not have good performance.

As mentioned above, one of the main drawbacks of FOSMC is chattering especially in the input of controllers that results in chattering in systems variables. Figure 9 shows the input control signals to the MSC and GSC in the FOSMC. As shown in Figure 9, the amplitude of chattering in controller signals is extremely high. To overcome this drawback, the STSMC is suitable solution. Figure 10 shows input control signals to the MSC and GSC in the super-twisting controller. As it is seen, the chattering in the input control signals is reduced.

The phase portraits corresponding to the reaching phase of the d-q axis sliding variables of the generator side and grid side controllers in the STSMC are depicted in Figure 11. These plots (Figures 11(a)11(e)) show the features illustrated by super-twisting controllers during the system convergence toward .

To show the compliance of the grid code requirements in reactive current injection by the designed system, the injected reactive current versus grid voltage drop is shown in Figure 12. As shown in Figure 9, the designed controller injects reactive current similar to Figure 1(b).

4.2. Operation with Asymmetrical Grid Fault

To evaluate the performance of proposed controller in asymmetrical voltage drop, a single-phase voltage drop is simulated in phase A. Figure 13 shows the PCC voltage in single-phase voltage drop.

Figure 14 shows DC-link voltage variation in single-phase voltage drop. As shown, there is low amplitude fluctuation in the DC-link voltage in asymmetrical voltage drop. However, the STSMC has good performance in this condition.

4.3. Operation in Normal Condition

To emulate the performance of the super-twisting controllers in normal condition, a wind velocity profile is defined as expressed in [22]:

The amplitude and frequency of wind velocity components are given in Table 2. The resultant wind velocity profile is as shown in Figure 15(a). Figure 15(b) shows wind turbine power (Pm) and injected active power (Pgrid) to PCC. Figure 15(c) shows the power coefficient of the wind turbine. The average of Cpis 0.4493 from t=5s to t=40s. Of course, when the wind speed increases from 12 m/s the turbine power reaches the nominal value. Hence, the pitch angle controller will be active to set turbine power in nominal power as shown in Figure 15(d). As a result, Cp is slightly less than 0.48. It is clear that with the STSMC, also, MPPT is perfectly achieved. As shown in Figure 15(e), the DC-link voltage is fixed in 1500 V. This shows the robust DC-link voltage against wind fluctuations.


1234567

0.70.10.70.50.40.30.2
0.02π0.31π1010π50

5. Conclusions

This paper presents a super-twisting sliding mode control for gearless PMSG-based variable speed wind turbines. First, the mathematical models of the wind turbine, PMSG, DC-link, and grid are represented. Then the super-twisting controllers are designed for generator and grid side converters. The proposed controller is compared with the proportional-integral (PI) controller and first-order sliding mode controller (FOSMC). By applying the proposed controller, the DC-link overvoltage is significantly reduced rather than the PI controller. Furthermore, chattering is reduced in input controller signals and DC-link voltage. Also, the reactive current injection is done according to modern grid codes. In order to evaluate the performance of the super-twisting controllers in unbalanced voltage drop, it has simulated a single-phase voltage drop. Furthermore, designed controller can greatly do the MPPT and regulates the DC-link voltage in reference value in normal condition. The robustness and effectiveness of the designed super-twisting controllers in different conditions are confirmed by simulation results.

Appendix

Nomenclature

Symbols
:Air density ()
:Radius of blade (m)
:Pitch angle (°)
:Tip-speed ratio
:Power coefficient
:Wind speed ()
:Total equivalent inertia ()
:Viscous friction coefficient )
:Number of poles
:Shaft speed ()
:Electrical angular velocity ()
:Resistance ()
:Inductance ()
:DC-link capacitance ()
k:Positive gain
:Grid side inverter voltage
:Voltage ()
:Current ()
:Flux ()
Subscripts
d, q:Direct and quadrature components
s:Stator of machine
f:Grid side
dc:DC-link
Superscripts
ref:Reference.

Data Availability

The MATLAB files data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that they have no conflict of interest.

References

  1. P. Michalak and J. Zimny, “Wind energy development in the world, Europe and Poland from 1995 to 2009; Current status and future perspectives,” Renewable & Sustainable Energy Reviews, vol. 15, no. 5, pp. 2330–2341, 2011. View at: Publisher Site | Google Scholar
  2. J. Chen, W. Yao, C. Zhang, Y. Ren, and L. Jiang, “Design of robust MPPT controller for grid-connected PMSG-Based wind turbine via perturbation observation based nonlinear adaptive control,” Journal of Renewable Energy, vol. 134, pp. 478–495, 2019. View at: Publisher Site | Google Scholar
  3. F. Blaabjerg and K. Ma, “Wind energy systems,” Proceedings of the IEEE, vol. 105, no. 11, pp. 2116–2131, 2017. View at: Publisher Site | Google Scholar
  4. F. Blaabjerg, Y. Yang, D. Yang, and X. Wang, “Distributed power-generation systems and protection,” Proceedings of the IEEE, vol. 105, no. 7, pp. 1311–1331, 2017. View at: Publisher Site | Google Scholar
  5. C.-N. Wang, W.-C. Lin, and X.-K. Le, “Modelling of a PMSG wind turbine with autonomous control,” Mathematical Problems in Engineering, vol. 2014, Article ID 856173, 9 pages, 2014. View at: Publisher Site | Google Scholar
  6. S. Li, T. A. Haskew, R. P. Swatloski, and W. Gathings, “Optimal and direct-current vector control of direct-driven PMSG wind turbines,” IEEE Transactions on Power Electronics, vol. 27, no. 5, pp. 2335–2337, 2012. View at: Publisher Site | Google Scholar
  7. P. Xing, L. Fu, G. Wang, Y. Wang, and Y. Zhang, “A compositive control method of low-voltage ride through for PMSG-based wind turbine generator system,” IET Generation, Transmission & Distribution, vol. 12, no. 1, pp. 117–125, 2017. View at: Google Scholar
  8. M. Nasiri, J. Milimonfared, and S. H. Fathi, “A review of low-voltage ride-through enhancement methods for permanent magnet synchronous generator based wind turbines,” Renewable & Sustainable Energy Reviews, vol. 47, pp. 399–415, 2015. View at: Publisher Site | Google Scholar
  9. A. D. Hansen and G. Michalke, “Multi-pole permanent magnet synchronous generator wind turbines' grid support capability in uninterrupted operation during grid faults,” IET Renewable Power Generation, vol. 3, no. 3, pp. 333–348, 2009. View at: Publisher Site | Google Scholar
  10. M. Nasiri and R. Mohammadi, “Peak current limitation for grid side inverter by limited active power in PMSG-based wind turbines during different grid faults,” IEEE Transactions on Sustainable Energy, vol. 8, no. 1, pp. 3–12, 2017. View at: Publisher Site | Google Scholar
  11. K.-H. Kim, Y.-C. Jeung, D.-C. Lee, and H.-G. Kim, “LVRT scheme of PMSG wind power systems based on feedback linearization,” IEEE Transactions on Power Electronics, vol. 27, no. 5, pp. 2376–2384, 2012. View at: Publisher Site | Google Scholar
  12. M. Nasiri, J. Milimonfared, and S. Fathi, “Efficient low-voltage ride-through nonlinear backstepping control strategy for PMSG-based wind turbine during the grid faults,” Journal of Operation and Automation in Power Engineering, vol. 6, no. 2, pp. 218–228, 2018. View at: Google Scholar
  13. M. Nasiri, J. Milimonfared, and S. H. Fathi, “Robust control of PMSG-based wind turbine under grid fault conditions,” Indian Journal of Science and Technology, vol. 8, no. 13, 2015. View at: Google Scholar
  14. H. Sun, Y. Han, and L. Zhang, “Maximum wind power tracking of doubly fed wind turbine system based on adaptive gain second-order sliding mode,” Journal of Control Science and Engineering, vol. 2018, Article ID 5342971, 11 pages, 2018. View at: Publisher Site | Google Scholar | MathSciNet
  15. M. Benbouzid, B. Beltran, Y. Amirat, G. Yao, J. Han, and H. Mangel, “Second-order sliding mode control for DFIG-based wind turbines fault ride-through capability enhancement,” ISA Transactions, vol. 53, no. 3, pp. 827–833, 2014. View at: Publisher Site | Google Scholar
  16. B. Beltran, M. E. H. Benbouzid, and T. Ahmed-Ali, “Second-order sliding mode control of a doubly fed induction generator driven wind turbine,” IEEE Transactions on Energy Conversion, vol. 27, no. 2, pp. 261–269, 2012. View at: Publisher Site | Google Scholar
  17. F. Z. Tria, K. Srairi, M. T. Benchouia, and M. E. H. Benbouzid, “An integral sliding mode controller with super-twisting algorithm for direct power control of wind generator based on a doubly fed induction generator,” International Journal of Systems Assurance Engineering and Management, vol. 8, no. 4, pp. 762–769, 2017. View at: Google Scholar
  18. Y. Errami, A. Obbadi, S. Sahnoun, and M. Benhmida, “Design and sliding mode control for PMSG based wind power system connected to a non-ideal grid voltages,” in Proceedings of the 2015 3rd International Renewable and Sustainable Energy Conference (IRSEC), IEEE, 2015. View at: Google Scholar
  19. A. Merabet, M. Aminul Islam, R. Beguenane, and H. Ibrahim, “Second-order sliding mode control for variable speed wind turbine experiment system,” Renewable Energy and Power Quality Journal, vol. 12, pp. 1–5, 2014. View at: Publisher Site | Google Scholar
  20. J. Hostettler and X. Wang, “Sliding mode control of a permanent magnet synchronous generator for variable speed wind energy conversion systems,” Systems Science & Control Engineering, vol. 3, no. 1, pp. 453–459, 2015. View at: Publisher Site | Google Scholar
  21. S. Heier, Grid Integration of Wind Energy: Onshore and Offshore Conversion Systems, John Wiley & Sons, 2014.
  22. M. Nasiri, J. Milimonfared, and S. H. Fathi, “Modeling, analysis and comparison of TSR and OTC methods for MPPT and power smoothing in permanent magnet synchronous generator-based wind turbines,” Energy Conversion and Management, vol. 86, pp. 892–900, 2014. View at: Publisher Site | Google Scholar
  23. W. Perruquetti and J. P. Barbot, Sliding Mode Control in Engineering, vol. 11, M. Dekker, 2002.
  24. Y. Shtessel, C. Edwards, L. Fridman, and A. Levant, “Conventional sliding mode observers,” in Sliding Mode Control and Observation, pp. 105–141, Springer, 2014. View at: Publisher Site | Google Scholar | MathSciNet
  25. M. A. Abdullah, A. H. M. Yatim, C. W. Tan, and R. Saidur, “A review of maximum power point tracking algorithms for wind energy systems,” Renewable & Sustainable Energy Reviews, vol. 16, no. 5, pp. 3220–3227, 2012. View at: Publisher Site | Google Scholar

Copyright © 2019 Mojtaba Nasiri et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.


More related articles

955 Views | 538 Downloads | 7 Citations
 PDF  Download Citation  Citation
 Download other formatsMore
 Order printed copiesOrder

Related articles

We are committed to sharing findings related to COVID-19 as quickly and safely as possible. Any author submitting a COVID-19 paper should notify us at help@hindawi.com to ensure their research is fast-tracked and made available on a preprint server as soon as possible. We will be providing unlimited waivers of publication charges for accepted articles related to COVID-19. Sign up here as a reviewer to help fast-track new submissions.