Next Article in Journal
Double Redundancy Electro-Hydrostatic Actuator Fault Diagnosis Method Based on Progressive Fault Diagnosis Method
Previous Article in Journal
A Backstepping Controller Based on a Model-Assisted Extended State Observer for a Slice Rotor Supported by Active Magnetic Bearings
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Composite Sliding Mode Control of High Precision Electromechanical Actuator Considering Friction Nonlinearity

1
Department of Aerospace Engineering, Harbin Engineering University, Harbin 150001, China
2
School of Electronics and Information Engineering, Zhongyuan University of Technology, Zhengzhou 450007, China
3
Shanghai Electro-Mechanical Engineering Institute, Shanghai 201109, China
*
Author to whom correspondence should be addressed.
Actuators 2022, 11(9), 265; https://doi.org/10.3390/act11090265
Submission received: 16 August 2022 / Revised: 7 September 2022 / Accepted: 8 September 2022 / Published: 13 September 2022
(This article belongs to the Section Control Systems)

Abstract

:
Friction nonlinearity, which is common in electromechanical actuator (EMA) systems, leads to undesired dynamic responses such as “flat top”, low-speed crawl, which brings challenges to high precision attitude control of flight vehicles. In order to improve the robustness of the actuator control system under friction nonlinearity, and suppress the chattering caused by high gain of sliding mode control (SMC), a composite SMC scheme based on modified extended state observer (MESO) is proposed. Nonlinear MESO is adopted for estimating the nonlinear friction dynamics, unmodeled disturbance, and external real-time load dynamics so as to compensate for their adverse effect. At the same time, in order to improve the robustness of EMA, and reduce the tracking error of the servo system, SMC is adopted to ensure the tracking error convergence in a finite time. The stability of the proposed method is proved, and the effectiveness is verified by simulations.

1. Introduction

For high-precision guided flight, the dynamic performance of EMA plays an important role in the attitude and trajectory control of flight vehicles. However, in the traditional attitude control system, EMA is generally simplified as a rigid model, and the nonlinear factors are ignored, then, the model is simplified into a first-order or second-order ideal linear system with a certain bandwidth, and the subsequent attitude controller design is carried out based on it. However, in real EMA and other electromechanical drive systems, nonlinearities, disturbances, and uncertainties are widespread. These factors may come from inside or outside of the system, such as unmodeled dynamics, parameters time-varying, friction [1], backlash [2,3,4], dead-zone [5,6,7,8,9], load disturbances, etc. Therefore, in high-precision situations, these disturbances should not be ignored. The traditional controller, such as PID, can not meet the technical requirements of the modern flight vehicle under these disturbances.
Friction, which exists extensively in EMA, is the main factor affecting the low-speed tracking performance of servo systems, it can cause position tracking errors and speed distortion [10], limit cycles oscillation (LCO), and so on. In order to improve the tracking performance, a corresponding control method must be undertaken to suppress it. Methods based on models [11,12,13,14] and model-independent [15] methods were proposed by researchers in recent years. Model-based compensation methods rely on the friction model and parameter identification of the model, scholars have proposed a variety of friction models. As the LuGre model can describe the friction’s dynamic and static characteristics more completely. Compensation methods based on the LuGre friction model have been widely researched [15,16,17]. However, the friction nonlinearity of EMA is uncertain, and the parameter identification accuracy is highly dependent on the speed signal’s quality, so real-time parameter identification is not easy to be realized. The compensation method based on model-independent has been widely developed because of its simple control algorithm and good real-time performance. Wang [15] proposed an Active Disturbance Rejection Control method for friction compensation and control, and the Extended State Observer (ESO) was used to estimate the system disturbance. The results show that the method can improve the control accuracy to some extent. However, its performance in high bandwidth has not been verified yet. Sliding Mode Control (SMC) has been widely used due to its strong disturbance rejection capability, insensitivity to parameter changes, simplicity, and easy implementation [18]. However, the discontinuous switching characteristics of traditional SMC in nature will cause chattering, and LCO problems will occur in conventional SMC scheme with switching mechanism in the presence of unmodeled dynamics, especially when the second or higher order unmodeled dynamics exists [19], therefore, it is necessary to improve the controller. Fallaha [19] used the exponential function to design nonlinear reaching law, experiments on the robot arm show the effectiveness of the method in chattering suppression, particularly in steady-state. Shepit [20] proposed a state feedback sliding mode controller to achieve a complex eigen-structure assignment, this can be advantageous in providing maximum flexibility in specifying closed loop dynamics. However, the chatter phenomenon still exists. Difonzo [21] characterizes the attractivity for a co-dimension two discontinuity manifold in piecewise smooth differential systems, where a new parameter is proposed to indicate how the integral curve evolves, sliding or spiraling, but how to use this method to design chattering-free controllers is not discussed. For some complex nonlinear dynamic systems, constructing generalized nonlinear Lyapunov functions to establish system stability might be simpler [22], however, constructing the proper global Lyapunov function would be quite tough.
External disturbances, nonlinearity, and unmodeled disturbances are the main sources of chattering in sliding mode control [23]. Aiming at the friction nonlinearity and disturbance of electromechanical systems, a new adaptive sliding mode control scheme is proposed by Wang [24] based on a neural network to estimate the nonlinear and unmodeled dynamics of friction. The effectiveness is verified, but the real-time performance of the intelligent method is relatively poor in the case of high dynamic response. The ESO is a novel observe algorithm, which is able to estimate internal and external disturbance online [25], Deng [25,26,27] integrates linear ESO into controllers to develop high performance control of the systems with nonlinear dynamics or disturbances and unmeasurable signals, the unmeasured velocity, model uncertainty, mismatched disturbance, and bounded disturbance are estimated by ESO, and the estimated information is used to design the controllers, and the control effect is verified. Ren [28] proposed a friction compensation method based on an ESO, the effectiveness and robustness were verified in compensating for the robot friction effects. Different from robot control, the aircraft steering gear only has a position detection device, and the hinge torque changes in real-time during the flight.
In this paper, the effect of EMA friction nonlinearity on the attitude control system of an axisymmetric flight vehicle is analyzed, and the influence of friction on EMA servo dynamic characteristics and its compensation method are studied. In order to compensate for the nonlinear dynamics of EMA, a composite control strategy based on MESO and SMC is proposed. In this strategy, MESO is used to estimate nonlinear dynamics, load disturbances, and unmodeled dynamics and compensate for them in real-time, and SMC based on reaching law is used to enhance the robustness.
This paper is organized as follows, modeling and analysis of EMA are presented in Section 2, in Section 3, a composite controller based on MESO and SMC is designed for the EMA model, and the convergence of the system is proved. In Section 4, the performance of the proposed controller is compared with traditional SMC and Active Disturbance Rejection Control. Finally, conclusions are drawn in Section 5.

2. Modeling and Analysis of EMA

2.1. Modeling of EMA

As shown in Figure 1 [29], the simplified EMA system of the flight vehicle consists of DC brushless servo motor, control surface, drive control device, mechanical transmission mechanism, etc.
In Figure 1, the mechanical transmission section of EMA contains two stages of reducer, J m , J 1 , J 2 , J L are the rotational inertia of the motor, the corresponding gear inertia, and load inertia, respectively. θ m and θ 1 are the output angles of motor and gear 1, i 1 and i 2 are the transmission ratios of two gears, and δ is the deflection angle of EMA.
For the brushless DC motor, ignoring viscous damping parameters, The mathematical model of the EMA system without a controller can be written as follows [10]:
The phase voltage equation is
U a ( t ) = R a I a ( t ) + L d I a ( t ) d t + E a ( t )
where U a ( t ) is the control voltage of the motor, R a is the winding phase resistance, I a ( t ) is phase current, L is equivalent phase winding inductance, and E a ( t ) is phase induced electromotive force (EMF).
Phase EMF equation is
E a ( t ) = K E ω m ( t ) = K E d θ m ( t ) d t
where ω m ( t ) is the angular velocity of the motor, θ m ( t ) is the angle of the motor, and K E is the back-EMF coefficient.
The output electromagnetic torque of the motor is
T m ( t ) = K T I a ( t )
where T m ( t ) is the output electromagnetic torque, and K T is the torque coefficient.
The torque balance equation of the motor can be written as
T m ( t ) = J m d ω m ( t ) d t + T L
where T L is the equivalent load torque.
Apply the Laplace transform to Equations (1)–(4), and after ignoring T L , the open-loop transfer function of the motor output angle is obtained
G U a θ m ( s ) = K t K e K t s + J m ( L s + R a ) s 2 = 1 / K e s ( τ m τ e s 2 + τ m s + 1 )
where τ m = J m R a K e K t is the mechanical time constant, and τ e = L R a is the electrical time constant.
If the elasticity of the mechanical transmission system is neglected, the rotational inertia of the load and gears can be reduced to the motor side, and the total transmission ratio is i = i 1 i 2 . Then, the model of EMA can be simplified to be a rigid one.
In order to study the effects of EMA’s friction nonlinearity on flight control, first of all, the friction torque to EMA should be analyzed. It is assumed that the load torque of the EMA includes only including hinge torque T h and friction torque T f , then δ and T L can be represented as
δ = θ m i T L = T h + T f i
The LuGre model [30] is adopted for friction modeling and is expressed as
d z d t = ω L ω L g ω L z σ 0 g ω L = T C + T S T C e ω L / ω s 2 T f = λ σ 0 z + σ 1 d z   d t + σ 2 ω L
where z is the average deformation of bristles, T C is coulomb friction torque, T S is maximum static friction torque, σ 0 , σ 1 , σ 2 are bristle stiffness coefficient, bristle damping coefficient and viscous friction coefficient, ω L = δ ˙ is the deflection angular velocity of the EMA output, ω s is the Stribeck velocity and λ is the variation coefficient of friction torque [31]. These parameters can be obtained by identification method [32]. From Equations (1)–(7), the model block diagram of EMA can be obtained as in Figure 2.
Based on the diagram and equations, an open-loop model of EMA with friction can be obtained as
δ ¨ = a δ ˙ + b u + d
where a = K e K t + σ 2 R a J m 2 R a i 2 , b = K t J m 2 i R a , d = T h + σ 0 z + σ 1 z ˙ J m 2 i 2 , J m 2 is the equivalent rotating inertia of EMA, u is the control voltage U a . For EMA, d is related to friction, hinge torque, etc., and it is bounded. Parameters related to EMA and friction identification parameters are shown in Table 1.

2.2. Influence of Friction Nonlinearity on Attitude Control of Flight Vehicle

In-flight attitude control, EMA is the inner loop, and the attitude control is the outer loop, the nonlinearity of the inner loop will lead to unexpected dynamics.

2.2.1. Friction Nonlinearity Effects on EMA Dynamics under PID Control

For the EMA model, take position error e p = δ c δ , PID controller can be described as
u = k p e p + k i e p + k d e ˙ p
where δ c is the position command value, and k p , k i , k d are the proportional, integral, and differential coefficients, respectively. Parameters were designed by ITAE criterion [33]: k p = 33, k i = 895, k d = 0.15, take λ = 4. Ignoring hinge torque, dynamic responses of position and speed under friction nonlinearity of EMA can be seen in Figure 3 and Figure 4.
The simulation results show that, when the position command is small and the frequency is low, friction nonlinearity of EMA may cause the “flat top” phenomenon of position tracking under PID control, meanwhile, speed curve distortion is generated.

2.2.2. Influence of EMA Friction Nonlinearity on Attitude of Flight Vehicle

  • Flight vehicle model and controller design
The flight vehicle in this paper has an axisymmetric cross structure, the length of the body is 4.76 m, the characteristic length is 0.273 m, and the booster time is 8.2 s, and the maximum range is 300 km. During the flight of the vehicle, fuel consumption, flight altitude, speed, trajectory, and other parameters are changed in real-time, flight attitude transfer function model between equivalent actuator to pitch angle is
G δ ϑ = a 25 s + a 25 a 34 a 33 a 35 a 24 s 3 + a 34 a 33 a 22 s 2 + a 22 a 33 a 34 a 24 s + a 33 a 24
where a 22 = M z ω z J z , a 24 = M z α J z , a 26 = 1 J z , a 33 = g sin θ V , a 34 = P + Y α m V , a 35 = Y δ z m V are corresponding dynamic coefficients, J z is the moment of inertia, m is the mass of the vehicle at the feature point, V is the velocity vector, g is gravitational acceleration, θ is trajectory inclination angle, P is thrust, Y α is derivative of lift generated by the angle of attack, M z α , M z ω z , M z δ z are the angle of attack moment coefficient, angular velocity moment coefficient, and actuator moment coefficient, respectively.
Take 10 s after launching as the feature point. At this point, the dynamic pressure is large and the aerodynamic nonlinear factor is strong, it is suitable to study the influence of EMA friction nonlinearity on the attitude of the flight vehicle. First of all, without considering the nonlinearity of EMA, design the PID attitude controller of the feature point. By substituting aerodynamic parameters, the model of the flight vehicle in 10 s is
G = 183.8 s + 102.9 s 3 + 0.5668 s 2 23.11 s + 0.1197
The diagram of the attitude control model is shown in Figure 5.
It is assumed that the EMA has ideal dynamic characteristics, then the PID parameters of the controller can be designed as k p f = 4.54, k i f = 6.44, k d f = 0.13, k w = 0.2. Figure 6 shows the open loop and closed loop Bode diagrams of the longitudinal attitude transfer function of the model.
It can be seen from Figure 6a that, at this point, the model is a non-minimum phase system and the system is statically unstable. After compensation, see Figure 6b, the system is stable and the control bandwidth is 18.8 rad/s.
2.
Analysis of EMA friction nonlinearity influence on attitude control
In the presence of friction nonlinearity, the dynamic characteristics of EMA are considered in Figure 5. The influence of EMA friction nonlinearity on attitude control is studied by simulation. The result is shown in Figure 7.
As can be seen from Figure 7, the pitching angle results in unexpected dynamics due to the friction nonlinearity of EMA. With the friction torque increasing, the overshot of the pitch angle response increases slightly, and when λ = 8, the dynamic response of the pitch angle appears LCO. Therefore, for high precision flight control, the existence of EMA friction will lead to undesirable attitude dynamic responses under PID control, which is unacceptable. So, it is necessary to compensate for the friction nonlinearity of EMA.

3. Composite SMC Based on MESO

In order to enhance the disturbance suppression and improve its robustness, SMC is often used [10,18,34] to achieve high dynamic performance of EMA. However, conventional SMC methods [35] always contain a discontinuous sign function and large switching gains, so control chattering is unavoidable. External disturbances, nonlinearity, and unmodeled disturbances are the main sources of chattering in sliding mode control [23], ESO is a good way to estimate these disturbances.
In this section, for the EMA system considering friction nonlinearity, in order to eliminate chattering and enhance the robustness, a composite SMC based on MESO is proposed.

3.1. Design of Composite SMC Based on MESO

In this section, the proposed composite SMC control strategy based on MESO is introduced in detail. The proposed control diagram can be seen in Figure 8.

3.2. Design of MESO

For the EMA model with friction, as shown in Figure 2 and Equation (8), T f is friction torque, T h is load hinge torque, take x 1 = δ , x 2 = δ ˙ as state variables, and treating the unknown dynamics, x 3 as an expand state, then, expanded EMA model can be rewritten as:
x ˙ 1 = x 2 x ˙ 2 = a x 2 + d + b u = x 3 + b u x ˙ 3 = w ( t ) y = x 1
where, a = K e K t + σ 2 R a J m 2 R a i 2 , b = K t J m 2 i R a , d = T h + σ 0 z + σ 1 z ˙ J m 2 i 2 , J m 2 is the equivalent rotational inertia of the EMA system. x 3 = a x 2 + d = f ( x 2 , d , t ) , that is the overall disturbance including friction nonlinearity which represent the real-time action of acceleration, it is assumed that the first derivative of x 3 exits. and w ( t ) is bounded. Since the LuGre friction model can be linearized around zero, and several dynamic friction models satisfy this assumption, such as the Dahl model, the assumption is reasonable [36].
Because nonlinear ESO can make full use of the characteristics of nonlinear functions to further improve the performance of the observer, nonlinear ESO is a better way to estimate disturbances. The traditional nonlinear state observer of the system (12) is given by
e 1 = z 1 y z ˙ 1 = z 2 β 1 e 1 z ˙ 2 = z 3 β 2 f a l ( e 1 , α 1 , δ 1 ) + b u z ˙ 3 = β 3 f a l ( e 1 , α 2 , δ 1 )
Equation (13) is ESO designed by Han [37]. It can be seen that ESO is only related to the information of u and y , and it’s independent of the disturbance models. In (13), β 1 , β 2 , β 3 are positive gain coefficients, 0 < α i < 1 , i = 1 , 2 , δ 1 > 0 is the width of the linear interval, and nonlinear error feedback function:
f a l ( e 1 , α i , δ 1 ) = e 1 α i s i g n ( e 1 ) , e 1 > δ 1 e 1 δ 1 1 α i , e 1 δ 1
It has been proven [38] that all states can be tracked well by designing reasonable parameters and ESO can be converged. However, in Equation (14), the nonlinear error feedback function is a piecewise function about e 1 , at the switching point δ 1 , there would be a value mutation of z 2 and z 3 , then the disturbance cannot be compensated effectively. To solve this issue, at δ 1 , differential linearization of the error feedback function is adopted. Taking the first quadrant as an example, from Equation (14), we can obtain
f a l ( e 1 , a i , δ 1 ) e 1 = a i e 1 a i , e 1 > δ 1 δ 1 a i 1 e 1 ,   0 < e 1 δ 1
That is, the gain on the right side of the switching point δ 1 is a i times that on the left side. Therefore, a continuous error feedback function [37] is introduced to solve this problem:
f a c e 1 , α i , λ = e 1 α i 2 π arctan λ e 1
Equation (16) is continuously differentiable, so, when the function of Equation (14) in Equation (13) changed into Equation (16), its dynamic gain will be continuous. Characteristic curves of these two error feedback functions f a c ( e 1 , 0.25 , 10 , 000 ) and f a l ( e 1 , 0.25 , 0.001 ) , is shown in Figure 9.
As seen in Figure 9, outside the linear interval δ 1 , the gain characteristics of these two functions are close to the same. However, when the error is inside the linear interval, f a c ( ) has a larger gain, which will lead to a better compensation effect.
Based on Equations (13) and (16), the MESO was designed as:
e 1 = z 1 y z ˙ 1 = z 2 β 1 e 1 z ˙ 2 = z 3 β 2 f a c ( e 1 , α 1 , λ 1 ) + b u z ˙ 3 = β 3 f a c ( e 1 , α 2 , λ 2 )
It can be seen from Equation (17), only u and y information of the plant is required in MESO.

3.3. Stability Proof of MESO

Let
f a c ( e 1 , α i , λ i ) = f a c ( e 1 , α i , λ i ) e 1 e 1 ,   i = 1 , 2
And take f i ( e 1 ) = f a c ( e 1 , α i , λ i ) e 1 , i = 1 , 2 , then, Equation (17) can be converted into matrix form:
z ˙ 1 z ˙ 2 z ˙ 3 = β 1 1 0 β 2 f 1 ( e 1 ) 0 1 β 3 f 2 ( e 1 ) 0 0 z 1 z 2 z 3 + β 1 β 2 f 1 ( e 1 ) β 3 f 2 ( e 1 ) δ + 0 b 0 u
Take the Laplace transform of Equation (19) and according to Klem’s rule, the transfer function of the extended state z 3 can be obtained as:
z 3 ( s ) = β 3 f 2 ( e 1 ) ( s 2 δ b u ) s 3 + β 1 s 2 + β 2 f 1 ( e 1 ) s + β 3 f 2 ( e 1 )
Its characteristic equation is s 3 + β 1 s 2 + β 2 f 1 ( e 1 ) s + β 3 f 2 ( e 1 ) = 0 . And according to Routh-Hurwitz Criterion, the conditions of its stability are:
(1). The coefficients of the characteristic equation are all positive.
(2). β 1 β 2 f 1 ( e 1 ) β 3 f 2 ( e 1 ) > 0 .
Firstly, let’s analyze the dynamic characteristics of f i ( e 1 ) . We take α 1 = α 2 = 0.5 , and λ 1 = λ 2 = 10 6 , then, curve of f ( e 1 ) can be seen as Figure 10.
As shown in Figure 10, f ( e 1 ) > 0, the larger the error is, the closer f ( e 1 ) goes to zero. Gain coefficients β i are all positive, so, the first condition is satisfied. If β i are properly designed and satisfied the inequation β 1 β 2 > β 3 , the second condition is satisfied, then, the system stability is satisfied.

3.4. Composite SMC Design

Since the nonlinear dynamics of friction, load disturbance, and unmodeled dynamics can be estimated by the extended state z 3 in real-time by MESO, assuming that the estimation is accurate, all dynamics involving friction and hinge torques are treated as disturbances, and there is no need to model these factors in SMC controller.
Take e p = δ δ c , and the switching function
s 1 = c e p + e ˙ p
Then, define the Lyapunov function as
V = 1 2 s 1 2
Obviously, V 0 , For system (8) with frictional nonlinearity, a composite SMC control law based on MESO is designed as
u = δ ¨ c c e ˙ p z 3 k s 1 b
V ˙ = s 1 ( x 3 z 3 ) k s 1 2
Assuming that the estimation error of x 3 is within a certain range, and does not exceed a certain value of Δ, i.e., x 3 z 3 Δ , then Equation (24) can lead toan inequality:
V ˙ s 1 Δ k s 1 2 2 k 1 V + 1 2 Δ 2
According to the lemma of solving the inequality in literature [39], and assuming that α = 2 k 1 , f = 1 2 Δ 2 , then, for any finite value α , the solution of Equation (25) is
V ( t ) e α t t 0 V t 0 + f t 0 t e α ( t τ ) d τ ,   t t 0 0
Integrate Equation (26), and substitute α and f into it, so we have
V ( t ) e 2 k 1 t t 0 V t 0 + 1 2 2 k 1 Δ 2 1 e α t t 0
Then, when k 1 2 , V ( t ) 0 , therefore, lim t V ( t ) = 1 2 2 k 1 Δ 2 , the rate of convergence is related to control gain k and MESO gain β 3 , according to the discussion above, control parameters β 1 , β 2 , β 3 and k can be designed to get satisfactory results.
Based on the above analysis, (17) and (23) constitute the composite controller proposed in this paper.

4. Numerical Simulation

In order to verify the effectiveness of the proposed controller against friction nonlinearity and external disturbance of EMA, ESO based nonlinear proportional derivative controller (ESO-PD) in literature [15] and conventional SMC [23] are compared with the controller proposed in this paper, and the effectiveness of the controller is simulated under different conditions. The friction torque is simulated by the LuGre friction model, and the EMA and friction parameters are listed in Table 1. In order to make the simulation close to the actual situation, the voltage saturation limit is added according to the rated voltage. The control parameters of these three controllers are shown in Table 2.

4.1. Low Frequency Tracking Performance

Simulation conditions: λ = 4 , no hinge torque, at a frequency of 0.15 Hz, amplitude of 0.1° sinusoidal command signal. Simulation results can be seen in Figure 11a,b.
It can be seen from Figure 11a that under the effect of friction nonlinearity, the speed curve under SMC has obvious chattering, while speed curves under ESO-PD and MESO-SMC are smooth. Figure 11b shows control voltage is chattering obviously under SMC, while the other two control methods have no chattering. Voltage chattering will cause speed chattering and undesirable vibration of mechanical equipment.
According to Figure 12a,b, under the action of friction nonlinearity, these three methods can track the position command, and the steady-state accuracy are all better than 5‰. It can be seen from Figure 12c that under different control, the friction torques are different due to the different velocity dynamics; under SMC control, the friction torques are slightly different from those under the other two control modes due to the speed chattering. Figure 12d shows the system disturbance can be estimated under both ESO and MESO. The estimated values of EMA are different in the initial stage but converge in the steady state.

4.2. Step Response

In the presence of friction nonlinearity and λ = 1 , the command signal is unit step, and the maximum hinge torque of 30 N. m is added from 5 s to the end, its dynamic response is shown in Figure 13.
Figure 13a shows that the three control methods have good tracking performance. The ESO-PD control has slight lag and overshoot when adding the hinge torque. Figure 13b shows that ESO-PD has the largest overshoot, the speed response performance under MESO-SMC is similar to that of SMC. Figure 13c shows SMC control voltage has a large chattering, and this phenomenon does not exist under the other two controllers. Figure 13d shows that the estimates of the two observers differ, but the trend is similar.

4.3. High-Frequency Dynamic Response

Considering the friction nonlinearity of EMA and λ = 1 , and the maximum hinge torque of 30 N. m is added from 5 s to the end, after 10 Hz sinusoidal signal commanding, its dynamic responses are shown in Figure 14, Figure 15 and Figure 16.
As shown in Figure 14a,b. Considering friction nonlinearity, hinge torque load has little influence on the control accuracy of the EMA system under simulation conditions, and both the MESO-SMC control method and SMC have strong robustness and tracking accuracy. But position tracking has a 14.4° phase lag and about 1% amplitude attenuation under ESO-PD control. According to Figure 14c,d, voltage chattering still exists under SMC, it is also observed that the control voltage can respond quickly without chattering under nonlinearity and hinge torque load, and dynamic performances are greatly improved by MESO-SMC.
As can be seen from Figure 15, due to phase lag of velocity, its friction torque also lags under ESO-PD control.
As shown in Figure 16a,b, under the simulation condition of 10 Hz, the proposed MESO can better estimate the system disturbance in real time. According to Figure 16c,d, it can be seen that there is certain tracking lag and distortion in the speed curve under the ESO-PD control method.

5. Conclusions

In this paper, to deal with the negative influence of the EMA friction nonlinearity of the flight vehicle on the flight attitude control, the friction compensation method of EMA is studied, the composite control strategy of SMC based on the MESO is proposed, and its theoretical analysis is given. By comparing the dynamic responses of the three different control strategies under different control commands, the superiority of the proposed control strategy is verified, and the following conclusions are drawn:
(1)
The MESO-SMC method in this paper has high position tracking accuracy, and the “flat top” and speed curve distortion phenomenon caused by friction are effectively compensated for.
(2)
In the case of high-frequency response, the position dynamics governed by the ESO-PD controller have a certain phase lag, and both MESO-SMC and SMC have good position tracking performance.
(3)
MESO-based controller has good estimation performance and can effectively estimate and compensate for internal and external disturbances.
(4)
Although SMC has good robustness, it will produce control chattering or even speed chattering. MESO-SMC can effectively suppress sliding mode chattering and achieve high-precision robust control.
In this paper, only the friction nonlinearity influence on the high-precision control of EMA is considered, and the modeling of EMA is simplified. However, there are still other nonlinearities such as backlash and dead zone in the actual EMA, these factors will be considered in further studies.

Author Contributions

Conceptualization, J.X.; methodology, B.F., H.Q. and Y.Y.; software, B.F.; validation, H.Q. and J.X.; writing—original draft preparation, B.F.; review and editing, Y.Y.; supervision, H.Q. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

We would like to show gratitude to Shaobo Wang and Caixu Shen for their help during this work.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Wei, R.; Dong, Q.; Zhang, X.; Li, Z. Friction Compensation Control of Electromechanical Actuator Based on Neural Network Adaptive Sliding Mode. Sensors 2021, 21, 1508. [Google Scholar]
  2. Merzouki, R.; Cadiou, J.C. Estimation of backlash phenomenon in the electromechanical actuator. Control Eng. Pract. 2005, 13, 973–983. [Google Scholar] [CrossRef]
  3. Merzouki, R.; Davila, J.A.; Fridman, L.; Cadiou, J.C. Backlash phenomenon observation and identification in electromechanical system. Control Eng. Pract. 2007, 15, 447–457. [Google Scholar] [CrossRef]
  4. Sun, G.; Zhao, J.; Chen, Q. Observer-based compensation control of servo systems with backlash. Asian J. Control 2019, 23, 499–512. [Google Scholar] [CrossRef]
  5. Elmic, R.R.; Lewis, F.L. Deadzone compensation in motion control systems using neural networks. IEEE Trans. Autom. Control 2000, 45, 602–613. [Google Scholar] [CrossRef]
  6. Zuo, Z.; Li, X.; Shi, Z. L1 adaptive control of uncertain gear transmission servo systems with deadzone nonlinearity. ISA Trans. 2015, 58, 67–75. [Google Scholar] [CrossRef]
  7. Zuo, Z.; Ju, X.; Ding, Z. Control of Gear Transmission Servo Systems With Asymmetric Deadzone Nonlinearity. IEEE Trans. Control Syst. Technol. 2016, 24, 1472–1479. [Google Scholar] [CrossRef]
  8. Fu, B.; Qi, H.; Xu, J.; Yang, Y.; Wang, S.; Gao, Q. Attitude Control in Ascent Phase of Missile Considering Actuator Non-Linearity and Wind Disturbance. Appl. Sci. 2019, 9, 5113. [Google Scholar] [CrossRef]
  9. Deng, W.; Yao, J.; Ma, D. Robust adaptive precision motion control of hydraulic actuators with valve dead-zone compensation. Isa Trans. 2017, 70, 269–278. [Google Scholar] [CrossRef]
  10. Man, Z.; Mao, D.; Zhang, M.; Guo, L.; Gong, M. A Hybrid Control with PID–Improved Sliding Mode for Flat-Top of Missile Electromechanical Actuator Systems. Sensors 2018, 18, 4449. [Google Scholar] [CrossRef]
  11. Keck, A.; Zimmermann, J.; Sawodny, O. Friction parameter identification and compensation using the ElastoPlastic friction model. Mechatronics 2017, 47, 168–182. [Google Scholar] [CrossRef]
  12. Khayati, K.; Bigras, P.; Dessaint, L.A. LuGre model-based friction compensation and positioning control for a pneumatic actuator using multi-objective output-feedback control via LMI optimization. Mechatronics 2009, 19, 535–547. [Google Scholar] [CrossRef]
  13. Capace, A.; Merola, A.; Cosentino, C.; Amato, F. A Multistate Friction Model for the Compensation of the Asymmetric Hysteresis in the Mechanical Response of Pneumatic Artificial Muscles. Actuators 2019, 8, 49. [Google Scholar] [CrossRef]
  14. Lrinc, M.; Béla, L. Identification and Model-based Compensation of Striebeck Friction. Acta Polytech. Hung. 2006, 3, 45–58. [Google Scholar]
  15. Wang, S.; Wu, Q.; Yang, X. Active disturbance rejection control of friction for optoelectronic telescopes. Int. J. Simul. Syst. 2016, 17, 35.31–35.35. [Google Scholar] [CrossRef]
  16. Freidovich, L.; Robertsson, A.; Shiriaev, A.; Johansson, R. LuGre-Model-Based Friction Compensation. IEEE Trans. Control Technol. 2009, 18, 194–200. [Google Scholar] [CrossRef]
  17. Yu, H.; Gao, H.; Deng, H.; Yuan, S.; Zhang, L. Synchronization Control with Adaptive Friction Compensation of Treadmill-based Testing Apparatus for Wheeled Planetary Rover. IEEE Trans. Ind. Electron. 2021, 69, 592–603. [Google Scholar] [CrossRef]
  18. Utkin, V.I.; Chang, H.C. Sliding mode control on electro-mechanical systems. Math. Probl. Eng. 2002, 8, 1–23. [Google Scholar] [CrossRef]
  19. Fallaha, C.J.; Saad, M.; Kanaan, H.Y.; Al-Haddad, K. Sliding-Mode Robot Control With Exponential Reaching Law. IEEE Trans. Ind. Electron. 2011, 58, 600–610. [Google Scholar] [CrossRef]
  20. Shepit, B.M.; Pieper, J.K. Sliding mode control design for complex valued sliding manifold. In Proceedings of the 40th IEEE Conference on Decision and Control (Cat. No.01CH37228), Orlando, FL, USA, 4–7 December 2001. [Google Scholar]
  21. Difonzo, F.V. A note on attractivity for the intersection of two discontinuity manifolds. Opusc. Math. 2020, 40, 685–702. [Google Scholar] [CrossRef]
  22. Haddad, W.M.; Chellaboina, V. Nonlinear dynamical systems and control. In Nonlinear Dynamical Systems and Control; Princeton University Press: Princeton, NJ, USA, 2011. [Google Scholar]
  23. Liu, J.; Sun, F. Research progress of sliding mode variable structure control theory and algorithm. Control Theory Appl. 2007, 24, 407–418. [Google Scholar] [CrossRef]
  24. Wang, S.; Chen, Q.; Ren, X.; Yu, H. Neural network-based adaptive funnel sliding mode control for servo mechanisms with friction compensation. Neurocomputing 2020, 377, 16–26. [Google Scholar] [CrossRef]
  25. Deng, W.; Yao, J. Extended-State-Observer-Based Adaptive Control of Electro-Hydraulic Servomechanisms without Velocity Measurement. IEEE ASME Trans. Mechatron. 2019, 25, 1151–1161. [Google Scholar] [CrossRef]
  26. Deng, W.; Yao, J.; Yaoyao, W.; Xiaowei, Y.; Jiuhui, C. Output feedback backstepping control of hydraulic actuators with valve dynamics compensation. Mech. Syst. Signal Processing 2021, 158, 107769. [Google Scholar] [CrossRef]
  27. Deng, W.; Yao, J.; Ma, D. Time-varying input delay compensation for nonlinear systems with additive disturbance: An output feedback approach. Int. J. Robust Nonlinear Control 2018, 28, 31–52. [Google Scholar] [CrossRef]
  28. Ren, C.; Li, X.; Yang, X.; Ma, S. Extended State Observer based Sliding Mode Control of an Omnidirectional Mobile Robot with Friction Compensation. IEEE Trans. Ind. Electron. 2019, 66, 9480–9489. [Google Scholar] [CrossRef]
  29. Shin, W.H.; Lee, S.J.; Lee, I.; Bae, J.S. Effects of actuator nonlinearity on aeroelastic characteristics of a control fin. J. Fluids Struct. 2007, 23, 1093–1105. [Google Scholar] [CrossRef]
  30. Canudas, d.W.C.; Olsson, H.; Astrom, K.J.; Lischinsky, P. A new model for control of systems with friction. IEEE Trans. Autom. Control 1995, 40, 419–425. [Google Scholar] [CrossRef]
  31. De Wit, C.C.; Lischinsky, P. Adaptive friction compensation with partially known dynamic friction model. Int. Comm. Radiat. Units Meas. 1997, 11, 65–80. [Google Scholar] [CrossRef]
  32. Wei, Y.U.; Jia-Guang, M.A.; Jin-Ying, L.I.; Xiao, J. Friction parameter identification and friction compensation for precision servo turning table. Opt. Precis. Eng. 2011, 19, 2736–2743. [Google Scholar] [CrossRef]
  33. Awouda, A.; Mamat, R.B. New PID Tuning Rule Using ITAE Criteria. In Proceedings of the 2010 the 2nd International Conference on Computer and Automation Engineering (ICCAE), Singapore, 26–28 February 2010; pp. 171–176. [Google Scholar]
  34. Utkin, V.I. Sliding mode control design principles and applications to electric drives. IEEE Trans. Ind. Electron. 2002, 40, 23–36. [Google Scholar] [CrossRef]
  35. Gao, W.; Wang, Y.; Homaifa. Discrete-time variable structure control systems. IEEE Trans. Ind. Electron. 1995, 42, 117–122. [Google Scholar]
  36. Olsson, H.; Strm, K.J.; Wit, C.; Gfvert, M.; Lischinsky, P. Friction Models and Friction Compensation. Eur. J. Control 1998, 4, 5517–5522. [Google Scholar] [CrossRef] [Green Version]
  37. Han, J. The Technique for Estimating and Compensating the Uncertainties: Active Disturbance Rejection Control Technique, 4th ed.; National Defense Industry Press: Beijing, China, 2008. [Google Scholar]
  38. Guo, B.Z.; Zhao, Z.L. On the convergence of an extended state observer for nonlinear systems with uncertainty. Syst. Control Lett. 2011, 60, 420–430. [Google Scholar] [CrossRef]
  39. Petros, A.; Ioannou, J.S. Robust Adaptive Control; Prentice Hall: Hoboken, NJ, USA, 1995. [Google Scholar]
Figure 1. Schematic diagram of EMA system.
Figure 1. Schematic diagram of EMA system.
Actuators 11 00265 g001
Figure 2. Diagram of simplified EMA model with friction.
Figure 2. Diagram of simplified EMA model with friction.
Actuators 11 00265 g002
Figure 3. Position response under friction.
Figure 3. Position response under friction.
Actuators 11 00265 g003
Figure 4. Speed response under friction.
Figure 4. Speed response under friction.
Actuators 11 00265 g004
Figure 5. Block diagram of longitudinal attitude control.
Figure 5. Block diagram of longitudinal attitude control.
Actuators 11 00265 g005
Figure 6. Bode diagrams of the model at the feature point: (a) Open loop bode diagram. (b) Close loop bode diagram under PID control.
Figure 6. Bode diagrams of the model at the feature point: (a) Open loop bode diagram. (b) Close loop bode diagram under PID control.
Actuators 11 00265 g006
Figure 7. Pitch response under friction.
Figure 7. Pitch response under friction.
Actuators 11 00265 g007
Figure 8. Composite SMC control diagram of EMA.
Figure 8. Composite SMC control diagram of EMA.
Actuators 11 00265 g008
Figure 9. Error feedback function gain curves comparison.
Figure 9. Error feedback function gain curves comparison.
Actuators 11 00265 g009
Figure 10. Curve of f ( e 1 ) .
Figure 10. Curve of f ( e 1 ) .
Actuators 11 00265 g010
Figure 11. Speed and control voltage comparisons. (a) Speed response. (b) Control voltage.
Figure 11. Speed and control voltage comparisons. (a) Speed response. (b) Control voltage.
Actuators 11 00265 g011
Figure 12. Comparisons of low-frequency simulation: (a) Position tracking. (b) Position tracking errors. (c) Friction torque curves. (d) Disturbance estimation.
Figure 12. Comparisons of low-frequency simulation: (a) Position tracking. (b) Position tracking errors. (c) Friction torque curves. (d) Disturbance estimation.
Actuators 11 00265 g012
Figure 13. Step response comparations. (a) Position tracking. (b) Speed response. (c) Control voltage. (d) Disturbance estimation.
Figure 13. Step response comparations. (a) Position tracking. (b) Speed response. (c) Control voltage. (d) Disturbance estimation.
Actuators 11 00265 g013
Figure 14. High-frequency dynamic simulation comparisons. (a) Position tracking. (b) Local curves of (a). (c) Control voltage. (d) Local curves of (c).
Figure 14. High-frequency dynamic simulation comparisons. (a) Position tracking. (b) Local curves of (a). (c) Control voltage. (d) Local curves of (c).
Actuators 11 00265 g014
Figure 15. Friction torque curves.
Figure 15. Friction torque curves.
Actuators 11 00265 g015
Figure 16. Disturbance estimation and velocity response. (a) Disturbance estimation curves. (b) Local curves of (a). (c) Velocity response curves. (d) Local curves of (c).
Figure 16. Disturbance estimation and velocity response. (a) Disturbance estimation curves. (b) Local curves of (a). (c) Velocity response curves. (d) Local curves of (c).
Actuators 11 00265 g016
Table 1. EMA parameters.
Table 1. EMA parameters.
SymbolValueSymbolValue
i 315 K t 0.056   ( N m / A )
J m 3.6   × 10 6   ( kg m 2 ) R a 3.15   ( Ω )
J L 5.5   × 10 3   ( kg m 2 ) U 28.0   ( V )
K e 0.056   ( V / rad / s ) L 3.2   ( mH )
σ 0 7800   ( N m rad 1 ) σ 1 37.5   ( N m s rad 1 )
σ 2 0.3   ( N m s rad 1 ) ω s 0.063   ( / s )
T C 1.4   ( N m ) T S 2.04   ( N m )
Table 2. Parameters of controllers.
Table 2. Parameters of controllers.
ControllerParametersValueParametersValue
MESO-SMC c 230 λ 1 ,   λ 2 1.0 × 10 6
k 3570 β 1 1.5 × 10 3
b 12.5 β 2 4.17 × 10 3
α 1 ,   α 2 0.5 β 3 4.0 × 10 6
SMC c 230 d −50
k 500 d 50
ε 0.5 b 12.5
a −220.5
ESO-PD ω c 500 b 12.5
ω o 1500 ξ 1
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Fu, B.; Qi, H.; Xu, J.; Yang, Y. Composite Sliding Mode Control of High Precision Electromechanical Actuator Considering Friction Nonlinearity. Actuators 2022, 11, 265. https://doi.org/10.3390/act11090265

AMA Style

Fu B, Qi H, Xu J, Yang Y. Composite Sliding Mode Control of High Precision Electromechanical Actuator Considering Friction Nonlinearity. Actuators. 2022; 11(9):265. https://doi.org/10.3390/act11090265

Chicago/Turabian Style

Fu, Bangsheng, Hui Qi, Jiangtao Xu, and Ya Yang. 2022. "Composite Sliding Mode Control of High Precision Electromechanical Actuator Considering Friction Nonlinearity" Actuators 11, no. 9: 265. https://doi.org/10.3390/act11090265

APA Style

Fu, B., Qi, H., Xu, J., & Yang, Y. (2022). Composite Sliding Mode Control of High Precision Electromechanical Actuator Considering Friction Nonlinearity. Actuators, 11(9), 265. https://doi.org/10.3390/act11090265

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

Article Metrics

Back to TopTop