Abstract

Friction is a nonlinear effect that occurs in all mechanical systems which may cause limit cycles, tracking errors, and other undesirable effects. Traditional static friction models cannot characterize all the friction situations. In recent years, neural network (NN) technique has been widely used to approximate the nonlinear function. In this paper, a new method which combines radical basis function neural network (BRFNN) with beetle antennae search (BAS) algorithm for modeling friction in a robotic joint is proposed. Velocity, load, and temperature are considered as the three factors that influence the static friction. It is shown that the proposed BAS-RBFNN possesses better performance in terms of faster convergence rate and higher accuracy.

1. Introduction

Friction can be defined as a kind of force resisting the relative motion of two surfaces in contact. It can be influenced by many factors, including greases, contact materials, temperature, relative velocity, wears, and load [1]. It may cause limit cycles, tracking errors, and other undesirable effects in control [2, 3]. Therefore, a precise friction model may considerably improve the efficiency for control purposes and other applications [47].

Generally, friction models can be divided into two categories: static models and dynamic models [8, 9]. In this paper, static friction is studied based on experiments in a robotic manipulator joint. Many static friction models have been proposed in the past, among which viscous and Coulomb are two typical ones describing the friction as a function of velocity:where is the friction force, is the relative velocity of the surfaces, is the Coulomb friction force, is the coefficient of the viscous friction, and is the sign function.

Stribeck [10] observed that the friction decreases as the velocity increases in the low-velocity Stribeck region and modeled the friction as below:where is the extra break-away force to start the movement and is the Stribeck velocity that defines the region in which such an effect is present. The Coulomb friction coefficient is linear with the load, whereas the viscous coefficient and Stribeck velocity keep steady when the load varies [11, 12]. Therefore, load should also be considered.

Then, the joint load extended friction model model was proposed:where is the load in the joint and and are the incremental parameters describing the changes in and , respectively.

In addition, recent research shows that temperature has great influence on friction. In [13], a new mathematical model which describes the velocity, temperature, and load torque dependent friction torque for collaborative industrial robot joints is built:

However, no model is more effective than the other since many factors can significantly affect the practical implementation and the performances of each scheme [6]. Therefore, a universal model is needed. In recent years, neural network (NN) has been widely employed to build the relationship, especially nonlinear relationship between inputs and outputs [14]. A friction compensation method using the fuzzy neural network for friction compensation is proposed in [15], where the simulation is conducted using a 3DOF planar robot manipulator. In [16], friction and the disturbance are identified by a neural network. The PID controller with neural network is designed. In [17], an adaptive control method for hybrid position/force control of robot manipulators based on neurofuzzy modeling is presented. Adaptive control with Artificial Neural Network- (ANN-) based compensation of dynamic uncertainties such as friction was proposed in [18]. Unlike many approaches for ANN-based control, parameter update equations for the ANN model and for the parametric adaptive model are driven by both the tracking error and the system identification error. In [19], an adaptive neural network friction compensator was presented for servocontrol of hard disk drives. The abovementioned control methods identify the friction and the disturbance simultaneously which leads to a decrease in the accuracy of friction models. Therefore, some researchers proposed a genetic algorithm-optimized BP neural network to model the relationship between velocity, joint load, and the static friction [20]. It is shown that the GA-BPNN model has a higher prediction precision and better generalization capability compared with the model. However, the effect of temperature is not considered in this model and the training speed of the BP network is slow and the convergence is difficult. Thus, we propose to use radical basis function neural network (RBFNN) which can increase the training speed by using a local transfer function so that only a few neurons have a nonzero response and become active to each input value. Numerical analysis has proved that compared with the BP neural network, RBFNN is superior in approximating continuous functions [21]. However, RBFNN also suffers from the problem that its performance depends heavily on the input data and initial values of the parameters. This paper employs the beetle antennae search (BAS) algorithm to optimize the parameters of RBFNN in order to minimize the required time and effort. Furthermore, temperature is also considered in this model.

The rest of the paper is organized as follows. In Section 2, the experimental method to obtain velocity, load, temperature, and friction is proposed. And the setting approaches of the velocity, load, and temperature are described in detail. The BAS-RBFNN model is presented in Section 3 and validated in Section 4. In Section 5, conclusions and future works are presented.

2. Measurement of Static Friction and Parameters Settings of Influence Factors

2.1. Measurement of Static Friction

For a manipulator, the dynamic model can be described aswhere is the inertial matrix, is the Coriolis and centrifugal terms, relates to the gravitational torque and is the friction torque in the joint, and is the input torque.

When the joint rotates in a constant velocity, the Coriolis force and inertial force can be ignored. Thus, the dynamic model can be expressed as follows:where is the moment obtained by the motor, is the angular acceleration of the motor, is the motor torque, is the friction torque, and is the gravitational torque which can be expressed aswhere is the position of a joint and and are the coefficients which can be obtained by the least-squared based method. Thus, equation (6) can be simplified as follows:

We make two trajectories in opposite directions with constant velocity and obtain the actual motor torque and . Since the gravitational torque is basically equal around the same , we have

Then, we have

Since the RV reducer is of high precision, friction forces in both directions are approximate. Dividing the two sides of equation (10) by 2 yields

2.2. Parameters Settings of Load, Velocity, and Temperature

In this paper, velocity, load, and temperature are considered as influence factors of friction torques. Thus, the data samples of friction should be obtained under different velocities, loads, and temperatures. Furthermore, the number of samples is of great significance for the training process and performance of the neural network. Too many or too few samples will lead to overfitting or underfitting.

For joint load samples, the gravitational torque is considered as the joint load and changes with different :

According to equation (7), the gravitational torque fitting curve in which and is in Figure 1. Then, we set 21 joint positions ranging from −60° to −20° with 2° step size to obtain the load samples.

Since the Stribeck phenomenon at a low speed is very prominent, the velocity sampling should be denser. In the viscous area, it can be distributed loosely [20]. Thus, the following 15 joint velocities are selected: 0.1, 0.3, 0.5, 1.0, 2.0, 4.0, 6.0, 9.0, 13, 18, 24, 31, 38, 44, and 50 degree/s.

For temperature samples, the range should be wide enough for current application. Therefore, it is set from 5°C to 45°C with 2°C step size.

The joint torque is obtained aswhere is the motor current and is a torque constant. The sampling frequency is 0.25 kHz.

In the sampling process, measurement noise, which may cause uncertainty and bias errors, is assumed to be an additive normally distributed zero-mean stochastic disturbance. Therefore, we set the joint move periodically for 10 times. All the data samples are obtained by averaging. Moreover, the current data in the joint motor is further filtered by median filtering.

3. The Design of BAS-RBFNN for Friction Modeling

3.1. Design of RBFNN for Modeling the Static Friction

BPNN is the most commonly used artificial neural network (ANN) due to good generalization capability and robustness. Compared with BPNN, RBFNN has the advantages of approximation, category, and convergence velocity.

The radial basis function neural network which is composed of three layers is a kind of feed-forward neural network. The input layer transfers the data to the hidden layer and the hidden layer is the basis function layer which is composed of a basis function. The most common basis function is a Gaussian function as below:where is a vector of n dimensions, is the center of the unit in the hidden layer, determines the extent of the basis function, and means the distance between x and . Only a small part of inputs is activated near .

The output can be expressed as the linear-weighted sum of the basic functions in the hidden layer, which is shown as below:

The error of neuron in the output layer iswhere and are the desired and actual value of the neuron in the output layer.

The total error of the output layer iswhere m is the number of neurons in the output layer.

The experiment takes load, velocity, and temperature as inputs and static friction torques as the outputs. The structure of the RNN is shown in Figure 2.

3.2. Beetle Antennas Search Algorithm

The beetle antennae search algorithm (BAS) is a metaheuristic algorithm that is inspired by the searching behavior of longhorn beetles [22, 23]. A vector at tth time instant () denotes the position of the beetle. At position x, is defined as the fitness function to represent the concentration of odour. Normally, we can use two rules to simplify the algorithm, including search and detection behavior. In an unknown environment, the beetle searches randomly. A normalized random unit vector is utilized to model the searching behavior as follows:where represents a random function and m is the number of dimensions of the position. and are defined to imitate the searching behavior of both the right and left antennae:where d is the sensing length of antennae.

An iterative form is proposed to model the detecting behavior:where is the step size and represents a sign function.

The antennae length and step size are updated according to the following formula:

3.3. The RBFNN Optimized by BAS

The validity and accuracy of RBFNN may be reduced if the weights and biases are improperly selected. The BAS algorithm is used to optimize the network parameters of RBFNN. In this section, the BAS-RBFNN method is proposed to overcome the shortcomings of low accuracy. The process is shown in Figure 3.

4. Validation and Discussion

All experiments were made on joint 2 of Robot SIASUN SR4C as seen in Figure 4(a). The D-H model is in Figure 4(b). The desired temperature was carefully held constant in a thermostatic laboratory, as shown in Figure 5. Before the experiments, the desired equilibrium joint temperature should be reached by running the second joint continuously.

4.1. The Effects of Velocity, Load, and Temperature on Static Friction

In this section, the effects of velocity, load, and temperature on the static friction are investigated, respectively. Based on the constant-speed tracking experiment, two combined factors dependence of friction is plotted in Figure 6.

As shown in Figure 6(a), the Stribeck phenomenon is very prominent and the viscous friction plays a heavy role at high speed [9]. The static friction decreases obviously with the increasing temperature, as shown in Figure 6(b). The reason is that the limiting shear stress decreases nonlinearly as the temperature rises [24]. According to [11], the friction is proportional to force between two surfaces in contact. The theory can be extended to a robot joint. Obviously, with the increase of gravity load, the value of friction increases, as shown in Figure 6(c).

The goodness-of-fit and errors between actual friction and model prediction for three kinds of models from equation (2) to equation (4) are presented in Table 1.

The normalized value of errors are presented in Figure 7, which indicates that both temperature and load have large influence on friction. Many researchers have only investigated two of these factors [20, 25, 26]. Therefore, it is of great significance to investigate all three factors.

4.2. Prediction Effects of BAS-RBFNN, RBFNN, and Velocity-Load-Temperature Model

To compare the prediction effects of three models, we take all the data as the train set and compare the prediction effects with the analytical model in equation (4), the related parameters of which are listed in Table 2.

Figure 8 shows that the mean error predicted using BAS-RBFNN is about 0.0011 Nm, which is lower than 0.0029 using RBFNN and 0.0058 Nm using the analytical model. Details are shown in Table 3. Compared with RBFNN, BAS-RBFNN has advantage in data fitting, especially for a higher accuracy prediction.

4.3. Comparison of Prediction Effects between GA-RBFNN and BAS-RBFNN

Genetic Algorithm (GA) is widely used for optimization. In this section, GA-RBFNN and BAS-RBFNN for friction modeling are compared. Table 4 shows that the accuracy of the two methods is nearly the same. In the training process, the mean squared error (MSE) is set to 0.0001. Figure 9 shows the training process of GA-RBFNN and BAS-RBFNN about the mean squared error (MSE). For GA-RBFNN, the target MSE is obtained at epoch 356, while for BAS-RBFNN, the target MSE is obtained at epoch 257. It can be seen that the proposed BAS-RBFNN possesses better performance in terms of faster convergence rate compared with GA-RBFNN.

4.4. Generalization Ability of BAS-RBFNN Model

In order to test the generalization ability of the BAS-RBFNN model, all the data are selected in unseen areas. The velocity, load, and temperature are set from 52°/s to 60°/s with 2° step, from −22° to −10° with 2° step, and from 45°C to 55°C with 2°C step, respectively. The error of static friction is plotted in Figure 10. It is shown that the maximum and average errors are 0.0063 Nm and 0.0012 Nm, which indicates that the proposed BAS-RBFNN has good generalization capability in unseen areas.

5. Conclusion and Future Work

In this paper, we propose a BAS-RBFNN method to describe the relationship between velocity, load, temperature, and friction. It is shown that the proposed BAS-RBFNN model has a better prediction precision and a higher convergence rate than RBFNN. In the future, some other influencing variables for static friction will be considered to improve the generalization capability. Furthermore, dynamic friction will be studied.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare no conflicts of interest.

Acknowledgments

This research was supported by the 2017 National Key R&D Program of China (no. 2017YFB1301400).