Next Article in Journal
Comparison of Optimization Methods for the Attitude Control of Satellites
Previous Article in Journal
LMANet: A Lightweight Asymmetric Semantic Segmentation Network Based on Multi-Scale Feature Extraction
Previous Article in Special Issue
Double-Loop Controller Design of a Single-Phase 3-Level Power Factor Correction Converter
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Robust PI-PD Controller Design: Industrial Simulation Case Studies and a Real-Time Application

1
Institute of Natural and Applied Sciences, Dicle University, 21280 Diyarbakır, Turkey
2
Electrical and Electronics Engineering Department, Dicle University, 21280 Diyarbakır, Turkey
3
ESTACA’Lab, Ecole Supérieure des Techniques Aéronautiques et de Construction Automobile, 53000 Laval, France
*
Author to whom correspondence should be addressed.
Electronics 2024, 13(17), 3362; https://doi.org/10.3390/electronics13173362
Submission received: 18 July 2024 / Revised: 18 August 2024 / Accepted: 21 August 2024 / Published: 24 August 2024
(This article belongs to the Special Issue Advanced Control Techniques for Power Converter and Drives)

Abstract

:
PI-PD controllers have superior performance compared to traditional PID controllers, especially for controlling unstable and integrating industrial processes with time delays. However, computing the four tuning parameters of this type of controller is not an easy task. Recently, there has been significant interest in determining the tuning rules for PI-PD controllers that utilize the stability region. Currently, most tuning rules for the PI-PD controller are presented graphically, which can be time-consuming and act as a barrier to their industrial application. There is a lack of analytical tuning guidelines in the literature to address this shortfall. However, the existing analytical tuning guidelines do not consider a rigorous design approach. This work proposes new robust analytical tuning criteria based on predefined gain and phase margin bounds, as well as the centroid of the stability region. The proposed method has been tested using various simulation studies related to a DC–DC buck converter, a DC motor, and a heat exchanger. The results indicate that the proposed tuning rules exhibit strong performance against parameter uncertainty with minimal overshoots. Furthermore, the suggested technique for simultaneous control of yaw and pitch angles has been tested in a real-time application using the twin rotor multi-input multi-output system (TRMS). Real-time results indicate that, compared to other methods under investigation, the suggested approach provides nearly minimal overshoots.

1. Introduction

The second order plus time delay (SOPTD) system representation is extensively utilized in the literature because it accurately represents the dynamics of numerous practical industrial applications such as temperature controllers and DC motors. However, long-time delays complicate such systems’ control problems as these delays cause uncertainty and slow responses with larger overshoots. Further, the presence of poles on the right-hand side or at the origin of the s-plane adds another complexity to the control design problem. Thus, a robust control approach is required to regulate such systems [1].
PI-PD controllers are strong candidates relative to PID controllers to solve the control mentioned above difficulties [2]. They use the PD part in an inner feedback loop to shift the poles of a process to a more desirable location so that the PI part can control it more effectively in the forward path [3]. Compared to typical PI/PID controllers, PI-PD controllers perform better, especially when it comes to controlling unstable and integrated processes. They also provide the designer more flexibility to enhance closed-loop performance because PI-PD incorporates four adjustable parameters [4]. However, because four controller parameters need to be adjusted, constructing a PI-PD controller is a difficult task [5]. Only a few applications have made use of the PI-PD controller due to the difficulty of calculating those four parameters [6]. In addition, compared to PID controllers, this controller’s design has received comparatively less study attention [7]. Therefore, more research on this type of controller is still necessary.
Utilizing the stability region approach to calculate all stabilizing controller gains is a widely researched topic (see [8,9,10,11] for examples). Since it provides straightforward and useful tuning guidelines and gets around the challenge of tuning the PI-PD controller, many researchers have recently explored using the centroid point of the stability region, which can be calculated by using the centroid of convex stability region (CCSR) [12] or the weighted geometrical center (WGC) approaches [13,14], to adjust the parameters of PI-PD controllers. Additionally, it has been suggested that the controller evaluated using the centroid point may produce robustness against parameter fluctuations, faster perturbation rejection, and improved set-point tracking [15,16]. WGC and CCSR have just lately been used to graphically adjust the PID controller’s gains, which are intended to control systems with non-integer order delay [17]. Additionally, a comparison study of the PI-PD controller tuning techniques employing CCSR and WGC has recently been conducted [14]. Because CCSR and WGC approaches have only been applied to specified transfer functions in the above-mentioned literature, they share the common drawback of requiring the design process to be repeated whenever the transfer function changes. This might be time-consuming for controller engineers and might require mathematical knowledge. The user-selected step size to be employed in the design phase causes a computational load, which is a unique drawback of the WGC approach [12]. Additionally, if the time lag is particularly significant, the WGC approach may yield an incorrect centroid location because of its sensitivity to the chosen step size [18]. The graphical way of determining the centroid location is a drawback unique to the CCSR approach [4]. Additionally, novel methods for determining a fractional order controller’s stability region’s centroid point have recently been introduced [19,20]. In addition, the centroid point has recently been used in the controller parameter design for DC–DC boost converters [21]. Nevertheless, the centroid point is still obtained graphically in all of the aforementioned stated approaches, which is laborious.
To remove the above-mentioned disadvantages of centroid point approaches, analytical tuning rules have been recently suggested [3,4,16,18,22]. However, these approaches have not considered incorporating a centroid point based on a robust design method. In fact, since the original work that formulated the stability region’s formulas for time delay systems and PI-PD controllers, as documented in Ref. [23], the majority of researchers have tuned these controllers using graphical approaches without taking analytical ones into account. Furthermore, robust design techniques as previously discussed have not taken into account analytical methodologies. Therefore, the main goal of this paper is to design a robust tuning rule based on the centroid point, which can specifically withstand the system parameter uncertainties and provide acceptable settling times and overshoots. The primary contributions made by this work are:
  • Relative to the centroid point approaches reported in Refs. [3,4,16,18,22], the proposed method incorporates a robust design approach based on predefined gain and phase margin boundaries, which gives the designer more flexibility to obtain the desired control performances.
  • Relative to the centroid point approaches reported in Refs. [12,14,19,20], the proposed methods are analytical and not graphical, which saves time and is easy to implement on the industrial level.
  • Compared to the centroid point derived based on a robust design technique for DC–DC converter control published in Ref. [24], the suggested approach is analytical and does not call for any graphical adjustments. Furthermore, any system that can be represented as a stable SOPTD can use the suggested approach.
To conclude, proposing an analytical robust design method for designing the gains of the PI-PD controller is the main contribution of this paper. This is the layout for the rest of the paper. The next part explains the PI-PD controller structure. The suggested methodology is expounded upon in Section 3. Section 4 describes the simulation and the real-time results. The last section deals with conclusions.

2. PI-PD Controller Structure

Every action taken by the PID regulator in a traditional closed-loop system occurs via the forward path. This can cause an unwanted occurrence known as a derivative kick. To solve this issue, the PD part of PI-PD regulators is moved to an internal feedback loop, which moves the poles of the plant transfer function to a more favorable location where the PI component operating in the onward path might more adequately control it [4,18]. The PI-PD regulator’s architecture is shown in Figure 1. The inner loop is composed of the transfer functions of the plant, G ( s ) , and the PD controller, C P D ( s ) , whereas the external loop is composed of the transfer functions of the inner loop and the PI regulator, C P I ( s ) . The transfer functions C P D ( s ) and C P I ( s ) of the PI-PD regulator are written as follows:
C P D ( s ) = k f + k d s
C P I ( s ) = k p + k i s
The behavior of many industrial plants can be described as a second-order transfer function with time delays as follows:
G ( s ) = K e τ s a s 2 + b s + c
In the above equation, the time delay, the process gain, and the coefficients of the transfer function are τ , K , a , b , and c , respectively.

3. Proposed Approach

The gain-phase margin tester, T ( s ) = A e ϕ s , is added to the forward path to find boundaries in the general stability boundary corresponding to the predetermined gain and phase margins. Figure 2 shows the structure of the PI-PD controller after adding the tester [23].
The following closed-loop characteristic equation can be used to characterize the inner feedback loop shown in Figure 2, which is made up of the PD controller transfer function, C P D ( s ) , and the process transfer functions, G ( s ) :
Δ ( s ) = 1 + C P D ( s ) G ( s ) Δ ( s ) = k f K e τ s + k d K s e τ s + a s 2 + b s + c
In the parameter space methodology, the root of a stable polynomial has three ways of crossing over the imaginary axis and becoming unstable. The real root boundary, the infinite root boundary, and the complex root boundary are the three boundaries that these ways define [4]. To find the real root boundary, we solve Equation (4) with s = 0 and Δ ( s ) = 0 . Therefore, k f = c / K defines the boundary. Since the controller parameters are absent from the coefficient of s 2 , the infinite root boundary cannot be part of the stability region. To find the complex root boundary, we enter s = j ω , e j τ ω = cos ( τ ω ) j sin ( τ ω ) , and j = 1 into Equation (4) as given below:
Δ ( s ) = a ω 2 + j b ω + c + k f K cos ( τ ω ) j k f K sin ( τ ω ) + j k d K ω cos ( τ ω ) + k d K ω sin ( τ ω ) = 0
The imaginary and real components in Equation (5) are then split and equaled to zero. Consequently, the following equations can be obtained:
k d ( ω cos ( τ ω ) ) k f ( sin ( τ ω ) ) = b K ω
k d ( ω sin ( τ ω ) ) + k f ( cos ( τ ω ) ) = a ω 2 c K
The expressions that follow can be obtained by solving Equations (6) and (7) for k d and k f :
k d = 1 K b cos ( τ ω ) + ( a ω c ω ) sin ( τ ω )
k f = 1 K ( a ω 2 c ) cos ( τ ω ) + b ω sin ( τ ω )
The inner feedback loop’s stability region is defined by Equations (8) and (9) and the line k f = c / K . By varying ω throughout the range of ε ω P D , Equations (8) and (9) are utilized to construct the stability region. ε is a small number that is used to avoid dividing by zero in Equation (8), while ω P D is found by taking the first root of the equation: ( a ω 2 c ) cos ( τ ω ) + b ω sin ( τ ω ) = c . The adjusting point of the PD controller can be found using the mathematical relations in the range of e f [18]:
x ¯ = 1 f e e f z ( x ) d x
y ¯ = 1 2 ( f e ) e f g ( x ) d x
The integrations provided in Equations (10) and (11) should be used to analytically calculate the settings of the PD controller. Also, the same mathematical relations will be used later for finding the settings of the PI controller. Considering f = ω P D , e = 0 , and z ( x ) = k d which is given in Equation (8), and g ( x ) = k f which is given in Equation (9), the gains of the PD controller are analytically computed using Equations (10) and (11). Thus, the following expressions are obtained and used to find the setting of this controller:
k d * = 1 K ω P D a sin τ ω P D τ ω P D cos τ ω P D τ 2 b sin τ ω P D τ c τ 6 ω P D 5 600 τ 4 ω P D 3 18 + τ 2 ω P D
k f * = 1 2 K ω P D b ( sin ( τ ω P D ) τ ω P D cos ( τ ω P D ) ) τ 2 + ( τ 2 ( a ω P D 2 c ) 2 a ) sin ( τ ω P D ) + 2 a τ ω P D cos ( τ ω P D ) τ 3
Here, it is important to illustrate that the first three terms of the Taylor series are used for approximation τ sin ( τ ω ) τ ω and consequently achieving the integration action in Equation (13). The first three terms of the Taylor series are selected as they can provide good computation accuracy without introducing mathematical complexity. To compute the controller parameters using Equations (12) and (13), one should first identify the parameters of the SOPTD system reported in Equation (3). After that, ω P D should be computed by taking the first root of the equation: ( a ω 2 c ) cos ( τ ω ) + b ω sin ( τ ω ) = c . This can be achieved by using a scientific package such as MATLAB.
To find the PI controller parameters k i * and k p * , the outer loop’s characteristic equation is given below:
Δ ( s ) = 1 + G P D ( s ) G ( s ) + A e j ϕ G P I ( s ) G ( s ) = a s 3 + b s 2 + c s + A K k p s e j ( τ ω + ϕ ) + A K k i e j ( τ ω + ϕ ) + K k f * s e τ s + K k d * s 2 e τ s
Afterward, the real root boundary, the infinite root boundary, and the complex root boundary are determined. For the outer loop, the real root boundary is found by solving Equation (14) for s = 0 and Δ ( s ) = 0 . Thus, this boundary is defined by k i = 0 . The infinite root boundary cannot be a component of the stability zone since the controller’s settings are missing from the coefficient of s 3 . The complex root boundary is found by entering s = j ω , e j τ ω = cos ( τ ω ) j sin ( τ ω ) , e j ϕ = cos ( ϕ ) j sin ( ϕ ) , and j = 1 into Equation (14). After that, imaginary and real components are split to obtain the following equations:
k p ( ω sin ( τ ω + ϕ ) ) + k i ( cos ( τ ω + ϕ ) ) = b A K ω 2 k f * A ω sin ( τ ω ) + k d * A ω 2 cos ( τ ω )
k p ( ω cos ( τ ω + ϕ ) ) k i ( sin ( τ ω + ϕ ) ) = a A K ω 3 c A K ω k f * A ω cos ( τ ω ) k d * A ω 2 sin ( τ ω )
The complex root boundary’s equations are given by solving Equations (15) and (16) for k p and k i as follows:
k p = b A K ω 2 cos ( τ ω + ϕ ) + b A K ω sin ( τ ω + ϕ ) k f * A + sin ( τ ω ) sin ( τ ω + ϕ ) k f * A cos ( τ ω ) cos ( τ ω + ϕ ) + k d * A ω cos ( τ ω ) sin ( τ ω + ϕ ) k d * A ω sin ( τ ω ) cos ( τ ω + ϕ ) c A K cos ( τ ω + ϕ )
k i = a A K ω 3 sin ( τ ω + ϕ ) + c A K ω sin ( τ ω + ϕ ) + b A K ω 2 cos ( τ ω + ϕ ) + k f * A ω cos ( τ ω ) sin ( τ ω + ϕ ) k f * A ω sin ( τ ω ) cos ( τ ω + ϕ ) + k d * A ω 2 sin ( τ ω ) sin ( τ ω + ϕ ) + k d * A ω 2 cos ( τ ω ) cos ( τ ω + ϕ )
The line k i = 0 and Equations (17) and (18) are used to determine the stability region of the outer loop. To plot Equations (17) and (18), ω is changed throughout the range of 0 ω P I , where ω P I is found by equating Equation (18) to zero and taking the first solution. After applying the formulae from Equations (10) and (11) to Equations (17) and (18), correspondingly, the settings of the PI controller are determined as follows:
k p * = 1 ω P I 2 a τ 2 ω P I 2 2 c τ 2 + 2 b τ 4 a sin τ ω P I + ϕ + 4 a τ 2 b τ 2 ω P I cos τ ω P I + ϕ 2 A K τ 3 + K k d τ 3 ω P I 2 sin ϕ 2 K k f τ 3 ω P I cos ϕ 2 A K τ 3 + 2 c τ 2 2 b τ + 4 a sin ϕ 2 A K τ 3
k i * = 1 2 ω P I 6 b τ 3 18 a τ 2 ω P I 2 + 6 c τ 2 12 b τ + 36 a sin τ ω P I + ϕ 6 A K τ 4 + 6 a τ 3 ω P I 3 + 6 c ω P I τ 3 + 12 b τ 2 36 a τ ω P I cos τ ω P I + ϕ 6 A K τ 4 + + 2 K k d τ 4 cos ϕ ω P I 3 + 3 K k f τ 4 sin ϕ ω P I 2 6 A K τ 4 6 c τ 2 12 b τ + 36 a sin ϕ 6 A K τ 4
One can find the controller parameters using the flowchart shown in Figure 3. Here, it is important to mention that no approximations have been used to find Equations (19) and (20).

4. Simulation and Real-Time Case Studies

Three simulation case studies and a real-time application will be shown in this section. The PD controller’s derivative component with the form C P D ( s ) = k f + k d s / ( t f s + 1 ) , where t f = ( k d ) / ( 10 k f ) , is used for noise rejection and practical application. To ensure a fair comparison with the methods in the literature, the derivative filter will also be introduced to all literature techniques under consideration. The examined approaches are numerically evaluated using the gain margin, the phase margin, the integral of absolute error (IAE), I A E = 0 e d t , and the total variation (TV), T V = 0 ( u ( t + 1 ) u ( t ) ) d t .
Case study 1: A DC–DC buck converter, which is shown in Figure 4, is studied here. The identified transfer function of this system is G ( s ) = 8 . 87 e 2 × 10 4 s 7 . 0341 × 10 7 s 2 + 0 . 0017 s + 1 [25]. The suggested method for different predefined gain and phase margin boundaries is applied to find out the controller settings using the procedures shown in Figure 3. Table 1 summarizes these settings’ values as well as the PID controller’s values reported in Reference [25].
Figure 5 and Figure 6 show the stability region with computed centroid points using the proposed method for inner and outer loops, respectively. Figure 6 shows the stability region for three predefined scenarios of the gain and phase margin boundaries, which are A = 1 ϕ = 0 , A = 1.5 ϕ = 0 , and A = 1.6 ϕ = π / 4 . The scenario corresponding to A = 1 ϕ = 0 means that no gain and phase margin boundaries are included in the design of the centroid point. First, the described controllers are assessed using the system parameters’ nominal values. Specifically, at t = 0 s, the system is given a unity step input. Additionally, at t = 0.02 s, a disturbance with an amplitude of 0.3 is introduced into the closed-loop system. Figure 7 illustrates the output voltage responses of the assumed scenarios along with the PID [25] and their control signals. From the figure, all the proposed scenarios give better results in terms of generating small overshoots relative to the PID controller [25], which appears to consume more energy with possible oscillation, as is clear in Figure 7b. İn fact, the overshoot and observed oscillations of the PID controller [25] might be a source of instability in practical implementations. On the other hand, the proposed method does not have an overshoot, which might make it more suitable for real-time applications. Additionally, Figure 7a shows that the proposed scenarios A = 1 ϕ = 0 relative to other scenarios and PID [25] quickly rejects the injected disturbances into the closed loop with an amplitude of 0.3 at t = 0.02 s. Thus, one might conclude that the proposed method, due to the freedom in the selection of the predefined values of the gain and phase margins, has the flexibility to be tuned by control engineers to generate the desired performances according to the application in hand. This can be seen numerically in Table 1, in which one of the proposed scenarios has at least the smallest values for IAE or TV. The optimal values of the gain and phase margins lie in the ranges 2 6 and 30 60 , respectively [4]. Therefore, it seems from Table 1 that the proposed method with all the investigated scenarios provides optimum values of the gain and phase margin relative to the PID controller [25], which has a gain margin of less than 2. The numerical observations, given in Table 1, illustrate again how the proposed method is flexible in terms of providing different desired control performances.
To evaluate the effectiveness of the approaches under evaluation in the event of a situation where the system parameters change, another test is carried out. This test is particularly important to show how the suggested method gives a robust performance against parameter uncertainties. To show the benefits of incorporating the predefined gain and phase margin boundaries into the proposed method, the system parameters τ , a , b , and c are changed by 500%. The system’s parameters are changed by 500% since the reported controllers start to show a noticeable influence. Figure 8 displays the findings of this experiment. It is evident that, in contrast to PID and scenario A = 1 ϕ = 0 , which display oscillatory behavior and the potential to lose stability, scenarios A = 1.5 ϕ = 0 and A = 1.6 ϕ = π / 4 maintain system stability with a very small overshoot. To sum up, the suggested scenarios examined in Figure 7 and Figure 8 demonstrate that the suggested approach gives the control engineer the flexibility and opportunity to choose the most suitable control performance, potentially weighing fast-tracking against robust performance in the face of system parameter fluctuations. This possibility is not available in the centroid tuning methods reported in Refs. [3,4,16,18,22], which forms the main contribution of this paper in addition to the analytical nature of the proposed method.
Case study 2: A position control problem for a DC motor is considered in this case study. The identified transfer function of the DC motor position is G ( s ) = 0 . 9843 e 0.02 s 0 . 0651 s 2 + s [26]. The parameters of the proposed method for the assumed values of A = 1.3 ϕ = 0 are computed using the steps given in Figure 3. The fractional order PI controller (FOPI) [26] and the PI-PD controller adjusted by CCSR [4] are used to compare the performance with the suggested approach. Table 2 contains a list of all the reported controller settings values together with the values of performance measures.
First of all, the investigated methods are tested using the nominal parameters and by applying a unity response and unity disturbance at t = 0 s and t = 3 s, respectively. Figure 9 displays the outcomes for this test. It is clear from the figure that FOPI [26] has poor performance with oscillatory response and long settling times. On the other hand, the proposed method offers a quick response with zero overshoot and a similar disturbance rejection capability relative to PI-PD [4], which has an oscillatory transient response. Also, the PI-PD [4] consumes more energy relative to the proposed controller as seen in the control signals given in Figure 9b.
The second test for the methods under investigation includes introducing parameter perturbations. By introducing 30% uncertainties in all system parameters, the benefits of incorporating the predetermined gain and phase margin boundaries into the design procedures of the proposed method are demonstrated in Figure 10. From the figure, the suggested strategy is the only one that can continue to function well even when there are fluctuations in the system parameters, while other approaches lose their stability entirely. This makes the proposed method suitable for controlling systems that have changing parameters.
The third test for the reported methods is performed using a variable input consisting of square and sinusoidal signals, as shown in Figure 11a. One can see from the figure that the proposed method has decent performance without introducing oscillations or overshoots. As can be observed from the IAE readings for the variable input given in Table 2, the PI-PD controller [4] provides the fewest tracking errors; but, as Figure 11a,b makes evident, it exhibits oscillatory transient responses. This might lead to instability in real-time scenarios. Also, the FOPI controller [26] has a response with overshoots and oscillations, as noticeable in Figure 11a. Thus, one can conclude that only the proposed method can give suitable performance without overshoots and oscillations and with reasonable tracking of errors.
The TV values, listed in Table 2, show that the proposed method has the smoothest control signal for servo or regulatory responses for both types of control signals. The IAE values presented in the same table demonstrate that, for either servo or regulatory responses, the suggested approach offers a trade-off performance between the strategies described in Refs. [4,26] for all the input types. However, different gain and phase margins might be chosen if control engineers want a different level of intended performance. This is because the suggested method has the flexibility to do so as shown in case study 1. Also, the proposed method has optimal values of the gain and phase margins relative to the other methods as clear from Table 2. To conclude, the above observations demonstrate, yet again, the flexibility, robustness, and importance of the suggested approach.
Case study 3: An industrial temperature control based on a heat exchanger, displayed in Figure 12, is considered under this case study. A shell-and-tube system is used to transfer the heat from a primary to a secondary flow. The process control goal is to maintain the secondary flow point’s temperature by regulating the primary flow. The system’s recognized transfer function is G ( s ) = 0 . 002 e 3 s s [2]. By following the methods shown in Figure 3, one can calculate the settings of the suggested technique for assumed scenarios of A = 1 ϕ = 0 and A = 1.3 ϕ = 0 . The values of these settings are summarized in Table 3, together with the reported values of the PI-PD controller plus the Smith predictor reported in reference [2]. PI-PD-SP is the abbreviation for the PI-PD controller plus the Smith predictor [2] in this case study. The goal of the comparison between the proposed method and PI-PD-SP is to show that the proposed method can perform well in the presence of long time delays, even though it has not incorporated another technique for handling the long time delay such as the Smith predictor.
Different tests are performed to test the proposed method. In the first test, a unity step input and disturbance with an amplitude of 3 are applied at t = 0 s and t = 80 s, respectively. The outcomes of this test are shown in Figure 13. The second test includes introducing variable input in the form of square input as shown in Figure 14. The third test is similar to the first one but after introducing 70% model uncertainty in the time delay. Its results are shown in Figure 15. The following observations are noted based on Figure 13, Figure 14 and Figure 15, as well as the performance metrics listed in Table 3:
  • As seen in Figure 13 and Figure 14, the suggested scenario, A = 1 ϕ = 0 , provides somewhat similar responses and results to the PI-PD-SP controller [2], even though it does have an additional mechanism to tackle long time delays.
  • In contrast to PI-PD-SP [2], which completely loses stability under 70% variations in the time delay as clear in Figure 15, the suggested scenarios A = 1 ϕ = 0 and A = 1.3 ϕ = 0 can keep the system stable with a satisfactory performance. Also, the robustness of the proposed method can be numerically seen from the real values of the gain and phase margin provided in Table 3.
  • As can be observed from the TV values listed in Table 3, the suggested scenarios provide smoother control signals. Furthermore, as shown in Figure 13b and Figure 14b, the suggested method uses less control power for the transient response than the PI-PD-SP [2].
Experimental results: Here, the proposed technique is tested using a real-time trial on a TRMS. The dynamics of the nonlinear TRMS system are comparable to those of a helicopter. As seen in Figure 16, it comprises two DC motors that drive the pitch and yaw propellers at the ends of a beam. The yaw rotor produces horizontal revolution on the yaw axis, whilst the pitch rotor is in charge of the upward motion on the pitch axis. There is a noticeable cross-coupling between the rotor’s actions, with each rotor affecting both angles. The yaw transfer function of yaw angle is identified as G ( s ) = 0 . 10 e 0 . 3227 s 0 . 95965 s 2 + s [4]. On the other hand, the transfer function of the pitch angle is identified using the input and output data as G ( s ) = 0 . 32 e s 2 . 347 s 2 + 0.1 s + 1 .
The suggested method in Figure 3 is used to compute the settings of the yaw and pitch PI-PD controllers using the assumed scenarios A = 1 ϕ = 0 and A = 1 ϕ = π / 3.7 , respectively. The proposed controller is implemented using the diagram shown in Figure 17. The PID controller described in Ref. [27] is contrasted with the suggested approach. Table 4 displays the IAE values together with the parameters for every method that has been reported. To ensure a fair comparison, it is vital to note that every controller evaluated is subjected to identical implementation conditions. An example of these conditions might include the absence of external disturbances such as wind.
The real-time outcomes for concurrently controlling the pitch and yaw angles are shown in Figure 18 and Figure 19. The simulation outcomes for the suggested method are also displayed in Figure 18a and Figure 19a. Interestingly, the proposed method yields somewhat similar real-time and simulation results, even though the TRMS is a highly nonlinear system with coupling effects that may be challenging to model using low-order transfer functions. The suggested approach exhibits rapid responses with shorter settling times and overshoot, particularly for the yaw response, as seen in Figure 18. Additionally, Figure 18b and Figure 19b demonstrate how much less energy the suggested controller uses than the PID controller [27]. Furthermore, even though the suggested approach is primarily suggested for single-input single-output systems, the suggested controller could continue to function well considering the existence of coupling influences between yaw and pitch angles. This is explained by the flexibility of the suggested approach, which enables the designer to use predetermined gain and phase margin boundaries to achieve the desired performance.

5. Conclusions

A robust centroid point design technique based on predefined phase and gain margins is presented in this paper. In comparison to the existing literature, the suggested approach is both analytical and robust. The control engineer has the flexibility to choose the required control performance according to the suggested approach. The analytical centroid point techniques described in the literature really lack this feature. The results of the simulations demonstrated that the suggested approach produces responses with reduced overshoots and settling times and performs well against parameter uncertainties. These findings are further supported by the results of the real-time test on TRMS, which shows that the proposed methodology can handle the coupling effects to larger levels when utilized for controlling multiple input and multiple output systems even if it is originally only designed to deal with single-input and single-output systems. Nonetheless, the coupling effects continue to influence the pitch angle control performance, albeit to a lesser extent. Future studies will therefore involve extending the suggested approach to handle the coupling effects in multiple-input and multiple-output systems.

Author Contributions

Methodology, F.A. and I.K.; software, F.A.; validation, F.A.; formal analysis, F.A.; investigation, F.A.; resources, F.A.; data curation, F.A.; writing—original draft, F.A.; writing—review and editing, I.K. and A.A.; visualization, F.A.; supervision, I.K; project administration, I.K.; funding acquisition, A.A. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

Data is contained within the article.

Acknowledgments

This study was implemented at Dicle University under the research project number: Mühendislik.22.002.

Conflicts of Interest

The authors declare no conflicts of interest.

List of Symbols

SymbolExplanation
K Process gain
τ Time delay
a , b , c Coefficients of the process transfer function
G ( s ) Process transfer function
G P D ( s ) PD controller’s transfer function
G P I ( s ) PI controller’s transfer function
k f The proportional gain of the PD controller
k d The derivative gain of the PD controller
k p The proportional gain of the PI controller
k i The integral gain of the PI controller
A Gain margin
θ Phase margin
ω P D The critical frequency of the inner loop
ω P I The critical frequency of the outer loop
e f The minimum and maximum integral limits
x ¯ The centroid coordinate on the x-axis
y ¯ The centroid coordinate on the y-axis
μ The degree of integral operator

References

  1. Alyoussef, F.; Kaya, I. Improved adaptive dynamic non-singular terminal sliding mode controller with fractional disturbance observer. Inf. Sci. 2023, 641, 119110. [Google Scholar] [CrossRef]
  2. Singha, P.; Das, D.; Chakraborty, S.; Lloyds Raja, G. Experimentally validated predictive PI-PD control strategy for delay-dominant chemical processes. Chem. Eng. Sci. 2024, 295, 120197. [Google Scholar] [CrossRef]
  3. Alyoussef, F.; Kaya, I. Novel Tuning Rules for Adjusting the Parameters of the PI-PD Controller for Controlling Unstable Processes. Adv. Transdiscipl. Eng. 2024, 46, 254–261. [Google Scholar]
  4. Alyoussef, F.; Kaya, I. Simple PI-PD tuning rules based on the centroid of the stability region for controlling unstable and integrating processes. ISA Trans. 2023, 134, 238–255. [Google Scholar] [CrossRef] [PubMed]
  5. Kaya, I. I-PD Controller Design for Integrating Time Delay Processes Based on Optimum Analytical Formulas. IFAC-PapersOnLine 2018, 51, 575–580. [Google Scholar] [CrossRef]
  6. Zou, H.; Li, H. Improved PI-PD control design using predictive functional optimization for temperature model of a fluidized catalytic cracking unit. ISA Trans. 2017, 67, 215–221. [Google Scholar] [CrossRef]
  7. Lloyds Raja, G.; Ali, A. New PI-PD Controller Design Strategy for Industrial Unstable and Integrating Processes with Dead Time and Inverse Response. J. Control Autom. Electr. Syst. 2021, 32, 266–280. [Google Scholar] [CrossRef]
  8. Mahapatro, S.R.; Subudhi, B. A New H∞ Weighted Sensitive Function-Based Robust Multi-Loop PID Controller for a Multi-Variable System. IEEE Trans. Circuits Syst. II Express Briefs 2024, 71, 1256–1260. [Google Scholar] [CrossRef]
  9. Jain, C.; Vaishnav, N.; Jain, A.K. MTPA Control of Double Inverter Fed Wound Rotor Induction Motor Drive. IEEE Trans. Energy Convers. 2023, 38, 2326–2343. [Google Scholar] [CrossRef]
  10. Vaishnav, N.; Jain, A. Stabilizing Sets of Current PI Controllers for IM Drives with and without LC Filter. IEEE Trans. Energy Convers. 2022, 37, 958–969. [Google Scholar] [CrossRef]
  11. Mehta, U.; Aryan, P.; Raja, G.L. Tri-Parametric Fractional-Order Controller Design for Integrating Systems with Time Delay. IEEE Trans. Circuits Syst. II Express Briefs 2023, 70, 4166–4170. [Google Scholar] [CrossRef]
  12. Onat, C. A new design method for PI–PD control of unstable processes with dead time. ISA Trans. 2019, 84, 69–81. [Google Scholar] [CrossRef] [PubMed]
  13. Ozyetkin, M.M.; Onat, C.; Tan, N. PI-PD controller design for time delay systems via the weighted geometrical center method. Asian J. Control 2020, 22, 1811–1826. [Google Scholar] [CrossRef]
  14. İrgan, H.; Menak, R.; Tan, N. A comparative study on PI-PD controller design using stability region centroid methods for unstable, integrating and resonant systems with time delay. Meas. Control 2024. [Google Scholar] [CrossRef]
  15. Imran Kalim, M.; Ali, A. New tuning strategy for series cascade control structure. IFAC-PapersOnLine 2020, 53, 195–200. [Google Scholar] [CrossRef]
  16. Alyoussef, F.; Kaya, I. Tuning proportional-integral controllers based on new analytical methods for finding centroid of stability locus for stable/unstable first-order plus dead-time processes. Proc. Inst. Mech. Eng. Part I J. Syst. Control Eng. 2022, 236, 818–831. [Google Scholar] [CrossRef]
  17. Ozyetkin, M.M. An approximation method and PID controller tuning for systems having integer order and non-integer order delay. Alex. Eng. J. 2022, 61, 11365–11375. [Google Scholar] [CrossRef]
  18. Alyoussef, F.; Kaya, I. Proportional–integral and proportional–derivative controller design based on analytically computed centroid point for controlling integrating processes. Proc. Inst. Mech. Eng. Part I J. Syst. Control Eng. 2023, 237, 1045–1065. [Google Scholar] [CrossRef]
  19. Yildirim, B.; Gheisarnejad, M.; Khooban, M.H. A Robust Non-Integer Controller Design for Load Frequency Control in Modern Marine Power Grids. IEEE Trans. Emerg. Top. Comput. Intell. 2021, 6, 852–866. [Google Scholar] [CrossRef]
  20. Yildirim, B.; Gheisarnejad, M.; Khooban, M.H. A New Parameter Tuning Technique for Noninteger Controllers in Low-Inertia Modern Power Grids. IEEE J. Emerg. Sel. Top. Ind. Electron. 2021, 3, 279–288. [Google Scholar] [CrossRef]
  21. Kalim, M.I.; Ali, A. A graphical approach for controller design with desired stability margins for a DC–DC boost converter. IET Power Electron. 2021, 14, 1323–1335. [Google Scholar] [CrossRef]
  22. Çökmez, E.; Kaya, İ. An analytical solution of fractional order PI controller design for stable/unstable/integrating processes with time delay. Turk. J. Electr. Eng. Comput. Sci. 2023, 31, 626–645. [Google Scholar] [CrossRef]
  23. Tan, N. Computation of stabilizing PI-PD controllers. Int. J. Control Autom. Syst. 2009, 7, 175–184. [Google Scholar] [CrossRef]
  24. Kumar, S.R.; Anwar, M.N.; Verma, P.; Ram, M.K.; Somanshu, S.R. H∞ Criterion Based PI Controller Design for DC-DC Step-Up Converters. IEEE J. Emerg. Sel. Top. Ind. Electron. 2023, 5, 402–413. [Google Scholar] [CrossRef]
  25. Ramana, K.V.; Majhi, S.; Gogoi, A.K. Modeling and Estimation of DC-DC Buck Converter Dynamics Using Relay Feedback Output with Performance Evaluation. IEEE Trans. Circuits Syst. II Express Briefs 2019, 66, 427–431. [Google Scholar]
  26. Lino, P.; Maione, G.; Stasi, S.; Padula, F.; Visioli, A. Synthesis of fractional-order PI controllers and fractional-order filters for industrial electrical drives. IEEE/CAA J. Autom. Sin. 2017, 4, 58–69. [Google Scholar] [CrossRef]
  27. Alyoussef, F.; Kaya, I. TRMS Experimental Results of New Nonlinear PID Tuned by de Algorithm. In Proceedings of the 2019 3rd International Conference on Applied Automation and Industrial Diagnostics, ICAAID 2019, Elazig, Turkey, 25–27 September 2019. [Google Scholar]
Figure 1. PI-PD controller.
Figure 1. PI-PD controller.
Electronics 13 03362 g001
Figure 2. PI-PD controller plus gain-phase margin tester.
Figure 2. PI-PD controller plus gain-phase margin tester.
Electronics 13 03362 g002
Figure 3. The flowchart for finding the PI-PD parameters using the proposed method.
Figure 3. The flowchart for finding the PI-PD parameters using the proposed method.
Electronics 13 03362 g003
Figure 4. DC–DC buck converter circuit [25].
Figure 4. DC–DC buck converter circuit [25].
Electronics 13 03362 g004
Figure 5. Stability region of the PD controller for Case study 1 for the nominal parameters.
Figure 5. Stability region of the PD controller for Case study 1 for the nominal parameters.
Electronics 13 03362 g005
Figure 6. Stability region of the PI controller for Case study 1 for the nominal parameters.
Figure 6. Stability region of the PI controller for Case study 1 for the nominal parameters.
Electronics 13 03362 g006
Figure 7. (a) Closed-loop responses for DC–DC buck converter under a 0.3 disturbance amplitude and for the nominal parameters; (b) corresponding control signals.
Figure 7. (a) Closed-loop responses for DC–DC buck converter under a 0.3 disturbance amplitude and for the nominal parameters; (b) corresponding control signals.
Electronics 13 03362 g007
Figure 8. (a) Closed-loop responses for DC–DC buck converter under a 0.3 disturbance amplitude and 500% model uncertainty; (b) corresponding control signals.
Figure 8. (a) Closed-loop responses for DC–DC buck converter under a 0.3 disturbance amplitude and 500% model uncertainty; (b) corresponding control signals.
Electronics 13 03362 g008
Figure 9. (a) Closed-loop responses for DC motor for a step input and under a unity disturbance amplitude for the nominal parameters; (b) corresponding control signals.
Figure 9. (a) Closed-loop responses for DC motor for a step input and under a unity disturbance amplitude for the nominal parameters; (b) corresponding control signals.
Electronics 13 03362 g009
Figure 10. (a) Closed-loop responses for DC motor for a step input and under a unity disturbance amplitude and 30% model uncertainty; (b) corresponding control signals.
Figure 10. (a) Closed-loop responses for DC motor for a step input and under a unity disturbance amplitude and 30% model uncertainty; (b) corresponding control signals.
Electronics 13 03362 g010
Figure 11. (a) Closed-loop responses for DC motor for a variable input and under the nominal parameters; (b) corresponding control signals.
Figure 11. (a) Closed-loop responses for DC motor for a variable input and under the nominal parameters; (b) corresponding control signals.
Electronics 13 03362 g011
Figure 12. Diagrammatic representation of temperature control using a heat exchanger.
Figure 12. Diagrammatic representation of temperature control using a heat exchanger.
Electronics 13 03362 g012
Figure 13. (a) Closed-loop responses for the temperature control for a step input and under the nominal parameters; (b) corresponding control signals.
Figure 13. (a) Closed-loop responses for the temperature control for a step input and under the nominal parameters; (b) corresponding control signals.
Electronics 13 03362 g013
Figure 14. (a) Closed-loop responses for the temperature control for a square input and under the nominal parameters; (b) corresponding control signals.
Figure 14. (a) Closed-loop responses for the temperature control for a square input and under the nominal parameters; (b) corresponding control signals.
Electronics 13 03362 g014
Figure 15. (a) Closed-loop responses for the temperature control under a unity disturbance amplitude and 70% time delay model uncertainty; (b) corresponding control signals.
Figure 15. (a) Closed-loop responses for the temperature control under a unity disturbance amplitude and 70% time delay model uncertainty; (b) corresponding control signals.
Electronics 13 03362 g015
Figure 16. TRMS setup.
Figure 16. TRMS setup.
Electronics 13 03362 g016
Figure 17. Diagram for implementing the PI-PD controller.
Figure 17. Diagram for implementing the PI-PD controller.
Electronics 13 03362 g017
Figure 18. (a) Real-time closed-loop yaw responses for TRMS and (b) corresponding control signals.
Figure 18. (a) Real-time closed-loop yaw responses for TRMS and (b) corresponding control signals.
Electronics 13 03362 g018
Figure 19. (a) Real-time closed-loop pitch responses for TRMS and (b) corresponding control signals.
Figure 19. (a) Real-time closed-loop pitch responses for TRMS and (b) corresponding control signals.
Electronics 13 03362 g019
Table 1. The controller parameters and performance measures for all the studied methods for case study 1.
Table 1. The controller parameters and performance measures for all the studied methods for case study 1.
MethodController ParametersReal Phase and
Gain Margins
IAETV
k f k d k p k i Gain
Margin
Phase
Margin
ServoReg.ServoReg.
Proposed
A = 1 ϕ = 0
0.5390.00020.104605.6192.69 43.48 0.00110.000510.000000610.0000011
Proposed
A = 1.5 ϕ = 0
0.5390.00020.069395.5782.74 49.23 0.00160.000740.000000360.0000010
Proposed
A = 1.6 ϕ = π / 4
0.5390.00020.170235.0892.58 45 . 18 0.00280.00130.000000710.0000012
PID [25]-0.00040.545626.4371.75 49 . 93 0.00130.000660.000000410.0000022
Table 2. The controller parameters and performance measures for all the studied methods for case study 2.
Table 2. The controller parameters and performance measures for all the studied methods for case study 2.
MethodController ParametersReal Phase and
Gain Margins
Step InputVariable Input
IAETVIAETV
k f k d k p k i μ Gain
Margin
Phase
Margin
ServoReg.ServoReg.ServoServo
Proposed
A = 1.3 ϕ = 0
38.8381.8668.175300.739-2.52 37.3 0.130.0030.0050.00050.130.0048
PI-PD [4]47.7152.70329.818823.207-1.67 22.9 0.100.0020.0580.00170.090.058
FOPI [26]--12.103123.7611.60.63 35.9 0.660.1150.0050.00070.660.0048
Table 3. The controller parameters and performance measures for Case study 3.
Table 3. The controller parameters and performance measures for Case study 3.
MethodController ParametersReal Phase and Gain MarginsStep InputSquare Input
IAETVIAETV
k f k d k p k i Gain Margin Phase
Margin
ServoReg.ServoReg.ServoServo
Proposed
A = 1 ϕ = 0
83.34022.7628.1132.18 32 . 63 10.2720.3700.0140.001330.8170.028
Proposed
A = 1.5 ϕ = 0
83.34015.1755.4092.45 39 . 87 15.4080.5550.0090.001146.2030.042
PI-PD-SP [2]130.15188.32626.01125.6011.93 27 . 85 5.1840.1200.0290.001815.5530.086
Table 4. The controller parameters and IAEs for all the studied methods for the real-time test.
Table 4. The controller parameters and IAEs for all the studied methods for the real-time test.
MethodController ParametersIAE
k f k d k p k i
Proposed Yaw   controller
A = 1 ϕ = 0
22.456316.60026.135214.05550.9711
Pitch   controller
A = 1 ϕ = π / 3.7
0.11222.65690.19591.13411.7484
PID [27]Yaw controller-5211.7774
Pitch controller-10561.2817
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Alyoussef, F.; Kaya, I.; Akrad, A. Robust PI-PD Controller Design: Industrial Simulation Case Studies and a Real-Time Application. Electronics 2024, 13, 3362. https://doi.org/10.3390/electronics13173362

AMA Style

Alyoussef F, Kaya I, Akrad A. Robust PI-PD Controller Design: Industrial Simulation Case Studies and a Real-Time Application. Electronics. 2024; 13(17):3362. https://doi.org/10.3390/electronics13173362

Chicago/Turabian Style

Alyoussef, Fadi, Ibrahim Kaya, and Ahmad Akrad. 2024. "Robust PI-PD Controller Design: Industrial Simulation Case Studies and a Real-Time Application" Electronics 13, no. 17: 3362. https://doi.org/10.3390/electronics13173362

APA Style

Alyoussef, F., Kaya, I., & Akrad, A. (2024). Robust PI-PD Controller Design: Industrial Simulation Case Studies and a Real-Time Application. Electronics, 13(17), 3362. https://doi.org/10.3390/electronics13173362

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop