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

: 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.


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 singleinput, 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 M s 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.

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: 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):

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: Using Laplace transformation, Equation (3) can derive the following transfer function: where α i and β i are arbitrary, real, and positive.

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.

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 (kgm 2 ) ω: 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: where 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, R a : armature resistance (Ohm); L a : armature inductance (H) K e : back emf constant (V/rad/s); K m : torque constant (Nm/A) J m : inertial moment of the motor shaft (kgm 2 ) ω: 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: where . 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: sin γ cos γ ω ω ω where

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 (kgm 2 ) ω: 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: where 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: Case 2: a rotational mechanism The block diagram, in this case, can be simplified as Figure 3 [6]: transfer function in this case is obtained: Case 2: a rotational mechanism The block diagram, in this case, can be simplified as Figure 3 [6]: where L J is the inertia moment of the rotational load; L D is viscous damping; and θ is the rotational angle output (rad).From Figure 3, the transfer function of the primary plant is derived as follows: where

The Control System Structure
The general cascade control system is shown in Figure 4, where  According to Section 2.4, a dynamic second-order system is used for the secondary model Equation (10).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: where

The Control System Structure
The general cascade control system is shown in Figure 4, where G p1 (s) and G p2 (s) are the transfer functions of the primary and secondary plants, respectively.Additionally, as mentioned above, G p2 (s) represents the transfer function of the motor Equation ( 10) and G p1 (s) describes the transfer function of the load in both cases (translational and rotational movement); G c1 (s) and G c2 (s) are the primary and secondary controllers; and d 1 and d 2 denote load disturbances affecting each control loop, respectively.

Case 2: a rotational mechanism
The block diagram, in this case, can be simplified as Figure 3 [6]: where L J is the inertia moment of the rotational load; L D is viscous damping; and θ is the rotational angle output (rad).From Figure 3, the transfer function of the primary plant is derived as follows: where

The Control System Structure
The general cascade control system is shown in Figure 4, where  According to Section 2.4, a dynamic second-order system is used for the secondary model Equation (10).According to Section 2.4, a dynamic second-order system is used for the secondary model Equation (10).
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:

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: From Equation ( 16), the function of the controller is derived as follows: where G p2 is the model of the system in the secondary loop, ideally G p2 = G p2 ; (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: where τ c2 is the desired time response of the controlled loop and τ c2 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 c2 is obtained: The complex form in frequency domain of the controller is derived by substituting s = jω into Equation ( 19): Comparing Equations ( 8) and ( 20), the control parameters of the secondary control loop are derived:

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 r 1 to y 1 ) is described as Equation ( 23): Fractal Fract.2024, 8, x FOR PEER REVIEW 7 of 18 The block diagram of the system is reduced in Figure 5 where ( ) y r 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): ) Case 1: The primary plant is a feed drive system, so its transfer function is: 1 ( ) 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: ( ) Similarly to Equation ( 17), the primary controller is also derived: ) Figure 5.The control structure design of the primary control loop.
Case 1: The primary plant is a feed drive system, so its transfer function is: 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: Similarly to Equation (17), the primary controller is also derived: K l s ; and, in this case, where λ is the fractional order; τ c1 is the desired response time of the primary output.
Replacing Equation ( 26) and G p into Equation ( 25), the primary controller is obtained: Substituting s = jω to convert G c1 into the frequency domain in a complex form: Comparing Equations ( 7) and ( 28), the control parameters of the primary control loop are obtained: Case 2: The primary plant is a rotational mechanism, G p1 (s) = K L s(τs+1) In this case, the equivalent transfer function of the primary plant is: Therefore, the primary controller could be calculated as follows: , and (y 1 /r 1 ) d is chosen as Equation (26) (case 1); therefore, the primary controller is rewritten as follows: where Similarly to case 1, substituting s = jω to Equation (33) and rewriting G c1 in a complex form: Comparing Equations ( 7) and (34), the control parameters of the primary control loop are obtained in this case: Fractal Fract.2024, 8, 244 8 of 17

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: 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 M s can be larger than 1, which degrades the system performance.Therefore, the peak value of M s can be used to measure the robustness of the controlled system.In general, for both robust stability and performance, the M s value should be close to 1. Therefore, in most of the literature, the typical M s 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): Compute K p2 , K i2 according to Equations ( 21) and ( 22) respectively 4: Calculate M s for each set of control parameters 5: ω = ω + 1 6: end while 7: Choose appropriate ω to get value of M s being closed to 1.2 8: end

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: where Π Oi denotes the set of output perturbed system models; G pi is the transfer function of the actual model including its nominal model G p1 , G p2 and lumped uncertainty at the output represented by E O1 , E O2 .From Equations ( 38) and (39), ∆ 1 , ∆ 2 are normalized perturbation with |∆ i (jω)| ≤ 1, ∀ω.
Fractal Fract.2024, 8, x FOR PEER REVIEW 10 of 18 The M-Δ structure using multiplicative output uncertainties of a cascade control system.
According to multiplicative relative uncertainties, the weighting scalar i W captures the variation of uncertainties over frequencies and has to be chosen to ensure [36]: To represent unmodelled dynamics, the common simple form of the multiplicative weight is as follows: where 0 r is the relative uncertainties at low frequencies and r ∞ is the magnitude of the weight at high frequencies.It is obvious that i W is considered as a high-pass filter and, therefore, w τ is chosen as a cut-off frequency (1/ w τ ) to ensure that the weight i W 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 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: Therefore, each element of M can be derived as follows: ( ) 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]: To represent unmodelled dynamics, the common simple form of the multiplicative weight is as follows: 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 ∆ = diag(∆ 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: Therefore, each element of M can be derived as follows: (43) Fractal Fract.2024, 8, 244 10 of 17 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].

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 p2 (s) = 33.12170.00001835s 2 + 0.0468s + 1 (47)

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 τ c2 = 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 (M s ) 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 M s with respect to ω as in Figure 7. From the figure, it can be seen that the M s value converges to 1.2 when increasing ω.However, when ω is greater than 230 (rad/s), M s becomes undetermined.Therefore, in this case, we choose ω = 200 (M s = 1.232), and as a result of that, the control parameters are obtained as follows:

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:

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: In this case, α = 0.6, τ c1 = 0.03, and from Equation (26), λ = 1.1 is chosen.From Equations ( 29) and (30), the control parameters are obtained: 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 d 1 for the primary output and d 2 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.Figures 8 and 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 (d 2 ) does not have any effect on the primary output.In contrast, at time t = 5, there is a variation in position due to d 1 ; 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.

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:


is chosen.From Equations ( 29) and (30), the control parameters are obtained: 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.Figures 8 and 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 Section 4. To obtain () 40), the nominal parameters of the inner tran function are changed in the range of ±50%.The corresponding relative err are shown as functions of frequency in Figure 11a, and it can be s that () I l  is 0.4 at low frequencies and 5 at high frequencies.W is obtained (from Equation ( 41)).
the outer loop (Figure 11b), it can be performed similarly, finally also obtaining transfer function of 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 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 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.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 W 1 and W 2 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.range, which ensures the robustness of the controlled system.35) and (36), the control parameters are obtained: 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).Figures 13 and  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.In this case, α = 0.9, τ c1 = 0.02, and from Equation (26), λ = 1.2 is chosen.From Equations ( 35) and ( 36), the control parameters are obtained: G c1 (s) = 8.8414 + 0.0115s 0.9 (54) 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).Figures 13 and 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.002s + 0.5 (0.002/1.1)s +  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.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.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.

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 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 µ values are practically the same as the ones of the nominal models.

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.

Figure 1 .
Figure 1.Block diagram of a DC servo motor.

Figure 2 .
Figure 2. The block diagram of the ball screw.

Figure 1 .
Figure 1.Block diagram of a DC servo motor.

Figure 2 .
Figure 2. The block diagram of the ball screw.Figure 2. The block diagram of the ball screw.

Figure 2 .
Figure 2. The block diagram of the ball screw.Figure 2. The block diagram of the ball screw.

and 2 d
the transfer functions of the primary and secondary plants, respectively.Additionally, as mentioned above, 2 ( ) p G s represents the transfer function of the motor Equation (10) and 1 ( ) p G s describes the transfer function of the load in both cases (translational and rotational movement); denote load disturbances affecting each control loop, respectively.

Figure 3 .
Figure 3.The block diagram of a rotational load.

Figure 4 .
Figure 4.The cascade control system structure.

Figure 5 .
Figure 5.The control structure design of the primary control loop.

Figure 6 .
Figure 6.The M-∆ structure using multiplicative output uncertainties of a cascade control system.

Figure 7 .
Figure 7.The relationship between maximum sensitivity peak and frequency.

Figure 7 .
Figure 7.The relationship between maximum sensitivity peak and frequency.

Figure 7 .
Figure 7.The relationship between maximum sensitivity peak and frequency.

Figure 8 .
Figure 8.The translation responses to the ramp input with disturbances.

Figure 8 . 18 Figure 9 .
Figure 8.The translation responses to the ramp input with disturbances.Fractal Fract.2024, 8, x FOR PEER REVIEW 13 of 18

Figure 9 .
Figure 9.The velocity responses of case 1 with and without disturbances.

Figure 9 .
Figure 9.The velocity responses of case 1 with and without disturbances.

Figure 10 .
Figure 10.The applied voltages in both cases.

Figure 9 .
Figure 9.The velocity responses of case 1 with and without disturbances.

Figure 10 . is chosen to guarantee that 2 W is large enough to satisfy 2 2 W 1 W
Figure 10.The applied voltages in both cases.Moreover, the robustness of the proposed controller will be shown, as mentioned in Section 4. To obtain () I l  in Equation (40), the nominal parameters of the inner transfer function are changed in the range of ±50%.The corresponding relative errors

Figure 11 .
Figure 11.(a,b) The scalar multiplicative weights of both control loops.

Figure 12 .
Figure 12.The SSV plot for robust stability of case 1.

Fractal 18 Figure 13 .
Figure 13.The rotational responses to the ramp input with disturbances.

Figure 14 .
Figure 14.The velocity responses of case 2 with and without disturbances.

Figure 13 .
Figure 13.The rotational responses to the ramp input with disturbances.

Figure 13 .
Figure 13.The rotational responses to the ramp input with disturbances.

Figure 14 .
Figure 14.The velocity responses of case 2 with and without disturbances.

Figure 13 .
Figure 13.The rotational responses to the ramp input with disturbances.

Figure 14 .
Figure 14.The velocity responses of case 2 with and without disturbances.

Figure 16 .
Figure 16.The scalar multiplicative weight of the primary loop of case 2.

Figure 16 .
Figure 16.The scalar multiplicative weight of the primary loop of case 2.

Figure 16 .
Figure 16.The scalar multiplicative weight of the primary loop of case 2.

Figure 17 .
Figure 17.The SSV plot for robust stability of case 2.

Figure 17 .
Figure 17.The SSV plot for robust stability of case 2.

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.