Backstepping Adaptive Neural Network Control for Electric Braking Systems of Aircrafts

This paper proposes an adaptive backstepping control algorithm for electric braking systems with electromechanical actuators (EMAs). First, the ideal mathematical model of the EMA is established, and the nonlinear factors are analyzed, such as the deformation of the reduction gear. Subsequently, the actual mathematical model of the EMA is rebuilt by combining the ideal model and the nonlinear factors. To realize high performance braking pressure control, the backstepping control method is adopted to address the mismatched uncertainties in the electric braking system, and a radial basis function (RBF) neural network is established to estimate the nonlinear functions in the control system. The experimental results indicate that the proposed braking pressure control strategy can improve the servo performance of the electric braking system. In addition, the hardware-in-loop (HIL) experimental results show that the proposed EMA controller can satisfy the requirements of the aircraft antilock braking systems.


Introduction
With the development of aviation technology, the performance of modern aircraft has also greatly improved and the requirement of braking system performance is more stringent [1][2][3][4][5]. The electric braking system is a new type of braking system that uses electromechanical actuators instead of traditional hydraulic actuators. Compared with the hydraulic braking system, the electric braking system has advantages of light weight, small volume, and high reliability which make it the development direction of aircraft braking systems in the future [6][7][8][9].
As a new type of actuator, the electromechanical actuator (EMA) is composed of a motor, reduction gear, and a ball screw. With the development of permanent magnet material technology, the brushless DC motor (BLDCM) has been widely used in many industrial and automotive applications. The BLDCM has the advantages of a high-power density, a compact structure, and high efficiency which make it suitable for driving the EMA [10][11][12]. From the analysis above, it can be noted that the electric braking system is a new kind of servo system with special applications. Therefore, an appropriate control strategy is needed to satisfy the high-performance requirement of the aircraft braking system.
The EMA is a complex mechatronic transmission system that contains friction, backlash, hysteresis, and other nonlinear factors. These characteristics make it difficult to accurately describe the dynamic mathematical model of the EMA. In addition, aircrafts usually work in harsh environments, such as large fluctuations in temperature, strenuous vibrations, and damp conditions. For these reasons, it is impossible to obtain an accurate mathematical model of the EMA and is therefore difficult for the EMA to achieve the high precision and fast response of the braking pressure control effect with a traditional linear controller. To design an appropriate controller, the mathematical model of the EMA must be reasonably simplified, and many studies have been working on this issue [13][14][15]. The 2 of 16 compound disturbance in the EMA is a combination of the internal deformation, external disturbance, time-varying parameters, and unmodeled dynamics, so an appropriate design idea is also essential to design the control law of the electric braking system [16,17].
Peng et al. proposed a sliding mode control method for the electromechanical braking system with a BLDCM, and the switching control law was adjusted by a fuzzy corrector [18]. Li et al. presented a dynamic surface-based control algorithm for solving the stability problem of slip ratio control in an electric braking system [19]. To improve the control ability and accuracy of the motor in a flexible-joint robot, a fuzzy PID control method was proposed in [20]. To overcome the parameter variability and unknown disturbance in the electric rudder servo system, Lv et al. proposed a backstepping power fast terminal sliding mode control algorithm in [21]. Liang et al. [22] designed an RBF neural network-based nonsingular fast terminal sliding mode control algorithm for braking systems with an EMA. The global stability of the control law was proved by using the Lyapunov function. In [23], Chen et al. proposed a sliding-mode extremum seeking controller for the all-electric active braking system in unmanned aerial vehicles and validated the servo performance of the designed controller through HIL experiments. Lin et al. [24] realized a non-mechanical antilock braking system (ABS) controller for the electric scooter, and the braking performance for the ABS action was further addressed via experiments.
Furthermore, most of the existing studies have focused on improving the performance of the EMA controller, and the braking pressure-motor angle models in the literature are generally equivalent to a first-order proportional linear model. However, as shown in Figure 2, the real curve of the EMA is not the case. Therefore, the traditional first-order mathematical model is not accurate enough to describe the actual EMA model. On the basis of [22], the deformation of the reduction gear under braking pressure is analyzed and combined with the ideal EMA mathematical model, and subsequently, the actual EMA mathematical model is rebuilt. The backstepping control method has shown its effectiveness in dealing with systems with multiple dynamics and mismatched uncertainties, but a major problem is that certain functions must be linear in the unknown parameters and tedious analysis is needed to determine the regression matrices. The RBF neural network has strong input and output mapping functions, which make it suitable to use a stable neural network controller to estimate certain nonlinear functions. With the help of RBF neural networks, the linear-in-the-parameters assumption of the nonlinear function and the determination of the regression matrices can be avoided. According to the characteristics of the actual EMA model, an adaptive backstepping control method is proposed and the RBF neural network is used to address the uncertain compound disturbance of the electric braking system. Compared with the routine control strategies, the servo performance of the braking pressure and the adaptive ability of the control system are improved significantly.
The remainder of this paper is organized as follows: Section 1 introduces the control problems of the electric braking system in aircrafts. Section 2 presents the working principle of the electric braking system and structure of the EMA. Then, the deformation of the reduction gear is analyzed in detail, the actual EMA mathematical model is established, and the corresponding control target is put forward. Section 3 presents the backstepping adaptive controller and the design method of the RBF neural network, and the stability analysis of the proposed controller is carried out using Lyapunov functions. Section 4 discusses the experimental results and the hardware in loop experiments. The conclusions are summarized in Section 5.
The main contributions of this paper are summarized as follows: first, the paper analyzed the ideal EMA mathematical model and the nonlinear factors, and then the actual EMA mathematical model is rebuilt by considering the deformation of the reduction gear. Second, the proposed adaptive nonlinear control method can overcome the nonlinear factors of the electric braking system and improve the response speed and control precision of the EMA. Third, the control effect of the newly designed EMA controller is fully validated through the HIL experiments.

Mathematical Model of the EMA
The electric braking system is composed of a braking controller, an EMA, a brake wheel, and a brake disc. The structure of the electric braking system in the aircraft is shown in Figure 1. The braking controller is the control and drive unit of the braking system. Its functions include braking signal receiving, data processing, and control algorithm implementation. Then, the corresponding motor driving signal is generated and used to drive the EMA to realize the braking actions. The braking pressure feedback, motor current feedback, and the rotor position feedback signals are also transformed to the braking controller, which will be used to achieve the braking pressure closed-loop control. The electric braking system is composed of a braking controller, an EMA, a brake wheel, and a brake disc. The structure of the electric braking system in the aircraft is shown in Figure 1. The braking controller is the control and drive unit of the braking system. Its functions include braking signal receiving, data processing, and control algorithm implementation. Then, the corresponding motor driving signal is generated and used to drive the EMA to realize the braking actions. The braking pressure feedback, motor current feedback, and the rotor position feedback signals are also transformed to the braking controller, which will be used to achieve the braking pressure closed-loop control.
At the same time, the aircraft speed and wheel speed are used to calculate the slip ratio, which is usually used in the antilock braking system to prevent wheels from locking. The slip ratio control is also finished in the braking controller by controlling the output braking pressure signal. No more expatiation is needed here as there are many studies.

Ideal EMA Mathematical Model
The brushless DC motor (BLDCM) is usually controlled through a variable armature voltage. To simplify the mathematical model of the EMA, the following assumptions are made: the current fluctuation caused by the commutations is ignored, the magnetic condition of the motor is unsaturated, and the hysteresis losses and the eddy-current losses are ignored.
The voltage equation and torque equation of the BLDCM are expressed as: where a U is the armature voltage, R is the stator resistance, i is the armature current, L is the stator inductance, a E is the back electromotive force, emf k is the back electromotive force constant, J is the inertia moment, m ω is the mechanical speed of the motor, v B is the viscous damping coefficient, e T is the electromagnetic torque, T k is the torque constant, and L T is the load torque The relationship between the vertical motion of the ball screw and the rotary motion of the motor is shown as: where ema x is the axis displacement of the EMA, 0 L is the lead of the ball screw, andη is the transmission ratio of the reduction gear. By considering the essential characteristic of the electric braking system, the output braking pressure of the EMA can be calculated as: At the same time, the aircraft speed and wheel speed are used to calculate the slip ratio, which is usually used in the antilock braking system to prevent wheels from locking. The slip ratio control is also finished in the braking controller by controlling the output braking pressure signal. No more expatiation is needed here as there are many studies.

Ideal EMA Mathematical Model
The brushless DC motor (BLDCM) is usually controlled through a variable armature voltage.
To simplify the mathematical model of the EMA, the following assumptions are made: the current fluctuation caused by the commutations is ignored, the magnetic condition of the motor is unsaturated, and the hysteresis losses and the eddy-current losses are ignored.
The voltage equation and torque equation of the BLDCM are expressed as: ( where U a is the armature voltage, R is the stator resistance, i is the armature current, L is the stator inductance, E a is the back electromotive force, k em f is the back electromotive force constant, J is the inertia moment, ω m is the mechanical speed of the motor, B v is the viscous damping coefficient, T e is the electromagnetic torque, k T is the torque constant, and T L is the load torque. The relationship between the vertical motion of the ball screw and the rotary motion of the motor is shown as: where x ema is the axis displacement of the EMA, L 0 is the lead of the ball screw, and η is the transmission ratio of the reduction gear. By considering the essential characteristic of the electric braking system, the output braking pressure of the EMA can be calculated as: where p is the braking pressure of the EMA, c b is the stiffness coefficient of the brake disc and x b is the lateral displacement of the brake disc. If the lateral displacement of the brake disc is assumed to be zero, Equation (3) can be simplified as: The load torque of the motor (T L ) can be divided into two parts: part of it is used to overcome the resistance torque of the ball screw, and the remainder is used to overcome the drive torque of the ball screw. Then, the load torque can be calculated using the following equation: where T rt is the resistance torque of the ball screw and T dt is the drive torque of the ball screw.
In the EMA of the electric braking system, the resistance torque T rt is assumed to be a constant value, and the drive torque T dt is related to the output braking pressure p. In most studies, T dt is equivalent to a linear function of the braking pressure as: However, this simplification has a significant effect on the control performance of the electric braking system. By analyzing the electric braking system, it can be found that when the changing rate of the braking pressure . p is positive, the drive torque is a resistance to the motor. When the changing rate of the braking pressure . p is negative, the drive torque is a motivation to the motor, as shown in Equation (7): To decrease the chattering of the braking pressure control effect, a saturation function is introduced instead of the sign function. According to Equations (2) and (4), the load of the motor can be expressed as: where the saturation function is defined as: where τ is the boundary coefficient of the saturation function. By taking consideration of Equations (1)- (9), selecting the state variables as x i = [p, ω m , i], the ideal mathematical model of the EMA can be expressed as:

Deformation of the Reduction Gear
As seen from Equation (10), the relationship between the braking pressure p and the rotation angle of motor θ m is established by a linear function. However, the ideal EMA mathematical model ignored many nonlinear factors, such as the deformation of the gear and gear clearance. As described in the literature [25,26], the reduction gear will exhibit deformation under a heavy load which can eventually affect the relationship between the braking pressure and the motor position, rendering it different from the ideal model. After many experiments, the actual curves of the braking pressure and the motor position are drawn and the comparison of the ideal and actual curves is shown in Figure 2.
Algorithms 2019, 12, x FOR PEER REVIEW 5 of 16 ignored many nonlinear factors, such as the deformation of the gear and gear clearance. As described in the literature [25,26], the reduction gear will exhibit deformation under a heavy load which can eventually affect the relationship between the braking pressure and the motor position, rendering it different from the ideal model. After many experiments, the actual curves of the braking pressure and the motor position are drawn and the comparison of the ideal and actual curves is shown in Figure 2.
Rotation Angle of Motor (rad) Braking Pressure (kg) It was assumed that the deformation of the reduction gear under the braking pressure p is g δ , and g δ is divided into two parts: the contact deformation 1 g δ and the bending deformation 2 g δ .
When analyzing the contact deformation of the reduction gear 1 g δ , it can be simplified as a contact force model of two cylinders with parallel axes. If the contact force is small, then the two cylinders can be thought of as touching along a line parallel to their axes. As the contact force increases, the contact line deforms into a contact surface along with the elastic deformation of the material. The interface width g w can be determined by the gear parameters, and the interface length g l can be determined from the Hertz theory [27]: where g α is the angle between the acting force g f and the y axis, g ρ is the radius of curvature of the contact surface, g μ and g E are the elastic constants of the reduction gear.
Then, the stress in the reduction gear can be expressed as: The contact deformation of the reduction gear can be calculated using Hooke's law: For a fixed reduction gear, all of the parameters are constants except for the acting force g f .
Then, Equation (14) can be simplified as: where 1 g c is a constant, which is determined by the reduction gear parameters.
When analyzing the bending deformation of the reduction gear 2 g δ , the research should focus on the deformation of the meshing point so the distributed force along the tooth surface is simplified It was assumed that the deformation of the reduction gear under the braking pressure p is δ g , and δ g is divided into two parts: the contact deformation δ g1 and the bending deformation δ g2 .
When analyzing the contact deformation of the reduction gear δ g1 , it can be simplified as a contact force model of two cylinders with parallel axes. If the contact force is small, then the two cylinders can be thought of as touching along a line parallel to their axes. As the contact force increases, the contact line deforms into a contact surface along with the elastic deformation of the material. The interface width w g can be determined by the gear parameters, and the interface length l g can be determined from the Hertz theory [27]: where α g is the angle between the acting force f g and the y axis,ρ g is the radius of curvature of the contact surface,µ g and E g are the elastic constants of the reduction gear. Then, the stress in the reduction gear can be expressed as: The contact deformation of the reduction gear can be calculated using Hooke's law: For a fixed reduction gear, all of the parameters are constants except for the acting force f g . Then, Equation (14) can be simplified as: where c g1 is a constant, which is determined by the reduction gear parameters. When analyzing the bending deformation of the reduction gear δ g2 , the research should focus on the deformation of the meshing point so the distributed force along the tooth surface is simplified into Algorithms 2019, 12, 215 6 of 16 a concentrated force, and the reduction gear is regarded as a cantilever of the elastic material. The force component in the y axis of the acting force f g is represented as f gy : The bending deformation of the reduction gear under f gy in the y axis is expressed as: where r g is the radius of the gear base,α g is the angle between f g and the y axis, E g is the elasticity modulus of the gear, I g is the section modulus of the gear and ∆x g is the arm length of the equivalent bending moment. Then, the bending deformation of the gear is generated as: Similarly, all of the parameters are constants except for the acting force f g . Equations (17) and (18) can be simplified as: where c g2 is a constant, which is determined by the reduction gear parameters.
After the deformation analysis, the relationship between the braking pressure p and the acting force f g is needed. It is expressed as: Then, the deformation of the reduction gear under the braking pressure is generated as:

Actual EMA Mathematical Model and the Control Target
If there is no deformation in the reduction gear, then it can be generated from Equations (2) and (4) that the ideal relationship between the braking pressure p and the motor position ∆θ m is: Considering the deformation in the reduction gear by combining Equations (21) and (22), the relationship between the braking pressure and the motor angle can be expressed by a second order equation [28][29][30]. Then, it can be obtained by using on a least square fit in Figure 2 as: where c 1 and c 2 are constants, and they are determined by the EMA parameters. Finally, in the EMA mathematical model, the quadratic item in Equation (23) is treated as a disturbance item. By taking the state variables x i = [p, ω m , i] into the equations, the actual state space model of the EMA is expressed as: The corresponding control target is proposed as follows: design a braking pressure controller of the EMA mathematical model described in Equation (24) so that the braking pressure p can track the desired braking pressure p * without errors in finite time.
It can be found that the electric braking system is essentially a high-performance pressure servo system which is driven by a BLDCM, and the controller design needs to consider the model uncertainty and the disturbance from the system. The unknown function comes from the deformation of the EMA and the nonlinear load torque of the motor, where the disturbance is mainly caused by parameter perturbation and the unmodeled dynamics of frictions. Therefore, the control algorithm to be designed needs to overcome these unknown disturbances to achieve high performance pressure control.

Control Strategy
The backstepping control method can decompose the complex nonlinear system into several sub systems, and the corresponding Lyapunov function of the whole system can be derived step by step [31,32]. According to the characteristics of the actual EMA mathematical model, this paper proposes an adaptive backstepping control method that uses an RBF neural network to approximate and adaptively cancel the unknown parts of the system.
To facilitate the design of the controller, the assumption that the reference braking pressure signal p * is continuous and its derivative exists is made.

Controller Design
Based on the backstepping theory, the high order system is divided into three sub systems, and the Lyapunov function of the system is built [33,34]. The error variables of system are defined as: where e 1 , e 2 , e 3 are the errors,α 1 and α 2 are the virtual control variables of the sub system.
Step 1: For the first sub system, the Lyapunov function is selected as: Obtain the derivative of Equation (26) as: The virtual control variable α 1 is designed as: Take Equation (28) into Equation (27) as: .
wheref 1 (x) is the designed neural network function, which will be used to approximate the unknown continuous function f 1 (x), and ∆ 1 (x) = f 1 (x) −f 1 (x) is the estimated error.
Step 2: For the second sub system, the item e 1 e 2 in (29) needs to be offset, and the Lyapunov function is selected as: Obtain the derivative of Equation (30) as: .
wheref 2 (x) is the designed neural network function, which will be used to approximate the unknown continuous function f 2 (x), and ∆ 2 (x) = f 2 (x) −f 2 (x) is the estimated error.
Step 3: For the third sub system, the Lyapunov function is selected as: Obtain the derivative of Equation (34) as: .
The control law of U a is designed as: Take Equation (36) into Equation (35) as: wheref 3 (x) is the designed neural network function, which will be used to approximate the unknown continuous function f 3 (x), and After the control variable U a is designed, we need to design the neural network functionf i (x) i = 1, 2, 3 and keep the control system stable.

The RBF Neural Network
The structure of an RBF neural network is shown in Figure 3. In Figure 3, the RBF neural network consists of the input layer, the basis function layer (hidden layer), and the output layer [35,36]. The input layer is used to complete the information transfer, the basis function layer is the respective field, and the output layer is the outputs of the RBF neural network. The numbers of the input layer and the output layer are usually determined by the characteristics of the control system. In this paper, there are three neurons in the input layer and one neuron in the output layer. A k-means clustering algorithm is used on selecting the number of the hidden layer and the corresponding control parameters of the neural network [37], and the number of neurons in the hidden layer is eventually determined to be five.
The Gaussian function is chosen as the radial basis function, as shown in Equation (38): where i c is the center of the basis function, i d is the scaling factor of the basis function, and ( ) i h  is the base Gaussian function. From Equation (37), it can be found that there are three unknown functions 1 f , 2 f and 3 f , so we use three RBFs to approximate them as follows (39): Define 1 f , 2 f and 3 f as follows (40): whereˆi W is the weight vector estimation.
Similarly, define Z ,Ẑ and Z  as follows (41): In Figure 3, the RBF neural network consists of the input layer, the basis function layer (hidden layer), and the output layer [35,36]. The input layer is used to complete the information transfer, the basis function layer is the respective field, and the output layer is the outputs of the RBF neural network. The numbers of the input layer and the output layer are usually determined by the characteristics of the control system. In this paper, there are three neurons in the input layer and one neuron in the output layer. A k-means clustering algorithm is used on selecting the number of the hidden layer and the corresponding control parameters of the neural network [37], and the number of neurons in the hidden layer is eventually determined to be five.
The Gaussian function is chosen as the radial basis function, as shown in Equation (38): where c i is the center of the basis function,d i is the scaling factor of the basis function, and h i (·) is the base Gaussian function. From Equation (37), it can be found that there are three unknown functions f 1 , f 2 and f 3 , so we use three RBFs to approximate them as follows (39): where W i is the ideal neural network weight vector and W i = [w 1i w 2i w 3i w 4i w 5i ] T , w Mi (M = 1, 2, 3, 4, 5) is the weight value, H i is the radial basis vector, and 2,3,4,5) is the Gaussian function, ε i is the approximation error,i = 1, 2, 3, and Definef 1 ,f 2 andf 3 as follows (40): whereŴ i is the weight vector estimation. Similarly, define Z,Ẑ and Z as follows (41): Design the Lyapunov function as: where V 3 = 1 2 ξ T ξ, Z F is the Forbenius norm of Z, Γ is a positive definite matrix with proper Design the RBF neural network weights adaptive law as follows: .
where W T i = W T i −Ŵ T i and i = 1, 2, 3. Then we can obtain: Z, by using the RBF network weights adaptive law in (43), .
V is expressed as: If K min is the minimum eigenvalue of K, then we can know that K min ξ 2 ≤ ξ T Kξ, according to the Schwarz inequality, it can be found that Take (47) into (46), we can obtain: V ≤ 0 will be true, and the control system will be stable. By adjusting the parameters ε N , n and K min properly, the dynamic performance of the controller will be improved.
Compared to the traditional backstepping control method, the RBF backstepping control method in this paper has strong robustness, benefiting from the universal approximation character of the RBF neural network.

Results and Analysis
To verify the servo performance of the proposed control method, the signal generator is taken as the braking pressure reference signal source, and a scope is used to record the braking pressure signals.
For practical purposes, the EMA is installed in the aircraft wheel as shown in Figure 4a, and the HIL test bench is shown in Figure 4b.

Results and Analysis
To verify the servo performance of the proposed control method, the signal generator is taken as the braking pressure reference signal source, and a scope is used to record the braking pressure signals. For practical purposes, the EMA is installed in the aircraft wheel as shown in Figure 4a, and the HIL test bench is shown in Figure 4b.  First, we need to design the control parameters of the neural networks, including the cluster centers of the basis function, the scaling factors of the basis function, and the weight value vectors of the output layer. As mentioned above, a k-means clustering algorithm is adopted to calculate the cluster centers and the scaling factors, and the weight value vectors can be obtained by a pseudo-inverse method. The workflow of the k-means is shown as follows: (1) Selecting M different vectors as the initial cluster center (4) Modifying the cluster centers by calculating the average values of the samples in each subset. It is calculated as Equation (51): where ( ) j U k is the j-th cluster domain, and i N is the samples number of ( ) j U k .
(5) Let k=k+1, and turns to step 2). Then repeat the process until the change in ( ) j c k is less than the preset threshold.
After the cluster centers are determined, the scaling factors can be obtained as: First, we need to design the control parameters of the neural networks, including the cluster centers of the basis function, the scaling factors of the basis function, and the weight value vectors of the output layer. As mentioned above, a k-means clustering algorithm is adopted to calculate the cluster centers and the scaling factors, and the weight value vectors can be obtained by a pseudo-inverse method. The workflow of the k-means is shown as follows: (1) Selecting M different vectors as the initial cluster center c 1 (0), c 2 (0), . . . , c M (0), where M is the number of the hidden nodes.
(3) Classifying the similar sample points into a class. When Equation (50) is satisfied, X p is classified as the j * class. Then, all of the samples are divided into M subsets as U 1 (k), U 2 (k), . . . , U M (k), and each of them is a clustering domain, which can be represented by the corresponding cluster center.
(4) Modifying the cluster centers by calculating the average values of the samples in each subset. It is calculated as Equation (51): where U j (k) is the j-th cluster domain, and N i is the samples number of U j (k). (5) Let k = k + 1, and turns to step (2). Then repeat the process until the change in c j (k) is less than the preset threshold.
After the cluster centers are determined, the scaling factors can be obtained as:

RBF Neural Networks Appropximate Performance
Then, the approximate performance of the designed RBF neural networks is verified by real-time experiments. As the neural networks proposed in this paper is used to approximate the unknown functions in the EMA model, we took a set of actual braking pressure data to verify it. The braking pressure ranges from 20 kg to 190 kg, the actual braking pressure and the ideal braking pressure are as shown in Figure 5a. It can be found that there are some differences between the actual pressure and the ideal model, due to the nonlinear factors such as the gear deformation. So, the proposed neural networks are used to approximate the unknown functions in the EMA model and its approximate performance of the braking pressure error is shown in Figure 5b. The experiments results indicate that the neural networks can approximate the braking pressure error effectively. This is also the basis of the high-performance pressure control.

RBF Neural Networks Appropximate Performance
Then, the approximate performance of the designed RBF neural networks is verified by real-time experiments. As the neural networks proposed in this paper is used to approximate the unknown functions in the EMA model, we took a set of actual braking pressure data to verify it. The braking pressure ranges from 20kg to 190kg, the actual braking pressure and the ideal braking pressure are as shown in Figure 5a. It can be found that there are some differences between the actual pressure and the ideal model, due to the nonlinear factors such as the gear deformation. So, the proposed neural networks are used to approximate the unknown functions in the EMA model and its approximate performance of the braking pressure error is shown in Figure 5b. The experiments results indicate that the neural networks can approximate the braking pressure error effectively. This is also the basis of the high-performance pressure control.

EMA Servo Performance
Subsequently, the servo performance of the EMA is verified by experiments. The sine wave and square wave are taken as the pressure reference signals respectively, and the reference signal and actual pressure signal are both recorded by the scope. To provide a sufficient comparison, a traditional proportional integral derivative (PID) controller is designed according to the requirements of the electric braking system, and the transfer function is: Due to the particularity of the electric braking system, the parameters of the PID controller is designed in two stages: when the braking pressure is increasing, the parameters is 18 , and when the braking pressure is decreasing, the parameters is 12 In the proposed braking controller, three RBFs are designed. Figure 6a,b shows that the frequency of the pressure signal is 3 Hz, the amplitude is 300 kg and the offset is 30 kg. It can be found that the proposed control algorithm in this paper has better tracking accuracy and smaller signal phase lag than the PID control algorithm. In Figure 6c,d, the reference signal is a square wave, which requires a more precise control algorithm. The frequency is

EMA Servo Performance
Subsequently, the servo performance of the EMA is verified by experiments. The sine wave and square wave are taken as the pressure reference signals respectively, and the reference signal and actual pressure signal are both recorded by the scope. To provide a sufficient comparison, a traditional proportional integral derivative (PID) controller is designed according to the requirements of the electric braking system, and the transfer function is: Due to the particularity of the electric braking system, the parameters of the PID controller is designed in two stages: when the braking pressure is increasing, the parameters is K p = 18, K i = 0.02, K d = 0.5, and when the braking pressure is decreasing, the parameters is K p = 12, In the proposed braking controller, three RBFs are designed. For each Gaussian function, the parameters of c i and d i are designed as [−8, 3, 0, 3, 8] and 20. The initial weight value of each neural net in the hidden layer is selected as 0.90. The parameters in the adaptive law is designed as µ 1 = µ 2 = µ 3 = 8, Γ 1 = diag{80, 80, 80}, Γ 2 = diag{15, 15, 15}, Γ 3 = diag{0.9, 0.9, 0.9}, and n = 0.2. Figure 6a,b shows that the frequency of the pressure signal is 3 Hz, the amplitude is 300 kg and the offset is 30 kg. It can be found that the proposed control algorithm in this paper has better tracking accuracy and smaller signal phase lag than the PID control algorithm. In Figure 6c,d, the reference signal is a square wave, which requires a more precise control algorithm. The frequency is 2 Hz, the amplitude is 300 kg, and the offset is 30 kg. Compared with the traditional PID control, the proposed method can increase the braking pressure much more quickly without overshot and the static error is much smaller than that in the PID control.

HIL Experimental Results
In practical work, the braking pressure signal is more irregular which requires higher control performance of the EMA. To verify the control effects of the proposed EMA controller in an actual working environment of the aircraft, an aircraft electric braking system HIL test bench was built as shown in Figure 4b. The aircraft landing model and the antilock controller were established by MATLAB-Simulink and downloaded to the target board. Then, the reference pressure signal was output to the EMA controller and the corresponding braking pressure is realized by EMA. The main structure of the HIL test bench was analyzed in research [14] and no more detail is provided here. The control target of the HIL experiment is to control the slip ratio around the optimal value to maintain the maximum adhesive coefficient between the aircraft wheel and the runway. The slip ratio λ and adhesive coefficient μ are defined as: where x V is the aircraft speed, w V is the wheel speed, and 1 c , 2 c , 3 c are constants.

HIL Experimental Results
In practical work, the braking pressure signal is more irregular which requires higher control performance of the EMA. To verify the control effects of the proposed EMA controller in an actual working environment of the aircraft, an aircraft electric braking system HIL test bench was built as shown in Figure 4b. The aircraft landing model and the antilock controller were established by MATLAB-Simulink and downloaded to the target board. Then, the reference pressure signal was output to the EMA controller and the corresponding braking pressure is realized by EMA. The main structure of the HIL test bench was analyzed in research [14] and no more detail is provided here. The control target of the HIL experiment is to control the slip ratio around the optimal value to maintain the maximum adhesive coefficient between the aircraft wheel and the runway. The slip ratio λ and adhesive coefficient µ are defined as: where V x is the aircraft speed,V w is the wheel speed, and c 1 , c 2 , c 3 are constants.
The initial conditions of the aircraft speed and the wheel speed are V x = 72 m/s and V w = 72 m/s, respectively. The limitation of the maximum braking pressure in the EMA controller is set at 225 kg. In the u-y function, the optimal wheel slip ratio is λ = 0.16 and the maximum adhesive coefficient is µ = 0.75. All of the experimental data are saved by the upper computer and drawn by MATLAB.
The HIL experimental results are shown in Figure 7a-d. It can be clearly observed that the aircraft speed and wheel speed descend steadily under the action of the braking pressure. After a short time transient, the optimal slip ratio is well tracked and the corresponding adhesive coefficient remains stable around its maximum value without significant overshoots, which means that the proposed EMA controller can meet the requirements of the aircraft braking system. The initial conditions of the aircraft speed and the wheel speed are 72 x V  m/s and 72 w V  m/s, respectively. The limitation of the maximum braking pressure in the EMA controller is set at 225 kg. In the u-y function, the optimal wheel slip ratio is 0.16   and the maximum adhesive coefficient is 0.75   . All of the experimental data are saved by the upper computer and drawn by MATLAB.
The HIL experimental results are shown in Figure 7a-d. It can be clearly observed that the aircraft speed and wheel speed descend steadily under the action of the braking pressure. After a short time transient, the optimal slip ratio is well tracked and the corresponding adhesive coefficient remains stable around its maximum value without significant overshoots, which means that the proposed EMA controller can meet the requirements of the aircraft braking system.

Conclusions
Based on the analysis of the electric braking system, the ideal mathematical model of the EMA is established and the deformation of the reduction gear is analyzed. Then, the actual mathematical model is obtained by reasonably simplifying the nonlinear factors. A new RBF neural network-based adaptive backstepping design principle is adopted to address the multiple dynamics in the electric braking system. The Lyapunov function analysis indicates that the proposed control strategy guaranteed that the tracking error of braking pressure is uniformly ultimately bounded. The experimental results show that the servo performance of the EMA is improved significantly with the proposed control strategy compared with the traditional controller. Furthermore, the experimental results of the HIL test show that the performance of the proposed EMA controller can meet the requirements of the aircraft antilock braking system.

Conclusions
Based on the analysis of the electric braking system, the ideal mathematical model of the EMA is established and the deformation of the reduction gear is analyzed. Then, the actual mathematical model is obtained by reasonably simplifying the nonlinear factors. A new RBF neural network-based adaptive backstepping design principle is adopted to address the multiple dynamics in the electric braking system. The Lyapunov function analysis indicates that the proposed control strategy guaranteed that the tracking error of braking pressure is uniformly ultimately bounded. The experimental results show that the servo performance of the EMA is improved significantly with the proposed control strategy compared with the traditional controller. Furthermore, the experimental results of the HIL test show that the performance of the proposed EMA controller can meet the requirements of the aircraft antilock braking system.