Performance and robustness analysis of V-Tiger PID controller for automatic voltage regulator

This paper presents a comprehensive study on the implementation and analysis of PID controllers in an automated voltage regulator (AVR) system. A novel tuning technique, Virtual Time response-based iterative gain evaluation and re-design (V-Tiger), is introduced to iteratively adjust PID gains for optimal control performance. The study begins with the development of a mathematical model for the AVR system and initialization of PID gains using the Pessen Integral Rule. Virtual time-response analysis is then conducted to evaluate system performance, followed by iterative gain adjustments using Particle Swarm Optimization (PSO) within the V-Tiger framework. MATLAB simulations are employed to implement various controllers, including the V-Tiger PID controller, and their performance is compared in terms of transient response, stability, and control signal generation. Robustness analysis is conducted to assess the system's stability under uncertainties, and worst-case gain analysis is performed to quantify robustness. The transient response of the AVR with the proposed PID controller is compared with other heuristic controllers such as the Flower Pollination Algorithm, Teaching–Learning-based Optimization, Pessen Integral Rule, and Zeigler-Nichols methods. By measuring the peak closed-loop gain of the AVR with the controller and adding uncertainty to the AVR's field exciter and amplifier, the robustness of proposed controller is determined. Plotting the performance degradation curves yields robust stability margins and the accompanying maximum uncertainty that the AVR can withstand without compromising its stability or performance. Based on the degradation curves, robust stability margin of the V-Tiger PID controller is estimated at 3.5. The worst-case peak gains are also estimated using the performance degradation curves. Future research directions include exploring novel optimization techniques for further enhancing control performance in various industrial applications.


Research gap and solution
Due to the system dynamics, nonlinear behavior of the AVR components, dynamics of the AC generator and operational constraints, the optimization of AVR is a challenging task [42][43][44] .The implementation of this innovative method may help to reduce the difficulties associated with AVR optimization 45,46 .V-Tiger is a newly introduced method for adjusting the gains of the PID controller.The advantages of the proposed tuning method are highlighted by comparing its performance with the other PID tuning methods.There are several methods for testing the dynamic system stability in control systems.Among these methods, a pole-zero map is frequently employed because of its ease of use.The Bode plot provides the plant stability in the frequency domain.The plant is more stable the larger the margins 47,48 .The discussion in 36,49 presents disk stability margins to examine the robustness of the controller and closed-loop stability.In comparison to traditional margin analysis, disk-based margin analysis offers a more robust stability guarantee.In particular research studies, the authors have employed plant parameter uncertainties of ± 25% and ± 50% to assess the robustness of their proposed controllers.However, the maximum uncertainty that a system can withstand without losing its stability is not discussed.This article uses performance degradation curves to illustrate the robustness of the controller and the maximum range of uncertainty that the AVR can withstand without losing its stability.

Research contributions
The main contributions of the research article are: 1.A novel tuning approach based on system virtual time response is proposed for tuning the gains of the PID controller.2. Demonstrated that the AVR equipped with the V-Tiger PID controller has a better transient response than the other tuning techniques listed for comparison.3. Illustrated the performance degradation curves of the AVR with various controllers.4. The robust stability margins and the corresponding maximum uncertainty, which a system can tolerate without compromising its stability or performance, are measured from the degradation curves.5.The degradation curves are also used to estimate the worst-case gains at the given uncertainty range.
The rest of this article is structured as follows: The closed-loop transfer function of the AVR is discussed in "Automatic voltage regulator model" section ; the V-Tiger concept and its application to modify the PID gains for the robust AVR is covered in "Evaluating PID controller gains" section; in the section "MATLAB simulation and analysis", the performance and robustness of the AVR utilizing various controllers is discussed.Lastly, this article concludes with research findings.

Mathematical model of AVR
A first-order mathematical model of a generator with an AVR can be represented using a simple transfer function.The first-order model typically captures the primary dynamic behavior of the system, specifically focusing on the response of the generator's voltage to changes in the reference voltage or load on the utility grid.In the context of a generator with AVR, this transfer function can be used to model the response of the generator's voltage to changes in the reference voltage.A simplified first-order model for the generator with AVR is Here V out (s) = AC generator output voltage; V f (s) = Generator field voltage; K g and τ g are the gain and time constants of the AC generator respectively.
In this model, a change in the field voltage V f (s) results in a dynamic response of the generator's output voltage V out (s).Note that this is a simplified representation, and actual generator systems may involve more complex models that consider additional dynamics, such as the electrical and mechanical dynamics of the generator 55,56 .
Generally, the exciter voltage output is the nonlinear function of the field voltage because of the saturation of the magnetic circuit.Hence there is no simple relationship between the generator output voltage and field exciter voltage.A linearized model of the field exciter in an AVR is obtained by ignoring the saturation effect and other nonlinearities.Practically, the field excitation system is used for adjusting the generator field current to control its output voltage.The simplified first-order model of the field excitation system is This model suggests that the field voltage V f (s) responds to changes in error voltage V e (s) (the difference between the reference voltage V ref (s) and the actual output voltage V oa (s)) with a gain K exc and a time constant τ exc .The time constant reflects the speed at which the field voltage adjusts to changes in the reference voltage.
The voltage feedback system is a crucial component in regulating the generator's output voltage 57,58 .The voltage feedback system in an AVR can be formulated as a first-order model by neglecting its saturation effect as an assumption.The first-order transfer function of the voltage feedback system is given below.
(1) where V out (s) = Generator output voltage; V oa (s) = Actual output voltage; K vf and τ vf are the gain and time constants of the voltage feedback sensor.The parameters K vf and τ vf are assumed to be constant, but in reality, these parameters might vary with operating conditions.
The excitation amplifier is an electronic device; it amplifies the voltage error signal before sending it to the excitation system 59 .A mathematical model of the amplifier can be derived by ignoring its nonlinearities and is represented below where V r = Amplifier output voltage; V e (s) = error between V ref (s) and V ao (s); K vf and τ vf are the gain and time constants of the voltage feedback sensor.The parameters K amp and τ amp are the gain and time constant of the amplifier respectively.
Utilizing the above models, the transfer function model of AVR with a power generation system is shown in Fig. 2a.
The closed-loop block diagram in Fig. 2a relating the generator output voltage V out (s) to the reference voltage V ref (s) is Equation ( 5) can be rewritten as For a unit step reference voltage V ref (s) = 1 s , the steady-state error (3) Hvf(s) Hvf(s) The steady-state error defined in (7) can be eliminated by adding the controller to AVR.The controller's role is not just about eliminating steady-state errors but also about maintaining voltage regulation, ensuring system stability, responding to dynamic changes, and preventing voltage violations in different operating conditions 60,61 .In this research article, a PID controller is used.A PID controller is a widely used feedback control system in various industrial applications.It regulates the response of the plant by adjusting the input based on the error generated by the required set-point and the actual response.
The control signal u(t) from the PID controller is mathematically represented as 62 .
where e(t) is the error at time t, K p is the proportional gain, K i is the integral gain and K d is the derivative gain.
Figure 2b shows the transfer function model of AVR with PID controller.The closed-loop transfer function of block diagram in Fig. 2b relating the generator output voltage V out (s) to the reference voltage V ref (s) is The nominal values and proposed uncertainty of the AVR system, presented in Fig. 2a, are given in Table 1 36,47 .By utilizing these nominal parameters, the closed-loop transfer function of the AVR is Also, the closed-loop transfer function of the AVR including the PID controller is given as

Evaluating PID controller gains
Tuning a PID controller is crucial for the proper functioning and performance of a control system 63,64 .Proper tuning ensures that the controller effectively regulates the system, minimizing errors and achieving the desired dynamic response 65 .This article proposes the Pessen Integral Rule (PIR) and Virtual Time response-based iterative gain evaluation and re-designs (V-Tiger) methods to tune the PID controller.The results of the proposed controller are compared with FPA, TLbO, PIR, and ZN-based PID controllers.

Pessen Integral Rule
The Pessen Integral Rule (PIR) is another tuning method for PID controllers.It was proposed by Leonard Pessen and is a modification of the ZN tuning method.This method emphasizes the integral action by setting the integral time constant to one-fourth of the ultimate period obtained from the Ziegler-Nichols method 66,67 .The Pessen Integral Rule aims to improve the transient response of the system by giving more weight to the integral term.Here are the steps to apply the Pessen Integral Rule: Step-I: Perform step response experiment www.nature.com/scientificreports/Using proportional control alone, do a step-response experiment to obtain the plant's ultimate gain (K u ) and ultimate period (P u ).Adjust the proportional gain (K p ) until the response exhibits the persistent oscillations as seen in Fig. 3.
Figure 3 states that the ultimate period (P u ) is the interval of time between two successive peaks, and the ultimate gain (K u ) is the proportional control gain at which persistent oscillations occur.

Step-II:
Calculate the integral time constant (T i ) Calculate the integral and derivative time constants by applying the Pessen integral rule i.e. integral time constant (T i ) = P u /4 and derivate time constant T d = 0.15P u .
Step-III: Set PID controller parameters Set the PID controller gains by using the Pessen Integral Rule as follows

Virtual time response based iterative gain evaluation and re-design (V-Tiger)
The V-Tiger approach to tune a PID (Proportional-Integral-Derivative) controller involves virtual time-response analysis and iterative gain adjustment to optimize the PID parameters for a given system 68,69 .Using this tuning procedure, the frequency characteristics of the controller are combined with the computed frequency components of the controlled object's input and output.The closed-loop system time response is obtained by first calculating the response in the frequency domain, and then using the inverse Fourier transformation.With this ideal, oneshot experimental data is used to compute the controlled plant's time response.The time response obtained by this method is known as Virtual time response 70 .
Let the controlled plant G(z) be a linear time invariant discrete-time SISO system.The input and output time series of the one-shot experimental data of G(z) are, respectively, u 0 (k) and y 0 (k).In this case, k is the sample number (= 1, 2, 3,…, n) and z is the shift operator.Let the time series of the reference signal (input), disturbance signal (output), and disturbance signal (input) to the closed-loop system be denoted by ref(k), δ(k), and δu(k).The discrete fourier transformations of u 0 (k), y 0 (k), ref(k), δ(k), and δ u (k) are u 0 (jω), y 0 (jω), ref(jω), δ(jω), and δ u (jω), respectively.The angular velocity (rad/sec) is denoted by ω, which may be expressed mathematically as follows in terms of sampling time (t s ): ω = 0, 2π/nt s , 4π/nt s . . .(n − 1)2π/nt s .The following presumptions were made to use this method: i.The input time series of one-shot experimental data u 0 (jω a ) ≠ 0 at a given angular velocity ω a if ref(jω a ), δ(jω a ), and δ u (jω a ) are not 0, and vice versa.ii.The closed-loop system should be stabilized by the controller C PID .If controller C PID (jω) ≠ 0, then there is no stability boundary pole for any value of ω. iii.The inverse of plant G -1 (jω a ) has a stability boundary pole at a given angular velocity ω a if C −1 PID (jω a ) ≠ 0, and vice versa.iv.The initial and final values of step response data u 00 and y 00 are considered as steady state.

Basic procedure
The under-dampened step response in Fig. 4a is copied, flipped, and coupled to the original step response as shown in Fig. 4b.In terms of u 00 (k) and y 00 (k), the input and output time series data (u 0 (k) and y 0 (k)) of G(z) are represented as The standard form of PID controller is u(s Since the frequency has only a steady-state response, the discrete-time frequency function of plant G and controller C PID can be written, respectively, as G(jω), C PID (jω).Now the time series of output y 0 (jω) can be represented as Figure 5 depicts the closed-loop plant when the controller (whose control performance is to be evaluated) is introduced.

Implementation of V-Tiger for controller gains
Before applying V-Tiger, model the system to be controlling.Develop a mathematical model that represents the dynamics of the plant under control.Initialize the PID controller with initial gain values.These initial values are based on any heuristic methods or values commonly used in similar systems.In this work, the initial PID gains are based on the Pessen Integral Rule.Perform the virtual time-response analysis using the modeled system and the initial PID controller settings.Simulate the closed-loop system and observe the system's response.Now quantify the performance metrics such as overshoot, settling time, and stability margin from the virtual timeresponse analysis.These metrics will be used to evaluate the performance of the system under the current PID controller settings.Based on this virtual time-response analysis, adjust the PID gains iteratively.To adjust the PID gains iteratively, Particle Swarm Optimization (PSO) is used in this research.Repeat the virtual time-response analysis with the updated PID gains.Measure the performance metrics again and assess whether the changes result in improved control performance.Continue the iterative process of adjusting PID gains, performing virtual experiments, and evaluating system performance until the desired control performance is achieved.The objective function is defined in (21) with the values of overshoot and stability margins as constraints.
To adjust the controller's behavior to fulfill the performance needs, satisfy the constraints, and enhance the robustness of the AVR, here are two weighting factors, ω 1 and ω 2 , selected as 4 and 2, respectively.Os = Overshoot.Figure 7 shows the flow chart of the V-Tiger.

MATLAB simulation and analysis
The suggested controllers are implemented in an AVR system that is simulated using 2021b MATLAB/Simulink.The nominal parameters and their associated uncertainties are shown in Table 1.By implementing the ZN ultimate sensitivity method to (10), the critical gain K u and the ultimate period P u are identified as 1.6053 and 1.121 s respectively.Using the Pessen Integral Rule (PIR), the PID controller gains are K p = 1.1235,K i = 2.509, and K d = 0.189.These gains are considered as initial PID controller gains for obtaining the optimum PID controller using Virtual Time response-based iterative gain evaluation and re-designs (V-Tiger) with the constraints (i) Overshoot ≤ 10%, (ii) Stability margins GM > 3 dB and PM > 20°.The search range for controller gains using V-Tiger is from zero to twice the PID gains of the Pessen Integral Rule.After carrying 30 simulation runs, the optimal gains of the V-Tiger PID controller are identified as K p = 0.8415, K i = 0.6163, and K d = 0.2838.

Transient response analysis
The transient response of the AVR without a controller possesses a percentage overshoot of 50.53%, settling time of 6.99 s with steady-state error of 0.0909 pu 36 .Figure 8 displays, for various controllers, the transient response (20) Os 100  According to Fig. 8 and Table 2, the AVR using the V-Tiger PID controller possesses superior transient response, more stable margin range, and damping ratio than the other controllers.Hence the AVR with V-Tiger PID controller is more stable and exhibits good performance.Figure 9 compares the control signals generated by the various PID controllers.Figure 10 illustrates how different PID controllers provide field excitation to AC generator.

Robustness analysis
Robustness analysis is a crucial aspect of system design and control engineering.It involves examining how well a system can maintain stability and desired performance levels in the presence of uncertainties or variations in its parameters.The goal is to identify the maximum amount of uncertainty that a system can tolerate without compromising its stability or performance.The trade-off curve, as shown in Fig. 11a, provides a visual representation of the relationship between the normalized amount of uncertainty in the system and its performance, specifically measured by the peak gain of the closed-loop transfer function using the Bode plot.In Fig. 11a, the value of normalized uncertainty (x) = 1 corresponds to the uncertainty ranges specified in the model.This is  considered the nominal level of uncertainty.When x = 2, it represents a system with twice as much uncertainty compared to the nominal system.x = 0 corresponds to the nominal system.Performance is represented on the y-axis by the peak gain of a closed-loop transfer function.From the graph, it can be inferred that the plant moves into an unstable zone when the normalized uncertainty rises with an increase in peak gain.
The performance degradation curves of the AVR with various controllers are shown in Fig. 11(b).This graph indicates that in the absence of a controller, the AVR system becomes unstable for even small amounts of system uncertainty.The AVR with V-Tiger PID controller is more stable than the FPA PID, TLbO PID, and ZN PID controllers in terms of system uncertainty as well.The peak gains and time delays of the AVR system with various PID controllers for the normalized uncertainty range of the amplifier and exciter are listed in Table 3.
According to Table 3, when the uncertainty of the amplifier and exciter are set to their nominal values (x = 0), the nominal system gain (peak) with V-Tiger-PID controller is around 0.025dB.The peak gain of the AVR with the V-Tiger PID controller is infinite at normalized uncertainty (x) of amplifier and exciter of AVR is 3.5 i.e. the AVR system becomes unstable when the uncertainty range of the amplifier and exciter of AVR is beyond 3.5 times the specified uncertainty of the AVR model.Therefore, the robust stability margin of the AVR with the V-Tiger PID controller is 3.5.A robust stability margin means it is the maximum value of uncertainty that the system can withstand and remains stable.Similarly, for the AVR using FPA-PID, TLbO-PID, PIR-PID, and ZN-PID controllers, respectively, the peak gains of the AVR become infinite beyond the normalized uncertainty x = 1.5, x = 1.5, x = 1, and x = 3 and the robust stability margins are identified as 1.5, 1.5, 1 and 3. Figure 12a and b show the AVR responses for different normalized uncertainty levels using V-Tiger and ZN PID controllers, respectively.4, for small amounts of uncertainty in the amplifier and exciter, the nominal AVR system without a controller becomes unstable.In contrast, the AVR with FPA-PID, TLbO-PID, PIR-PID, and ZN-PID controllers becomes unstable when the normalized uncertainty in the amplifier and exciter is x = 1.5, x = 1.5, x = 1, and x = 3 correspondingly.With the V-Tiger PID controller, the permitted normalized uncertainty is measured as 3.5.This investigation makes it evident that the AVR with V-Tiger PID controller shows superior robustness performance than the other controllers.

Worst-case gain
The worst-case gain (K wc ) is the peak gain at the specific uncertainty range.This value is the counterpart of the robust performance margin.Figure 13 shows the performance degradation curve of the AVR with the V-Tiger PID controller and it rises monotonically in proportion to the degree of uncertainty.The worst-case gain for the AVR with the V-Tiger-PID controller at the given uncertainty (x = 1) is around 7.95 dB, according to Fig. 13.When the specified uncertainty of the AVR using the V-Tiger-PID controller doubles (x = 2), the worst-cage gain hits 15.4 dB.Similarly, the AVR with FPA-PID, TLbO-PID, PIR-PID, and ZN-PID controllers has worst-case gains of about 15 dB, 15.1 dB, 47 dB, and 8.98 dB for the specified uncertainty.The AVR with FPA-PID, TLbO-PID, and PIR-PID controllers is unstable when the stated uncertainty is doubled; however, the AVR with ZN-PID controller has the worst-case gain of 17.6 dB.

Conclusions and future research directions
In this study, a novel tuning technique known as Virtual Time response based iterative gain evaluation and redesign (V-Tiger) is proposed, and its use in adjusting PID controller gains for automated voltage regulator (AVR) systems is examined.The AVR system with V-Tiger tuning exhibits better transient response than other tuning techniques, as shown in Fig. 8 and Table 2.In particular, it performs better than the ZN, PIR, FPA, and TLbO PID controllers in terms of low peak overshoot, short settling time, and higher damping factor.Furthermore, by adding uncertainty to the field exciter and amplifier, the transient response of the AVR with the V-Tiger and ZN PID controllers is also shown.The performance degradation curves of AVR with different controllers are plotted in Fig. 11b.From the curves, the robust stability margins and worst-case stability gains are measured.The robust stability margins of the AVR system using the ZN PID and V-Tiger PID controllers are 3 and 3.5, respectively, for the maximum allowable uncertainty of the amplifier and field exciter.Moreover, the worst-case stability gain with a ZN PID controller is 8.67 dB, and with a V-Tiger PID controller, it is 7.95 dB.Future research will use novel optimization tuning techniques, such as the honey badger algorithm and the dandelion optimizer, to examine the transient responsiveness of the generator output and the robustness of the AVR with the controller.Also, the proposed methodologies will be used for Automatic generation control, Power System Stabilizer, speed control of industrial DC drives, and distributed network voltage profile improvement.
Future research in this area shows potential for further advancements in voltage regulation and control strategies, with an emphasis on improving system performance and robustness.One possible avenue of investigation is the use of novel optimization techniques, such as the honey badger algorithm and the dandelion optimizer, to fine-tune control parameters and improve system responsiveness.Furthermore, there exists a potential to broaden the scope of the proposed approaches beyond the context of voltage regulation in power generation plants.Future research may look into the effectiveness of these technologies in alternative areas, including Automatic Generation Control (AGC), Power System Stabilization (PSS), and speed regulation of industrial DC drives.Moreover, research studies may explore the integration of these control strategies into distributed energy systems with the aim of improving voltage profile management and ensuring grid stability.Furthermore, the development of better predictive modeling tools may provide more precise forecasting of voltage variations, allowing for proactive mitigation strategies.Integration with machine learning algorithms could improve system adaptability and efficiency by allowing for real-time decision-making depending on changing grid conditions.
In conclusion, future research should focus on enhancing the current level of knowledge in voltage regulation and control, exploring applications that involve various disciplines, and addressing new challenges related to sustainability and resilience in power systems.

Figure 1 .
Figure 1.Typical arrangement of power generation with AVR.

Figure 2 .
Figure 2. (a) Transfer function model of AVR with the power generation system.(b) Transfer function model of AVR with PID controller.

Figure 9 .
Figure 9.Comparison of control signals by different PID controllers.

Figure 12 .
Figure 12.(a) Response of the AVR with V-Tiger-PID controller for normalized uncertainty.(b) Response of the AVR with ZN-PID controller for normalized uncertainty.

PerformanceFigure 13 .
Figure 13.Performance degradation curve of AVR with V-Tiger PID controller.

Table 1 .
Nominal values and proposed uncertainty of AVR components.

Table 2 .
Comparison of time response, stability performance, and closed-loop poles with damping ratio for different controllers.T r = Rise time; T s = Settling time; %Os = Percentage Overshoot; GM = Gain margin at rad/ sec; PM = Phase margin at rad/sec. of the generator output voltage (pu).Table2lists the step response specifications, stability performance, and closed loop poles with the damping ratio of the AVR utilizing various controllers.

Table 3 .
Peak gains and time delays of AVR for the normalized uncertainty of amplifier and exciter.K Pg = Peak gain (Closed loop), dB; t d = Delay margin, sec.

Table 4
shows the comparison of AVR's transient response for normalized uncertainty of amplifier and exciter.From Table

Table 4 .
AVR's transient response for a range of normalized uncertainty of amplifier and exciter.