Establishment of an Eleven-Freedom-Degree Coupling Dynamic Model of Heavy Vehicle-Pavement

: Considering the actual situation of moving vehicles acting on road surfaces, a dynamic model of a heavy vehicle–road coupling system was established based on the traditional vehicle–road vibration model. Firstly, a seven-degree-of-freedom vehicle model was established, and the vibration characteristics of the road subsystem were considered part of the whole system. The excitation effect of road roughness on the vehicle model was considered, and the dynamic model of the coupling system was ﬁnally obtained by combining the displacement compatibility principle of the contact relationship between the wheels and the road surface. The results show that the maximum bending value of the surface course reaches 2.37 mm. The maximum shear stress in the middle part of the surface course is 43,858 Pa. The vertical dynamic stress in the middle part of the surface course is larger, reaching 119,373 Pa, while the value of the vertical dynamic stress in the subgrade is much smaller, reaching 5824 Pa. The coupling dynamic model can reﬂect the relationship between the moving vehicle and the road and the dynamic performance, which not only provides theoretical support for the design parameters of heavy-duty vehicles but also provides a reference for the design of road durability. is much smaller, and the vertical dynamic stress in the middle of the subgrade is 5824 Pa. In the vertical direction of the road, the vertical dynamic stress of each layer of the road structure decreases rapidly. In contrast, the dynamic stress on the surface of the asphalt base course is 0.602 times the dynamic stress on the surface course, and when it reaches the bottom of the subbase, the dynamic stress decreases to 0.038 times the surface stress.


Introduction
In recent years, heavy-duty vehicle transportation has become important all around the world [1,2]. The concept of dynamic interactions between the vehicle and the road surface has been widely studied, and a series of problems caused by the coupling vibrations of the vehicle structure and road have become very prominent. Thanks to the in-depth study of this issue by scholars, research on the dynamic interaction between vehicles and road systems [3,4], and the vehicle-road dynamics model of coupled systems [5][6][7][8][9][10][11], our understanding in this area is continuously developing. The key problems that hinder the long-term development of heavy haul transportation have been partially solved, including the impact coefficient of vehicles, the vibration of vehicles under the uneven road surface, and the instantaneous response of road vibrations.
As mentioned above, the dynamic characteristics of the roadbed are regarded as one of the important components of the system. At present, a vehicle-road system can be divided into a vehicle vibration model and a road subsystem model. The vehicle is the moving part, which is the main excitation part that causes the vehicle and road vibration. The pavement and the following parts are the road subsystem. Due to the long-term bearing of the dynamic load of the upper vehicle, including the repeated load of heavy vehicle movement, subgrade diseases form easily. In particular, the deformation of flexible asphalt pavement is often affected by many factors such as traffic congestion, joints, load capacity, and repeated vehicle loads [12,13]. All these factors affect the normal service of the road and make the road surface irregular, seriously affecting the normal operation of the road and even damaging it. The above analysis shows that the damage of roads not only involves damage to the road surface but also the excessive deformation of the roadbed. Therefore, when analyzing the dynamic interaction between vehicle and road, the vehicle vibration excitation is the key problem, which is, therefore, one of the main concerns of this work.
Scholars have studied the dynamic problems of vehicle and highway subgrades from different perspectives [14,15]. However, few of the existing models consider the vehicle-road system as a whole, and most models analyze the road dynamic performance under the dynamic load of vehicles. The literature [16] describes the lateral dynamics model of vehicles under road environment excitation. The authors of [17,18] studied the dynamic response of the vehicle-pavement coupling system based on the nonlinear Timoshenko beam method and the vehicle-road dynamic response of the multi-degree-of-freedom vehicle with a double-layer, rectangular, thin plate. Zhang et al. [19] studied the road deformation and crack propagation path under the action of a quarter of vehicles using a modified two-parameter foundation plate. The above analysis fails to fully reflect the interaction of a vehicle-subgrade system in the process of vehicle moving, and this is another main objective of the present study.
In the present study, we establish a seven-degree-of-freedom vehicle dynamic model and use the harmonic superposition method to simulate the process of road random vibration. The dynamic interaction between vehicle and road is analyzed theoretically. Then, the validity of the model is verified by numerical analysis, and the dynamic response is analyzed. The organization of the paper is as follows. In Section 2, the dynamic model and equilibrium equation of heavy load vehicles are established. In Section 3, the coupling dynamic equation of heavy load vehicle-road is established. Then, the vehicle-road system is analyzed from multiple angles through numerical simulation, and the validity of the established equation is verified in Section 4.

Simplified Vehicle Road Dynamic Model
According to different research needs, traditional vehicle models can come in three different styles. They are the quarter two-degree-of-freedom vehicle model, the half fourdegree-of-freedom vehicle model, and the whole vehicle model, respectively. The first two models are widely used due to their simple calculation, but the parameters of vibration are not fully considered and the calculation accuracy is not high. Based on the traditional model, a 7-DOF vehicle model is established in this paper, as shown in Figure 1. In this model, the car body and wheels are treated as rigid bodies. A vehicle consisting of a single body and several wheels passes through the line at a uniform speed; for the body itself, the freedom of nod vibration, deflection vibration, and vertical vibration will be considered [20], along with the vibration of the four wheels, for a total of seven degrees of freedom. Furthermore, it is assumed that the wheels of the vehicle are in close contact with the road surface. In this paper, we mainly study the vibration of the vehicle-road system caused by the excitation of road roughness. In particular, this paper studies the effects on vehicle ride performance, including comfort, safety, and irregularity. To simplify the analysis process, it was simplified into the 7-DOF vehicle structure model shown in Figure 1.

Dynamic Balance of Vehicle-Road System
Based on the above description, combined with the vehicle vibration balance equation and the road vibration balance equation, the dynamic balance equation of the vehicleroad system is established. The stress analysis of the vehicle model and the body are shown in Figures 2 and 3, respectively.

Dynamic Balance of Vehicle-Road System
Based on the above description, combined with the vehicle vibration balance equation and the road vibration balance equation, the dynamic balance equation of the vehicleroad system is established. The stress analysis of the vehicle model and the body are shown in Figures 2 and 3, respectively.

Dynamic Balance of Vehicle-Road System
Based on the above description, combined with the vehicle vibration balance equation and the road vibration balance equation, the dynamic balance equation of the vehicleroad system is established. The stress analysis of the vehicle model and the body are shown in Figures 2 and 3, respectively.
and F tc1 , F ts1 , F tc2 , F ts2 , F tc3 , F ts3 , F tc4 , and F ts4 represent the interaction force between the unsprung mass part of the vehicle and the body. According to the force balance conditions in Figures 2 and 3, the balance equations of vertical vibration, pitch vibration, and roll vibration of the body are as follows [21][22][23]: The force analysis of the vehicle's unsprung mass model is shown in Figure 4.
The force analysis of the vehicle's unsprung mass model is shown in Figure 4. Where i = 1, 2, 3, and 4 represent the right front unsprung mass part, the rig unsprung mass part, the left front unsprung mass part, and the left rear unsprung part, respectively. The force balance equations of the four unsprung mass parts are lows: The stress analysis of the wheel is drawn according to the stress of the vehicle as shown in Figure 5. Where i = 1, 2, 3, and 4 represent the right front unsprung mass part, the right rear unsprung mass part, the left front unsprung mass part, and the left rear unsprung mass part, respectively. The force balance equations of the four unsprung mass parts are as follows: The stress analysis of the wheel is drawn according to the stress of the vehicle wheel, as shown in Figure 5. Where, i = 1, 2, 3, and 4 represent the four wheels, respectively. According to th balance, the balance equations are obtained: Where, i = 1, 2, 3, and 4 represent the four wheels, respectively. According to the force balance, the balance equations are obtained: ..

Z l3
(10) In the above equations and pictures (Equations (1)- (11) and Figures 1-5), m s represents the body mass; θ is the body nod displacement angle, ψ represents the body turnover displacement angle; m ti represents the unsprung mass of vehicle front and rear; m li represents wheel mass; C ti represents the damping coefficient of the corresponding wheel before and after the vehicle; C si represents the damping coefficient of the corresponding suspension system before and after the vehicle; K ti represents the stiffness coefficients of the corresponding wheels before and after the vehicle; K si represent the stiffness coefficients of the front and rear suspension of the vehicle; Z s represents the vertical vibration displacement of the body; Z ti represent the vertical vibration displacements of the four unsprung mass parts corresponding to the front and rear of the vehicle; Z li represents the vertical vibration displacement of the front and rear wheels of the vehicle; J X represents the rotational inertia of the body around the x-axis; J Y represents the rotational inertia of the body around the y-axis; and f ti represents the dynamic loads of the four wheels.
Based on the balance equation of each part of the vehicle system mentioned above, the vector expression can be obtained as: where M represents the mass matrix of the vehicle; C represents the damping matrix of the vehicle; K represents the stiffness matrix of the vehicle; Z represents the displacement matrix of the vehicle; F t represents the dynamic load vector acted on the road surface by the vehicle system. The above matrix expressions are as follows: The boundary constraint conditions of the finite element model of the road multi-layer system need to be simplified according to the force and loading conditions of the road under the actual vehicle traveling conditions. When the vehicle is driving in the center of the road, the response at the far side of the road is negligible. Because of the thick soil foundation, the response at the bottom is too small and negligible. Therefore, it can be known from Saint-Venant's principle that longitudinal constraints are imposed on both sides of the principle and fixed displacement are imposed on the bottom to obtain the finite element model of the road [24][25][26], as shown in Figure 6. road under the actual vehicle traveling conditions. When the vehicle is driving in the cen ter of the road, the response at the far side of the road is negligible. Because of the thick soil foundation, the response at the bottom is too small and negligible. Therefore, it can be known from Saint-Venant's principle that longitudinal constraints are imposed on both sides of the principle and fixed displacement are imposed on the bottom to obtain the finite element model of the road [24][25][26], as shown in Figure 6.

Pavement Vibration Model
The roughness of the road surface is called unevenness in highway engineering. I originates from the accidental factors that cannot be controlled artificially in the process of construction and the irregular vertical undulation of the ground caused by repeated loads caused by vehicles. The unevenness of the road will cause the vibration of the vehicle Based on the establishment of the road finite element model in the previous section, accord ing to the forced vibration equation of the damped system with multiple degrees of free dom, the vibration balance equation of the road structure model can be obtained as follows where

Pavement Vibration Model
The roughness of the road surface is called unevenness in highway engineering. It originates from the accidental factors that cannot be controlled artificially in the process of construction and the irregular vertical undulation of the ground caused by repeated loads caused by vehicles. The unevenness of the road will cause the vibration of the vehicle. Based on the establishment of the road finite element model in the previous section, according to the forced vibration equation of the damped system with multiple degrees of freedom, the vibration balance equation of the road structure model can be obtained as follows: δ} represents the acceleration vector, and { f } represents the vertical load vector, including the dynamic load and static load.
To solve the common vibration problem caused by vehicles acting on the road, the vehicle vibration balance equation and road vibration balance equation in Equations (12) and (13) are needed to calculate the dynamic response of the whole system. To simplify the calculation without losing accuracy and reliability, the high-order term of vibration mode superposition is discarded, and only the contact between the vehicle and the road is considered; that is, the node between the wheel and the road. The vertical vibration displacement at any cross-section x position of the road pavement is obtained by superposition of the vibration mode functions of each order of the road: where A n represents the corresponding generalized coordinates and φ n (x i ) is the vibration mode function at the corresponding road horizontal position x and the nth order.

Model of Road Roughness
The methods of road unevenness can be divided into static cross-section measurement, dynamic cross-section measurement, reaction level measurement system, and subjective evaluation method. The road smoothness is expressed as the superposition of some sine or cosine waves with random phases. For the superposition of N similar sine waves, the smoothness of a random road surface can be expressed as [24]: where, Z(x) represents the unevenness value of the random pavement, x represents the displacement of the pavement along the horizontal direction, θ represents the random number of [0-2π], and n mid_i represents the intermediate value of the spatial frequency of the pavement flatness between each cell. Assuming the vehicle speed is v, let x = vt, and the first-order and second-order derivatives of the above formula can obtain the speed and acceleration of a random road, respectively. . ..

Wheel-Road Displacement Coupled Vibration
It has been assumed that the vehicle is always in contact with the pavement during driving, and the vertical displacement of the pavement is divided into pavement flatness z p (x i ) and pavement vibration vertical displacement z r (x i ); then, the vertical displacement of the wheel can be expressed as: where z p (x i ) is the road surface flatness value of tire i (i represents the wheel position) at the horizontal direction x of the road. Substituting Equation (14) into Equation (18), we get: The contact between vehicle wheel and the road surface can be simplified according to Hertz's contact theory. Consider the coupling relationship between the wheel and the road surface to be a nonlinear contact force, as shown in Figure 7. The contact between vehicle wheel and the road surface can be simplified according to Hertz's contact theory. Consider the coupling relationship between the wheel and the road surface to be a nonlinear contact force, as shown in Figure 7. According to the contact relation in Hertz's law, when further simplified the vertical contact force between the wheel and the road surface is expressed as:

Establishment of Vehicle-Road Coupling Dynamic Analysis Model
As mentioned above, the vehicle vibration balance equation and road vibration balance equation are combined to obtain the vibration coupling relationship model under the condition of displacement compatibility: According to the contact relation in Hertz's law, when further simplified the vertical contact force between the wheel and the road surface is expressed as:

Establishment of Vehicle-Road Coupling Dynamic Analysis Model
As mentioned above, the vehicle vibration balance equation and road vibration balance equation are combined to obtain the vibration coupling relationship model under the condition of displacement compatibility: ..
Symmetry 2022, 14, 250 11 of 23 In the above equation, the total number of equations is N + 7, where N represents the modal equations of the road subsystem and generally takes the low-order vibration mode. The dynamic equilibrium equation will be solved by the Newmark-Beta method. The vibration equation of the standard vehicle-road system can be obtained by rearranging Equation (21): where [M] represents the total mass matrix, [C] represents the total damping matrix, [K] represents the total stiffness matrix, and [F] represents the total load matrix. The matrix expressions are as follows: ..

Z t4
It is difficult to see from the above matrices that the total mass matrix [M], the total damping matrix [C], and the total stiffness matrix [K] change with the road position x (x 1 , x 2 , x 3 , and x 4 represent the positions of the four wheels). Therefore, the vehicle road system coupling dynamic balance equations can be regarded as a highly second-order nonlinear differential equation. By solving the equation, the vertical vibration displacement Z s , vertical velocity . Z s , and vertical acceleration .. Z s of the vehicle body; the vertical vibration displacement Z t1 , Z t2 , Z t3 , and Z t4 of each unsprung mass part of the vehicle body, the

ψ.
The vehicle is regarded as a multi-stiffness system. According to D'Alembert's principle, considering the dynamic balance conditions of each rigid body in the system, the road vibration balance equation and the vehicle dynamic balance equation are listed according to the coupling model under the contact force condition in Figure 8. The dynamic balance equation of the vehicle-road system is obtained by com Equation (20). The total number of equations is N + 7, where N represents the num modal equations of the path sub-model.
The Newmark-β and Park integral methods are used to solve the above equ and the responses of the vehicle and road in the time domain, such as the force, d ment, and vibration velocity, can be obtained [27].

Theoretical Analysis of Vehicle-Road System
According to the vehicle-road coupling dynamic equation and the dynamic equation listed above, the following vehicle driving mass and road design parame be obtained using the Newmark-β method in a self-programmed program. It is as that the heavy vehicle has a mass of 20 t and runs at a speed of 80 km/h on a class The road structure adopts planar 8-node PLANE82 unit, and the road model take along the road longitudes. The calculated step time is 0.02. See Appendixes A an detailed parameters.

Maximum Deflection of Pavement
According to the relevant parameters, an instantaneous dynamics calculation ried out in this paper. Figure 9 shows the curve of deflection value change at a n road surface under a moving load. Table 1 shows the deflection values of each laye road structure in the vertical direction. It can be seen that the bending value of the The dynamic balance equation of the vehicle-road system is obtained by combining Equation (20). The total number of equations is N + 7, where N represents the number of modal equations of the path sub-model.
The Newmark-β and Park integral methods are used to solve the above equations, and the responses of the vehicle and road in the time domain, such as the force, displacement, and vibration velocity, can be obtained [27].

Theoretical Analysis of Vehicle-Road System
According to the vehicle-road coupling dynamic equation and the dynamic balance equation listed above, the following vehicle driving mass and road design parameters can be obtained using the Newmark-β method in a self-programmed program. It is assumed that the heavy vehicle has a mass of 20 t and runs at a speed of 80 km/h on a class B road. The road structure adopts planar 8-node PLANE82 unit, and the road model takes 100 m along the road longitudes. The calculated step time is 0.02. See Appendices A and B for detailed parameters.

Maximum Deflection of Pavement
According to the relevant parameters, an instantaneous dynamics calculation is carried out in this paper. Figure 9 shows the curve of deflection value change at a node of road surface under a moving load. Table 1 shows the deflection values of each layer of the road structure in the vertical direction. It can be seen that the bending value of the surface layer is the largest, and the bending value of the surface layer and surface layer changes relatively little. When the subgrade is used, the bending settlement value changes greatly, and the fundamental reason is that the material properties of subgrade soil differ greatly from the surface layer, base layer, and cushion layer. The smaller the elastic modulus of the material is, the larger the vertical deformation is. As the bottom layer of the soil foundation is set as a constraint in the calculation, its displacement is zero.

Tensile Stress Analysis
The surface of the asphalt pavement undergoes deflection changes unde dynamic load of the vehicle. At the same time, the road structure layer al bending and tensile stress. Under repeated vehicle loads, this bending and will cause bending and tensile fatigue damage to the asphalt road structure Figure 10 shows the time-history curve of the maximum tensile stress of each structural layer on the road. It can be seen from the figure that the ten the road is sometimes positive and sometimes negative, but it is the maximu surface course and the base course. Figure 11 shows the horizontal stress nephogram at the position of t tensile stress at each node of the road structure layer. From this figure, it can the maximum tensile stress mainly occurs in the road surface course at the  The surface of the asphalt pavement undergoes deflection changes under the coupled dynamic load of the vehicle. At the same time, the road structure layer also undergoes bending and tensile stress. Under repeated vehicle loads, this bending and tensile stress will cause bending and tensile fatigue damage to the asphalt road structure. Figure 10 shows the time-history curve of the maximum tensile stress at each node of each structural layer on the road. It can be seen from the figure that the tensile stress on the road is sometimes positive and sometimes negative, but it is the maximum at the road surface course and the base course. Figure 11 shows the horizontal stress nephogram at the position of the maximum tensile stress at each node of the road structure layer. From this figure, it can be seen that the maximum tensile stress mainly occurs in the road surface course at the vehicle load and the combination part of the road surface course, bed course, and subgrade at the load, while the combination part of the road surface layer and the cushion is the largest. This shows that each structural layer of the road is relatively easy to damage in these parts. For the tensile stress of the road structure, it is easy to cause the cracking of road pavement on the surface course. However, there is a large tensile stress at the joint between the structural layers, which will increase the wear between the structural layers and directly lead to the damage of the road.

Shear Stress Analysis
According to the existing service conditions of the highway, the asphalt pavement is mainly damaged by shear stress. In asphalt pavement, it is easy to cause the shear failure of asphalt structural layers, especially under high temperatures. Shear failure is one of the most common problems with asphalt pavements.

Shear Stress Analysis
According to the existing service conditions of the highway, the asphalt pavement is mainly damaged by shear stress. In asphalt pavement, it is easy to cause the shear failure of asphalt structural layers, especially under high temperatures. Shear failure is one of the most common problems with asphalt pavements.

Shear Stress Analysis
According to the existing service conditions of the highway, the asphalt pavement is mainly damaged by shear stress. In asphalt pavement, it is easy to cause the shear failure of asphalt structural layers, especially under high temperatures. Shear failure is one of the most common problems with asphalt pavements. Figure 12 shows the time-history variation diagram of the horizontal shear stress of each node in the pavement structure layer under the action of the coupled dynamic load. It can be seen from the figure that the horizontal shear stress of the road surface is the most obvious, while the horizontal shear stress of the base and subgrade surface is relatively small, indicating that the horizontal sheer force of the asphalt road structure mainly has an effect on the asphalt road surface, which reminds us to pay special attention to the shear condition of the road surface in the design and construction process. At the same time, it can be seen that the horizontal shear stress can be positive or negative within the range of influence, and its value is approximately symmetric.
14, x FOR PEER REVIEW 17 of 24 Figure 12 shows the time-history variation diagram of the horizontal shear stress of each node in the pavement structure layer under the action of the coupled dynamic load. It can be seen from the figure that the horizontal shear stress of the road surface is the most obvious, while the horizontal shear stress of the base and subgrade surface is relatively small, indicating that the horizontal sheer force of the asphalt road structure mainly has an effect on the asphalt road surface, which reminds us to pay special attention to the shear condition of the road surface in the design and construction process. At the same time, it can be seen that the horizontal shear stress can be positive or negative within the range of influence, and its value is approximately symmetric.  Figure 13 shows the variation in the shear stress of nodes in the pavement structure layer at the midpoint of the road, along with the depth of the road section. According to the data in the figure, the proportion of cushion and roadbed with horizontal shear stress below 1 m is almost zero. The shear stress of asphalt pavement mainly occurs in the road surface and the base layer, the maximum shear stress is in the middle of the surface layer, and the value of shear stress is 43,858 Pa.  Figure 13 shows the variation in the shear stress of nodes in the pavement structure layer at the midpoint of the road, along with the depth of the road section. According to the data in the figure, the proportion of cushion and roadbed with horizontal shear stress below 1 m is almost zero. The shear stress of asphalt pavement mainly occurs in the road surface and the base layer, the maximum shear stress is in the middle of the surface layer, and the value of shear stress is 43,858 Pa.
layer at the midpoint of the road, along with the depth of the road section. the data in the figure, the proportion of cushion and roadbed with horizonta below 1 m is almost zero. The shear stress of asphalt pavement mainly occur surface and the base layer, the maximum shear stress is in the middle of the s and the value of shear stress is 43,858 Pa.   Figure 14 shows the variation of road vertical dynamic stress with depth when the vehicle travels to about 50 m. It can be seen from the figure that, like the horizontal shear stress, the value of the vertical dynamic stress in the number of road surface layers and the base layer is larger, and the dynamic stress in the middle of the surface layer is 119,373 Pa. The value of vertical dynamic stress in the road subgrade is much smaller, and the vertical dynamic stress in the middle of the subgrade is 5824 Pa. In the vertical direction of the road, the vertical dynamic stress of each layer of the road structure decreases rapidly. In contrast, the dynamic stress on the surface of the asphalt base course is 0.602 times the dynamic stress on the surface course, and when it reaches the bottom of the subbase, the dynamic stress decreases to 0.038 times the surface stress.

Vertical Dynamic Stress Analysis
Symmetry 2022, 14, x FOR PEER REVIEW Figure 14 shows the variation of road vertical dynamic stress with depth wh vehicle travels to about 50 m. It can be seen from the figure that, like the horizonta stress, the value of the vertical dynamic stress in the number of road surface lay the base layer is larger, and the dynamic stress in the middle of the surface layer is Pa. The value of vertical dynamic stress in the road subgrade is much smaller, a vertical dynamic stress in the middle of the subgrade is 5824 Pa. In the vertical di of the road, the vertical dynamic stress of each layer of the road structure decreas idly. In contrast, the dynamic stress on the surface of the asphalt base course is 0.60 the dynamic stress on the surface course, and when it reaches the bottom of the su the dynamic stress decreases to 0.038 times the surface stress.    Figure 15 shows the time-history variation diagram of the vertical dynamic stress of each layer. It can be seen from the diagram that the vertical dynamic stress of road is sometimes positive and sometimes negative, and the positive value shows tension while the negative value shows compression, and each value shows the maximum value when the asphalt road surface course. Figure 15 shows the time-history variation diagram of the vertical dynamic s each layer. It can be seen from the diagram that the vertical dynamic stress of sometimes positive and sometimes negative, and the positive value shows tensio the negative value shows compression, and each value shows the maximum valu the asphalt road surface course.   Figure 16 shows the distribution of vertical dynamic stress on asphalt roads. It can be seen from the figure that vertical dynamic stress is mainly concentrated in asphalt road surface and base, within a relatively small range of the whole road length. Due to the alternating action of vertical dynamic stress tension and compression, the damage of road surface is accelerated.  Figure 16 shows the distribution of vertical dynamic stress on asphalt roads. It can be seen from the figure that vertical dynamic stress is mainly concentrated in asphalt road surface and base, within a relatively small range of the whole road length. Due to the alternating action of vertical dynamic stress tension and compression, the damage of road surface is accelerated.

Stress Analysis of Pavement Considering Tire Starting Force and Braking Force
When the vehicle is running on the road, in addition to the dynamic load in the vertical direction of the vehicle, it also receives the starting force and braking force in the horizontal direction when the vehicle starts and brakes. In fact, when the vehicle is driving on the road, the horizontal starting and braking force are the reasons for the large shear stress on the road surface. China's urban road design code also requires us to consider the horizontal force applied to the road when the vehicle's braking starts.
The horizontal shear stress change mileage curve of the road under the action of ve-

Stress Analysis of Pavement Considering Tire Starting Force and Braking Force
When the vehicle is running on the road, in addition to the dynamic load in the vertical direction of the vehicle, it also receives the starting force and braking force in the horizontal direction when the vehicle starts and brakes. In fact, when the vehicle is driving on the road, the horizontal starting and braking force are the reasons for the large shear stress on the road surface. China's urban road design code also requires us to consider the horizontal force applied to the road when the vehicle's braking starts.
The horizontal shear stress change mileage curve of the road under the action of vehicle vertical coupling dynamic load, horizontal starting force, and braking force is shown in Figures 17 and 18. The horizontal shear stress of surface course and the base course is obvious during vehicle starting and braking. Compared with normal conditions, the horizontal shear stress of the road during starting and braking is larger. This is not only the main factor affecting vehicle starting and braking but also a factor affecting road service life.  In the process of vehicle starting and braking, the road surface course and the base course are stressed the most, and the shear stress of the base course is larger than that of the surface course, which is more obvious when the vehicle is braking. From the perspective of the vertical road displacement, the vehicle can affect the front of the road when it is started, and the vehicle can still affect the back road when it is braked, as shown in  In the process of vehicle starting and braking, the road surface course and the base course are stressed the most, and the shear stress of the base course is larger than that of the surface course, which is more obvious when the vehicle is braking. From the perspective of the vertical road displacement, the vehicle can affect the front of the road when it is started, and the vehicle can still affect the back road when it is braked, as shown in In the process of vehicle starting and braking, the road surface course and the base course are stressed the most, and the shear stress of the base course is larger than that of the surface course, which is more obvious when the vehicle is braking. From the perspective of the vertical road displacement, the vehicle can affect the front of the road when it is started, and the vehicle can still affect the back road when it is braked, as shown in Figures 19 and 20.

Conclusions
The heavy-duty vehicle road coupling model not only adopts the vibration coupling relationship under the condition of displacement compatibility between vehicle and road but also tries to adopt the vibration coupling relationship under the condition of contact force between vehicle and road. Further research is conducted on the coupling relationship, which provides a new idea for the future development of a vehicle-road coupling vibration model. This paper studies the vehicle-road coupling dynamic analysis model and vehicle-road interaction analysis under heavy load conditions, which will play a sig-

Conclusions
The heavy-duty vehicle road coupling model not only adopts the vibration coupling relationship under the condition of displacement compatibility between vehicle and road but also tries to adopt the vibration coupling relationship under the condition of contact force between vehicle and road. Further research is conducted on the coupling relationship, which provides a new idea for the future development of a vehicle-road coupling vibration model. This paper studies the vehicle-road coupling dynamic analysis model and vehicle-road interaction analysis under heavy load conditions, which will play a sig-

Conclusions
The heavy-duty vehicle road coupling model not only adopts the vibration coupling relationship under the condition of displacement compatibility between vehicle and road but also tries to adopt the vibration coupling relationship under the condition of contact force between vehicle and road. Further research is conducted on the coupling relationship, which provides a new idea for the future development of a vehicle-road coupling vibration model. This paper studies the vehicle-road coupling dynamic analysis model and vehicle-road interaction analysis under heavy load conditions, which will play a significant role in the whole design method and theory of vehicles and roads and will also have very important application value to future engineering practice.
In this paper, the subgrade is regarded as part of the vehicle road dynamic system, and the vehicle road coupling dynamic model under heavy load is established. The results show that the equation can better reflect the dynamic performance of a vehicle-road system; Additionally, it better reflects the coupling vibration relationship and interaction between the vehicle and the road. The model will provide a theoretical basis for a more accurate analysis of the interaction between heavy-duty vehicles and roads and then provide a theoretical basis for important parameters in heavy-duty vehicle design and road durability design. The coupling dynamic equation established in this paper can analyze the instantaneous dynamic performance of the coupling system by the finite element method.
The two-dimensional finite element model is used in the analysis of road coupled dynamic load in this paper, and the three-dimensional finite element model can be used to simulate the spatial dynamic response of the road. This will be the focus of future research.

Symbols of Vehicle Parameters Value
Rear-wheel damping coefficient (C t2 /kN·m −1 ·s −1 ) 24 Front-wheel damping coefficient (C t1 /kN·m −1 ·s −1 ) 6 The distance from the center of gravity of the car to the rear axle (λ 2 /m) 1.5 The distance from the center of gravity of the car to the front axle (λ 1 /m) 2.5 The distance from the center of gravity of the car to the right wheel (λ 4 /m) 0.9 The distance from the center of gravity of the car to the left wheel (λ 3 /m) 0.9