Next Article in Journal
Bringing the Social Back into Sustainability: Why Integrative Negotiation Matters
Next Article in Special Issue
An Evaluation Model for the Comfort of Vehicle Intelligent Cockpits Based on Passenger Experience
Previous Article in Journal
Perception of the Impacts of Tourism by the Administrations of Protected Areas and Sustainable Tourism (Un)Development in Slovakia
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Prediction of Clearance Vibration for Intelligent Vehicles Motion Control

1
School of Mechanical Engineering, Beijing Institute of Technology, Beijing 100081, China
2
School of Mechanical and Electrical Engineering, Changchun University of Science and Technology, Changchun 130013, China
3
The 55 Research Institute of China North Industries Group Corporation Limited, Changchun 130012, China
*
Author to whom correspondence should be addressed.
Sustainability 2022, 14(11), 6698; https://doi.org/10.3390/su14116698
Submission received: 20 April 2022 / Revised: 27 May 2022 / Accepted: 27 May 2022 / Published: 30 May 2022
(This article belongs to the Special Issue Smart Transportation and Intelligent and Connected Driving)

Abstract

:
Motion control analysis should consider the system’s uncertainty to ensure the intelligent vehicle’s autonomy. The clearance structure of the transmission shaft is modeled as a cantilever beam with double clearance to predict the clearance vibration for mitigating the nonlinearity. Based on the Kelvin–Voigt collision model, a clearance model was developed using time-varying parameters identified by the wavelet transform. Comparing the frequency response functions (FRF) of the initial model with constant parameters and the updated model with time-varying parameters, the experimental results from the updated model indicate that the modal assurance criterion (MAC) is increased by 42.92%, 31.08%, 38.97%, and 50.74% in the first-four order. Cross-signature assurance criteria (CSAC) and cross-signature scale factor (CSF) have been increased by 6.55% and 12.37%. The control method based on the clearance model has been verified. In the case of 120 km/h, compared with model-predictive control (MPC) and sliding mode control (SMC), the peak of the lateral position error was reduced by 35.7% and 14.3%, and the peak of the heading error was reduced by 50% and 15.6%.

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, l 1 , l 2 , l f represent the distance between the fixed end O to the clearance 1, clearance 2, and the excitation. E , A , I , ρ , l are Young’s modulus, cross-sectional area, moment of inertia, density, and length of the cantilever beam. The clearance values: d = d 1 = d 2 = d l d b / 2 . d b and d l are the outer diameter of the beam and the inner diameter of the bearing bush. k d 1 , k d 2 and c d 1 , c d 2 are stiffness and damping at clearance 1 and 2. x L , t is the vertical displacement,
x L , t = 1 n Y i L p i t
Y i x , p i t , i are modal shape function, coordinates, and the mode order. The boundary conditions of the cantilever beam are
2 x L , t L 2 L = 0 = 0 ;   3 x L , t L 3 L = 0 = 0 ;   2 x L , t L 2 L = l = 0 ;   3 x L , t L 3 L = l = 0
Based on the Hamilton principle and assumed modal method, a dynamic model of the cantilever beam with double clearance is established:
M p i p ¨ i + C p i p ˙ i + K p i p i + Y i l 1 f 1 t + Y i l 2 f 2 t + V i l 3 ϕ B W t = Y i l f f t
M p i , K p i , C p i represent the i-order modal mass, stiffness, and damping. f d 1 t , f d 2 t , f t represent the nonlinear forces at positions 1, 2, l f . The following formula can calculate each parameter:
M p i = 0 l ρ A Y i 2 d x
K p i = 0 l E I Y i 2 d x = 0 l Y i E I Y i d x
C p i = η 1 M p i + η 2 K p i
Y i x = cos β i x c h β i x + sin λ i s h λ i cos λ i c h λ i sin β i x s h β i x
η 1 and η 2 represent Rayleigh damping parameters. β i is the order of frequency equation. The characteristic roots λ i of the frequency equation satisfy Equation (8), and λ i = 2 i + 1 π / 2 .
cos λ i c h λ i = 1
The hysteresis model is composed of a linear spring-damper k B W c b w 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:
ϕ B W ( x , x ˙ ) = λ k x + ( 1 λ ) k x B W x ˙ B W = A x ˙ ζ x ˙ x B W n B W 1 x B W γ x ˙ x B W n B W
The restoring force is constituted by linear force λ k x and hysteresis force ( 1 λ ) k x B W . λ determines the weight of the nonlinear and linear. A , ζ , γ control the shape of the hysteresis loop. The order is usually 2.
Based on the nonlinear clearance force f d t 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:
f d t = c d t g ˙ t + k d g t
g t are the displacement and velocity functions at the clearance position. The expression is as follows:
g t = x l , t d x l , t d 0 d < x l , t < d x l , t + d x l , t d
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 g t at a certain point of the cantilever beam with double clearance under excitation. The wavelet transform is W a , b . The ridges of the wavelet transform a , b 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: x t = e ξ t ω t cos ω t + φ , the instantaneous damping ratio is
ξ = d I n A t ω d t
In combination with the wavelet center frequency ω 0 , scale factor a , translation factor b , wavelet coefficient W a , b , and instantaneous amplitude A t , the instantaneous frequency ω is represented as follows:
ω = ω 0 a ;   A t = 2 W a , b a
Assuming that the mass is constant and known, the time-varying stiffness obtained by the identification system is
k ( t ) = m ω 2
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
c ( t ) = α 0 m + α 1 k ( t ) + α 3 j k ( t )
m and k ( t ) 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
α 0 = 2 ω ξ α 1 = 2 ξ / ω α 3 = 2 ξ
As a consequence, the nonlinear damping model is expressed as follows:
c ( t ) = 4 m ω ξ t + 2 ξ j m ω 2
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: c ( t ) = 0.02 , t < 2 0.04 , t 2 , k ( t ) = 8000 , t < 2 5000 , t 2 .
Working condition 2: Continuous function value: c ( t ) = 0.02 0.02 t , k ( t ) = 8000 1000 t .
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: c ( t ) = 0.06 0.0002 t , 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 2 5 , 2 8 has three discrete scale sequences. Translation parameters: b = a 16 k + x 8 ; variable parameters: k = 0 : 0.005 : T / 4 a .
  • a = 2 5 , 2 6 , 2 7 , 2 8 ;
  • a = 2 6 , 2 6.25 , 2 6.5 , 2 6.75 , 2 7 ;
  • a = 2 6.5 , 2 6.6 , 2 6.7 , 2 6.8 , 2 6.9 , 2 7 .
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 a = 2 5 . The discrimination accuracy is improved obviously with a = 2 6 and a = 2 7 . However, the identification results are dramatically reduced with a = 2 8 . Therefore, the analysis scale should be selected between a = 2 6 and a = 2 7 . In sequence 2, the accuracy of identification results is better than other scales with a = 2 6.5 , a = 2 6.5 , a = 2 7 . 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 a = 2 6.8 , 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.

5. Conclusions

Vibration reduces the transmission shaft control accuracy and the response time resulting from the clearance. The intelligent shaft is abstracted to a cantilever beam with double clearance to forecast vibrations for error compensation by introducing time-varying parameters. The results show the following:
  • With the Hamilton principle and the assumed modal method, a dynamic model of the cantilever beam with double clearance is constructed. A clearance modeling method is proposed to describe clearance dynamics with time-varying parameters. A wavelet transform was used to identify the time-varying stiffness of clearance by estimating the natural frequency and damping ratio. The time-varying damping of clearance was constructed using Rayleigh and complex damping. The clearance model accurately describes the dynamical clearance variation using time-varying damping and stiffness.
  • An experimental model was constructed to verify the proposed modeling method. Compared to the constant model, the maximum natural frequency error of the time-varying model was reduced from 5.38% to 1.29%. The maximum MAC also increased from 0.4173 to 0.8492. The MRE and RMSE of each measurement point decreased by 27.67% and 63.41%. The FRAC increased from 0.7971 to 0.9465. CSAC and CSF increased by 6.55% and 12.37%, which is over 90%, particularly near the natural frequency position. Therefore, the above results demonstrate the validity and accuracy of the modeling method.
  • The control method based on the clearance model was verified. At 72 km/h, compared with SMC, the tracking error peak of SMC-CM decreased by 21.2% and the reliability increased by 12.7%. Compared to MPC and MSC, the lateral position error peak of SMC-CM decreased by 50% and 14.3%, the heading error peak decreased by 35.7% and 15.6%, and the yaw velocity peak decreased by 23.5% and 3.7%. Peak lateral accelerations decreased by 27.8% and 18.8%. The effectiveness of the control method is verified.
This paper aims to examine the influence of mechanical transmission system errors on control system accuracy. In the future, tracking effects and handling reliability will be improved in extreme conditions. The influence of clearance factors in other vehicle mechanical systems on handling reliability and safety will also be considered, along with optimizing the control system.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/su14116698/s1, The data: The data of FRF.

Author Contributions

Writing—original draft, Y.Z.; Writing—review & editing, F.Z., W.W., F.M., D.Z. and H.W. All authors have read and agreed to the published version of the manuscript.

Funding

This research is supported by the State Department project of China (Grant No. JCKY2017208A001, MKF20200022) and the National Science Fund of China (Grant No. 51805502). The authors gratefully acknowledge the facilities provided by the Industrial and Intelligent System Engineering Laboratory (IISEL) at the Beijing Institute of Technology.

Data Availability Statement

The data is in the Supplementary Materials.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Wang, W.; Bengler, K.; Jiang, X. Green Intelligent Transportation Systems. In Proceedings of the 7th International Conference on Green Intelligent Transportation System and Safety, Nanjing, China, 1–4 July 2016; Wang, W., Bengler, K., Jiang, X., Eds.; Lecture Notes in Electrical Engineering. Springer: Singapore, 2018; Volume 419. [Google Scholar]
  2. Fu, T.; Yao, C.; Long, M.; Gu, M.; Liu, Z. Overview of Longitudinal and Lateral Control for Intelligent Vehicle Path Tracking. In Proceedings of the Chinese Intelligent Automation Conference, Jiangsu, China, 20–22 September 2019. [Google Scholar]
  3. Colbaugh, R. Adaptive Stabilization of Mobile Manipulators. J. Robot. Syst. 1998, 15, 511–523. [Google Scholar] [CrossRef]
  4. Huang, L. Speed Control of Differentially Driven Wheeled Mobile Robots—Model-Based Adaptive Approach. J. Robot. Syst. 2005, 22, 323–332. [Google Scholar] [CrossRef]
  5. Pu, S.T.; Li, S.W.; Fan, R.C.; Ter-Feng, W. Point Stabilization Control of a Car-like Mobil Robot in Hierarchical Skew Symmetry Chained Form. In Proceedings of the IEEE International Conference on Networking, Sensing and Control, Taipei, Taiwan, 21–23 March 2004; IEEE: Taipei, Taiwan, 2004; Volume 2, pp. 1346–1351. [Google Scholar]
  6. Fierro, R.; Lewis, F.L. Control of a Nonholomic Mobile Robot: Backstepping Kinematics into Dynamics. J. Robot. Syst. 1997, 14, 149–163. [Google Scholar] [CrossRef]
  7. Dixon, W.E.; Dawson, D.M.; Zhang, F.; Zergeroglu, E. Global Exponential Tracking Control of a Mobile Robot System via a PE Condition. IEEE Trans. Syst. Man Cybern. Part B 2000, 30, 129–142. [Google Scholar] [CrossRef]
  8. Tayebi, A.; Rachid, A. Adaptive Controller for Non-Holonomic Mobile Robots with Matched Uncertainties. Adv. Robot. 2000, 14, 105–118. [Google Scholar] [CrossRef] [Green Version]
  9. Lihan, S.; Jie, M.; Baoqing, Y. Fuzzy PID Design of Vehicle Attitude Control Systems. In Proceedings of the 2020 Chinese Control and Decision Conference (CCDC), Hefei, China, 22–24 August 2020; pp. 1826–1830. [Google Scholar]
  10. Majid, M.; Arshad, M. A Fuzzy Self-Adaptive PID Tracking Control of Autonomous Surface Vehicle. In Proceedings of the 2015 IEEE International Conference on Control System, Computing and Engineering (ICCSCE), Batu Ferringhi, Penang, Malaysia, 27–29 November 2015; pp. 458–463. [Google Scholar]
  11. Shuguang, W.; Xiaojun, M.; Qinghan, Z.; Chunguang, L. Adaptive Sliding Mode Wheel Slip Control of Electric Vehicles Based on Feedback Linearization. Fire Control Command Control 2016, 41, 23–27. [Google Scholar]
  12. Temiz, O.; Cakmakci, M.; Yildiz, Y. A Fault Tolerant Vehicle Stability Control Using Adaptive Control Allocation. In Proceedings of the Dynamic Systems and Control Conference, Atlanta, GA, USA, 30 September–3 October 2018; Volume 51890, pp. 1–9. [Google Scholar]
  13. Lucibello, P.; Oriolo, G. Robust Stabilization via Iterative State Steering with an Application to Chained-Form Systems. Automatica 2001, 37, 71–79. [Google Scholar] [CrossRef]
  14. Yongoug, C.; Chongkug, P.; Harashima, F. A Position Control Differential Drive Wheeled Mobile Robot. IEEE Trans. Ind. Electron. 2001, 48, 853–863. [Google Scholar] [CrossRef]
  15. Hamel, T.; Meizel, D. On Robustness and Precision of Mobile Robots Missions. Automatica 2001, 37, 437–444. [Google Scholar] [CrossRef]
  16. Lee, T.C. Practical Stabilization for Nonholonomic Chained Systems with Fast Convergence, Pole-Placement and Robustness. In Proceedings of the 2002 IEEE International Conference on Robotics and Automation, Washington, DC, USA, 11–15 May 2002; Volume 4, pp. 3534–3539. [Google Scholar]
  17. Ying, M.J. Robust Control with Decoupling Performance for Steering and Traction of 4WS Vehicles under Velocity-Varying Motion. IEEE Trans. Contr. Syst. Technol. 2000, 8, 554–569. [Google Scholar]
  18. Xu, X.; Mi, J.; Wen, C.Y.; Wang, F.; Ma, S.D.; Tao, T. Robust H∞ Optimal Control for the Lateral Stability of Semitrailer Trucks. Mod. Manuf. Eng. 2019, 45, 1283–1293. [Google Scholar]
  19. Li, P.; Li, P.; Zhao, J.; Zhang, B. Robust Gain-Scheduling Static Output-Feedback H∞ Control of Vehicle Lateral Stability with Heuristic Approach. Inf. Sci. 2021, 546, 220–233. [Google Scholar] [CrossRef]
  20. Chen, T.; Chen, D. Lateral Control of Intelligent Vehicle Based on Neural Networks Sliding Mode. Transducer Microsyst. Technol. 2017, 36, 63–67. [Google Scholar]
  21. Chang, Z.; Hao, L.; Yan, Q.; Ye, T. Research on Manipulator Tracking Control Algorithm Based on RBF Neural Network. J. Phys. Conf. Ser. 2021, 1802, 032072. [Google Scholar] [CrossRef]
  22. Jia, L.; Zhu, Z. Improved Fractional-Order Integral Sliding Mode Control for AUV Based on RBF Neural Network. In Proceedings of the 2019 Chinese Automation Congress (CAC), Hangzhou, China, 22–24 November 2019; pp. 4809–4814. [Google Scholar]
  23. Henaff, P.; Chocron, O. Adaptive Learning Control in Evolutionary Design of Mobile Robots. In Proceedings of the IEEE International Conference on Systems, Man and Cybernetics, Yasmine Hammamet, Tunisia, 6–9 October 2022; Volume 3, p. 5. [Google Scholar]
  24. Kim, S.; Kim, T.I.; Jang, K.S.; Jung, S.; Kim, S.S. Control Experiment of a Wheeled Drive Mobile Pendulum Using Neural Network. In Proceedings of the 30th Annual Conference of IEEE Industrial Electronics Society, Busan, Korea, 2–6 November 2004; Volume 3, pp. 2234–2239. [Google Scholar]
  25. Rusu, P.; Petriu, E.M.; Whalen, T.E.; Cornell, A.; Spoelder, H.J. Behavior-Based Neuro-Fuzzy Controller for Mobile Robot Navigation. IEEE Trans. Instrum. Meas. 2003, 52, 1335–1340. [Google Scholar] [CrossRef]
  26. Wang, Z.; Ge, S.; Lee, T. Adaptive Neural Network Control of a Wheeled Mobile Robot Violating the Pure Nonholonomic Constraint. In Proceedings of the 2004 43rd IEEE Conference on Decision and Control, Atlantis, Paradise Island, Bahamas, 14–17 December 2004; Volume 5, pp. 5198–5203. [Google Scholar]
  27. Fierro, R.; Lewis, F.L. Control of a Nonholonomic Mobile Robot Using Neural Networks. IEEE Trans. Neural Netw. 1998, 9, 589–600. [Google Scholar] [CrossRef] [Green Version]
  28. Choomuang, R.; Afzulpurkar, N. Hybrid Kalman Filter/Fuzzy Logic Based Position Control of Autonomous Mobile Robot. Int. J. Adv. Robot. Syst. 2005, 2, 20. [Google Scholar] [CrossRef] [Green Version]
  29. Maalouf, E.; Saad, M.; Saliah, H. A Higher Level Path Tracking Controller for a Four-Wheel Differentially Steered Mobile Robot. Robot. Auton. Syst. 2006, 54, 23–33. [Google Scholar] [CrossRef]
  30. Butcher, M.; Karimi, A.; Longchamp, R. A Statistical Analysis of Certain Iterative Learning Control Algorithms. Int. J. Control 2008, 81, 156–166. [Google Scholar] [CrossRef]
  31. Utkin, V.; Poznyak, A.; Orlov, Y.; Polyakov, A. Conventional and High Order Sliding Mode Control. J. Frankl. Inst. 2020, 357, 10244–10261. [Google Scholar] [CrossRef]
  32. Tian, Q.; Xiao, Q.; Sun, Y.; Hu, H.; Liu, H.; Flores, P. Coupling Dynamics of a Geared Multibody System Supported by ElastoHydroDynamic Lubricated Cylindrical Joints. Multibody Syst. Dyn. 2015, 33, 259–284. [Google Scholar] [CrossRef] [Green Version]
  33. Liu, J.; Li, B. Theoretical and Experimental Identification of Clearance Nonlinearities for a Continuum Structure. J. Comput. Nonlinear Dyn. 2016, 11, 041019. [Google Scholar] [CrossRef]
  34. Yang, G.; Yang, J.; Qiang, C.; Ge, J.; Chen, Q. Natural Frequencies of a Cantilever Beam and Block System with Clearance While Block Staying on given Position. J. Vib. Control 2013, 19, 262–275. [Google Scholar] [CrossRef]
  35. Bai, Z.F.; Zhao, Y. Dynamic Behaviour Analysis of Planar Mechanical Systems with Clearance in Revolute Joints Using a New Hybrid Contact Force Model. Int. J. Mech. Sci. 2012, 54, 190–205. [Google Scholar] [CrossRef]
  36. Khatiwada, S.; Chouw, N.; Butterworth, J.W. A Generic Structural Pounding Model Using Numerically Exact Displacement Proportional Damping. Eng. Struct. 2014, 62–63, 33–41. [Google Scholar] [CrossRef]
  37. Marques, F.; Isaac, F.; Dourado, N.; Flores, P. An Enhanced Formulation to Model Spatial Revolute Joints with Radial and Axial Clearances. Mech. Mach. Theory 2017, 116, 123–144. [Google Scholar] [CrossRef]
  38. Liu, J.; Li, B.; Jin, W.; Han, L.; Quan, S. Experiments on Clearance Identification in Cantilever Beams Reduced from Artillery Mechanism. J. Mech. Eng. Sci. 2017, 231, 1010–1032. [Google Scholar] [CrossRef]
  39. Mayergoyz, I.D. Mathematical Models of Hysteresis and Their Applications; Academic Press: New York, NY, USA, 2003. [Google Scholar]
  40. Goldsmith, W.; Frasier, J.T. Impact: The Theory and Physical Behavior of Colliding Solids. J. Appl. Mech. 1961, 28, 639. [Google Scholar] [CrossRef]
  41. Xu, C.; Wang, C.; Gao, J. Modal Identification of Linear Time-Varying Systems Using Continuous Wavelet Transform. J. Harbin Inst. Technol. 2015, 22, 30–36. [Google Scholar]
  42. Hu, G. Modern Signal Processing Course; Tsinghua University Press: Beijing, China, 2004. [Google Scholar]
  43. Rayleigh, J.; Lindsay, R.B. The Theory of Sound; Dover Publications: New York, NY, USA, 1945; Volume 2. [Google Scholar]
  44. Myklestad, N.O. The Concept of Complex Damping. J. Appl. Mech. 1952, 19, 284–286. [Google Scholar] [CrossRef]
  45. Dong, X.Y.; Peng, Z.R.; Yin, H.; Dong, H.T. Distance Coefficient-Fisher Information Criterion for Optimal Sensor Placement. CAAI Trans. Intell. Syst. 2017, 12, 32–37. [Google Scholar]
  46. Wen, H.B.; Wang, G.Z. Updating Finite Element Model of Cabin with Floating Raft through Sensitivity Analysis of FEF Correlation Coefficient. J. Jiangsu Univ. Sci. Technol. 2005, 19, 75–78. [Google Scholar]
Figure 1. Simplified model of cantilever beam with double clearance.
Figure 1. Simplified model of cantilever beam with double clearance.
Sustainability 14 06698 g001
Figure 2. Identification results: (a) is step damping identification results; (b) is step stiffness identification results.
Figure 2. Identification results: (a) is step damping identification results; (b) is step stiffness identification results.
Sustainability 14 06698 g002
Figure 3. Identification results: (a) is continuous damping identification results; (b) is continuous stiffness identification results.
Figure 3. Identification results: (a) is continuous damping identification results; (b) is continuous stiffness identification results.
Sustainability 14 06698 g003
Figure 4. Parameter identification based on different sequences: (ac) represent the comparison of identification results at scales 1, 2, and 3, respectively.
Figure 4. Parameter identification based on different sequences: (ac) represent the comparison of identification results at scales 1, 2, and 3, respectively.
Sustainability 14 06698 g004
Figure 5. Structure model.
Figure 5. Structure model.
Sustainability 14 06698 g005
Figure 6. Test model.
Figure 6. Test model.
Sustainability 14 06698 g006
Figure 7. FRF of clearance position 1 and 2: (a) is position 1 for sensor 6; (b) is position 2 for sensor 7.
Figure 7. FRF of clearance position 1 and 2: (a) is position 1 for sensor 6; (b) is position 2 for sensor 7.
Sustainability 14 06698 g007
Figure 8. FRF of other position: (af) represent the comparison of FRF at sensor measuring points 1, 2, 3, 4, 5, 8.
Figure 8. FRF of other position: (af) represent the comparison of FRF at sensor measuring points 1, 2, 3, 4, 5, 8.
Sustainability 14 06698 g008
Figure 9. Comparison of CSAC: (ah) represent the comparison of CSAC at sensor measuring points 1, 2, 3, 4, 5, 6, 7, 8.
Figure 9. Comparison of CSAC: (ah) represent the comparison of CSAC at sensor measuring points 1, 2, 3, 4, 5, 6, 7, 8.
Sustainability 14 06698 g009aSustainability 14 06698 g009b
Figure 10. Comparison of CSF: (ah) represent the comparison of CSF at sensor measuring points 1, 2, 3, 4, 5, 6, 7, 8.
Figure 10. Comparison of CSF: (ah) represent the comparison of CSF at sensor measuring points 1, 2, 3, 4, 5, 6, 7, 8.
Sustainability 14 06698 g010
Figure 11. Error results at 72 km/h speed: (a) is lateral position error; (b) is yaw velocity error.
Figure 11. Error results at 72 km/h speed: (a) is lateral position error; (b) is yaw velocity error.
Sustainability 14 06698 g011
Figure 12. Error results at high speed: (a) is lateral position error; (b) is yaw velocity error.
Figure 12. Error results at high speed: (a) is lateral position error; (b) is yaw velocity error.
Sustainability 14 06698 g012
Figure 13. Error results: (a) is lateral acceleration error; (b) is yaw velocity error.
Figure 13. Error results: (a) is lateral acceleration error; (b) is yaw velocity error.
Sustainability 14 06698 g013
Table 1. Comparison of natural frequency errors.
Table 1. Comparison of natural frequency errors.
OrderFrequency (Hz)Initial Model Based on Constant DampingUpdated Model Based on Time-Varying Damping
Frequency (Hz)RE (%)MACFrequency (Hz)RE (%)MAC
19.259.151.0811 0.5378 9.20.5405 0.9422
239.7538.353.5220 0.6283 40.251.2579 0.9117
392.8597.855.3850 0.5265 91.651.2924 0.8627
4188.85191.351.3238 0.4183 187.250.8472 0.8492
Table 2. FRF errors in measuring points.
Table 2. FRF errors in measuring points.
PositionInitial Model Based on Constant DampingUpdated Model Based on Time-Varying Damping
MRERMSEFRACCSACCSFMRERMSEFRACCSACCSF
10.8013 1.4181 0.8643 0.88910.84100.5763 0.7481 0.9443 0.96330.9619
20.8078 1.3410 0.7708 0.89210.84850.5827 0.5902 0.9507 0.95310.9608
30.7963 1.5589 0.8593 0.88470.83630.5785 0.6026 0.9465 0.96770.9656
40.8050 1.3795 0.8680 0.88120.84590.5810 0.5946 0.9490 0.94870.9652
50.7970 1.6391 0.6600 0.89730.84310.5799 0.6574 0.9479 0.93750.9599
60.7938 1.9037 0.8568 0.89560.83920.5777 0.0976 0.9457 0.95670.9568
70.8007 1.0062 0.6373 0.88700.83710.5778 0.2419 0.9458 0.95320.9521
80.7969 1.3674 0.8599 0.89120.84570.5743 0.7174 0.9423 0.93720.9657
Means0.7999 1.4517 0.7971 0.8898 0.8421 0.5785 0.5312 0.9465 0.9522 0.9610
Table 3. Comparison of different control methods.
Table 3. Comparison of different control methods.
ParameterSMC-CMSMC
Lateral position error (m)0.0250.033
Yaw velocity (deg/s)2.833.24
Table 4. Error comparison of different control methods.
Table 4. Error comparison of different control methods.
ParameterSMC-CMSMCMPC
Lateral position error (m)0.060.070.12
Yaw velocity (deg/s)0.270.320.42
Table 5. Peak reliability comparison of different control methods.
Table 5. Peak reliability comparison of different control methods.
ParameterSMC-CMSMCMPC
Lateral acceleration (deg/s2) −5.2−5.4−6.8
Yaw velocity (deg/s) −2.6−3.2−3.6
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Zhang, Y.; Zhang, F.; Wang, W.; Meng, F.; Zhang, D.; Wang, H. Prediction of Clearance Vibration for Intelligent Vehicles Motion Control. Sustainability 2022, 14, 6698. https://doi.org/10.3390/su14116698

AMA Style

Zhang Y, Zhang F, Wang W, Meng F, Zhang D, Wang H. Prediction of Clearance Vibration for Intelligent Vehicles Motion Control. Sustainability. 2022; 14(11):6698. https://doi.org/10.3390/su14116698

Chicago/Turabian Style

Zhang, Yunhe, Faping Zhang, Wuhong Wang, Fanjun Meng, Dashun Zhang, and Haixun Wang. 2022. "Prediction of Clearance Vibration for Intelligent Vehicles Motion Control" Sustainability 14, no. 11: 6698. https://doi.org/10.3390/su14116698

APA Style

Zhang, Y., Zhang, F., Wang, W., Meng, F., Zhang, D., & Wang, H. (2022). Prediction of Clearance Vibration for Intelligent Vehicles Motion Control. Sustainability, 14(11), 6698. https://doi.org/10.3390/su14116698

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