Next Article in Journal
Heterogeneity of Micro- and Nanopore Structure of Lacustrine Shales with Complex Lamina Structure
Next Article in Special Issue
Hover Flight Improvement of a Quadrotor Unmanned Aerial Vehicle Using PID Controllers with an Integral Effect Based on the Riemann–Liouville Fractional-Order Operator: A Deterministic Approach
Previous Article in Journal
Fixed Point Theorems: Exploring Applications in Fractional Differential Equations for Economic Growth
Previous Article in Special Issue
Fractional-Order Control Method Based on Twin-Delayed Deep Deterministic Policy Gradient Algorithm
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Robust Fractional-Order PI/PD Controllers for a Cascade Control Structure of Servo Systems

Faculty of Mechanical Engineering, Ho Chi Minh City University of Technology and Education, 01 Vo Van Ngan St, Thu Duc City, Ho Chi Minh City 71300, Vietnam
*
Author to whom correspondence should be addressed.
Fractal Fract. 2024, 8(4), 244; https://doi.org/10.3390/fractalfract8040244
Submission received: 17 March 2024 / Revised: 14 April 2024 / Accepted: 17 April 2024 / Published: 22 April 2024
(This article belongs to the Special Issue Applications of Fractional-Order Systems to Automatic Control)

Abstract

:
In this paper, a cascade control structure is suggested to control servo systems that normally include a servo motor in coupling with two kinds of mechanism elements, a translational or rotational movement. These kinds of systems have high demands for performance in terms of fastest response and no overshoot/oscillation to a ramp function input. The fractional-order proportional integral (FOPI) and proportional derivative (FOPD) controllers are addressed to deal with those control problems due to their flexibility in tuning rules and robustness. The tuning rules are designed in the frequency domain based on the concept of the direct synthesis method and also ensure the robust stability of controlled systems by using the maximum sensitivity function. The M-Δ structure, using multiplicative output uncertainties for both control loops simultaneously, is addressed to justify the robustness of the controlled systems. Simulation studies are considered for two kinds of plants that prove the effectiveness of the proposed method, with good tracking of the ramp function input under the effects of the disturbances. In addition, the robustness of the controlled system is illustrated by a structured singular value (µ) plot in which its value is less than 1 over the frequency range.

1. Introduction

In general, a servo system is an electromechanical drive that includes a servo motor coupled with a ball screw or a rotational load as a mechanism element. Servo systems have been widely used in industrial machines such as computerized numerical controls (CNCs), robots, and so on, with strict requirements for speed and position control, which are essentially fast responses, and no overshoots or oscillations to a ramp function input. The well-known control structure for these systems consists of two control loops in a cascade scheme [1], also called cascade control structure (CCS), which was first introduced in chemical process control [2] and gradually became common in industry applications. This structure consists of two control loops: the inner or secondary loop and the outer or primary loop. The main advantage of this scheme is to reduce the disturbance of the inner loop and improve the servo response of the outer one. This means that if a disturbance affects the inner loop, it will be attenuated before having an impact on the primary output. Therefore, currently, most of the commercial controllers for servo systems utilize this configuration with a proportional integral controller (PI) for the velocity control in the inner loop and a proportional (P) or a proportional derivative (PD) controller for the position control in the outer loop [3,4,5,6,7,8]. However, it is still difficult to derive analytical tuning rules to meet the increasingly high requirements of motion control and also ensure robustness to system uncertainties [9].
In the literature, for cascade control systems, many studies have been conducted to improve the control performances. Lee et al. [10] proposed a method for tuning the PID controller based on the internal model control (IMC) structure for both secondary and primary loops. Many researchers deploy this approach to solve the tuning problems for a variety of plants such as first-order plus time delay, second-order plus time delay, integrating systems, and unstable systems [11,12,13,14,15,16]. However, all mentioned methods are only suitable for the processes with sluggish responses and overshoot at the output being accepted. To improve system performance, various works have been reported to use some advanced control techniques such as cascaded sliding mode control [17,18,19], robust adaptive control [20,21], adaptive fuzzy control [9,22], and so on. These approaches have shown better performances compared to traditional controllers. However, it is hard for them to find their way into industrial applications due to their complexity and the requirement of enhanced modeling parameters. In addition, some nonlinear strategies, especially sliding mode control, normally have an issue in manipulated variables such as discontinuous control signals or ripple effects, which could cause harm to actuators.
Recently, fractional calculus has attracted researchers’ attention in the control field in terms of combination with some conventional controllers such as fractional-order PID, fractional-order nonlinear control [18,19], and fractional-order fuzzy control [23]. Due to the most popular PID controller in industrial applications, the generalization of the PID controller, which is called a fractional-order PID controller (FOPID) or PIλDµ [24], for which, λ and μ are the fractional-order integral and derivative, respectively, deserves to be studied thoroughly for future controllers. The FOPID controller affords more flexibility in tuning rules due to having two extra tuning parameters and more robustness than the integer-order one. In the last two decades, there have been many works that have been proposed to tune the FOPID controller. Normally, they could be categorized into two approaches: the frequency-based method and the time domain-based heuristic algorithm. In the frequency domain, the internal model control (IMC)-based scheme or Bode’s ideal transfer function is used to reduce the number of tuning parameters; and finally, some criteria are addressed to ensure robust performance, such as using constraints on phase margin and gain crossover frequency [25,26,27,28,29,30]. However, most are only available for single-input, single-output (SISO) systems. Therefore, expanding some existing design methods to a complex structure is necessary and deserves to attract more attention from researchers. In the second approach, some evolutionary algorithms such as the genetic algorithm (GA) and particle swarm optimization (PSO) are used to solve the control design problems [31,32]. In [32], the authors proposed the multi-objective PSO to design a robust FOPID controller for first-order plus delay time processes. The objective function considered both system performances in terms of set-point tracking as well as disturbance rejection and system robustness using the maximum sensitivity function. However, this is only verified by SISO systems. Recently, some modern metaheuristic methods [23] such as marine predators algorithm (MPA) [33] or a combination with reinforcement learning for online tuning [34] also have been proposed. From the review, most FOPID controllers are used for SISO or non-linear systems, there are only a few works that use fractional-order controllers for cascade structures in the literature. In [35], the authors proposed the FOPI controller based on internal model control for the inner loop and a PD controller for the outer loop of a parallel cascade control. However, the robust stability is only investigated by simulation studies without a specific systematic method.
In this work, a unified approach to tuning rules for both closed loops in the cascade control system is proposed. In the secondary loop, the FOPI controller is suggested, and its tuning rules are based on the direct synthesis method in the frequency domain. To guarantee the robustness of the controlled system, the maximum sensitivity function is adopted to find out the range of frequency where the Ms value is as close to 1.2 as possible [15,27,36]. Whereas, in the primary loop, the FOPD controller will be adopted and the corresponding tuning rules are similar to those of the inner loop to meet the high requirements of a position servo system, i.e., fast response with no overshoot and oscillation. To justify the robust stability of the whole control system (CCS), the M-Δ structure with multiplicative output uncertainties is considered [36]. Moreover, in this study, the parametric uncertainties of both control loops are investigated simultaneously. Simulation studies will be conducted in this work to illustrate the performance of the controlled system.
This paper is organized as follows. Section 2 briefly introduces fractional calculus and fractional-order controllers including FOPD and FOPI in the frequency domain that are used in this study. In addition, the mechanism systems are also introduced in this section. The proposed method is mentioned in Section 3 including the control structure and also the tuning rules for both control loops. The robust stability is analyzed in Section 4. The simulation studies with two different mechanisms will be discussed in Section 5. Finally, conclusions are given in Section 6.

2. Preliminaries

2.1. Fractional Calculus

Fractional calculus is a generalization of ordinary calculus. It develops a functional operator, D, associated to the order, v   v , that generalizes usual derivatives (for positive v) and integrals (for negative v). The most commonly used is the Riemann–Liouville definition [37], which is generalized by the following equation:
D t v a f t = 1 Γ n v d n d t n a t f τ t τ v n + 1 d t ,   n 1 < v < n
where Γ denotes Euler’s gamma function. a and t are the limits.
Note that the Laplace transform of the fractional derivative/integral in Equation (1) follows the rule for zero initial condition for order v   0 < v < 1 :
L D t ± v a f t = s ± v F s

2.2. Fractional Linear Model

For a single-input, single-output linear time invariant system where both input (u(t)) and output (y(t)) signals are relaxed at t = 0, the fractional-order differential equation can be expressed as:
i = 0 n a i D 0 α i y t = j = 0 m b j D 0 β j u t
Using Laplace transformation, Equation (3) can derive the following transfer function:
G s = Y s U s = b m s β m + b m 1 s β m 1 + + b 0 s β 0 a n s α n + a n 1 s α n 1 + + a 0 s α 0
where α i and β i are arbitrary, real, and positive.

2.3. The FOPD/FOPI Controller in Frequency Domain

Based on the fractional-order PID controller (PIλDµ) proposed by Podlubny [24], two fractional-order PD and PI controllers are derived as follows, respectively:
G c 1 s = K p 1 + K d 1 s α   0 < α < 2
G c 2 s = K p 2 + K i 2 s β   0 < β < 2
By substituting s = j ω into (5) and (6), the FOPD and FOPI controllers are presented in the frequency domain in Equations (7) and (8), respectively:
G c 1 j ω = K p 1 + K d 1 ω α cos γ d + j K d 1 ω α sin γ d
G c 2 j ω = K p 2 + K i 2 cos γ i ω β j K i 2 sin γ i ω β
where γ d = π α 2 and γ i = π β 2 .

2.4. System Modeling

In general, a single-axis servo system includes a servo motor that is used to drive a ball screw for translational movement or a rotational mechanism. In this work, a DC servo motor is considered due to its simplicity in modeling and effectiveness in positioning control. Normally, the control system has both the position and velocity control loops where the velocity and position information are estimated based on feedback signals of an encoder and the parameters of the mechanism. The mathematical model of a DC servo motor is available in numerous studies [4,5,6] and is also described as a block diagram in Figure 1.
where,
Ra: armature resistance (Ohm); La: armature inductance (H)
Ke: back emf constant (V/rad/s); Km: torque constant (Nm/A)
Jm: inertial moment of the motor shaft (kgm2)
ω: angular velocity of the motor shaft (rad/s)
u: applied voltage to the motor
From Figure 1, the equivalent transfer function of the DC servo motor is obtained as follows:
G M ( s ) = ω ( s ) U ( s ) = K m L a J m s 2 + ( L a b + R a J m ) s + R a b + K m K e
G M ( s ) = K 2 τ 2 s 2 + τ 1 s + 1
where K 2 = K m K m K e + R a b , τ 2 = L a J m K m K e + R a b , τ 1 = L a b + R a J m K m K e + R a b .
  • Case 1: a ball screw system, also called a feed drive system
The block diagram of the mechanical element [6,8] is obtained, as seen in Figure 2:
where ω is the angular velocity (rad/s) of the motor; K l = l 2 π , and l is the lead of the ball screw (mm); and y is the position of the sliding stable (mm). Therefore, the primary transfer function in this case is obtained:
G l ( s ) = K l s
  • Case 2: a rotational mechanism
The block diagram, in this case, can be simplified as Figure 3 [6]:
where J L is the inertia moment of the rotational load; D L is viscous damping; and θ is the rotational angle output (rad). From Figure 3, the transfer function of the primary plant is derived as follows:
G L ( s ) = K L s ( τ s + 1 )
where K L = 1 D L ; τ = J L D L .

3. Analytical Design of FOPI/FOPD Controllers for the Cascade Scheme

3.1. The Control System Structure

The general cascade control system is shown in Figure 4, where G p 1 ( s ) and G p 2 ( s ) are the transfer functions of the primary and secondary plants, respectively. Additionally, as mentioned above, G p 2 ( s ) represents the transfer function of the motor Equation (10) and G p 1 ( s ) describes the transfer function of the load in both cases (translational and rotational movement); G c 1 ( s ) and G c 2 ( s ) are the primary and secondary controllers; and d 1 and d 2 denote load disturbances affecting each control loop, respectively.
According to Section 2.4, a dynamic second-order system is used for the secondary model Equation (10).
G p 2 ( s ) = G M ( s ) = K 2 τ 1 s 2 + τ 2 s + 1
Whereas, for the primary loop, an integrating system or an integrating first-order stable system is adopted as Equation (11) or Equation (12) for both cases of the load, respectively:
G p 1 ( s ) = G l ( s ) = K l s
G p 1 ( s ) = G L ( s ) = K L s ( τ s + 1 )

3.2. Design of Secondary Controller-Based Direct Synthesis Method

Considering the secondary control loop where its input and output are r 2 and y 2 , respectively, the closed-loop transfer function is obtained:
y 2 r 2 = G c 2 G p 2 1 + G c 2 G p 2
From Equation (16), the function of the controller is derived as follows:
G c 2 = 1 G ˜ p 2 y 2 / r 2 d 1 y 2 / r 2 d
where G ˜ p 2 is the model of the system in the secondary loop, ideally G ˜ p 2 = G p 2 ; y 2 / r 2 d represents the desired response of the controlled system. In this paper, a first-order transfer function is chosen for the inner loop; therefore, its transfer function has the following form:
y 2 r 2 d = 1 τ c 2 s + 1
where τ c 2 is the desired time response of the controlled loop and τ c 2 is chosen to ensure the fastest response with no overshoot of the controlled system.
Replacing (13) and (18) into (17), the general form of G c 2 is obtained:
G c 2 = τ 1 s 2 + τ 2 s + 1 K 2 1 τ c 2 s
The complex form in frequency domain of the controller is derived by substituting s = j ω into Equation (19):
G c 2 = τ 2 K 2 τ c 2 j 1 τ 1 ω 2 K 2 τ c 2 ω
Comparing Equations (8) and (20), the control parameters of the secondary control loop are derived:
K i 2 = ω β 1 τ 1 ω 2 K 2 τ c 2 ω sin γ i
K p 2 = τ 2 K 2 τ c 2 ω β 1 τ 1 ω 2 cos γ i K 2 τ c 2 ω β + 1 sin γ i

3.3. FOPD Controller Design for the Primary Control Loop

The block diagram of the system is reduced in Figure 5 where y 2 / r 2 d the desired closed-loop transfer function of the inner loop as mentioned above. The closed-loop transfer function of the outer loop (from r1 to y1) is described as Equation (23):
Y 1 R 1 = G c 1 ( s ) y 2 / r 2 d G p 1 ( s ) 1 + G c 1 ( s ) y 2 / r 2 d G p 1 ( s )
  • Case 1: The primary plant is a feed drive system, so its transfer function is:   G p 1 ( s ) = K l s
The equivalent of the primary system is considered as the desired transfer function of the inner loop in series with the primary plant. Therefore, the primary transfer function is obtained as follows:
G p ( s ) = y 2 / r 2 d G p 1 ( s ) = 1 τ c 2 s + 1 K l s
Similarly to Equation (17), the primary controller is also derived:
G c 1 = 1 G p y 1 / r 1 1 y 1 / r 1 1 G ˜ p y 1 / r 1 d 1 y 1 / r 1 d
Ideally, G ˜ p ( s ) = G p ( s ) = 1 τ c 2 s + 1 K l s ; and, in this case,
y 1 r 1 d = 1 τ c 1 s λ + 1   ( 1 < λ < 2 )
where λ is the fractional order; τ c 1 is the desired response time of the primary output.
Replacing Equation (26) and G ˜ p into Equation (25), the primary controller is obtained:
G c 1 = s ( τ c 2 s + 1 ) K l 1 τ c 1 s λ = ( τ c 2 s + 1 ) K l 1 τ c 1 s λ 1 ,   where   λ 1 = λ 1   ( 0 < λ 1 < 1 )
Substituting s = j ω to convert G c 1 into the frequency domain in a complex form:
G c 1 = τ c 2 ω sin γ 1 + cos γ 1 τ c 1 K l ω λ 1 + j τ c 2 ω cos γ 1 sin γ 1 τ c 1 K l ω λ 1
Comparing Equations (7) and (28), the control parameters of the primary control loop are obtained:
K d 1 = τ c 2 ω cos γ 1 sin γ 1 τ c 1 K l ω λ 1 + α sin γ d
K p 1 = sin ( γ 1 + γ d ) τ c 2 ω cos ( γ 1 + γ d ) K l τ c 1 ω λ 1 sin γ d
  • Case 2: The primary plant is a rotational mechanism,   G p 1 ( s ) = K L s ( τ s + 1 )
In this case, the equivalent transfer function of the primary plant is:
G p ( s ) = y 2 / r 2 d   G p 1 ( s ) = 1 τ c 2 s + 1 K L s ( τ s + 1 )
Therefore, the primary controller could be calculated as follows:
G c 1 ( s ) = 1 G p y 1 / r 1 1 y 1 / r 1 = 1 G ˜ p y 1 / r 1 d 1 y 1 / r 1 d
Ideally, G ˜ p ( s ) = G p ( s ) = 1 τ c 2 s + 1 K L s ( τ s + 1 ) , and y 1 / r 1 d is chosen as Equation (26) (case 1); therefore, the primary controller is rewritten as follows:
G c 1 ( s ) = s ( τ c 2 s + 1 ) ( τ s + 1 ) K L 1 τ c 1 s λ = ( τ c 2 s + 1 ) ( τ s + 1 ) K L τ c 1 s λ 1
where λ 1 = λ 1 ( 0 < λ 1 < 1 ).
Similarly to case 1, substituting s = j ω to Equation (33) and rewriting G c 1 in a complex form:
G c 1 = ( 1 τ τ c 2 ω 2 ) cos γ 1 + ω ( τ + τ c 2 ) sin γ 1 K L τ c 1 ω λ 1 + j ( τ + τ c 2 ) ω cos γ 1 ( 1 τ τ c 2 ω 2 ) sin γ 1 K L τ c 1 ω λ 1
Comparing Equations (7) and (34), the control parameters of the primary control loop are obtained in this case:
K d 1 = ( τ + τ c 2 ) ω cos γ 1 ( 1 τ τ c 2 ω 2 ) sin γ 1 K L τ c 1 ω λ 1 + α sin γ d
K p 1 = ( 1 τ τ c 2 ω 2 ) sin ( γ 1 + γ d ) ( τ + τ c 2 ) ω cos ( γ 1 + γ d ) K L τ c 1 ω λ 1 sin γ d

3.4. Tuning Procedure Using Maximum Sensitivity Function

From the previous section, it can be seen that the controller parameters are tuned based on some parameters that play an important role in achieving better performance as well as maintaining robustness to model uncertainties. In the proposed method, the maximum sensitivity function is adopted to guarantee the robust stability of each control loop. For a classical feedback control system, the maximum sensitivity function is defined as the following equation:
M s = max ω S ( j ω )
where S = 1 + L 1 , and L is an open-loop transfer function of the system, which normally includes the controller and the system model. Normally, S is small at low frequencies and reaches 1 at high frequencies. However, at some intermediate frequencies, in practice, a peak value of Ms can be larger than 1, which degrades the system performance. Therefore, the peak value of Ms can be used to measure the robustness of the controlled system. In general, for both robust stability and performance, the Ms value should be close to 1. Therefore, in most of the literature, the typical Ms range is chosen from 1.2 to 2 to ensure the robust stability of control systems [36].
For the cascade control scheme in Figure 1, the maximum sensitivity function is suggested for both control loops (inner and outer) as a tuning criterion and its value is assigned as close as possible to 1.2. In this work, this value is used to choose the appropriate frequency to calculate the control parameters of the inner loop ( ω in Equations (21) and (22)), and this frequency is addressed for the outer loop as well. The algorithm to obtain ω is simple and is described as follows (Algorithm 1):
Algorithm 1: The guideline for tuning parameters based on Ms value
1:  Initialization
   ω = 1 ; choose  τ c 2  and  β
2:  while  ω < 1000  do
3:    Compute  K p 2 ,   K i 2  according to Equations (21) and (22) respectively
4:    Calculate Ms for each set of control parameters
5:     ω = ω + 1
6:  end while
7:  Choose appropriate  ω  to get value of Ms being closed to 1.2
8:  end

4. Robustness Analysis

In this work, the proposed tuning rules are based on the models of actual systems, which are usually represented by nominal models. The differences between the actual system and its model are considered model mismatch or model uncertainty and may degrade some performance indices in terms of servomechanism and regulator problems. Therefore, it is important to analyze the robust stability of model-based methods for control systems in the presence of uncertainties. The multiple sources of parametric and/or unmodelled dynamics uncertainty are commonly grouped into a single lumped perturbation of a specific structure [36].
In this paper, multiplicative output uncertainty is suggested for each control loop of the cascade structure. Therefore, robust stability analysis is carried out herein by considering the multiplicative output uncertainties of each system parameters simultaneously, as shown in Figure 6.
It can be expressed as follows:
Π O 1 :   G ¯ p 1 = ( 1 + E O 1 ) G p 1 ;   E O 1 = W 1 Δ 1
Π O 2 :   G ¯ p 2 = ( 1 + E O 2 ) G p 2 ;   E O 2 = W 2 Δ 2
where Π O i denotes the set of output perturbed system models; G ¯ p i is the transfer function of the actual model including its nominal model G p 1 , G p 2 and lumped uncertainty at the output represented by E O 1 , E O 2 . From Equations (38) and (39), Δ 1 , Δ 2 are normalized perturbation with Δ i ( j ω ) 1 ,   ω .
Figure 6. The M-Δ structure using multiplicative output uncertainties of a cascade control system.
Figure 6. The M-Δ structure using multiplicative output uncertainties of a cascade control system.
Fractalfract 08 00244 g006
According to multiplicative relative uncertainties, the weighting scalar W i captures the variation of uncertainties over frequencies and has to be chosen to ensure [36]:
W i ( j ω ) l I ( ω ) = max G ¯ p i Π o i G ¯ p i ( j ω ) G p i ( j ω ) G p i ( j ω ) ,   ω
To represent unmodelled dynamics, the common simple form of the multiplicative weight is as follows:
W i = τ w s + r 0 ( τ w / r ) s + 1
where r 0 is the relative uncertainties at low frequencies and r is the magnitude of the weight at high frequencies. It is obvious that W i is considered as a high-pass filter and, therefore, τ w is chosen as a cut-off frequency (1/ τ w ) to ensure that the weight W i has a larger magnitude than the largest magnitude of the relative uncertainties in Equation (40).
The structure singular value (denoted SSV, mu, or µ) is suggested to measure the robustness of the control system [36]. In Figure 6, the cascade control scheme with the multiplicative output uncertainty of each loop is rearranged in the M-Δ structure where Δ = d i a g ( Δ 1 , Δ 2 ) and M is a 2 × 2 matrix which includes all remaining blocks such as plants, controllers, and weight factors. To derive M, the relationships between y Δ and u Δ are obtained first:
y Δ 1 = W 1 G p 1 1 + G p 2 G c 2 + G p 1 G p 2 G c 1 G c 2 1 u Δ 2 G p 2 G c 1 G c 2 u Δ 1 y Δ 2 = 1 + G p 2 G c 2 + G p 1 G p 2 G c 1 G c 2 1 G c 1 G c 2 G p 2 W 2 u Δ 1 G p 2 G c 2 + G p 1 G p 2 G c 1 G c 2 W 2 u Δ 2
Therefore, each element of M can be derived as follows:
M 11 = W 1 G p 1 G p 2 G c 1 G c 2 1 + G p 2 G c 2 + G p 1 G p 2 G c 1 G c 2 1
M 12 = W 1 G p 1 1 + G p 2 G c 2 + G p 1 G p 2 G c 1 G c 2 1
M 21 = W 2 G p 2 G c 1 G c 2 1 + G p 2 G c 2 + G p 1 G p 2 G c 1 G c 2 1
M 22 = W 2 G p 2 G c 2 + G c 1 G c 2 G p 1 1 + G p 2 G c 2 + G p 1 G p 2 G c 1 G c 2 1
The µ-synthesis states that the control system is stable for all allowed perturbations with Δ i ( j ω ) 1 ,   ω if and only if μ M ( j ω ) < 1 ,   ω [36].

5. Simulation Study

In this paper, the performance indices including the integral absolute error (IAE), the integral time-weighted absolute error (ITAE), and the total variant (TV) are addressed and the smaller values of those indicate better performance. We use the motor and load parameters as seen in Table 1 [38].
From the parameters in Table 1 and Equation (10), the transfer function of the motor is derived:
G p 2 ( s ) = 33.1217 0.00001835 s 2 + 0.0468 s + 1

5.1. The Secondary Loop Control Design (Velocity Control)

From Equations (21) and (22), two controller parameters of the inner loop could be obtained. The desired time constant is chosen as τ c 2 = 0.001 and the fractional order of the integral term β = 1.2 . In this work, to guarantee the robustness of each single loop, the maximum sensitivity function (Ms) is adopted to play a part in tuning control parameters. The proposed algorithm (Section 3.4) is used in this case to obtain values of Ms with respect to ω as in Figure 7. From the figure, it can be seen that the Ms value converges to 1.2 when increasing ω . However, when ω is greater than 230 (rad/s), Ms becomes undetermined. Therefore, in this case, we choose ω = 200 (Ms = 1.232), and as a result of that, the control parameters are obtained as follows:
G c 2 s = 1.426 + 24.365 s 1.2

5.2. The Primary Loop Control Design (Position Control)

For the outer loop, in this work, two kinds of plants are considered as mentioned above, a sliding stable or feed drive system and a rotational mechanism. Using the parameters in Table 1, the transfer function of each case is obtained as follows:
Case   1 :   G p 1 ( s ) = 0.0016 s
In this case, α = 0.6 , τ c 1 = 0.03 , and from Equation (26), λ = 1.1 is chosen. From Equations (29) and (30), the control parameters are obtained:
G c 1 s = 12 , 196 + 26.0769 s 0.6
To justify the effectiveness of the proposed method, the controller settings of both control loops are simulated by providing a ramp function input. In addition, to investigate the effects of disturbances in terms of d1 for the primary output and d2 for the secondary output as in Figure 4, two step changes of these are inserted at times t = 1 (s) and t = 5 (s), respectively. Figure 8 and Figure 9 illustrate the closed-loop responses at the primary output (position) and secondary output (velocity) in this case. From the figures, it can be seen that the tracking control completely meets the servo performance requirements (fast response and no overshoot or oscillation). The disturbance of the secondary loop (d2) does not have any effect on the primary output. In contrast, at time t = 5, there is a variation in position due to d1; however, the primary controller maintains the performance in an excellent way. The control signals are shown in Figure 10 to prove that the proposed method has smooth control signals, which is the disadvantage of some advanced control techniques.
Moreover, the robustness of the proposed controller will be shown, as mentioned in Section 4. To obtain l I ( ω ) in Equation (40), the nominal parameters of the inner transfer function are changed in the range of ±50%. The corresponding relative errors ( G ¯ p 1 G p 1 ) / G p 1 are shown as functions of frequency in Figure 11a, and it can be seen that l I ( ω ) is 0.4 at low frequencies and 5 at high frequencies. τ w is chosen to guarantee that W 2 is large enough to satisfy W 2 > l I ( ω ) at all frequencies. From the figure, τ w = 0.0667 ; therefore, the transfer function of W 2 is obtained (from Equation (41)). For the outer loop (Figure 11b), it can be performed similarly, finally also obtaining the transfer function of W 1 in this case with τ w = 0.01 ,   r 0 = 0.4 , and r = 1.5 .
W 1 = 0.01 s + 0.4 ( 0.01 / 1.5 ) s + 1
W 2 = 0.0667 s + 0.4 ( 0.0667 / 5 ) s + 1
Figure 12 illustrates the structured singular value (SSV) of the M-Δ structure using From Equations (43)–(46), when two sources of parametric uncertainties of both control loops affect simultaneously, two scalar weight transfer functions W1 and W2 are as above. As mentioned in Section 4, the condition for robust stability is that the µ value has to be always less than 1. From the figure, it is obvious that the peak of µ is 0.4 over the frequency range, which ensures the robustness of the controlled system.
Case   2 :   G p 1 ( s ) = 2 s ( 0.0014 s + 1 )
In this case, α = 0.9 , τ c 1 = 0.02 , and from Equation (26), λ = 1.2 is chosen. From Equations (35) and (36), the control parameters are obtained:
G c 1 s = 8.8414 + 0.0115 s 0.9
In this case, the position output is the rotational angle (rad), and the ramp function input is adopted to verify the position control problem. Two disturbances described in Figure 4 are illustrated by adding step signals at time t = 1 (s) and t = 5 (s). Figure 13 and Figure 14 show the closed-loop responses at the primary output (position) and secondary output (velocity) in two cases. From the figure, it can be seen that the tracking control completely satisfies the servo performance requirements (fast response and no overshoot or oscillation). The disturbance in the inner loop has no effect on the primary output. This means that the effectiveness of the cascade scheme is clarified in this case. The primary controller plays an important part in maintaining the position in a servo manner. The control signals are shown in Figure 15. It can be seen that, in this case, the proposed method also has a smooth control signal, and only are there sudden changes in the manipulated variable due to the disturbances.
This is similar to case 1, to obtain l I ( ω ) in Equation (40), the nominal parameters of the inner transfer function are changed in the range of ±50%, approximately. Figure 16 illustrates the relative errors and it can be seen that l I ( ω ) is 0.5 at low frequencies and 1.1 at high frequencies. τ w is chosen as 0.002 to satisfy W 1 > l I ( ω ) at all frequencies. Therefore, the transfer function of W 1 is as Equation (55). Note that the transfer function of W 2 is still the same as the one in case 1.
W 1 = 0.002 s + 0.5 ( 0.002 / 1.1 ) s + 1
Figure 17 illustrates the SSV plot of the M-Δ structure in this case. From the figure, it can be seen that the peak of µ is 0.52 over the frequency range, which ensures the robustness of the controlled system.
The robust performance of the controlled systems is also evaluated by perturbing ±20% to all plant parameters. The performance indices are tabulated in Table 2 of nominal models and perturbed ones. From the data, it is obvious that the controlled systems still keep good tracking control in the presence of uncertainties, and the µ values are practically the same as the ones of the nominal models.

6. Conclusions

The fractional-order control designs for servo systems are studied in this work. To achieve the strict requirements of control performance in terms of fast response and no overshot/oscillation, the cascade control structure including two control loops, the secondary loop (inner loop) for velocity control and the primary loop (outer loop) for position control, is proposed. The FOPI and FOPD are suggested for the inner and outer, respectively, and the analytical design methods in the frequency domain are also proposed to derive the tuning rules for both controllers. In the inner loop, the maximum sensitivity function is adopted to choose the appropriate frequency to ensure the robustness of this control loop. The M-Δ structure is also addressed to justify the robustness of the controlled system using multiplicative output uncertainties, for which parametric uncertainties of both control loops are considered simultaneously. The simulation results indicate that the proposed method consistently satisfies the system requirements with a fast and no overshoot closed-loop time response to ramp function inputs. The performance indexes also clarify the effectiveness of the proposed method. The robust stability is also verified by the µ value, which has a peak value of less than 1 over the frequency range in all cases.

Author Contributions

Conceptualization, T.N.L.V. and L.H.G.; methodology, V.L.C. and T.N.L.V.; software, N.H.N.; validation, V.L.C. and T.N.L.V.; formal analysis, V.L.C.; investigation, V.L.C.; writing—original draft preparation, V.L.C.; writing—review and editing, V.L.C. and T.N.L.V.; visualization, V.L.C. and N.H.N.; supervision, T.N.L.V.; project administration, T.N.L.V.; funding acquisition, T.N.L.V. and L.H.G. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Ho Chi Minh City University of Technology and Education, grant number T2023-01NTĐ.

Data Availability Statement

Data sharing is not applicable to this article.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Lim, H.; Seo, J.W.; Choi, C.H. Position control of XY table in CNC machining center with non-rigid ball screw. In Proceedings of the 2000 American Control Conference, ACC (IEEE Cat. No. 00CH36334), Chicago, IL, USA, 28–30 June 2000; Volume 3. [Google Scholar]
  2. Franks, R.G.; Worley, C.W. Quantitative analysis of cascade control. Ind. Eng. Chem. 1956, 48, 1074–1079. [Google Scholar] [CrossRef]
  3. Chen, J.S.; Huang, Y.; Cheng, K. Mechanical model and contouring analysis of high-speed ball-screw drive systems with compliance effect. Int. J. Adv. Manuf. Technol. 2004, 24, 241–250. [Google Scholar] [CrossRef]
  4. Altintas, Y.; Verl, A.; Brecher, C.; Uriarte, L.; Pritschow, G. Machine feed drives. CIRP Ann. Manuf. Technol. 2011, 60, 779–796. [Google Scholar] [CrossRef]
  5. Hwan, S.S.; Kang, S.K.; Chung, D.H.; Stroud, I. Theory and Design of CNC Systems; Springer: Cham, Switzerland, 2008; pp. 157–185. [Google Scholar]
  6. Nakamura, M.; Goto, S.; Kyura, N. Mechatronic Servo System Control; Springer: Cham, Switzerland, 2004. [Google Scholar]
  7. Sato, R.; Tsutsumi, M. Mathematical model of feed drive systems consisting of AC servo motor and linear ball guide. J. Jpn. Soc. Precis. Eng. 2005, 71, 633–638. [Google Scholar] [CrossRef]
  8. Sato, R.; Tsutsumi, M. Modeling and controller tuning techniques for feed drive systems. In Proceedings of the ASME 2005 International Mechanical Engineering Congress and Exposition, Orlando, FL, USA, 5–11 November 2005; pp. 669–679. [Google Scholar]
  9. Liu, Y.; Wang, Z.Z.; Wang, Y.F.; Wang, D.H.; Xu, J.F. Cascade tracking control of servo motor with robust adaptive fuzzy compensation. Inf. Sci. 2021, 569, 450–468. [Google Scholar] [CrossRef]
  10. Lee, Y.; Park, S.; Lee, M. PID controller tuning to obtain desired closed loop responses for cascade control systems. Ind. Eng. Chem. Res. 1998, 37, 1859–1865. [Google Scholar] [CrossRef]
  11. Rao, A.S.; Seethaladevi, S.; Uma, S.; Chidambaram, M. Enhancing the performance of parallel cascade control using Smith predictor. ISA Trans. 2009, 48, 220–227. [Google Scholar] [PubMed]
  12. Uma, S.; Chidambaram, M.; Rao, A.S.; Yoo, C.K. Enhanced control of integrating cascade processes with time delays using modified Smith predictor. Chem. Eng. Sci. 2010, 65, 1065–1075. [Google Scholar] [CrossRef]
  13. Padhan, D.G.; Majhi, S. An improved parallel cascade control structure for processes with time delay. J. Process Control 2012, 22, 884–898. [Google Scholar] [CrossRef]
  14. Santosh, S.; Chidambaram, M. A simple method of tuning parallel cascade controllers for unstable FOPTD systems. ISA Trans. 2016, 65, 475–486. [Google Scholar] [CrossRef]
  15. Raja, G.L.; Ali, A. Modified parallel cascade control strategy for stable, unstable and integrating processes. ISA Trans. 2016, 65, 394–406. [Google Scholar] [CrossRef] [PubMed]
  16. Raja, G.L.; Ali, A. Smith predictor based parallel cascade control strategy for unstable and integrating processes with large time delay. J. Process Control 2017, 52, 57–65. [Google Scholar] [CrossRef]
  17. Neubauer, M.; Brenner, F.; Hinze, C.; Verl, A. Cascaded sliding mode position control (SMC-PI) for an improved dynamic behavior of elastic feed drives. Int. J. Mach. Tools Manuf. 2021, 169, 103796. [Google Scholar] [CrossRef]
  18. Sun, G.; Wu, L.; Kuang, Z.; Ma, Z.; Liu, J. Practical tracking control of linear motor via fractional-order sliding mode. Automatica 2018, 94, 221–235. [Google Scholar] [CrossRef]
  19. Wang, J.; Shao, C.; Chen, Y.Q. Fractional order sliding mode control via disturbance observer for a class of fractional order systems with mismatched disturbance. Mechatronics 2018, 53, 8–19. [Google Scholar] [CrossRef]
  20. Wang, S.; Yu, H.; Yu, J. Robust adaptive tracking control for servo mechanisms with continuous friction compensation. Control Eng. Pract. 2019, 87, 76–82. [Google Scholar] [CrossRef]
  21. Xu, H.; Hinze, C.; Lechler, A.; Verl, A. Robust µ parameterization with low tuning complexity of cascaded control for feed drives. Control Eng. Pract. 2023, 138, 105607. [Google Scholar] [CrossRef]
  22. Zou, M.; Yu, J.; Ma, Y.; Zhao, L.; Lin, C. Command filtering-based adaptive fuzzy control for permanent magnet synchronous motors with full-state constraints. Inf. Sci. 2020, 546, 1–12. [Google Scholar] [CrossRef]
  23. Nassef, A.M.; Abdelkareem, M.A.; Maghrabie, H.M.; Baroutaji, A. Metaheuristic-based algorithms for optimizing fractional-order controllers—A recent, systematic, and comprehensive preview. Fractal Fract. 2023, 7, 553. [Google Scholar] [CrossRef]
  24. Podlubny, I. Fractional-Order Systems and PIλDμ-Controllers. IEEE Trans. Autom. Control 1999, 44, 208–214. [Google Scholar] [CrossRef]
  25. Padula, F.; Visioli, A. Tuning rules for optimal PID and fractional-order PID controllers. J. Process Control 2011, 21, 69–81. [Google Scholar] [CrossRef]
  26. Vu, T.N.L.; Lee, M. Analytical design of fractional-order proportional-integral controllers for time-delay processes. ISA Trans. 2013, 52, 583–591. [Google Scholar] [CrossRef] [PubMed]
  27. Beschi, M.; Padula, F.; Visioli, A. Fractional robust PID control of a solar furnace. Control Eng. Pract. 2016, 56, 190–199. [Google Scholar] [CrossRef]
  28. Li, D.; Liu, L.; Jin, Q.; Hirasawa, K. Maximum sensitivity based fractional IMC-PID controller design for non-integer order system with time delay. J. Process Control 2015, 31, 17–29. [Google Scholar] [CrossRef]
  29. Vilanova, R.; Arrieta, O.; Ponsa, P. Robust PI/PID controllers for load disturbance based on direct synthesis. ISA Trans. 2018, 81, 177–196. [Google Scholar] [CrossRef] [PubMed]
  30. Yumuk, E.; Guzelkaya, M.; Eksin, I. Analytical fractional PID controller design based on Bode’s ideal transfer function plus time delay. ISA Trans. 2019, 91, 196–206. [Google Scholar] [CrossRef]
  31. Moradi, M. A genetic-multivariable fractional order PID control to multi-input multi-output processes. J. Process Control 2014, 24, 336–343. [Google Scholar] [CrossRef]
  32. Sánchez, H.S.; Padula, F.; Visioli, A.; Vilanova, R. Tuning rules for robust FOPID controllers based on multi-objective optimization with FOPDT models. ISA Trans. 2017, 66, 334–361. [Google Scholar] [CrossRef] [PubMed]
  33. Tumari, M.Z.M.; Ahmad, M.A.; Suid, M.H.; Hao, M.R. An improved marine predators algorithm-tuned fractional-order PID controller for automatic voltage regulator systems. Fractal Fract. 2023, 7, 571. [Google Scholar] [CrossRef]
  34. Shalaby, R.; El-Hossainy, M.; Abo-Zalam, B.; Mahmoud, T.A. Optimal fractional-order PID controller based on fractional-order actor-critic algorithm. Neural Comput. Appl. 2022, 35, 2347–2380. [Google Scholar] [CrossRef]
  35. Vu, T.N.L.; Chuong, V.L.; Truong, N.T.N.; Jung, J.H. Analytical Design of Fractional-Order PI Controller for Parallel Cascade Control Systems. Appl. Sci. 2022, 12, 2222. [Google Scholar] [CrossRef]
  36. Skogestad, S.; Postlethwaithe, I. Multivariable Feedback Control Analysis and Design, 2nd ed.; John Wiley & Sons: Hoboken, NJ, USA, 2001; pp. 253–351. [Google Scholar]
  37. Monje, C.A.; Chen, Y.; Vinagre, B.M.; Xue, D.; Feliu-Batlle, V. Fractional-Order Systems and Controls, Fundamentals and Applications; Springer: London, UK, 2010. [Google Scholar]
  38. Choudhary, K.; Qureshi, M.S.; Singh, B. Position Control of DC Motor using Improved Sliding Mode Control Techniques. Int. J. Appl. Eng. Res. 2018, 13, 15–19. [Google Scholar]
Figure 1. Block diagram of a DC servo motor.
Figure 1. Block diagram of a DC servo motor.
Fractalfract 08 00244 g001
Figure 2. The block diagram of the ball screw.
Figure 2. The block diagram of the ball screw.
Fractalfract 08 00244 g002
Figure 3. The block diagram of a rotational load.
Figure 3. The block diagram of a rotational load.
Fractalfract 08 00244 g003
Figure 4. The cascade control system structure.
Figure 4. The cascade control system structure.
Fractalfract 08 00244 g004
Figure 5. The control structure design of the primary control loop.
Figure 5. The control structure design of the primary control loop.
Fractalfract 08 00244 g005
Figure 7. The relationship between maximum sensitivity peak and frequency.
Figure 7. The relationship between maximum sensitivity peak and frequency.
Fractalfract 08 00244 g007
Figure 8. The translation responses to the ramp input with disturbances.
Figure 8. The translation responses to the ramp input with disturbances.
Fractalfract 08 00244 g008
Figure 9. The velocity responses of case 1 with and without disturbances.
Figure 9. The velocity responses of case 1 with and without disturbances.
Fractalfract 08 00244 g009
Figure 10. The applied voltages in both cases.
Figure 10. The applied voltages in both cases.
Fractalfract 08 00244 g010
Figure 11. (a,b) The scalar multiplicative weights of both control loops.
Figure 11. (a,b) The scalar multiplicative weights of both control loops.
Fractalfract 08 00244 g011
Figure 12. The SSV plot for robust stability of case 1.
Figure 12. The SSV plot for robust stability of case 1.
Fractalfract 08 00244 g012
Figure 13. The rotational responses to the ramp input with disturbances.
Figure 13. The rotational responses to the ramp input with disturbances.
Fractalfract 08 00244 g013
Figure 14. The velocity responses of case 2 with and without disturbances.
Figure 14. The velocity responses of case 2 with and without disturbances.
Fractalfract 08 00244 g014
Figure 15. Applied voltages for both cases.
Figure 15. Applied voltages for both cases.
Fractalfract 08 00244 g015
Figure 16. The scalar multiplicative weight of the primary loop of case 2.
Figure 16. The scalar multiplicative weight of the primary loop of case 2.
Fractalfract 08 00244 g016
Figure 17. The SSV plot for robust stability of case 2.
Figure 17. The SSV plot for robust stability of case 2.
Fractalfract 08 00244 g017
Table 1. System parameters.
Table 1. System parameters.
SymbolsValuesUnits
La8 × 10−5H
Ra0.316Ohm
Jm1.34 × 10−5kgm2
b1.82 × 10−5Nm/rad/s
Km0.03Nm/A
Ke0.03V/rad/s
l0.01m
JL6.75 × 10−5kgm2
DL1.82 × 10−4Nm/rad/s
Table 2. The performance indices in nominal and perturbed models for both cases.
Table 2. The performance indices in nominal and perturbed models for both cases.
Case 1Case 2
IAEITAETVµ (M)IAEITAETVµ (M)
Nominal0.05150.054321.060.405.6555.9711.9490.524
Perturbed (±20%)0.09620.369318.630.428.80129.8831.5430.5243
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

Chuong, V.L.; Nam, N.H.; Giang, L.H.; Vu, T.N.L. Robust Fractional-Order PI/PD Controllers for a Cascade Control Structure of Servo Systems. Fractal Fract. 2024, 8, 244. https://doi.org/10.3390/fractalfract8040244

AMA Style

Chuong VL, Nam NH, Giang LH, Vu TNL. Robust Fractional-Order PI/PD Controllers for a Cascade Control Structure of Servo Systems. Fractal and Fractional. 2024; 8(4):244. https://doi.org/10.3390/fractalfract8040244

Chicago/Turabian Style

Chuong, Vo Lam, Ngo Hong Nam, Le Hieu Giang, and Truong Nguyen Luan Vu. 2024. "Robust Fractional-Order PI/PD Controllers for a Cascade Control Structure of Servo Systems" Fractal and Fractional 8, no. 4: 244. https://doi.org/10.3390/fractalfract8040244

APA Style

Chuong, V. L., Nam, N. H., Giang, L. H., & Vu, T. N. L. (2024). Robust Fractional-Order PI/PD Controllers for a Cascade Control Structure of Servo Systems. Fractal and Fractional, 8(4), 244. https://doi.org/10.3390/fractalfract8040244

Article Metrics

Back to TopTop