An Optimized Fractional-Order PID Horizontal Vibration Control Approach for a High-Speed Elevator

Abstract

Due to factors such as uneven guide rails and airflow disturbance in the hoistway, high-speed elevators may experience significant vibrations during operation. This paper proposes an optimized fractional-order PID (FOPID) method to suppress vibrations of high-speed elevators. First, an accurate horizontal vibration model is established for the elevator car, in which the car frame and body are separate. Then, taking the control cost and the system performance as objective functions, we obtained an optimized FOPID controller based on multi-objective genetic algorithm optimization. Finally, the effectiveness of the controller in reducing elevator vibration was verified through numerical simulation. The results indicate that the horizontal acceleration controlled by the FOPID controller is reduced by about 68% compared to the case without a controller and about 25% compared to the conventional PID controller.

1. Introduction

As the development of urban construction advances, towers have increased rapidly, and high-speed elevators increase correspondingly. The speed rise of elevators causes the horizontal vibration to be more obvious. Large horizontal vibrations can affect the comfort and safety of elevators, leading to a decrease in elevator lifespan. Therefore, the vibration suppression of elevators is a key issue in the field of elevator research, which is significant for the development of the elevator industry.

The main factors causing the vibration are the defects of the guide rail or the non-standard installation and the influence of the airflow in the shaft. Passive damping and active damping are common methods for suppressing the horizontal vibration of elevators [1]. The passive damping method generally suppresses vibration by optimizing the structure or the car parameters. However, due to the complexity of the elevator operation and multiple factors, the effect of passive damping to suppress vibration is difficult to meet the demand. The active damping applies the active control force through the corresponding actuator and suppresses vibration through control algorithms [2]. Active damping is better than passive damping in suppressing vibration and meeting the demand, but it will increase the cost. Due to low speed and small vibration, low-speed elevators usually do not use active damping methods. However, high-speed elevators cannot ignore the vibration caused by speed, so an active damping method is needed. The usual executive structure of the active damping method is the active guide shoe, which plays an important role in reducing vibration and improving elevator ride comfort. It can reduce vibration by the force generated by the actuator installed on the passive guide shoe opposite to the vibration direction. Due to their good effect and strong adaptability, guide shoes have become an effective means of suppressing elevator vibration.

Cao et al. [3] established a 6-degree of freedom (DOF) horizontal vibration model for an elevator car and proposed an H2/generalized H2 hybrid control strategy. The results prove that the control strategy is effective. Feng et al. [4,5] simplified the guide shoe as a spring-mass damping system, treating the car body and car frame as a rigid body. Zhang et al. [6] designed a new type of shock absorber using a linear motor and designed a PID controller that combines a linear prediction model with a backpropagation neural network. From the results, it can be seen that the proposed control method is superior to traditional proportional integral differential controllers in suppressing vibration. Chen et al. [7] designed an optimal guaranteed performance state feedback controller based on a linear convex optimization method to minimize the H2 performance index and achieve the specified H∞ performance level. This control strategy is effective in vibration suppression and improving elevator comfort. Although the above scholars have tried to suppress the horizontal vibration, they have ignored the impact of treating elevators as rigid bodies on vibration response. He et al. [8] designed a pneumatic active guide shoe and a fuzzy neural network intelligent vibration reduction controller based on the Mamdani model, which can effectively suppress the horizontal vibration of high-speed elevators. Santo et al. [9] established a 3-DOF elevator car model, derived a dynamic equation for the impact of rail deformation on it, and proposed an SDRE control strategy that can effectively control system effects and help maintain the integrity of elevator components. Wang et al. [10] established a more practical elevator car horizontal vibration control model, proposed a predictive sliding mode controller based on the adaptive fuzzy-based (PSMC-AF) and designed a fuzzy logic system (FLS) to approach the compound uncertainty disturbance term (CUDT) online. Zhang [11] used the FA-BP algorithm to train fuzzy neural network (FNN) auxiliary controllers, and the elevator experiments have demonstrated the effectiveness. Tussett et al. [12] adopted a linear quadratic regulator (LQR) control strategy to control the magnetorheological damper, and the results validate the effectiveness. Zhao et al. [13] proposed an acceleration feedback-based active control method to suppress the vibration of elevators. Numerical simulations have demonstrated the effectiveness of the method. Wang et al. [14] proposed a semi-active high-speed elevator car horizontal vibration reduction method based on the improved particle swarm algorithm, and it has been verified through experiments. Due to the complexity of the controller designed using the above method, problems such as high cost or difficult implementation conditions may arise in practical use. The idea of combining fractional calculus theory with PID controllers has brought a new type of PID controller, known as the FOPID controller [15]. As a product of the combination of fractional calculus and PID controllers, FOPID not only has the advantage of simple control structure, strong adaptability, and independence of the controlled object model but also has better performance due to the introduction of more parameters. However, the design difficulty of FOPID controllers increases [16,17]. Some analytical and heuristic methods have been used in the design of FOPID controllers in different fields [18,19]. Some scholars use different algorithms for the design of FOPID controllers. Ataşlar [20] optimized the parameters of the FOPID controller using the GWO PSO algorithm. Idir [21] proposed the improved artificial bee colony (IABC) low-order approximation (LOA)-based fractional order PID (IABC/LOA-FOPID). Silaa [22] used the gray wolf optimizer (GWO) and extended gray wolf optimizer (EGWO) to obtain the appropriate parameters of the FOPID controller. Mok [23] proposes a modified smoothed function algorithm (MSFA)-based method to tune the FOPID controller of the AVR system. Ahmeda [24] designed a FOPID controller based on ILC (iterative learning control). The results demonstrate the effectiveness of the FOPID controller. Currently, no scholars have applied FOPID to elevator vibration suppression.

Regarding the above issues, this paper establishes a model that treats elevators not as rigid bodies but with the car body and frame separate, and designs a FOPID controller. Selecting system performance and control cost as objective functions, a multi-objective genetic algorithm (MOGA) is used to optimize the controller parameters. Through MATLAB/Simulink simulation verification, the FOPID controller can effectively suppress elevator vibration, which has a better effect than PID controllers.

The main contributions and innovations of this paper are as follows:

• An optimization design model for elevator vibration suppression is established, and an optimized FOPID controller is designed to suppress horizontal vibration effectively.
• Selecting control cost and system performance as the objective functions of the optimization algorithm, MOGA is used to optimize the parameters of both the PID controller and FOPID controller.
• Through multiple iterations, the proposed optimized FOPID controller can achieve the expected system performance at a limited control cost, which can reduce the horizontal acceleration by about 68% compared to the case without a controller and about 25% compared to PID control.

The remainder of the article is arranged as follows. In Section 2, a dynamic model of the horizontal vibration of high-speed elevators is established. Section 3 proposes an optimized FOPID control method based on a multi-objective genetic algorithm. Section 4 proposes numerical results and validation. Section 5 provides concluding remarks.

2. Dynamic Model of Horizontal Vibration of High-Speed Elevator

The elevator car is composed of a car body and a car frame. The car body is connected to the car frame through support rubber, and the four guide shoes guide the elevator to move. The guide shoes and the support rubber work together to reduce vibration. The elevator car is a relatively complex multi-degree of freedom vibration system. Referring to the experimentally validated model [14], we optimized the complex model [25] to obtain an elevator model that meets our requirements. The car system is simplified to the following lateral vibration dynamics model, as shown in Figure 1.

In the model, $m_c$ and $I_c$ are the car body’s mass and rotational inertia, $m_f$ and $I_f$ are the car frame’s mass and rotational inertia. $k_1$ and $c_1$ are the stiffness and damping coefficient of the active guide shoes, while $k_2$ and $c_2$ are the stiffness and damping coefficients of the supporting rubber. The active control force is $F_i(i=1\cdots 4)$; The guide shoe displacement excitation is $d_i(i=1\cdots 4)$;

$$X_a=\left[\begin{array}{c}{x}_1\\ x_2\end{array}\right]\hspace{0.17em}{X}_b=\left[\begin{array}{c}{x}_3\\ x_4\end{array}\right]\hspace{0.17em}{X}_c=\left[\begin{array}{c}{x}_c\\ {\theta }_c\end{array}\right]\hspace{0.17em}{X}_f=\left[\begin{array}{c}{x}_f\\ {\theta }_f\end{array}\right]$$

(1)

$x_i(i=1\cdots 4)$ is the horizontal displacement of corresponding points. $x_c$, $x_f$ are the displacement of the car body’s barycenter and frame’s barycenter, respectively. ${\theta }_c$, ${\theta }_f$ are the rotation angle of the car body’s barycenter and frame’s barycenter, respectively.

$$\begin{gathered} x_1=x_c+l_1{\theta }_c \\ {x}_2=x_c-l_2{\theta }_c \\ {x}_3=x_f+l_3{\theta }_f \\ {x}_4=x_f-l_4{\theta }_f \end{gathered}$$

(2)

$x_2=l_1$,$l_2$ is the distance between the rubber and the car body’s barycenter. $l_3$,$l_4$ are distance between active guide shoes and the car frame’s barycenter.

To facilitate writing, the following transformations are made:

$$X_a=H_c{X}_c\hspace{0.17em}{X}_b=H_f{X}_f$$

(3)

Including:

$$H_c=\left[\begin{array}{cc}1& l_1\\ 1& -l_2\end{array}\right],H_f=\left[\begin{array}{cc}1& l_3\\ 1& -l_4\end{array}\right]$$

(4)

Define transformation matrix:

$$H_P=\left[\begin{array}{cc}1& l_5\\ 1& -l_6\end{array}\right]$$

(5)

$l_5$ and $l_6$ represents the distance from the rubber to the car frame’s barycenter.

The dynamic equations for the car body and car frame in the horizontal direction are as follows:

$$M_c{\ddot{X}}_c={2H}_c^T(-K_2{X}_a-C_2{\dot{X}}_a+K_2{H}_P{X}_f+C_2{H}_P{\dot{X}}_f)$$

(6)

$$\begin{gathered} M_f{\ddot{X}}_f={2H}_f^T\left(-K_1{X}_b-C_1\dot{X_b}+K_1{D}_l+C_1\dot{D_l+}{K}_1{D}_r+C_1\dot{D_r}+F_l-F_r\right) \\ -{2H}_P^T(-K_2{X}_a-C_2{\dot{X}}_a+K_2{H}_P{X}_f+C_2{H}_P{\dot{X}}_f) \end{gathered}$$

(7)

In Equations (6) and (7), $M_c=diag(m_c,I_c)$, $M_f=diag(m_f,I_f)$, $K_1=diag(k_1,k_1)$, $K_2=diag\left(k_2,k_2\right),C_1=diag(c_1,c_1)$, $C_2=diag(c_2,c_2)$, $F_l=\left[\begin{array}{c}{F}_1\\ F_3\end{array}\right]$, $F_r=\left[\begin{array}{c}{F}_2\\ F_4\end{array}\right]$, $D_l=\left[\begin{array}{c}{d}_1\\ d_3\end{array}\right]$, $D_r=\left[\begin{array}{c}{d}_2\\ d_4\end{array}\right]$.

By selecting appropriate state variables, the above motion equation can be transformed into a system state equation. In this paper, the state equation is obtained by fusing the differential terms of rail displacement excitation into the state equation.

The lateral vibration state equation can be expressed as:

$$\dot{X}=AX+BU+EW$$

(8)

$$Y=CX+DU+E_{1}W$$

(9)

The state variable of the system is: $X=\left[X_b{X}_a{\dot{X}}_b-H_f{M}_f{H}_f^T{C}_1\left(D_l+D_r\right)\dot{X_a}\right]$; the active control of the system is: $U={\left[F_l{F}_r\right]}^T$; the disturbance of the system is: $W={\left[D_l{D}_r\right]}^T$; the output of the system is: $Y=\left[\begin{array}{ccc}\ddot{X_a}& X_d& X_e\end{array}\right]$; where:

$$X_d=X_b-\frac{(D_l+D_r)}{2}{X}_e=X_a-X_b$$

(10)

In Equation (10), $\ddot{X_a}$ represents the acceleration of the car body, $X_d$ represents the offset from the centerline, $X_e$ represents the offset between the car body and frame.

In Equations (8) and (9):

$$\begin{gathered} A=\left[\begin{array}{cccc}0& 0& I_{2\times 2}& 0\\ 0& 0& 0& I_{2\times 2}\\ -2H_f{M}_f^{-1}(H_f^T{K}_1+H_p^T{K}_2{H}_p{H}_f^{-1})& 2H_f{M}_f^{-1}{H}_p^T{K}_2& -2H_f{M}_f^{-1}(H_f^T{C}_1+H_p^T{C}_2{H}_p{H}_f^{-1})& 2H_f{M}_f^{-1}{H}_p^T{C}_2\\ 2H_c{M}_c^{-1}{H}_c^T{K}_2{H}_p{H}_f^{-1}& -2H_c{M}_c^{-1}{H}_c^T{K}_2& 2H_c{M}_c^{-1}{H}_c^T{C}_2{H}_p{H}_f^{-1}& -2H_c{M}_c^{-1}{H}_c^T{C}_2\end{array}\right] \\ B=\left[\begin{array}{cc}0& 0\\ 0& 0\\ H_f{M}_f^{-1}{H}_f^T& {-H}_f{M}_f^{-1}{H}_f^T\\ 0& 0\end{array}\right] \\ c=\left[\begin{array}{cccc}2H_c{M}_c^{-1}{H}_c^T{K}_2{H}_p{H}_f^{-1}& -2H_c{M}_c^{-1}{H}_c^T{K}_2& 2H_c{M}_c^{-1}{H}_c^T{C}_2{H}_p{H}_f^{-1}& -2H_c{M}_c^{-1}{H}_c^T{C}_2\\ \mathrm{I}& 0& 0& 0\\ -\mathrm{I}& \mathrm{I}& 0& 0\end{array}\right] \\ D=\left[\begin{array}{cc}0& 0\\ 0& 0\\ 0& 0\end{array}\right]\hspace{0.17em}{\mathrm{E}}_1=\left[\begin{array}{cc}2{\mathrm{H}}_{\mathrm{c}}{\mathrm{M}}_{\mathrm{c}}^{-1}{\mathrm{H}}_{\mathrm{c}}^{\mathrm{T}}{\mathrm{C}}_2{\mathrm{H}}_{\mathrm{p}}{\mathrm{M}}_{\mathrm{f}}^{-1}{\mathrm{H}}_{\mathrm{f}}^{\mathrm{T}}{\mathrm{C}}_1& 2{\mathrm{H}}_{\mathrm{c}}{\mathrm{M}}_{\mathrm{c}}^{-1}{\mathrm{H}}_{\mathrm{c}}^{\mathrm{T}}{\mathrm{C}}_2{\mathrm{H}}_{\mathrm{p}}{\mathrm{M}}_{\mathrm{f}}^{-1}{\mathrm{H}}_{\mathrm{f}}^{\mathrm{T}}{\mathrm{C}}_1\\ -\frac{\mathrm{I}}{2}& \frac{\mathrm{I}}{2}\\ 0& 0\end{array}\right] \\ E=\left[\begin{array}{cc}{H}_f{M}_f^{-1}{H}_f^T{C}_1& H_f{M}_f^{-1}{H}_f^T{C}_1\\ 0& 0\\ H_f{M}_f^{-1}\left(H_f^T\left(K_1-2C_1{H}_f{M}_f^{-1}{H}_f^T{C}_1\right)-2H_P^T{C}_2{H}_P{M}_f^{-1}{H}_f^T{C}_1\right)& H_f{M}_f^{-1}\left(H_f^T\left(K_1-2C_1{H}_f{M}_f^{-1}{H}_f^T{C}_1\right)-2H_P^T{C}_2{H}_P{M}_f^{-1}{H}_f^T{C}_1\right)\\ 2H_c{M}_c^{-1}{H}_c^T{C}_2{H}_p{M}_f^{-1}{H}_f^T{C}_1& 2H_c{M}_c^{-1}{H}_c^T{C}_2{H}_p{M}_f^{-1}{H}_f^T{C}_1\end{array}\right] \end{gathered}$$

3. Optimized FOPID Control Method Based on MOGA

3.1. Control Problem Description

The horizontal vibration suppression model of the elevator is shown in Figure 2:

In Figure 2, G (s) is the elevator vibration dynamics model established in the previous section, C (s) is the designed FOPID controller, and W (s) is the disturbance caused by rail irregularities.

The sensitivity of the human body to horizontal vibration is higher than that of vertical vibration[14]. Therefore, when suppressing vibration, it is sufficient to use the horizontal vibration as the control target.

The RMS of acceleration at the car body’s bottom and the relative displacement is selected as $J_1$. The RMS of the control force output by the actuator is taken as the control cost $J_2$.

$$J_1={‖w_1\begin{array}{ccc}\ddot{X_a}& w_2{X}_d& w_3{X}_e\end{array}‖}_2\hspace{0.17em}{J}_2={‖w_{4}U‖}_2$$

(11)

where $w_1$, $w_2$, $w_3$, $w_4$ are the regularization coefficients. We can adjust $w_1$, $w_2$, $w_3$ to change the weight of the corresponding output in the objective function $J_1$. $w_4$ is used to adjust $J_2$ for subsequent analysis.

The multi-objective integrated design model for elevator lateral vibration control is as follows:

$$\begin{gathered} min:J_1\hspace{0.17em}and\hspace{0.17em}{J}_2 \\ subjectto:\hspace{0.17em}p=max\left(Re\left(eig\left(LFT\left(G,C\right)\right)\right)\right)<0 \\ q\in \left(q_0,q_1\right) \\ Designvariable:q \end{gathered}$$

(12)

where $q$ represents the corresponding parameter of C (s), including Kp, Ki, Kd, the integral order λ and differential order $\mu$, constraint $p$ represents the maximum real part of the eigenvalue of $LFT\left(G,C\right)$ to ensure system stability; $q_0$ and $q_1$ represent the corresponding parameters’ upper and lower boundaries.

3.2. Fractional Order PID Controller

There are two main definitions of fractional calculus: Riemann-Liouville and Grunwald-Letnikov definitions.

Riemann-Liouville definition:

$$a{D}_t^{\alpha }f\left(t\right)=\frac{1}{\Gamma \left(m-\alpha \right)}{\left(\frac{d}{dt}\right)}^m{\int }_a^t\frac{f\left(\tau \right)}{{\left(t-\tau \right)}^{\alpha -m+1}}d\tau$$

(13)

In Equation (13), α is the order of the calculus, and a and t are the upper and lower limits of the integration. For $m-1$ < α < m, m∈ N, where $\Gamma \left(·\right)$ is Euler’s gamma function.

The Laplace transform of the RL fractional derivative is expressed as follows:

$$L\left\{a{D}_t^{\alpha }\right\}={\int }_0^{\infty }{e}^{-st}{D}_t^{\alpha }f\left(t\right)dt=s^{\alpha }F\left(s\right)-{\sum }_{k=0}^{m-1}{s}^{\alpha -k-1}{D}^{k}f\left(0\right)$$

(14)

When the initial condition is 0, Equation (12) can be expressed as:

$$L\left\{0D_t^{\alpha }\right\}=s^{\alpha }F\left(s\right)$$

(15)

FOPID’s transfer function is as follows:

$$G(s)=K_p+\frac{K_i}{s^{\lambda }}+K_d{s}^u$$

(16)

The integral order λ and differential order $\mu$ can take any value from 0 to 1. Conventional PID controllers are special cases where both integral and differential orders are taken as 1. The two order parameters bring more freedom to the controller’s adjustment. The parameter plane of the FOPID controller is shown in Figure 3.

3.3. Multi-Objective Genetic Algorithm

The vibration suppression problem of elevators is transformed into a multi-objective optimization problem shown in Equation (16), which achieves optimal system performance with minimal control costs. FOPID controllers introduce more parameters, which brings more degrees of freedom, but also increases the difficulty of controller design. Traditional PID controller design methods are challenging to achieve the desired results for FOPID controllers. Therefore, MOGA is selected to optimize the parameters of FOPID controllers.

Definition 1.

Pareto solution.

Select the objective function as $F\left(X\right)=\left\{f_1\left(x\right),f_2\left(x\right),f_3\left(x\right),\cdots \cdots f_k\left(x\right)\right\}$, $k\ge 2$; For two solution vectors $x_0$ and $x_1$ located in the solution space, $\exists i\in \left[1,k\right],f_i\left(x_0\right)>f_i\left(x_1\right)$, and $\forall i\in \left[1,k\right]$, $f_i\left(x_0\right)\ge f_i\left(x_1\right)$, It can be considered that the solution $x_0$ Dominant Solution $x_1$, or solve $x_0$ dominates solution $x_1$. If there is no other solution in the solution space that is superior to this solution, the solution is called a Pareto solution. A Pareto optimal solution set is represented as follows:

$$PS=\left\{u\in \mathsf{\Omega }|\nexists v\in \mathsf{\Omega },\mathrm{v}\succ \mathrm{u}\right\}$$

(17)

where $u$ is the solution vector and $\mathsf{\Omega }$ is the solution space.

Definition 2.

Noninferior solution sorting.

Noninferior solution ranking is important in MOGA, which determines the search direction of optimization algorithms. The specific process is as follows: 1. Let all solution spaces be S and compare them to find a noninferior solution set $S_1$. 2. Make S = S −$S_1$. Repeat the comparison process to find the non-inferior solution set $S_2$. 3. Repeat process 2 until S is an empty set. $S_1$ is called the first order noninferior solution, $S_2$ is called the second level noninferior solution, and so on. Obviously, the solution in $S_1$ is the best individual.

The process of MOGA is shown in Figure 4:

As mentioned in Section 3.1, the vibration suppression problem of elevators can be simplified into two control objectives: system performance and control cost. The parameters of the controller (Kp, Ki, Kd, etc.) are the variables that need to be optimized. Take two control objectives as the objective functions of a multi-objective genetic algorithm, and iteratively optimize all controller parameters as the population of the algorithm, ultimately obtaining the optimal controller parameters.

4. Numerical Results and Validation

To verify the effectiveness, a model was built using MATLAB/Simulink for simulation verification. The system is built in Simulink based on the state space equations given in Equations (8) and (9), and the simulation is completed by using the FOMCON toolbox [26] to import the FOPID controller.

4.1. Model Parameters

The specific parameters are shown in

Table 1. Parameters of the elevator model.

Parameter	Unit	Value
$m_c$	$\mathrm{k}\mathrm{g}$	1100
$m_f$	$\mathrm{k}\mathrm{g}$	2400
$I_c$	$\mathrm{k}\mathrm{g}·{\mathrm{m}}^2$	1600
$I_f$	$\mathrm{k}\mathrm{g}·{\mathrm{m}}^2$	8600
$k_1$	$\mathrm{N}/\mathrm{m}$	120,000
$k_2$	$\mathrm{N}/\mathrm{m}$	170,000
$c_1$	$\mathrm{N}·\mathrm{s}/\mathrm{m}$	2000
$c_2$	$\mathrm{N}·\mathrm{s}/\mathrm{m}$	2300
$l_1$	$\mathrm{m}$	1.6
$l_2$	$\mathrm{m}$	1.4
$l_3$	$\mathrm{m}$	3.3
$l_4$	$\mathrm{m}$	4.5
$l_5$	$\mathrm{m}$	1.25
$l_6$	$\mathrm{m}$	1.75

.

Most of the energy excited by the guide rail is located at low frequencies. The white noise signal passing through the low-pass filter is selected to simulate actual rail excitation. To approach the actual excitation signal of the guide rail more closely, select the appropriate amplitude of white noise. The white noise signal is output through the model in Simulink. Set the noise power to $1\times {10}^7$, the sample time to $1\times {10}^3$, and the seed to 23,341. The transfer function of the low-pass filter [26] is shown in Equation (18):

$$w_d=\frac{2\pi \times 10}{s+2\pi \times 10}$$

(18)

In practical applications, both the magnitude of the control force that activates the actuator and the travel of the mechanism should be considered while suppressing elevator vibration. Therefore, control costs should be limited so that the designed controller can meet actual needs. The optimal design model obtained after limiting the control costs is as follows:

$$\begin{gathered} min:J_{1}and{J}_2 \\ subjectto:p_1=max\left(Re\left(eig\left(LFT\left(G,C\right)\right)\right)\right)<0 \\ {p}_2=J_2<J_{2max} \\ q\in \left(q_0,q_1\right) \\ Designvariable:q \end{gathered}$$

(19)

The corresponding regularization coefficients are: ${\mathrm{w}}_1$ = 1000, ${\mathrm{w}}_2=10$, ${\mathrm{w}}_3$ = 0.02, $w_4$ = 0.001.

When using the MOGA algorithm to optimize the parameters, it is necessary to determine the search interval of the controller parameters first. Due to the elevator system being a complex MIMO system, the root locus method or the Routh Hurwitz method used to determine the search interval of the SISO system is no longer applicable. The initial interval of MIMO system parameters is generally determined using random area testing or decision statistics.

The FOPID controller designed in this paper is similar in structure to a PID controller. All parameters of the controller are nonnegative, and the lower bound of the initial search interval does not need to be searched. However, the designed controller has many parameters; it is hard or even unable to find a solution when all parameters are tested in a random area together. Referring to the parameter tuning method of conventional PID controllers, the parameter range of the designed controller is determined step by step, and the parameter range is gradually optimized in the order of proportional parameters, integral parameters, and differential parameters, as shown in Figure 5.

4.2. Results and Discussion

According to the above process, the parameters of the controller are optimized by using MOGA. Kp, Ki, Kd, the integral order λ, and differential order $\mu$ are chosen as variables which are all nonnegative. The upper bounds for λ and $\mu$ are 1. The upper bounds for Kp, Ki, and Kd change in [1$e^3$,1$e^4$]. The population size is 400, and the maximum algebra is 50. The function tolerance is $1\times {10}^{-4}$. The frontier of the Pareto solution obtained through the optimization is shown in Figure 6. The horizontal and vertical axes represent the two objective functions of the optimization algorithm. The horizontal axis represents system performance, while the vertical axis represents control cost. The curve in the Figure 6 consists of all feasible solutions obtained through optimization, namely the Pareto solution set, with each point representing a set of feasible controller parameters. Select three points $C_1$, $C_2$, $C_3$ on this curve. Three sets of controller parameters are brought into Simulink for simulation, and their differences in time domain response can be analyzed.

Figure 7 shows a comparison of the bottom horizontal acceleration of the car body with different controllers ($C_1$, $C_2$, $C_3$) and without control. According to the curve in Figure 7, compared to the case without control, the three groups of controllers can reduce the horizontal acceleration response. In the curve shown in Figure 6, the point near the vertical axis has better system performance and greater control costs. The controller parameters can be reasonably selected based on actual needs.

Figure 8 compares the relative centerline displacement of the car frame bottom under three sets of controllers and the situation without a controller. The image curve shows that the controller has almost no effect on the relative displacement because the relative displacement occupies a small weight in the system performance $J_1$. The controller does not increase the displacement, which meets the control expectations.

Figure 9 compares the bottom offset between the car body and frame under the control of three sets of controllers and without control. The curve results show that under the control of three sets of controllers (C1/C2/C3), the relative offset of the system decreases compared to the case without controllers.

Figure 10 compares the magnitude of vibration without a controller and controlled by optimized FOPID controller C2. The vibration of the elevator is mainly caused by the excitation of the guide rail, and the energy of the guide rail excitation is mainly concentrated in the low-frequency region. Figure 10 shows that the optimized FOPID controller can effectively reduce the vibration amplitude of the elevator compared to no controller. Figure 11 compares the output forces of three group controllers (C1/C2/C3). The output force corresponding to the controller is proportional to the control cost.

As shown in

Table 2. Comparison of corresponding indicators.

	Uncontrolled	C1	C2	C3
$J_1$	4.817 × 104	9.616 × 103	1.519 × 104	2.784 × 104
$J_2$	0	106.3202	77.7123	47.9846
Relative displacement of car frame bottom $\left(X_d\right)$	0.0020	0.0019	0.0019	0.0019
Relative offset ($X_e$)	8.5845 × 10−4	1.4406 × 10−4	2.3772 × 10−4	4.5487 × 10−4
Horizontal acceleration of the car body ($\ddot{x_2}$)	0.3933	0.0785	0.1240	0.2273
Control force ($F_1$)	0	647.8452	464.9079	283.5612
Control force ($F_2$)	0	75.9603	10.6733	19.4114
Control force ($F_3$)	0	573.9356	431.7551	271.3025
Control force ($F_4$)	0	14.4725	10.3735	19.6114

, the RMS values of system performance indicators, control cost, horizontal acceleration, and other parameters are compared. Taking controller C2 as an example, compared with those without a controller $J_1$ has been reduced by 68.46%, and the horizontal acceleration of the main control target car bottom has been reduced by 68.47%. The results show that the designed FOPID controller effectively suppresses elevator vibration compared to the uncontrolled situation, and the inhibition effect is related to the control cost.

Using the same optimization algorithm to optimize the conventional PID controller, a Pareto solution set like the FOPID controller is obtained (Figure 12). The two curves are shown below. Using the controller C2 designed above as the comparison object, the control effects between PID and FOPID are compared and analyzed.

The control cost of the designed FOPID controller C2 is 77. Select a point on the PID curve where the control cost is similar and bring the corresponding parameters and PID controller into the model for simulation analysis. Figure 13 compares the horizontal acceleration under the control of two types of controllers, respectively. From the curve in Figure 13, compared to the PID controller, the horizontal acceleration controlled by the FOPID controller is smaller. Figure 14 compares the displacement at the car frame’s bottom relative to the centerline under two control situations. The two controllers have almost the same control effect on the bottom displacement relative to the centerline. Figure 15 shows the bottom offset between the car body and the car frame under two control situations. From Figure 15, the relative deviation of FOPID is slightly smaller than that of PID. Figure 16 compares the output forces of the active mechanism under the two controllers. From Figure 16, the output forces of the two controllers are almost the same because they have similar control costs.

To sum up, both controllers can suppress elevator vibration, but when the output force is basically equal, the parameters such as elevator acceleration and car offset under optimized FOPID controller control are superior to optimized PID control. As shown in

Table 3. Comparison of corresponding indicators between PID and FOPID(C2).

	C2	PID
$J_1$	1.519 × 104	2.0229 × 104
$J_2$	77.7123	80.2407
Relative displacement of car frame bottom ($X_d$)	0.0019	0.0022
Relative deviation ($X_e$)	2.3772 × 10−4	3.2882 × 10−4
Horizontal acceleration of the car body ($\ddot{x_2}$)	0.1240	0.1652
Control force ($F_1$)	464.9079	489.1896
Control force ($F_2$)	10.6733	126.8436
Control force ($F_3$)	431.7551	429.4039
Control force ($F_4$)	10.3735	40.8524

, comparing the RMS of system performance indicators, control cost, horizontal acceleration, and other parameters under the two control conditions, compared to the PID controller, $J_1$ of FOPID controller decreased by 24.9%, the horizontal acceleration at the bottom of the car body decreased by 24.97%, the relative displacement at the bottom of the car frame decreased by 13.63%, and the relative bottom displacement between the car body and frame decreased by 27.71%. It can be found that FOPID is superior to PID control.

5. Conclusions

To suppress elevator vibration caused by rail irregularities, a FOPID controller optimized based on MOGA is proposed for elevator car systems using ARG active guide shoes. The main conclusions are as follows:

• Considering elevators’ vibration characteristics during actual operation, a dynamic model of the elevator with the car body and frame separated is established. White noise with a low-pass filter is used to simulate rail irregularities, and the simulation experiment is more accurate.
• An optimization design model for elevator vibration suppression is established, and MOGA is used to optimize the FOPID controller to limit the control cost and ensure the feasibility of the designed controller.
• From the results, the horizontal acceleration of the system has been reduced by 68% under optimized FOPID control compared to no control and reduced by about 25% compared to the optimized PID control, fully demonstrating the effectiveness of the designed controller.