Interval Recognition Algorithm of the Pavement Surface Condition Based on Lagrange Interpolation Method

.e accurate recognition of road condition is one of the important factors that influence vehicle safety performance. .is paper comes up with an original mathematical method of an interval recognition algorithm of the pavement surface condition based on Lagrange interpolation. .e ordinate of the peak point is solved by the Lagrange interpolation method, and the pavement surface condition is deduced by the interval identification algorithm. .e simulation results from six typical roads and the varied pavement surface show that besides the cobblestone pavement which is not common in the daily road, the estimation error of the initial tire-road friction coefficient by the Lagrange interpolation method is less than 2%, the pavement surface condition can be identified by interval recognition algorithm quickly and accurately, and the response time is less than 0.2 seconds.


Introduction
In order to reduce the occurrence of traffic accidents, major automobile companies and research institutes began to research and develop vehicle active safety control technology, such as ABS, TCS, EBD, ESP, and ACC. e essence of these active control strategies is to adjust the force acting on the tire and pavement surface so as to achieve the purpose of braking (driving) control and vehicle stability control, which are restricted by the adhesion conditions between the tire and the pavement surface [1,2]. For example, the maximum brake deceleration can be provided by a dry asphalt pavement with the adhesion coefficient being about 1, 20 times than that by an ice pavement with the adhesion coefficient being about 0.05.
is means that on different pavement surface conditions, the braking distance of the car may be 20 times different, which is crucial to the safe operation of the car.
Besides the quality of active safety control strategy, the significant automated driving depends to a great extent on whether the current road adhesion coefficient can be fully utilized.
e key to improving the utilization ratio of road adhesion coefficient is to obtain the maximum tire-road friction coefficient in real time. e tire-road friction coefficient μ is defined formally as follows [3]: where F z , F x , and F y are the normal forces, longitudinal forces, and side forces acting on the tire, respectively. For this wheel, maximum tire-road friction coefficient μ max is the maximum achievable value of |μ|. If we only consider the longitudinal motion of the vehicle, then the side force F y can be neglected:

Review of Previous Results
e road recognition method is traditionally examined with "cause-based" approaches and "effect-based" approaches according to the principles of the recognition method [4,5], as shown in Figure 1.
e "cause-based" strategies try to measure factors that lead to changes in friction and then attempt to predict μ max via some additional sensors. e aim of the "cause-based" approach is to find a correlation between the sensor data and tire-road friction-related parameters such as roughness and lubrication of the pavement surface. e most commonly used sensors are optical sensors, cameras [6,7], and acoustic sensors [8].
e main advantage of the "cause-based" method is that it can identify the pavement surface condition before the vehicle passes, which is of great significance to future pilotless. Besides, this method is not limited by the current road adhesion coefficient and can predict the maximum tireroad friction coefficient precisely under the condition that the friction between the tire roads is low.
However, this method exists two main disadvantages: first, the experiment-based approach considers only the road level and not the tire wear of the vehicle itself. Even on the road with good tire-road friction coefficients, tire with severe tread wear will not get good adhesion, which is very harmful to the control of vehicles after identifying the maximum tireroad friction coefficient. e second disadvantage is that the experiment-based approach usually requires additional sensors, which are usually expensive and vulnerable to external interference. When the sensor is polluted, it is harmful for the safe driving of the car. e "effect-based" approach, on the other hand, estimates the tire-road friction coefficient by testing the motion response generated by the change in the road adhesion coefficient on the vehicle body or tire.
It can be divided into tire-based response and vehicle dynamics response to identify the tire-road friction coefficient.
e tire-based response can be further divided into two methods of measuring the tire noise response [9] and measuring the tread deformation [10]. Both of these two methods do not get rid of the drawbacks of requiring additional sensors. e vehicle dynamics response approaches try to estimate the friction using mathematical models. e most commonly mathematical models are wheel dynamic model [11,12], roll dynamic model [13], quarter car model [14], four-wheel vehicle model [3], bicycle model [15,16], and tire model such as Pacejka tire model, Dugoff tire model, and Burckhardt tire model [17]. e vehicle dynamics response approaches do not require additional sensors and the accuracy of the recognition is within the allowable range.
But it still exists disadvantages that usually need to the tire close to locked state, which is harmful to the control of the vehicle, to estimate the maximum tire-road friction coefficient.
Besides the quality of active safety control strategy, proper maintenance of a road network needs to evaluate the pavement surface condition [18].
In order to avoid the problem of high cost and easy to be polluted of the sensor method, and the disadvantage of the model-based method that can identify the maximum tireroad friction coefficient only when the tire is locked in the critical state, an original mathematical method was proposed. e Lagrange interpolation method in mathematics for reference is used to determine the maximum friction coefficient, and an interval recognition algorithm is used to determine the final pavement surface condition. e organization of this paper is as follows. In Section 3, the system model is explained briefly, which discusses the wheel dynamics model and tire model. Section 4 expounds the overall approach to tire-road friction estimation. In this section, the acquisition of real-time road adhesion coefficient, the calculation of braking torque, and the application of the Lagrange interpolation method are introduced in detail. In Section 5, the simulation results of different roads and the varied pavement surface are given to verify the effectiveness of the algorithm. Finally, the paper ends with conclusion in Section 6.

Wheel Dynamics Model.
Suppose the vehicle runs straight on the horizontal road, ignoring the air resistance and rolling resistance, the single-wheel model is shown as in Figure 2.
is model assumes that the normal load is considered to be small and ignorable, and the dynamic differential equations are described as follows: where J is the wheel inertia, ω is the wheel angular velocity, T b is the braking torque, F x is the friction force, R is the wheel radius, m is the 1/4 vehicle mass, V is the wheel center speed, and F z is the vehicle normal load. Moreover, μ(s) expresses that the road adhesion coefficient μ is the equation of the slip ratio S, which is specifically introduced in the next section. Consider only the case of braking, the slip ratio S can be written as follows [19]: e calculation of the slip ratio S needs to know the wheel center speed V and the wheel angular velocity ω,  which can be obtained by equation (3) after knowing the braking torque T b and the real-time tire-road friction coefficient μ. Detailed calculation processes of the braking torque T b and road adhesion coefficient μ are presented in Section 4.

Tire Model.
Compared with Pacejka tire model, Dugoff tire model, and so on, Burckhardt tire model has the advantages of simple formula complexity, fewer parameters, and accuracy within the permissible range [20]. erefore, in this paper, the Burckhardt tire model which expresses the relationship between seven typical tire-road friction coefficients and slip ratio is adopted. e expression can be written as follows: where c i (i � 1, . . . , 4) are fitted by experimental data corresponding to different roads. e −c 4 v reflects the change of tire-road friction coefficient caused by the change of the wheel center speed V.
If the influence of velocity change is neglected, the model expression can be simplified as follows: where c 1 , c 2 , and c 3 for different roads are listed in Table 1 [18,20]. Figure 3 is drawn according to the data provided in Table 1.

Description of Approach.
From the observation of Figure 3, we can find that μ is an increasing function of S until reaching a critical slip value, where μ reaches μ max and then decreases. Moreover, we can also find that the peak points of other pavement curves are approximately in a straight line with a slope of "K." If we ignore the "Dry Cobblestone" road, which is not common, the trend of different road curves is more consistent.
erefore, we assume that the road curve equation to be estimated has the same trend of change.
Based on this idea, we hope to find the relationship between peak points of various known typical road curves first. en, an original mathematical method is proposed as is shown in Figure 4: First, we get the straight line whose slope is "K," which is fitted by the peak point of typical road curve ("Oblique line 1" in Figure 4). Suppose "Curve 1" is the road curve which we need to estimate, and point "b0" is the known point on "Curve 1." en, we can find a straight line ("Oblique line 2" in Figure 4) whose slope is "K" across "b0." Furthermore, we can easily find the intersection "b1 ∼ b6" of "Oblique line 2" with six other known typical road curves. Using the Lagrange interpolation method, we can express the relationship between "b0" and b i (i � 1, . . . , 6), which is the relationship between "Curve 1" and six known typical road curves.
Since the variation trend of each road curve is same, we can reverse the ordinate value of the peak point "a0" of "Curve 1" when this relationship and the coordinates of the peak points of the six typical roads are known.
Finally, based on the interval in which the ordinate value of "a0" located, the type of pavement surface that needs to be estimated can be determined.
After knowing the type of pavement surface, we can easily get the maximum tire-road friction coefficient μ max through the Burckhardt tire model.

Mathematical Problems in Engineering
It should be noted that "Oblique line 1" in Figure 4 is fitted by the peak coordinates of known typical road curves. In this paper, the coordinate values of "a1 ∼ a6," which are the intersection of "Oblique line 1" and the other six typical road curves, are approximate to the coordinate values of six typical road curves at the peak point. e overall approach used in this paper for the estimation of pavement surface condition consists of the following four steps: Firstly, the linear equation whose slope is "K" is obtained by fitting data points, which are the coordinates of peak points calculated from different typical road curves. In the second step, the observer and the wheel dynamics model are used to solve real-time road friction coefficient and real-time slip ratio, respectively. In step 3, the Lagrange interpolation method is used to obtain the ordinate value of the peak point "a0." Finally, based on the pavement surface recognition interval in which the ordinate value of "a0" located, the pavement surface condition is determined. e detailed introduction is as follows.

Solution of the Linear Equation.
As described in 4.1, the line with a slope of "K" is fitted by the vertices of a typical road curve. In order to find out the linear equation, we need to find out the coordinate of the peak point of the typical road curve first.
By using the method of finding the extremism, the derivative of equation (6) for S can be obtained: e maximum tire-road friction coefficient μ max and the optimal slip ratio S 0 which is the coordinate value of the peak point of the typical road curve are obtained, respectively, as shown in Table 2.
"Oblique line 1" can be obtained by fitting the points in Table 2. e diagram of "Oblique line 1" is displayed in Figure 5. e expression of the fitting equation "Oblique line 1" is In this paper, the intersection of "Oblique line 1" and the other six typical road curves are approximate to the coordinate values of six typical road curves at the peak point.

Solution of the Real-Time Tire-Road Friction Coefficient.
In order to achieve real-time road adhesion coefficient, the state observer is designed. Via equation (3), we can get It is a first-order nonlinear system of the wheel brake, in which the braking torque T b is the control input and the wheel speed ω is the system output. e item with the tireroad friction coefficient (mgR/J)μ can be considered as additional interference f, which is the expansion state variable of the system.
Replace the following equation as x 1 � ω, Equation (9) can be transformed into the following form: where x 1 and x 2 are the observed values of x 1 and x 2 , respectively, u and x are the input, x 1 , x 2 are the output value of the state observer, ε > 0, and α 1 and α 2 are both positive real numbers. Once the braking torque and angular velocity are known, the real-time tire-road friction coefficient can be obtained by the designed observer. e value of angular velocity is measured by the vehicle angular velocity sensor which is well equipped, and the calculation method of braking torque is introduced in next section.

Solution of the Real-Time Slip Ratio.
rough the introduction of the tire dynamics model in Section 3, the slip ratio S can be calculated as long as the braking torque T b and the real-time tire-road friction coefficient μ are known.
Since μ has been introduced in Section 4.3, the following focuses on the solution of the T b .
In this paper, the slip ratio S is taken as the control objective, and the sliding mode variable structure control method is selected to control the braking torque. When the control system moves on the sliding surface, the real-time slip ratio reaches the target slip ratio while keeping the slip ratio constant. en, _ S � 0. Using equation (4), the expression of _ S is obtained as Let _ S � 0 and bring equation (9) into equation (12) to get the braking torque T b1 on the sliding surface as follows: When the control system is outside the sliding surface, it should be controlled to the sliding surface as soon as possible, and the jitter problem should be reduced when the sliding surface is reached. e expression of the exponential approach rate is selected as [21] _ s � ε · sgn(s) − ks. (14) By introducing equation (14) into equation (12), the braking torque T b2 outside the sliding mode surface can be obtained as To sum up, we get the expression of braking torque as After knowing S and μ which is equal to the coordinate value of "b0," the expression of "Oblique line 2" can be obtained by the mathematical method. en, the intersection coordinates of "Oblique line 2" and six other typical road curves are obtained. e following is the introduction of the Lagrange interpolation method for obtaining the ordinate value of the peak point "a0" of "Curve 1."

Solution of the Ordinate Value of "a0" by Lagrange
Interpolation Method. Just as described in Section 4.1, Lagrange interpolation is mainly used to describe the relationship between curves. Comparing the coordinate value of "b0" with each intersection coordinates of "Oblique line 2" and six other typical road curves, we can get the relationship between "Curve 1" and other six typical road curves by the Lagrange interpolation method.
Let λ i (i � 1, . . . , 6) represent the relationship between "Curve 1" and six other standard road curves, respectively. e expression of the basis function of the Lagrange interpolation λ i (i � 1, . . . , 6) can be represented as where μ bi (i � 0, . . . , 6) represent the ordinate value of the point "b i " (i � 0, . . . , 6). According to λ i (i � 0, . . . , 6) and the ordinate values of the peak point of six topical road curves (as shown in Table 2), we can get the ordinate value of the point "a0," which is the peak point of "Curve 1." Define that μ ai (i � 0, . . . , 6) represent the ordinate value of point "a i " (i � 0, . . . , 6). en, the expression of μ a0 is as follows: μ a0 � λ 1 μ a1 + λ 2 μ a2 + λ 3 μ a3 + λ 4 μ a4 + λ 5 μ a5 + λ 6 μ a6 . (18)   Mathematical Problems in Engineering 4.6. Determination of the Maximum Tire-Road Friction Coefficient. After obtaining the ordinate value of the peak point "a0," we also need to know the pavement surface recognition interval in which the ordinate value of "a0" located. en, the next is a detailed introduction of interval selection. In order to improve the stability of interval recognition, the average value of the theoretical maximum tire-road friction coefficient of adjacent roads is selected as the boundary of the recognition interval. e calculation equation is as follows: where L i represents the boundary value of the recognition interval.    According to equation (19) and Table 2, the recognition interval of the typical pavement surface can be obtained as in Table 3. For the convenience of description, we numbered the pavement surface, respectively, as follows.

Simulation and Analysis
Simulation is conducted in the MATLAB/Simulink environment to show the result of the algorithm. e initial vehicle velocity is set at 100 km/h, and the main parameters of the wheel model are shown in Table 4. Figure 6 shows the state responses of simulation. e ordinate value of "a0" and road numbers obtained by the recognition algorithm are shown in Table 5: It is observed from Table 5 that excluding the error of wet cobblestone is close to 5%, the estimated error of the "a0"coordinate by Lagrange interpolation is less than 2%. e method of solving the ordinate value of "a0" by the Lagrange interpolation method can be proved to be effective and has the advantage of requiring fewer data points. en, the pavement surface condition can be deduced by the interval identification algorithm based on the recognition interval in which the ordinate value of "a0" is located. After identifying the pavement surface condition by the interval algorithm, the type of pavement can be accurately determined.
In order to further verify the effectiveness of the method, the more varied pavement surface was selected for simulation. e selected pavement surface was changed from the asphalt road with an adhesion coefficient of 0.8 to the dry asphalt road with an adhesion coefficient of 1.17. e simulation results are as follows.
According to Figure 7, it can be seen that the pavement surface condition can be accurately identified by the interval recognition algorithm based on Lagrange interpolation, and the response time is less than 0.2 seconds.

Conclusion
An original mathematical method of the interval recognition algorithm of road adhesion coefficient based on Lagrange interpolation is presented. e ordinate of the peak point of the road curve to be estimated is solved by the Lagrange interpolation method, and the pavement surface condition is deduced by the interval identification algorithm based on the recognition interval in which the ordinate value of "a0" located. e simulation results from six typical roads show that besides the cobblestone pavement which is not common in the daily road, the estimation error of the ordinate of the peak point by the Lagrange interpolation method is less than 2%, and the pavement surface condition can be accurately identified by the interval recognition algorithm. In addition, the simulation results of the varied pavement surface show that the method can also identify the road type quickly and accurately, and the response time is less than 0.2 seconds. erefore, the proposed novel mathematical method can not only maintain the accuracy of pavement surface condition estimation but also meet the requirements of the rapid response and have the advantage of requiring fewer data points.

Data Availability
e data used to support the findings of this study are included within the article.