1. Introduction
Intelligent vehicles provide an alternative solution to the growing problems of pollution, traffic accidents, and energy shortages [
1]. Motion control is key to achieving autonomous movement. Based on the information provided by sensors, the intelligent vehicle is controlled in real time to follow the desired trajectory [
2]. The control methods include adaptive control, robust control, and intelligent control.
For adaptive control, Colbaugh et al. combined adaptive control and homogeneous system theory to address the motion control for uncertain mobile manipulators [
3]. Huang et al. were able to control the motion of wheeled mobile robots by adaptive control. Compared to a model-free PID controller, the tracking effect of track was faster and more accurate [
4]. Tasi designed a new adaptive controller to provide progressive stabilization of intelligent parking spaces [
5]. Fukao, Dixon, and Tayebi et al. developed a controller for tracing uncertain mobile robots using the backstepping method [
6,
7,
8]. Sun provided the vehicle attitude control system with a fuzzy PID. Based on the fuzzy relationship between three PID control parameters and deviation signals, the online PID control parameters were adjusted to reduce the deviation of the vehicle attitude control system. As a result, the random interference and uncertainty of vehicle motion were reduced [
9]. For the ASV and the sway, surge, and yaw directions, Maajid used three parallel fuzzy adaptive PID controllers [
10]. In order to overcome the high uncertainties in acceleration and braking, Wei proposed an adaptive sliding mode control method based on feedback linearization [
11]. To ensure that the vehicle follows the trajectory created by the driver, Temiz proposed a parameter integration fault-tolerant vehicle monitoring algorithm based on online updates [
12]. In the process, adaptive control is suitable for solving nonlinear, time-varying, and hysteretic problems due to its high robustness. However, adaptive control fails to provide a general solution. Control accuracy and parameter estimation are inconsistent in some controllers. They do not meet the requirements for dynamic performance during the startup and transition process.
For robust control, Lucibello et al. designed robust stabilization controllers for nonlinear chains to ensure robustness [
13]. In order to achieve precise path tracking of a two-degree-of-freedom mobile robot (WMR), Chung et al. developed a double loop (speed loop and position loop) path tracking controller based on robust control methods [
14]. In order to track wheeled robots with bounded errors in the position feedback information, Hamel et al. developed a robust tracking controller [
15]. Lee used the switching controller method to stabilize nonholonomic chain systems, ensuring the system’s robustness and fast convergence of errors [
16]. Jia investigated the robust control of four-wheel-drive vehicles with uncertain velocity, moment of inertia, mass, and friction coefficient. Based on the established vehicle dynamics model, a robust control scheme was developed to limit the influence of changes in vehicle parameters on steering performance [
17]. To improve the lateral stability of tractors and semitrailers, Liu proposed a robust H∞optimal control scheme [
18]. Li proposed a method to improve lateral vehicle stability using output feedback control [
19]. Robust control systems are particularly effective for ensuring stability and reliability as the primary objectives. The dynamic process characteristics can be predicted without establishing a formal process model. Therefore, robust control is particularly suitable for systems with large uncertainty ranges and small stability margins. However, due to the robust control system generally failing to operate in the optimal state, steady-state accuracy is relatively low.
For intelligent control, Chen proposed the RBF neural network sliding mode control method to improve control accuracy and robustness due to many external disturbances during driving [
20]. Chang proposed controlling trajectories to enhance efficiency and accuracy by a neural network [
21]. Jia developed a robust adaptive control algorithm using RBF neural networks and improved the fractional integral sliding mode control. In complex underwater conditions, accuracy and stability were enhanced in the underwater docking process for autonomous underwater vehicles (AUVs) [
22]. Henaff et al. developed an intelligent vehicle controller based on neural networks [
23]. In a dynamic model of the inverted pendulum system, Kim compensated for parameter uncertainty errors using a neural network control strategy [
24]. Rusu and Wang et al. combined neural network control and backstepping to create a tracking controller for mobile robots [
25,
26]. By combining neural network and backstepping control methods, Fierro et al. developed a tracking controller for the intelligent vehicle to ensure asymptotic stability [
27]. Choomuang developed an obstacle avoidance/tracking controller for the intelligent vehicle using the Kalman filter and fuzzy control method [
28]. Because intelligent vehicles operate in complex environments containing unpredictable disturbances, the traditional dynamic model based on precise motion control was no longer appropriate. Butcher and Eli developed a fuzzy control system to ensure accurate tracking and robustness [
29,
30]. Intelligent control enables systems to run efficiently, solve generalized problems, and offer fault tolerance.
Moreover, it is capable of compensating and judging. Intelligent control does not clearly distinguish between controller and object. Compared with the traditional automatic control, the intelligent control objects have substantial uncertainty. The intelligent control is also limited to the particular requirements and conditions of the algorithm.
Furthermore, a sliding mode control method is also an exhaustive control algorithm. Due to its ability to overcome the system uncertainty, sliding mode control exhibits strong robustness to interference and unmodeled dynamics in the nonlinear control systems [
31]. The clearance of the transmission shaft is introduced into the control method to improve the accuracy and stability further.
However, the system nonlinearity is usually ignored in motion control or merely mitigated through simple error compensation. The previous study focused on the uncertainty in motion control based on simplifying the dynamics model. However, the simplified model cannot reflect the actual situation when the speed is high, and the mass is large. Intelligent vehicles are complex systems with multiple inputs and outputs, strong coupling, and significant nonlinearities. The mass, inertia, and tire friction coefficient parameters have measurement noise and uncertainty.
Additionally, the actuator exhibits nonlinearity (clearance, friction, and dead zone). Due to the clearance nonlinearity, the control precision is limited by the qualitative analysis, not quantitative analysis. Thus, it is imperative to analyze the quantitative analysis of clearance nonlinearity to eliminate the influence and effectively track the trajectory.
For clearance analysis, Tian developed various dynamic models to describe the dynamic clearance characteristics, establishing a theoretical basis for the subsequent research on clearance dynamics [
32]. Using the clearance model abstracted into the spring-damping model, Liu constructed a dynamic cantilever beam model with double clearance. Dynamic information determines the clearance value to control the nonlinear clearance error within a reasonable range [
33]. Yang carried out an experimental and numerical analysis for the free vibration modes of the clearance cantilever beam, derived the mathematical expression of the simplified clearance stiffness, and verified the calculated frequency with the experimental data. It was found that the measured frequency was in good agreement with the predicted frequency [
34]. To study the dynamic behavior of planar mechanical systems with joint clearance, Bai and Zhao developed a hybrid contact force model in conjunction with Lankarani–Nikravesh and Winkler elastic foundation models [
35]. Khantiwada proposed a relationship between force and deformation based on the Hunt–Crossley collision model and accurately calculated damping constants. Comparisons of linear and nonlinear Hunt–Crossley models with other models were performed to verify the model validity [
36]. A rotary joint model with radial and axial clearances was developed by Marques [
37].
The desired trajectory was tracked by controlling the transmission mechanism under a control algorithm, taking the clearance factor into account. Control precision affects the safety and stability of intelligent vehicles. The vehicle’s speed is usually set to a constant value subjectively in intelligent vehicle lateral control research on narrow curved roads. The control strategy can track the desired trajectory accurately. However, the transmission system accuracy may influence the vehicle tracking control with the preview system for roads with large curves. With high speed and large curvatures, a small transmission error can adversely affect the performance and stability of the intelligent vehicle tracking control. Furthermore, the speed is set subjectively as the expected speed lacks the ability to adapt to the vehicle in real time. As a result, the vehicle’s state cannot be accurately adapted for velocity changes in real time, resulting in sideslip or even roll-over, which directly impacts the vehicle’s control effect and local stability. In order to improve the motion control performance of an intelligent vehicle, the nonlinearity influence is estimated by predicting the clearance vibration. Stability and rapidity are improved to meet the driving requirements. This paper uses the time-varying parameters to describe clearance dynamics in contact moments.
The paper is organized into three sections. First, a dynamic clearance model is presented based on the Hamilton principle and the modal hypothesis. The hysteresis model is introduced to describe the lubrication impact on clearance further. Second, the wavelet transform determines the time-varying stiffness and damping. Finally, the effectiveness of the modeling method is verified by establishing an experimental model and comparing the relevant error indicators of the FRF. The clearance model is introduced into the sliding mode control, and a control method considering the gap nonlinearity is proposed, verified by simulation. Finally, the effectiveness of the modeling method is examined by establishing an experimental model and comparing the relevant error indicators. A control method based on clearance nonlinearity is proposed and validated by simulation.
2. Nonlinear Modeling for Clearance
The intelligent vehicle actuator consists of a connecting rod transmission mechanism, and the output shaft is connected to the front wheel to control the wheel deflection. A typical nonlinearity in the transfer system is the shaft clearance. In the transmission link of the servo system, the nonlinear position error is caused by the clearance. As a result, the dynamic performance and steady-state accuracy are negatively affected. In contact collision between the axis and bearing for a long time, clearance value will also increase, aggravating the nonlinear vibration of the transmission shaft. At the same time, the change in clearance size also impacts the control precision.
As the demand for high-precision control systems increases, more research is being conducted on clearance nonlinearity. Because of clearance in multiple systems, damping and stiffness exhibit different nonlinearities. Transmission system clearances directly impact the error compensation strategy and the control precision.
Therefore, the paper analyzes vibration generated by driveshaft clearance to improve the vibration prediction accuracy and ensure precise error compensation. The shaft is abstracted as a cantilever beam based on the dynamic characteristics. As illustrated in
Figure 1, based on the clearance model [
38], the bearing bush lengths and the grease hysteresis in clearance are introduced to construct a cantilever beam with double clearance. The left end of the beam is fixed, while the right end can be moved.
In
Figure 1,
,
,
represent the distance between the fixed end
to the clearance 1, clearance 2, and the excitation.
,
,
,
,
are Young’s modulus, cross-sectional area, moment of inertia, density, and length of the cantilever beam. The clearance values:
.
and
are the outer diameter of the beam and the inner diameter of the bearing bush.
,
and
,
are stiffness and damping at clearance 1 and 2.
is the vertical displacement,
,
,
are modal shape function, coordinates, and the mode order. The boundary conditions of the cantilever beam are
Based on the Hamilton principle and assumed modal method, a dynamic model of the cantilever beam with double clearance is established:
,
,
represent the i-order modal mass, stiffness, and damping.
,
,
represent the nonlinear forces at positions 1, 2,
. The following formula can calculate each parameter:
and
represent Rayleigh damping parameters.
is the order of frequency equation. The characteristic roots
of the frequency equation satisfy Equation (8), and
.
The hysteresis model is composed of a linear spring-damper
and a nonlinear B-W element. The B-W element provides the hysteresis force by the Coulomb friction block. The mechanical model has two degrees of freedom: particle displacement and hysteresis displacement [
39]. The dynamic equation is expressed as follows:
The restoring force is constituted by linear force and hysteresis force . determines the weight of the nonlinear and linear. , , control the shape of the hysteresis loop. The order is usually 2.
Based on the nonlinear clearance force
and the damping loss in the collision process, the Kelvin–Voigt model (a set of parallel linear spring-damping units) is adopted to describe the collision characteristics [
40], expressed as follows:
are the displacement and velocity functions at the clearance position. The expression is as follows:
A dynamic model of the double-clearance cantilever beam was developed to describe the dynamic characteristics of the drive shaft. The grease facilitated a smooth movement during collisions and recovery at the clearance. The grease forms change with collision, resulting in a complex interactive relationship. However, the constant damping in the Kelvin–Voigt model cannot accurately replicate the actual vibration laws. In order to overcome the limitations of the Kelvin–Voigt model, a spring-damping model using time-varying parameters is developed to model the dynamic process. Energy loss is nonlinear and time-varying at the re-contact moment after passing the contact point.
The time-varying parameters of the collision model are identified using wavelet transform theory. The validity of the identification algorithm is examined in the following section.
3. Parameter Identification Based on Wavelet Transform
This section aims to establish a time-varying parameters clearance model that can accurately predict the vibration caused by the clearance. The instantaneous damping ratio and frequency are calculated with wavelet ridges and skeletons. Due to the nonlinearity of clearance and lubrication, it is impossible to identify the contact force using the same scale accurately. Wavelet analysis can provide a variable time-frequency window to obtain more local information in some methods for parameter identification. Based on the analysis of high-frequency signals, the adaptive window size narrows automatically and grows wider for low-frequency signals [
41]. For non-stationary response signals of the output system, time-frequency analysis provides an effective method for identifying time-varying parameters.
3.1. Time-Varying Parameter of Clearance Structures
The test signal is assumed to be
at a certain point of the cantilever beam with double clearance under excitation. The wavelet transform is
. The ridges of the wavelet transform
represent a series of points. A wavelet ridge’s envelope coefficient is defined as the skeleton. A wavelet ridge and skeleton model can estimate the system’s natural frequency and damping ratio [
42].
For the response signal of the cantilever beam system:
, the instantaneous damping ratio is
In combination with the wavelet center frequency
, scale factor
, translation factor
, wavelet coefficient
, and instantaneous amplitude
, the instantaneous frequency
is represented as follows:
Assuming that the mass is constant and known, the time-varying stiffness obtained by the identification system is
As a result of the internal friction of solid materials and the effects of decoupling, the clearance vibration can be represented by the classical Rayleigh damping model and the complex damping model [
43,
44]. The dynamics model may be decoupled while maintaining the modal information. It can be written as
and
represent the equivalent mass and time-varying stiffness of the cantilever beam at the clearance position. Furthermore, the coefficients of the Rayleigh damping model and complex damping model are expressed as
As a consequence, the nonlinear damping model is expressed as follows:
The instantaneous frequency and damping ratio are estimated by analyzing the wavelet ridge and skeleton using the variable time-frequency window. Clearance models are built with time-varying stiffness and damping. Next, the wavelet transform algorithm is demonstrated numerically.
3.2. Parameter Identification and Verification
This section uses step and continuous function parameter models to verify identification accuracy. In order to analyze the performance of the identification algorithm against noise, SNR of 20 dB and 100 dB are applied. Other mechanical parameters did not change over time.
Working condition 1: Step function value: , .
Working condition 2: Continuous function value: , .
The wavelet transform is used to identify the damping and stiffness. In
Figure 2 and
Figure 3, damping and stiffness identification results for SNR = 20 and SNR = 100 are presented. For SNR = 100, the identified and theoretical values were kept within a low error range. Compared to SNR = 100, the relative error between the identified and theoretical values is greater with SNR = 20. In this case, it can be seen from the comparison that the identification effect of damping parameters is well within the stable range. However, in abrupt changes, the identification effect will fluctuate. As a whole, the identification result is tracked accurately, and it fluctuates around the theoretical value. Therefore, the method can accurately identify the theoretical value change at the mutation moment.
The effectiveness of the identification method has been demonstrated in numerous cases. Noise and scale factors affect the identification accuracy. According to Formula (13), the scale factor is a significant determinant of identification accuracy, influencing the relationship between accuracy and efficiency. The scale factor reduces as the translation factor increases. As part of maintaining an exact similarity, the scale factor is used to “extend” or “compress” the wavelet. It is necessary to select the scale parameter precisely to improve the time-varying parameter model.
In order to reduce the error and improve the identification accuracy, it is necessary to analyze the design parameters in the identification process. In the case of the continuous damping function: , multi-scale identification is conducted to compare and analyze the identification results over a range of scales. The selection and dispersion of scale factors in the recognition process is described in computational efficiency and parameter recognition accuracy.
Setting the range has three discrete scale sequences. Translation parameters: ; variable parameters: .
;
;
.
Based on the method proposed in this section, three different analytical scale sequences are used, as shown in
Figure 4. The solid line represents the theoretical value, and the dotted line represents the identification result at different scales.
In sequence 1, identification results were the worst with . The discrimination accuracy is improved obviously with and . However, the identification results are dramatically reduced with . Therefore, the analysis scale should be selected between and . In sequence 2, the accuracy of identification results is better than other scales with , , . However, as the scale factor value becomes smaller, the calculation amount increases. The number of corresponding translation parameters is 3200, 3800, 4400, 5200, and 6400 in five factors. In the smaller scale interval of sequence 3, the dispersion is more intensive. The identification result is the best at , but the amount of calculation and analysis increases significantly.
The smaller the analysis scale between the discrete scales of the three sequences, the more intensive the corresponding translation parameters. However, the more accurate the analysis, the more computational time and resources are required. In addition, estimating the scale factor affects identification accuracy and calculation efficiency. The wide intervals and long discrete methods should be selected based on the above analysis.
Based on the identified time-varying stiffness and damping, a time-varying parameter clearance model is derived to describe the vibration caused by clearance. The credibility of the identification algorithm is verified in the case of continuous and step functions. An experimental model based on time-varying parameters is used to determine the validity of the cantilever beam with double clearance.
4. Results and Discussion
In order to accurately predict the transmission shaft vibration caused by the clearance, the transmission shaft is abstracted as a cantilever beam with double clearance. A dynamic model of the cantilever beam double clearance has been developed using the wavelet transform to identify time-varying parameters in clearance. Next, the accuracy of vibration prediction is demonstrated by constructing a test platform.
4.1. Test Model
Figure 5 shows the structure model, while
Figure 6 shows the test model. Sensors 1–8 are installed on the beam. A Fisher information criterion is used to optimize the arrangement [
45]. Shakers applied random loads to the bottom plate, and sensors measured the dynamic response of each position. The signal is processed by the data acquisition system and transmitted by the power amplifier to the computer.
4.2. Validity of Clearance Model
FRF is calculated by testing the measured position. In
Figure 7, the FRF comparison is performed for the clearance positions (sensors 6 and 7). The blue, black, and red lines indicate the initial model based on constant values, the updated model based on parameter identification, and the test model. Initial and updated models agree well with experimental results at first-order and second-order frequencies. FRF error increases significantly between the initial and test models at third-order and fourth-order frequencies. By comparing the amplitude of FRF, the updated model is closer to the test model.
Figure 8 illustrates the comparison of FRF at the measurement points 1, 2, 3, 4, 5, and 8. First-order and second-order models may better reflect the test model with the initial and updated models. However, the difference between the initial and test models increases after the second order. In contrast, the updated model can reflect the changes in the FRF of the test model. Nevertheless, the updated model will contain some errors at the fourth-order position. The noise increased with identification times, which led to reduced accuracy. However, the updated model can better reflect the dynamic characteristics of the test model.
A correlation analysis is performed using relative error and MAC to examine the differences between the initial and updated models in natural frequencies, as shown in
Table 1. Compared to the original model, relative errors of the updated model are reduced by 0.54%, 2.26%, 4.09%, and 0.47%. The MAC grows by 42.923%, 31.083%, 38.97%, and 50.743%.
In addition to comparing errors at natural frequencies, the error at other frequencies may also reflect the correlation degree between the models. Next, the MRE, RMSE, FRAC, CSAC, and CSF are examined as evaluation indicators. CSF and CSAC indicate the FRF correlation degree between the numerical and test models at each frequency. This value ranges from 0 to 1, for which a value closer to 1 indicates greater accuracy [
46]. The comparison is presented in
Table 2.
Following the comparison, the relative error and RMSE are reduced by 27.67% and 63.41% for the updated model based on parameter identification. The FRAC has improved by 15.79%, while the value exceeds 94%. As shown in
Figure 9 and
Figure 10, CSAC and CSF have increased by 6.65% and 12.37%, over 90%. In the updated model, the CSAC and CSF have increased significantly at natural frequencies close to one. It is confirmed that the updated model can reflect the frequency domain characteristics of the test model.
In this section, experimental verification is carried out using the established experimental platform. In comparing the FRF of the initial and updated model, the model was evaluated by the relative error and MAC to prove the validity of the updated model. In addition, by comparing the MRE, RMSE, FRAC, CSAC, and CSF, the updated model would well reflect the dynamic characteristics of the test model.
Experimental verification is carried out in this section using the established experimental platform. Using the FRF of the initial and updated model, the model was evaluated by the relative error and MAC to determine the validity of the updated model. Furthermore, MRE, RMSE, FRAC, CSAC, and CSF are compared. The updated model would reflect the dynamic characteristics of the test model.
4.3. Validity of Motion Control
Based on the above analysis, a dynamic prediction model is formulated to consider clearance influence. Matlab/Simulink and Carsim are used to establish the sliding mode control model based on the clearance model (SMC-CM) for performance evaluation by introducing the dynamic prediction model into the sliding mode control system. To assess the feasibility of the SMC-CM, the sliding mode control (SMC) was compared to the SMC-CM at a speed of 72 km/h. Then, SMC-CM was compared with MPC and SMC to verify its superiority at 120 km/h.
As can be seen from the comparison results between SMC-CM and SMC at 72 km/h speed, both lateral position error and yaw velocity can be maintained within a tight range in
Figure 11.
Table 3 illustrates the lateral position error and yaw velocity peak pairs for different control methods. In comparison with SMC, the SMC-CM tracking error peak has decreased by approximately 21.2%, and reliability has improved by approximately 12.7%. The results indicate that SMC-CM is a viable solution under low-speed conditions.
In order to determine the control effect of the proposed control method in medium and high-speed driving, MPC, SMC, and SMC-CM results were compared to verify the superiority of SMC-CM at high speed.
Figure 12 illustrates the results of the lateral position error and heading angle error. The maximum errors of the three controllable methods all occur near the maximum curvature. These errors are shown in
Table 4. Compared with the MPC, the lateral position error peak and heading error peak of SMC-CM are reduced by 50% and 35.7%. In comparison with SMC, the lateral position error peak value and the heading error peak of SMC-CM are reduced by 14.3% and 15.6%. The accuracy of SMC-CM is superior to MPC and has superior improvement over SMC.
Figure 13 depicts a comparison of yaw velocity and lateral acceleration.
Table 5 illustrates the peak reliability of MPC, SMC, and SMC-CM. In comparison with MPC, the SMC-CM has a 23.5% increase in peak yaw velocity and a 27.8% increase in peak lateral acceleration. The peaks are reduced by 3.7% and 18.8% compared to SMC. The results indicate that the lateral reliability of SMC-CM is superior to MPC.
In simulation experiments, SMC-CM further improved the vehicle tracking effect compared to SMC and MPC. The tracking effect is optimized further for variable speeds and curvatures. On a real device, the radial displacement sensor is arranged vertically on the transmission shaft bench to measure the radial vibration. In other words, the two measuring points are 90 degrees apart. Sensors for X and Y directions are located on either side of the vertical centerline and are positioned as close as possible to the shaft. In addition to a wide range of interference factors in the measurement process, noise reduction of vibration data is also a major concern. In order to adjust the driving state of the vehicle, a motion control system predicts the vibration law of the drive shaft using the analyzed data.