# Research on Road Roughness Based on NARX Neural Network

Authors: Yingjie Liu; Dawei Cui
Journal: Mathematical Problems in Engineering
 (2021)

---

## Abstract

In order to solve the problem of road roughness identification, a study on the nonlinear autoregressive with exogenous inputs (NARX) neural network identification method was carried out in the paper. Firstly, a 7-DOF plane model of vehicle vibration system was established to obtain the vertical acceleration and elevation acceleration of the body, which were set as ideal input samples for the neural network. Then, based on the plane model, with common speed, the road roughness was solved as the ideal output sample of the NARX neural network, and the road roughness of B-level and C-level was identified. The results show that the proposed method has ideal identification accuracy and strong antinoise ability. The relative error of C-level road roughness is larger than that of B-level road roughness. The identified road roughness can provide a theoretical basis for analyzing the dynamic response of expressway roads.


---

## Body

## 1. Introduction

Road roughness is usually used to describe the degree of undulation of the road surface. It is the main incentive in the driving process of a vehicle and affects various aspects such as riding comfort, handling stability, component fatigue life, transportation efficiency, and fuel consumption. Therefore, it is of great significance for automotive engineering technicians to study and analyze road roughness [1].Road roughness is an important evaluation indicator of road surface, which not only affects the driving comfort of drivers and passengers but also is related to vehicle vibration, driving speed, tire friction and wear, and vehicle operating costs. When the vehicle is driving, the unevenness of the road surface will arouse vibration of the vehicle. When the vibration reaches a certain level, it will make the passengers feel uncomfortable and fatigued. The roughness of the road, which is the main vibration excitation source of the vehicle, plays a vital role in the research and development of the durability, handling stability, reliability, NVH (noise, vibration, and harshness), and several parts and systems of vehicle [2].The ride comfort of a vehicle is a comprehensive reflection of the vehicle body’s excitation on the road surface. The vibration caused by the excitation of the uneven road surface not only affects the driving smoothness and safety of the vehicle but also affects the fatigue life of the vehicle parts [3].In the analysis of vehicle ride comfort and handling stability, whether the road roughness can accurately reflect the actual road surface has a fundamental impact on the accuracy of the analysis.Therefore, the determination of a reasonable road roughness is one of the problems to be solved in the research of vehicle ride comfort and handling stability. However, the early road roughness was obtained by measuring and recording the unevenness of the road surface by a special road measurement and data processing as well as calculation. Because the road measurement is expensive, it is difficult for some inspection departments to equip it. With the development of computer technology and other related disciplines, people have begun to model road surfaces and simulate them in order to analyze and predict vehicle ride comfort. At present, scholars from various countries have conducted a lot of research on the digital simulation of road roughness, including pseudowhite noise method and harmonic superposition algorithm. However, the simulation accuracy of these methods is not very high. The paper applies load identification technology to identify road roughness [4, 5].The road roughness identification problem has been studied in the literature. A brief review is presented in what follows.Intelligent control methods have been used to study the problem of road roughness detection. For example, Gabrielli et al. [6] enhanced state of the art by introducing a Siamese Convolutional Neural Network architecture able to achieve improved results for the classification of the road surface roughness. Liu et al. [7] proposed a road roughness acquisition system design and an improved restricted Boltzmann machine deep neural network algorithm based on Adaboost Backward Propagation algorithm for road roughness detection to realize classification of different pavements. Liu et al. [8] developed an intelligent shift control method based on the identification of road roughness level.In recent years, vehicle stability improvement by considering road roughness has been a subject undergoing intense studies. For example, Liu et al. [9] proposed a new method of constructing road roughness in left and right wheel paths based on Specified Power Spectral Density (PSD) and coherence function. Gorges et al. [10] presented an experimental validation of a road roughness classification method and proposed an impact detection strategy for two-wheeled vehicles, including a classification of service loads, mild special events, and severe special events. Múčka et al. [11] presented original new results on measured whole-body vibration (WBV) on passenger’s seat surface and seat base as a function of measured road unevenness in the longitudinal direction. Mastinu et al. [12] developed a new and relatively simple analytical model. Ma et al. [13] proposed a slip ratio controller based on model predictive control (MPC) considering road roughness for electric vehicle (EV) to improve the driving/braking performance. Múčka [14] analyzed measured vibration dose value (VDV) in passenger cars as a function of longitudinal road roughness.Different estimation approaches have attracted the attention of research professionals. For example, Kang et al. [15] proposed a new road roughness estimation method by using the discrete Kalman filter with unknown input (DKF-UI). Kansake and Frimpong [16] presented mathematical models for estimating tire dynamic forces on haul roads. Kim et al. [17] presented an improved discrete Kalman filter for simultaneously estimating both all state variables and the unknown road roughness input for a vehicle suspension control system that played a key role in the ride quality and handling performance while driving the vehicle. Hamersma and Els [18] addressed a priori knowledge needed by predicting what the suspension forces would be before changing the suspension setting.Many researchers have focused on improving the comfort and vibrational response of the vehicle. For example, Liu and Cui [19] proposed a collaborative simulation algorithm for the performances between ride comfort and handling stability of a vehicle considering the connection and mutual influences. Jiang et al. [20] investigated the problem of interactive road roughness identification and proposed a novel inverse algorithm based on the knowledge of a vehicle dynamic characteristics and dynamic responses. Godbole et al. [21] aimed at understanding the likelihood of occurrence of vertical separation between the cargo and the trailer of a truck resulting from the amplified vertical vibrational response of the vehicle.In the literature reviewed, most of the multi-DOF vehicle models were simplified. Some models ignored the effect of kinematic changes on the tire. A small number of researchers had focused on the nonlinear model, and they had usually neglected the effect on solving methods from dynamic simulation and low-amplitude vibration ride comfort. For some researchers, the response of the selected vehicles was different. For example, some of them did not consider the rationality of each vehicle response as an input, and the others did not consider the appropriate combination which was conducive to application in actual vehicles of vehicle responses as an input scheme.This paper is dedicated to solving the problem of road roughness identification of different road level, providing a theoretical and methodological basis for road unevenness identification with the NARX neural network. It is different from static neural networks such as backpropagation (BP) neural network. Through feedback and memory, the NARX neural network can retain the data at the previous moment and add it to the calculation at the next moment. And also the NARX neural network can make the network be with dynamic characteristics and the system information be more complete greatly, improving the accuracy of road unevenness recognition and prediction.In the paper, based on the 7-DOF vehicle vibration model, the NARX neural network is used to establish identification vehicle responses which can provide a theoretical and methodological basis for practical application of the NARX neural network to identify road roughness. The remainder of the paper is organized as follows. Section2 presents the description of the NARX neural network. Section 3 provides the vehicle vibration model, and Section 4 presents the neural network design. In Section 5, the road roughness is identified based on NARX neural network. Finally, Section 6 concludes the paper and provides future research directions.


## 2. NARX Neural Network

### 2.1. Structure of the NARX Neural Network

NARX neural network is a classic dynamic neural network, which is composed of TDL delay layer, input layer, output layer, and hidden layer. The TDL delay layer delays the output of the neural network and feeds it back to the input layer to improve the prediction effect. The input layer concentrates the input and establishes weight and threshold operations on the input and then transfers the operation results to the hidden layer. The hidden layer can be one or more layers, which receives data, then calculates weights and thresholds, and finally transfers them to the output layer. The output layer compares the transferred calculation output with the expected output to obtain the error value and then corrects the weight, threshold, and the output data. The classic NARX neural network structure is shown in Figure1.Figure 1
Typical NARX neural network structure.The TDL performs c-order delay processing on inputxi and output yk. The number of input neurons is n; the hidden layer is 1; the number of the neurons in the hidden layer is l; the neuron in the output layer is m; the threshold of the hidden layer is aj; the threshold of the output layer is bk; and the expected output is ok; the error between the expected output and the calculated output is ek. The weight between the input layer and the hidden layer is ωij; the weight between the output delay layer and the hidden layer is ωsj; the weight between the hidden layer and the output layer is ωjk.


### 2.2. NARX Neural Network Structure Design

The number of neurons in the time delay layer, input layer, and output layer of the NARX neural network is determined by the input and output of the actual problem. Generally, the output is only delayed once.In practical applications, the hidden layer is generally taken as one layer, and the numberl of neuron nodes in the hidden layer is [22](1)l=m+n+a,where a is a constant between 0 and 10; m is the neuron in the output layer; and n is the number of input neurons.Generally, the transfer functionf of the hidden layer is set as logsig, and the transfer function g of the output layer is set as purelin:(2)fx=11+e−x,gx=x.The weights and thresholds of the neural network are adjusted by the training algorithm. The NARX neural network has 14 training algorithms. The default training algorithm is trainlm, which has fast training speed and high identification accuracy [22].When evaluating the performance of the NARX neural network for road roughness identification results, the paper provides two commonly used evaluation indicators in neural network research [23], which are the correlation coefficient (R) and the root mean square error (RMSE). R refers to the degree of linear correlation between the identified road roughness and the ideal road roughness when the neural network input is given. The higher R, the better the identification effect. The range of R is (−1, 1). If R is positive, it is a positive correlation; otherwise, it is a negative correlation. According to the R standard, the absolute value of R can be divided into 0.0∼0.2, 0.2∼0.4, 0.4∼0.6, 0.6∼0.8, and 0.8∼1.0 levels, which are, respectively, corresponding to extremely weak correlation or irrelevance, weak correlation, medium degree correlation, strong correlation, and extremely strong correlation.The calculation formula ofR is(3)R=∑i=1nyi−y¯ai−a¯∑i=1nyi−y¯2∑i=1nai−a¯2,where yi and ai are the ith points of the identified and expected road roughness, respectively; a¯ are the average values of the identified and expected road roughness respectively; and n is the number of sample points.RMSE that is known as the standard error is the arithmetic square root of variance. In the process of data training and debugging, the programmed neural network will produce certain random errors. RMSE can reflect the stability of the network identification result. The smaller the value, the more stable and reliable the identification result.The calculation formula of RMSE is(4)RMSE=∑i=1nai−yi2n.


## 2.1. Structure of the NARX Neural Network

NARX neural network is a classic dynamic neural network, which is composed of TDL delay layer, input layer, output layer, and hidden layer. The TDL delay layer delays the output of the neural network and feeds it back to the input layer to improve the prediction effect. The input layer concentrates the input and establishes weight and threshold operations on the input and then transfers the operation results to the hidden layer. The hidden layer can be one or more layers, which receives data, then calculates weights and thresholds, and finally transfers them to the output layer. The output layer compares the transferred calculation output with the expected output to obtain the error value and then corrects the weight, threshold, and the output data. The classic NARX neural network structure is shown in Figure1.Figure 1
Typical NARX neural network structure.The TDL performs c-order delay processing on inputxi and output yk. The number of input neurons is n; the hidden layer is 1; the number of the neurons in the hidden layer is l; the neuron in the output layer is m; the threshold of the hidden layer is aj; the threshold of the output layer is bk; and the expected output is ok; the error between the expected output and the calculated output is ek. The weight between the input layer and the hidden layer is ωij; the weight between the output delay layer and the hidden layer is ωsj; the weight between the hidden layer and the output layer is ωjk.


## 2.2. NARX Neural Network Structure Design

The number of neurons in the time delay layer, input layer, and output layer of the NARX neural network is determined by the input and output of the actual problem. Generally, the output is only delayed once.In practical applications, the hidden layer is generally taken as one layer, and the numberl of neuron nodes in the hidden layer is [22](1)l=m+n+a,where a is a constant between 0 and 10; m is the neuron in the output layer; and n is the number of input neurons.Generally, the transfer functionf of the hidden layer is set as logsig, and the transfer function g of the output layer is set as purelin:(2)fx=11+e−x,gx=x.The weights and thresholds of the neural network are adjusted by the training algorithm. The NARX neural network has 14 training algorithms. The default training algorithm is trainlm, which has fast training speed and high identification accuracy [22].When evaluating the performance of the NARX neural network for road roughness identification results, the paper provides two commonly used evaluation indicators in neural network research [23], which are the correlation coefficient (R) and the root mean square error (RMSE). R refers to the degree of linear correlation between the identified road roughness and the ideal road roughness when the neural network input is given. The higher R, the better the identification effect. The range of R is (−1, 1). If R is positive, it is a positive correlation; otherwise, it is a negative correlation. According to the R standard, the absolute value of R can be divided into 0.0∼0.2, 0.2∼0.4, 0.4∼0.6, 0.6∼0.8, and 0.8∼1.0 levels, which are, respectively, corresponding to extremely weak correlation or irrelevance, weak correlation, medium degree correlation, strong correlation, and extremely strong correlation.The calculation formula ofR is(3)R=∑i=1nyi−y¯ai−a¯∑i=1nyi−y¯2∑i=1nai−a¯2,where yi and ai are the ith points of the identified and expected road roughness, respectively; a¯ are the average values of the identified and expected road roughness respectively; and n is the number of sample points.RMSE that is known as the standard error is the arithmetic square root of variance. In the process of data training and debugging, the programmed neural network will produce certain random errors. RMSE can reflect the stability of the network identification result. The smaller the value, the more stable and reliable the identification result.The calculation formula of RMSE is(4)RMSE=∑i=1nai−yi2n.


## 3. Vehicle Vibration Model

A simplified 7-DOF planar model including the vertical, roll, and pitch movement of the body is shown in Figure1. And the parameters are shown in Table 1.Table 1
Parameters.
ParameterSymbolBody massmFront wheel massm1 and m2Rear wheel massm3 and m4Displacement of the body mass in the vertical directionzcYaw moment of inertia around thex-axisIxYaw moment of inertia around they-axisIyDistance of center of gravity to front axleaDistance of center of gravity to rear axlebDisplacement excitation of the road affecting the left front tirez01Displacement excitation of the road affecting the right front tirez02Displacement excitation of the road affecting the left rear tirez03Displacement excitation of the road affecting the right rear tirez04Synthesized stiffness of front tirek1 and k3Synthesized stiffness of rear tirek5 and k7Stiffness of front suspensionk2 and k4Stiffness of rear suspensionk6 and k8Track base of the vehicleLDisplacement of the left front wheels in the vertical directionz1Displacement of the right front wheels in the vertical directionz2Displacement of the left rear wheels in the vertical directionz3Displacement of the right rear wheels in the vertical directionz4Vertical displacement of body centroidz5Roll angle displacement of the bodyz6Pitch angle displacement of the bodyz7Equivalent damping coefficient of front tirec1 and c3Equivalent damping coefficient of rear tirec5 and c7Equivalent damping coefficient of front suspensionc2 and c4Equivalent damping coefficient of rear suspensionc6 and c8To establish the 7-DOF vehicle dynamic model, the following assumptions must be made firstly:(i)
The body is regarded as a rigid body with concentrated mass; only the influence of vertical, roll, and pitch vibrations of the body on the ride comfort is considered; the influence of torsional vibration of the body on the ride comfort is ignored(ii)
It is assumed that the body is a linear system(iii)
The vehicle is symmetrical along the longitudinal centerline and moves in a straight line at a constant speed(iv)
The road surface is isotropic(v)
Other vibration sources except the road surface are ignoredAccording to the Lagrange equation, the vibration equation of the vehicle can be obtained as follows (Figure2):(5)MZ¨+CZ˙+KZ=KtQ+CtQ˙,where M is the mass matrix of the vehicle; C is the damping matrix of the vehicle; K is the stiffness matrix of the vehicle; Z is the displacement vector corresponding to the vehicle’s respective degrees; Q is the pavement input vector; Kt is the stiffness matrix of the tire; Ct is the damping matrix of the tire:(6)M=diagm1m2m3m4mIxIy,Z=z1z2z3z4z5z6z7T,Q=z01z02z03z04T,Kt=k10000k30000k50000k7000000000000,Ct=c10000c30000c50000c7000000000000,K=k1+k2000−k2−L2k2ak20k3+k400−k4L2k4ak400k5+k60−k6−L2k6−bk6000k7+k8−k8L2k8−bk8−k2−k4−k6−k8k2+k4+k6+k8L2k2−k4+k6−k8−ak2−ak4+bk6+bk8−L2k2L2k4−L2k6L2k8L2k2−k4+k6−k8L24k2+k4+k6+k8L2−ak2+ak4+bk6−bk8ak2ak4−bk6−bk8−ak2−ak4+bk6+bk8L2−ak2+ak4+bk6−bk8a2k2+k4+b2k6+k8,C=c1+c2000−c2−L2c2ac20c3+c400−c4L2c4ac400c5+c60−c6−L2c6−bc6000c7+c8−c8L2c8−bc8−c2−c4−c6−c8c2+c4+c6+c8L2c2−c4+c6−c8−ac2−ac4+bc6+bc8−L2c2L2c4−L2c6L2c8L2c2−c4+c6−c8L24c2+c4+c6+c8L2−ac2+ac4+bc6−bc8ac2ac4−bc6−bc8−ac2−ac4+bc6+bc8L2−ac2+ac4+bc6−bc8a2c2+c4+b2c6+c8.Figure 2
7-DOF vehicle vibration model.


## 4. Neural Network Design

It is assumed that the input white noise of one wheel trackx is Wx, and the input white noise of the other wheel track y is Wy. According to the random vibration theory, for a linear system, when the input is a stable process, the output is also a stable process. In this way, the problem can be virtualized as follows: giving one input of the wheel track and finding the output of another wheel track. And the input and output are connected by the coherence function cohxy (f). That is, the transfer relationship between the left and right wheels is established, and another random process of the wheel track is derived from the given random process.The front wheel is set as the reference point, and the transfer function of the white noise inputWx and Wy between the left and right wheels is(7)Gxy=WysWxs=a0+a1s+⋯+ansnb0+b1s+⋯+bnsn,where Wxs and Wys are the road white noise input under the left and right wheels.Then, equation (8) can be obtained by taking the 2-order approximation:(8)Gxy=a0+a1s+a2s2b0+b1s+b2s2,where the coefficients are obtained by fitting and transforming the measured coherence functions of different roads.Equation (9) can be obtained by introducing the intermediate variable Ms:(9)Gxy=WysWxs=a0s−2+a1s−1+a2Msb0s−2+b1s−1+b2Ms,where Ms is an intermediate variable.Then, equations (10) and (11) can be obtained from the numerator and denominator, respectively:(10)WxsMs=b0s−2+b1s−1+b2Ms,(11)Wys=a0s−2+a1s−1+a2Ms.It is set thatx1=L−1s−1Ms and x2=L−1s−2Ms. At the same time, because Lx˙2=s⋅Lx2=s⋅s−2Ms=s−1Ms=Lx1, x˙2=x1.The correlation state equation of the two wheels is obtained by taking inverse Laplace transform on equation (10) and combining with x˙2=x1:(12)x˙1x˙2=−b1b2−b0b210x1x2+1b20Wx.Similarly, the output equation can be obtained by performing inverse Laplace transform on equation (11).(13)Wy=a2x˙1+a0x2+a1x1.Then, equation (14) can be obtained by substituting x˙1 into equation (13)(14)Wy=a1−a2b1b2a0−a2b0b2x1x2+a2b2Wx.When the road surface is generated, it is possible to obtain the white noise inputWy of another wheel track through the intermediate state variables x1 and x2 of equation (12) while generating one wheel track Wx. And then, according to equation (13), the space correlation time-domain model of the first and second wheel tracks can be obtained.The road surface roughness, which is the input of vehicle vibration, mainly uses the road surface power spectrum density to describe its statistical characteristics. According to ISO/TC108/SC2N67, the classification standard of different road level is shown in Table2.Table 2
Classification standard of road level.
Road levelGqn0/10−6m3Lower limitGeometric meanUpper limitB3264128C128256512Gqn0/10−6m3 is the power spectrum density.For the simulation, the initial values of the weightsωij, ωjk, and ωsj and thresholds aj and bk are set as a random number between [0, 1]. And the maximum number of iteration is set as 1000.The NARX neural network is programmed by MATLAB software. First, the input and target output of the training set are imported into the neural network program. The hidden layer is a single layer, the number of neurons in the input layer is 7, and the number of neurons in the output layer is 1. The simulation time is 12.5 s, and the sampling interval is 0.01 s. The number of sampling points is 1250. The front wheel road surface unevenness of 850 points is used as the output of the training set, and the other front wheel road surface unevenness of 400 points is used as the output of the test set to determine the effect of neural network recognition.In order to verify that the algorithm is valid, the solution is calculated by MATLAB software in a computer whose CPU is 2.8 GHz/Pentium IV and the operating system is Window XP.And the parameters arem = 920 kg, Iy = 948 kg m2, Ix = 239 kg m2, a = 1.109 m, b = 1.09 m, L = 1.22 m, m1 = m2 = m3 = m4 = 25 kg, c1 = c3 = c5 = c7 = 40 N·s/m, k2 = k4 = 16.7 kN/m, k6 = k8 = 21.1 kN/m, c6 = c8 = 3340 N·s/m, c2 = c4 = 2845 N·s/m, and k1 = k3 = k5 = k7 = 147 kN/m.The unevenness excitations of a vehicle on the B- level and C-level roads are shown in Figures3–6. It can be seen from the figures that the unevenness of the C-level road is larger than that of the B-level road. And there is a time delay in the unevenness of the front and rear wheels.Figure 3
Road roughness of the left front and rear wheels (B-level).Figure 4
Road roughness of the right front and rear wheels (B-level).Figure 5
Road roughness of the left front and rear wheels (C-level).Figure 6
Road roughness of the right front and rear wheels (C-level).


## 5. Road Roughness Identification Based on NARX Neural Network

In this paper, based on the 7-DOF vehicle dynamic model, the road roughness identification technology is studied. Among the dynamic responses of the vehicle, the vehicle response that is reasonable and has a significant impact on the identification result is selected, and the repeated vehicle responses are removed to establish a nonlinear relationship between input and output data in the NARX neural network. The NARX neural network is programmed by MATLAB software. Firstly, the input and target output of the training set are imported into the neural network program, and then the input and output are set as 1 : 2.


### 5.1. Simulation Results of B-Level Road Roughness

According to the 7-DOF vehicle data, the MATLAB software is used for programming, and the vertical acceleration and elevation angular acceleration of the body of the B-level road surface, which are substituted to the trained network for calculation, are simulated as shown in Figures7 and 8, respectively.Figure 7
Vertical acceleration of the body of B-level road surface.Figure 8
Elevation angular acceleration of B-level road surface.Figure9 shows the result of the NARX network identifying the road surface excitation of the left front wheel on the B-level road. Figure 10 shows the relative error between the identification value and the fitted value of the filtered white noise. Figures 11 and 12 are the simulation results after adding noise. It can be seen that the relative error is small, the antinoise ability is strong, and the identification accuracy is higher.Figure 9
Identification and fitted values of road roughness (without noise).Figure 10
Relative error of road roughness (without noise).Figure 11
Identification and fitted values of road roughness (with noise).Figure 12
Relative error of road roughness (with noise).


### 5.2. Simulation Results of C-Level Road Roughness

According to the 7-DOF vehicle data, the MATLAB software is used for programming, and the vertical acceleration and elevation angular acceleration of the body of the C-level road surface, which are substituted to the trained network for calculation, are simulated as shown in Figures13 and 14, respectively.Figure 13
Vertical acceleration of the body of C-level road surface.Figure 14
Elevation angular acceleration of C-level road surface.Figure15 shows the result of the NARX network identifying the road surface excitation of the left front wheel on the C-level road. Figure 16 shows the relative error between the identification value and the fitted value of the filtered white noise. Figures 17 and 18 are the simulation results after adding noise. From the figures, it can be seen that the identification accuracy is higher.Figure 15
Identification and fitted values of road roughness (without noise).Figure 16
Relative error of road roughness (without noise).Figure 17
Identification and fitted values of road roughness (with noise).Figure 18
Relative error of road roughness (with noise).The evaluation index values for different road levels are shown in Table3.Table 3
Evaluation index values for different road levels.
Road levelR/%RMSEB95.6900.0085C92.3800.0268It can be seen from Table3 that the NARX neural network has an ideal identification effect on different levels of roads. It can be seen from the table that R has reached more than 90%, indicating that the actual road roughness identification is very consistent with the expected road roughness.And also, compared with the traditional method (BP neural network), the results indicate that the NARX neural network has better recognition accuracy than that of the BP neural network for road roughness identification.


### 5.3. Test Verification Based on ADAMS/View

The virtual test system not only can be used as the preparatory work for the real test but also can replace the traditional test to a certain extent. Compared with traditional tests, it has the following advantages:(i)
It can greatly reduce the number of prototype manufacturing tests and shorten the new product test cycle. At the same time, the cost of actual testing is reduced(ii)
The application of virtual test technology in the development of complex products can realize the mutual feedback of information between designers and product users in the design phase so that designers can absorb and adopt suggestions for new products in all aspects(iii)
The virtual test technology replaces the actual test, and the test process can be repeatedDue to the limitations of real vehicle test conditions and the fact that the ADAMS virtual test has the above advantages, the paper uses the ADAMS/View virtual simulation test instead of the real vehicle test to verify the neural network identification method.A vehicle is a complex vibration system. Generally, different types of simplified modeling of the vehicle system are required according to the research objectives and actual effects when analyzing the driving comfort and handling stability of the vehicle. The vehicle model mainly includes vehicle body model, double-wishbone front independent suspension model, steering mechanism model, inclined arm rear suspension model, tire model, and road surface spectrum.In the process of vehicle modeling, the following simplifications and assumptions are made for each component model:(i)
The body part does not need to be modeled according to the actual structure and is integrated into a rigid body(ii)
The sprung mass is integrated into the body rigid body(iii)
The tires, dampers, springs, rubber components, and so on are deformable, and the other parts are regarded as rigid bodies(iv)
The friction between the motion pairs is ignoredThe body is simplified into a rigid body, and the mass is concentrated in the center of mass. Since the paper mainly researches the road roughness identification of the vehicle, the specific structure of other parts such as the engine and brake assembly is ignored. It adopts a double-wishbone independent suspension model, which is mainly composed of kingpin, upper and lower wishbone, pull arm, steering knuckle, steering tie rod, and wheels. Each spherical hinge provides 3 constraint equations, each rotary hinge provides 4 constraint equations, and each movable hinge provides 5 constraint equations. The steering mechanism model is mainly composed of steering tie rod, steering rocker arm, empty rocker arm, steering wheel, steering drive shaft, and so on. It adopts an inclined arm type rear suspension model, which is mainly composed of a left inclined arm and a right inclined arm.ADAMS/View provides five types of tire models: Delft tire model, Fiala tire model, Smithers tire model, UA tire model, and User Defined tire model. Among them, the Fiala tire model, UA tire model, and User Defined tire model are analytical models, and the Delft tire model and Smithers tire model are test models. UA tire is the most used tire model, and the paper uses the UA tire model.Because the input of the vehicle ride comfort analysis is a random road roughness input, a random road surface should be appropriately constructed as the input of the virtual prototype model established above. The generation of random roads must not only conform to a certain random distribution law but also meet the requirements of tires. Random roads must meet the tire requirements: firstly, the position of the road should be below the tire; secondly, the upward direction of the road should point to the side where the tire is located; thirdly, the size of the road should be determined according to the needs of simulation.In ADAMS, an uneven road is a three-dimensional surface composed of a series of triangular plane elements. The triangular element is composed of a series of nodes. Thex and y coordinates of these nodes must meet certain rules, and the z coordinate only represents the width of the road. These nodes form a pavement unit according to a certain rule. And then, the static and dynamic friction coefficient in the pavement unit is set to simulate the real road surface.A complete vehicle model is established by connecting each subsystem model established above according to the corresponding constraints shown in Figure19.Figure 19
Vehicle model based on ADAMS/View.The built vehicle model and the random road surface generated by the road generation software are selected. And the simulation time and step length are set in the ADAMS software for simulation so that the vehicle can drive at a constant speed of 60 km/h. When the vehicle is simulated for a period of time on random roads of grades A, B, D, and E, the measured vertical acceleration and elevation acceleration of the body are set as the neural network training samples.The vertical acceleration and elevation angular acceleration of the body obtained by ADAMS/View on the C-level road, as shown in Figures20 and 21, are chosen, respectively, and substituted into the trained network for calculation. The identification value and virtual test value are shown in Figure 22. The relative error of identification is shown in Figure 23. It can be seen that the identification value is not much different from the ADAMS/View virtual test value, which verifies the effectiveness and correctness of using the NARX neural network to identify road roughness.Figure 20
Vertical acceleration of the body of C-level road surface.Figure 21
Elevation angular acceleration of C-level road surface.Figure 22
Identification and test values of C-level road roughness.Figure 23
Relative error of C-level road roughness.


## 5.1. Simulation Results of B-Level Road Roughness

According to the 7-DOF vehicle data, the MATLAB software is used for programming, and the vertical acceleration and elevation angular acceleration of the body of the B-level road surface, which are substituted to the trained network for calculation, are simulated as shown in Figures7 and 8, respectively.Figure 7
Vertical acceleration of the body of B-level road surface.Figure 8
Elevation angular acceleration of B-level road surface.Figure9 shows the result of the NARX network identifying the road surface excitation of the left front wheel on the B-level road. Figure 10 shows the relative error between the identification value and the fitted value of the filtered white noise. Figures 11 and 12 are the simulation results after adding noise. It can be seen that the relative error is small, the antinoise ability is strong, and the identification accuracy is higher.Figure 9
Identification and fitted values of road roughness (without noise).Figure 10
Relative error of road roughness (without noise).Figure 11
Identification and fitted values of road roughness (with noise).Figure 12
Relative error of road roughness (with noise).


## 5.2. Simulation Results of C-Level Road Roughness

According to the 7-DOF vehicle data, the MATLAB software is used for programming, and the vertical acceleration and elevation angular acceleration of the body of the C-level road surface, which are substituted to the trained network for calculation, are simulated as shown in Figures13 and 14, respectively.Figure 13
Vertical acceleration of the body of C-level road surface.Figure 14
Elevation angular acceleration of C-level road surface.Figure15 shows the result of the NARX network identifying the road surface excitation of the left front wheel on the C-level road. Figure 16 shows the relative error between the identification value and the fitted value of the filtered white noise. Figures 17 and 18 are the simulation results after adding noise. From the figures, it can be seen that the identification accuracy is higher.Figure 15
Identification and fitted values of road roughness (without noise).Figure 16
Relative error of road roughness (without noise).Figure 17
Identification and fitted values of road roughness (with noise).Figure 18
Relative error of road roughness (with noise).The evaluation index values for different road levels are shown in Table3.Table 3
Evaluation index values for different road levels.
Road levelR/%RMSEB95.6900.0085C92.3800.0268It can be seen from Table3 that the NARX neural network has an ideal identification effect on different levels of roads. It can be seen from the table that R has reached more than 90%, indicating that the actual road roughness identification is very consistent with the expected road roughness.And also, compared with the traditional method (BP neural network), the results indicate that the NARX neural network has better recognition accuracy than that of the BP neural network for road roughness identification.


## 5.3. Test Verification Based on ADAMS/View

The virtual test system not only can be used as the preparatory work for the real test but also can replace the traditional test to a certain extent. Compared with traditional tests, it has the following advantages:(i)
It can greatly reduce the number of prototype manufacturing tests and shorten the new product test cycle. At the same time, the cost of actual testing is reduced(ii)
The application of virtual test technology in the development of complex products can realize the mutual feedback of information between designers and product users in the design phase so that designers can absorb and adopt suggestions for new products in all aspects(iii)
The virtual test technology replaces the actual test, and the test process can be repeatedDue to the limitations of real vehicle test conditions and the fact that the ADAMS virtual test has the above advantages, the paper uses the ADAMS/View virtual simulation test instead of the real vehicle test to verify the neural network identification method.A vehicle is a complex vibration system. Generally, different types of simplified modeling of the vehicle system are required according to the research objectives and actual effects when analyzing the driving comfort and handling stability of the vehicle. The vehicle model mainly includes vehicle body model, double-wishbone front independent suspension model, steering mechanism model, inclined arm rear suspension model, tire model, and road surface spectrum.In the process of vehicle modeling, the following simplifications and assumptions are made for each component model:(i)
The body part does not need to be modeled according to the actual structure and is integrated into a rigid body(ii)
The sprung mass is integrated into the body rigid body(iii)
The tires, dampers, springs, rubber components, and so on are deformable, and the other parts are regarded as rigid bodies(iv)
The friction between the motion pairs is ignoredThe body is simplified into a rigid body, and the mass is concentrated in the center of mass. Since the paper mainly researches the road roughness identification of the vehicle, the specific structure of other parts such as the engine and brake assembly is ignored. It adopts a double-wishbone independent suspension model, which is mainly composed of kingpin, upper and lower wishbone, pull arm, steering knuckle, steering tie rod, and wheels. Each spherical hinge provides 3 constraint equations, each rotary hinge provides 4 constraint equations, and each movable hinge provides 5 constraint equations. The steering mechanism model is mainly composed of steering tie rod, steering rocker arm, empty rocker arm, steering wheel, steering drive shaft, and so on. It adopts an inclined arm type rear suspension model, which is mainly composed of a left inclined arm and a right inclined arm.ADAMS/View provides five types of tire models: Delft tire model, Fiala tire model, Smithers tire model, UA tire model, and User Defined tire model. Among them, the Fiala tire model, UA tire model, and User Defined tire model are analytical models, and the Delft tire model and Smithers tire model are test models. UA tire is the most used tire model, and the paper uses the UA tire model.Because the input of the vehicle ride comfort analysis is a random road roughness input, a random road surface should be appropriately constructed as the input of the virtual prototype model established above. The generation of random roads must not only conform to a certain random distribution law but also meet the requirements of tires. Random roads must meet the tire requirements: firstly, the position of the road should be below the tire; secondly, the upward direction of the road should point to the side where the tire is located; thirdly, the size of the road should be determined according to the needs of simulation.In ADAMS, an uneven road is a three-dimensional surface composed of a series of triangular plane elements. The triangular element is composed of a series of nodes. Thex and y coordinates of these nodes must meet certain rules, and the z coordinate only represents the width of the road. These nodes form a pavement unit according to a certain rule. And then, the static and dynamic friction coefficient in the pavement unit is set to simulate the real road surface.A complete vehicle model is established by connecting each subsystem model established above according to the corresponding constraints shown in Figure19.Figure 19
Vehicle model based on ADAMS/View.The built vehicle model and the random road surface generated by the road generation software are selected. And the simulation time and step length are set in the ADAMS software for simulation so that the vehicle can drive at a constant speed of 60 km/h. When the vehicle is simulated for a period of time on random roads of grades A, B, D, and E, the measured vertical acceleration and elevation acceleration of the body are set as the neural network training samples.The vertical acceleration and elevation angular acceleration of the body obtained by ADAMS/View on the C-level road, as shown in Figures20 and 21, are chosen, respectively, and substituted into the trained network for calculation. The identification value and virtual test value are shown in Figure 22. The relative error of identification is shown in Figure 23. It can be seen that the identification value is not much different from the ADAMS/View virtual test value, which verifies the effectiveness and correctness of using the NARX neural network to identify road roughness.Figure 20
Vertical acceleration of the body of C-level road surface.Figure 21
Elevation angular acceleration of C-level road surface.Figure 22
Identification and test values of C-level road roughness.Figure 23
Relative error of C-level road roughness.


## 6. Conclusion

A 7-DOF ride comfort model is established, and the NARX neural network and its training process and structure design are studied and summarized. The input and output of the NARX neural network are obtained through simulation. The correlation coefficient and root mean square error are used as the evaluation indicators of the identification effect of the NARX neural network. The results show that the obtained neural network has anR of 95.690% and an RMSE of 0.0085 at a common speed of 60 km/h on B-level road and an R of 92.380% and an RMSE of 0.0268 at a common speed of 60 km/h on C-level road, which has ideal identification accuracy. And also the maximum value of the relative error of the identified B-level and C-level road surface is 0.22% and 0.41%. The maximum value of the relative error of the identified B-level and C-level road surface after adding noise is 3.7% and 4.65%. Therefore, using the method to identify road roughness has a small relative error and strong antinoise ability. The identified road roughness will have an important reference value for in-depth analysis of the dynamic response of the highway pavement structure.The research in this paper is currently limited to the overall level of the road without considering the impact of the road slope. In fact, road slope recognition is also a very important aspect of vehicle dynamics research. Combining road slope recognition and slope road surface unevenness recognition will become a direction for future research. In addition, how to establish real road test verification is also another direction for future research. The advancement of these two directions requires more comprehensive research and analysis of road roughness recognition methods based on neural networks in theory and methods.



---

*Source: 9173870-2021-12-14.xml*