Dynamic Analysis and Seat Selection of Bus Driving Comfort under Different Road Conditions

Featured Application: By giving speciﬁc road excitation, the differential equations and analytical solutions derived in this paper can be used to quickly estimate the vertical displacement and acceleration of the bus passing through this section of the road. The theoretical results and formulae in this paper can provide some guidelines for bus design and bus driving. Abstract: The comfort of a bus running on different road conditions is a matter of public concern. In this paper, the differential equations of motion are established for a bus running on different road conditions and the whole driving process is mechanically analyzed. Firstly, the bump degree at different positions is quantitatively analyzed and it is found that the rear row is bumpier on different roads. Then, the relationship between the speed of the bus and the vertical displacement and acceleration is quantitatively described. Regardless of the speed, a similar displacement and acceleration will be eventually achieved, but the speed is higher, and the duration of maximum displacement and acceleration is longer. When the speed is 8 m/s, resonance occurs on the bus during road condition II . Finally, the change in vertical displacement and acceleration under the action of different spring stiffness coefﬁcient ratios of the front and rear wheels is quantitatively analyzed. High stiffness ratios mean less displacement and acceleration. By establishing an actual excitation road surface, the differential equations and analytical solutions in this paper can be used to roughly analyze the mechanical response of a traveling bus. These results can provide some guidance for the design and driving of buses.


Introduction
According to the Statistical Bulletin on the Development of the Transportation Industry in 2021 issued by the Chinese government in May 2022, the total mileage of roads in China reached 5.2807 million kilometers, and the number of urban buses and train reached 709,400.Through the study in Malaysia, Kuys et al. [1] found that all countries in the world are facing the problems of increasing urban roads, increasing numbers of buses, and increasing road congestion.With the continuous increase in highway mileage, the comfort of buses has gradually attracted attention.Previous studies found that bus comfort feelings mainly come from noise, temperature, vibration, and acceleration, as listed in Table 1.Among these four factors, vibration and acceleration are the most popular two.
Vertical vibration is the most common form of vibration experienced by passengers in a bus.In 1962, Coermann [2] abstracted the human body into a mechanical impedance structure with multiple degrees of freedom and discussed the biomechanical responses of human body sitting posture to vertical vibration.Gao et al. [3] established a two-degreeof-freedom vertical vibration model of the human body and discussed the influence of body parameters on driver comfort.Bazil et al. [4] discussed the influence of the angle between the fixed backrest and seat on human comfort under different vibration situations through experiments.Their results show that the angle cannot be expressed by a single frequency weighting, rather a multi-frequency weighting function with strong nonlinearity.Xie et al. [5] tested the vibration of a multi-floor three-dimensional passenger station onsite and found that the passenger bus causes floor vertical vibration and thus causes the passenger discomfort in the passenger station.Faster bus speed induces stronger discomfort.The passengers on a train have similar discomfort when the train speed is fast.Zhou et al. [6] verified experimentally that the frequency dependence of the vibration discomfort depends on the acceleration and force at the seat interface with the body, which further verifies that the human body is less affected on the dynamic force than the acceleration.Even a small vertical acceleration can cause discomfort.Vertical acceleration is the main influencing factor in passenger comfort.affects the comfort of passengers and poses a certain health risk to passengers.In addition to the vertical vibration caused by the bus itself, the uneven ground affects the vertical acceleration.The vibration amplitude is larger, the acceleration is stronger, and the discomfort is stronger.Zhang et al. [36] established the random excitation of the road surface by using the filtered white noise method and analyzed the vertical acceleration caused by the uneven road surface through simulation.They evaluated the ride comfort.Tang et al. [37] used the neural network algorithm to analyze the random excitation of the road surface, thus establishing an annoyance rate model to quantitatively evaluate the comfort of passengers.Bogsjö et al. [38] proposed three new road profile models to quantitatively study passenger comfort and vehicle damage.Chen et al. [39] used experiments to prove the impact of deceleration strips on human comfort.Their results show that the characteristic vibration dose value can be used as a method to evaluate passenger comfort.Agostinacchio et al. [40] proved that an irregular road surface generates a considerable instantaneous dynamic load, which does not contribute to the total response but does affect the human comfort of riding.For the bumpy analysis of the front and rear seats in a bus, Wu and Hao [41] analyzed and demonstrated that the bumpier rear seat is due to the location of the bus gravity and the wheel from the perspective of the differential equation of fixed axle rotation.Song [42] performed a qualitative analysis for the safety requirements of automobiles and reached a similar conclusion.Their common shortcoming is the lack of vehicle dynamics analysis, and their conclusion that the rear row of the bus is bumpier is worthy of further investigation.
The bus speed is also a major factor affecting the comfort of passengers, especially on uneven roads.Sekulić et al. [19] took the intercity bus IK-301 as an oscillation model with ten degrees of freedom, used the power spectral density of asphalt concrete pavement roughness to generate a total road excitation, and analyzed the impact of vibration on passenger comfort.Bae et al. [33] discussed the impact of speed on passenger comfort, and then proposed a time-optimal speed planning method, which can improve passenger comfort by reasonably planning the driving path.Barabino et al. [43] showed that speed has a great impact on passenger comfort.Nguyen et al. [44] established three relatively simple mathematical models, focusing on the actual driving behavior of urban buses, and studied the impact of acceleration and speed changes on human comfort.In summary, the speed of a bus has a certain impact on the comfort of passengers.
In order to reduce the net displacement of the rear wheel, engineers often assign different values to the spring stiffness of the front wheel and the rear wheel.Zhang et al. [45] used MATLAB to carry out the bus driving matrix and set the front and rear wheel stiffness as 155,400 N/m and 250,000 N/m, respectively, in the calculation.The stability of the bus was analyzed through the derivation of a rigid body dynamics equation.Kong et al. [46] evaluated the quality of the suspension system of the bus through experiments.In the experiment, P-2 parabolic spring was used.The stiffness of the rear wheel was 1.28 times that of the front wheel on average.Long et al. [47] used simulink toolbox to study the influence of vehicle suspension stiffness on bus ride comfort.Their study found that the stiffness of the rear wheel was positively correlated with passenger comfort.Olmeda et al. [48] also explained the influence of the front and rear wheel stiffness coefficient on bus comfort when carrying out a bus lateral dynamics simulation.However, in the equation of its calculation, there are many state variables that have interference and are treated as noise when processing data.It can be seen that the spring stiffness coefficients of the front and rear wheels of the passenger car will also affect the comfort of passengers.
Previous studies focused on vibration and acceleration analysis, but they were too focused on a vehicle analysis of the response of the whole vehicle or the interaction between the vehicle and the uneven ground.Few studies have focused on the coupling analysis of vehicle parts and uneven pavement, and the study of vehicle parts is more conducive to bus design.There are still few studies on the impact of bus speed and the stiffness of front and rear wheel springs on passenger comfort.Therefore, this paper will start with the above problems and conduct quantitative research on the whole process of bus driving based on mathematical and physical methods through self-defining the uneven road surface.This paper is organized as follows: In Section 2, the force of the front and rear wheels of the bus when driving on a flat road is analyzed.In Section 3, the differential equation of the whole process of bus driving is established based on the user-defined uneven road surface and the local displacement characteristics of the front and rear seats of the bus are analyzed.In Section 4, the mechanical responses of the bus under three different road conditions and different driving speeds are analyzed by using the differential equation of motion.In Section 5, the influence of front and rear wheel spring stiffness on passenger car comfort is analyzed.Conclusions are drawn in Section 6.

Mechanical Analysis of Bus Running on Horizontal Road
For a bus running on uneven roads, the displacement of the rear seats is significantly greater than that of the front seats.This phenomenon can be clearly observed with the naked eye.Most buses are driven by rear wheels.This means that the engine drives the rear wheels to rotate and move forward depending on the friction between the rear wheels and the ground.In order to make the rear wheel of bus have greater friction with the ground, it is necessary to move the center of gravity of the bus as far back as possible, so as to ensure that the support force of the ground to the rear wheel is greater than that of the front wheel.The center of gravity of the actual running bus is often placed between the front and rear wheels near 1/3 of the rear wheels.This can ensure that the rear wheels bear 2/3 of the vehicle weight and the front wheels bear 1/3 of the vehicle weight.This section analyzes the force of a bus when the bus is normally running on a horizontal road.All bodies of the bus are taken as the rigid bodies to establish the mechanical model bus as shown in Figure 1.The bus is assumed to move in a horizontal straight line with acceleration a (m/s 2 ), where the total mass of the bus is m (kg), the height of the center of mass G from the ground is h (m), and the distance from the front and rear axles to the vertical line of the center of mass is equal to b (m) and c (m), respectively.
ground, it is necessary to move the center of gravity of the bus as far back as possibl as to ensure that the support force of the ground to the rear wheel is greater than th the front wheel.The center of gravity of the actual running bus is often placed betw the front and rear wheels near 1/3 of the rear wheels.This can ensure that the rear wh bear 2/3 of the vehicle weight and the front wheels bear 1/3 of the vehicle weight.
This section analyzes the force of a bus when the bus is normally running on a h zontal road.All bodies of the bus are taken as the rigid bodies to establish the mecha model bus as shown in Figure 1.The bus is assumed to move in a horizontal straight with acceleration a (m/s 2 ), where the total mass of the bus is m (kg), the height of the ce of mass G from the ground is h (m), and the distance from the front and rear axles to vertical line of the center of mass is equal to b (m) and c (m), respectively.Using D'Alembert's principle [49], the added virtual inertial force is g F ma = The bus gravity is

P mg =
where Fg is the inertial force, N; P is the gravity, N; g is the acceleration of gravity, m Moment analysis at point A of the front wheel obtains Similarly, the moment analysis at point B of the front wheel is expressed as where FNA (N) is the ground support force on the front wheels and FNB (N) is the gro support force on the rear wheels.
According to Equations ( 1)-( 4), one can obtain Taking the Yutong ZK6105HNGS1 bus [50] as an example, the whole-body ma the bus is 12,000 kg, the body length is 10 m, the height of the gravity center is 1.5 m Using D'Alembert's principle [49], the added virtual inertial force is The bus gravity is where F g is the inertial force, N; P is the gravity, N; g is the acceleration of gravity, m/s 2 .Moment analysis at point A of the front wheel obtains Similarly, the moment analysis at point B of the front wheel is expressed as where F NA (N) is the ground support force on the front wheels and F NB (N) is the ground support force on the rear wheels.
According to Equations ( 1)-( 4), one can obtain Taking the Yutong ZK6105HNGS1 bus [50] as an example, the whole-body mass of the bus is 12,000 kg, the body length is 10 m, the height of the gravity center is 1.5 m, the wheelbase is 9 m, the front wheel is 1/3 from the center of gravity, and the rear wheel is 2/3 from the center of gravity.The relationship between the support force and the acceleration of the front and rear wheels on the ground is shown in Figure 2, which is based on Equation (5).The normal maximum acceleration of the bus is 8 m/s 2 [51], and the braking acceleration is −10 m/s 2 [52].The critical takeoff means that the support force F N is zero.Figure 2 shows that the critical takeoff condition is met with difficulty when the bus is running on a horizontal road.When the bus acceleration a > −10 m/s 2 , the supporting force of the ground to the rear wheel is always greater than that of the front wheel.According to the isolation method, it is not difficult to infer that the passengers in the rear row bear greater acceleration and greater acceleration changes.
eration of the front and rear wheels on the ground is shown in Figure 2, w Equation (5).The normal maximum acceleration of the bus is 8 m/s 2 [51] acceleration is −10 m/s 2 [52].The critical takeoff means that the suppor Figure 2 shows that the critical takeoff condition is met with difficulty running on a horizontal road.When the bus acceleration a > −10 m/s 2 , the of the ground to the rear wheel is always greater than that of the front w to the isolation method, it is not difficult to infer that the passengers in greater acceleration and greater acceleration changes.

Mechanical Analysis of a Bus Running on an Uneven Road
The smoothness of the road surface seriously affects the safety of bu tuation of ISO-B pavement is 0.1 m [54], and some vertical acceleratio which is more serious in rural regions [55].Zhang et al. [56] used the po sity formula of road roughness in the Chinese national standard GB/T703 the road roughness.They obtained the maximum difference of 0.2 m an of 2.5 m/s 2 .When the human body feels a vertical acceleration above 0.5 feel obvious discomfort [36,57].In order to facilitate mechanical analys sonably simplified as a plane problem and the displacement of the bus road surface is assumed to be y = hcos (2πx/s) and x = vt.The carriage a are connected by the spring coefficient k.It is worth noting that the effec be not considered in this paper, which means that energy does not decay of the bus.All conclusions of the analysis are based on the bus limit state cive to a better explanation of the problem.The vertical displacement o the center of gravity, and the rear wheel are defined as uA, uG, and uB, res difference between two adjacent wave peaks with a convex periodic bou rameters are consistent with Figure 1.Three road conditions will be an bus passes through the uneven road section.The initial state includes two bears the bus gravity and produces a certain displacement.Second, the tion of the bus has just reached the starting point of the uneven road secti model is shown in Figure 3.In this section, if there is no special explanat ment is considered as vertical displacement, and the acceleration is al vertical acceleration.Downhill means that the bus is in road condition I (0

Mechanical Analysis of a Bus Running on an Uneven Road
The smoothness of the road surface seriously affects the safety of buses [53].The fluctuation of ISO-B pavement is 0.1 m [54], and some vertical acceleration reaches 2 m/s 2 , which is more serious in rural regions [55].Zhang et al. [56] used the power spectral density formula of road roughness in the Chinese national standard GB/T7031-2005 to inverse the road roughness.They obtained the maximum difference of 0.2 m and an acceleration of 2.5 m/s 2 .When the human body feels a vertical acceleration above 0.5 m/s 2 , it begins to feel obvious discomfort [36,57].In order to facilitate mechanical analysis, the bus is reasonably simplified as a plane problem and the displacement of the bus is analyzed.The road surface is assumed to be y = hcos (2πx/s) and x = vt.The carriage and tire of the bus are connected by the spring coefficient k.It is worth noting that the effect of dampers will be not considered in this paper, which means that energy does not decay with the progress of the bus.All conclusions of the analysis are based on the bus limit state, which is conducive to a better explanation of the problem.The vertical displacement of the front wheel, the center of gravity, and the rear wheel are defined as u A , u G , and u B , respectively.s is the difference between two adjacent wave peaks with a convex periodic boundary.Other parameters are consistent with Figure 1.Three road conditions will be analyzed when the bus passes through the uneven road section.The initial state includes two: First, the spring bears the bus gravity and produces a certain displacement.Second, the front wheel position of the bus has just reached the starting point of the uneven road section.The analytical model is shown in Figure 3.In this section, if there is no special explanation, the displacement is considered as vertical displacement, and the acceleration is also considered as vertical acceleration.Downhill means that the bus is in road condition I (0 ≤ x ≤ (b + c)) and climbing means that the bus is in road condition III (2ns

Road Condition I
The bus moves forward.The front wheel drives into the uneven road section, but the rear wheel is still in the flat road section; thus, the bus travel distance is 0 ≤ x ≤ (b + c).The vertical displacement increment of two points A and B of the bus is denoted as δA and δB.Since the front wheel has driven into the pit, there is a geometric relationship: where yA and yB are the vertical coordinates of the front and rear wheels of the bus.

cos
All bodies of the bus are regarded as rigid bodies.The vertical displacement coordinate uG at the center of gravity can be solved by using a similar triangle method: The upper part of the bus is analyzed by the free-body method.The force analysis is shown in Figure 4:

The equation of motion is
The equilibrium equation for calculating the moment of force at point G is Combining Equations ( 8)-( 10) obtains the differential equation of motion: A simplified mechanical analysis of uneven pavement.

Road Condition I
The bus moves forward.The front wheel drives into the uneven road section, but the rear wheel is still in the flat road section; thus, the bus travel distance is 0 ≤ x ≤ (b + c).The vertical displacement increment of two points A and B of the bus is denoted as δ A and δ B .Since the front wheel has driven into the pit, there is a geometric relationship: where y A and y B are the vertical coordinates of the front and rear wheels of the bus.
All bodies of the bus are regarded as rigid bodies.The vertical displacement coordinate u G at the center of gravity can be solved by using a similar triangle method: The upper part of the bus is analyzed by the free-body method.The force analysis is shown in Figure 4:

Road Condition I
The bus moves forward.The front wheel drives into the uneven road section, but the rear wheel is still in the flat road section; thus, the bus travel distance is 0 ≤ x ≤ (b + c).The vertical displacement increment of two points A and B of the bus is denoted as δA and δB.Since the front wheel has driven into the pit, there is a geometric relationship: where yA and yB are the vertical coordinates of the front and rear wheels of the bus.

cos
All bodies of the bus are regarded as rigid bodies.The vertical displacement coordinate uG at the center of gravity can be solved by using a similar triangle method: The upper part of the bus is analyzed by the free-body method.The force analysis is shown in Figure 4:

The equation of motion is
The equilibrium equation for calculating the moment of force at point G is Combining Equations ( 8)-( 10) obtains the differential equation of motion: The equation of motion is The equilibrium equation for calculating the moment of force at point G is Combining Equations ( 8)-( 10) obtains the differential equation of motion: Its general solution is obtained with where C 1 and C 2 are the coefficients to be determined by boundary conditions.Combining Equations ( 6) and ( 10) obtains u B : When x = 0 and the front wheel moves to point O, the boundary condition is Therefore, the coefficients in Equation ( 12) can be determined as The vertical displacement of bus points A and B is then expressed as The last subscript 1 in Equations ( 15) and ( 16) denotes road condition I.This paper uses MATLAB to solve Equation (16). Figure 5 shows the result of only net displacement, which means the change in displacement and acceleration of the front and rear wheels in the first stage.The maximum net displacement difference is 5.93 cm at the front wheel and 11.86 cm at the rear wheel.When the front wheel drives on the uneven road and the rear wheel is still on the flat road, the displacement near the rear wheel position is greater.The net vertical acceleration of the front and rear wheels is more than 3 m/s 2 , the peak net acceleration of the rear wheels is 5.31 m/s 2 , and the peak acceleration of the front wheels is 3.32 m/s 2 .The rear wheel acceleration is generally higher than the front wheel acceleration.Thus, the rear wheels are bumpier in the road condition I.

Road Condition II
Both the front and rear wheels drive into an uneven road section, the bus position is (b + c) ≤ x ≤ 2ns (n ∈ N*).N* is a positive integer.The general differential equation of bus motion is still Equation (11).Its boundary conditions at this stage become the following: when x = b + c, the motion state of road condition II should satisfy Equation (16) in road condition I.The continuity of bus displacement is The subscript 2 represents road condition II.

Road Condition II
Both the front and rear wheels drive into an uneven road section, the bus position is (b + c) ≤ x ≤ 2ns (n ∈ N*).N* is a positive integer.The general differential equation of bus motion is still Equation (11).Its boundary conditions at this stage become the following: when x = b + c, the motion state of road condition II should satisfy Equation ( 16) in road condition I.The continuity of bus displacement is The subscript 2 represents road condition II.When x = 2ns (n ∈ N*), the front and rear wheels are all located on uneven road sections.The instantaneous state of the rear wheel entering the uneven road section is The status of the front and rear wheels in this process is ( ) The coefficients obtained by Equations ( 17)- (19) and Equation ( 12) are  When x = 2ns (n ∈ N*), the front and rear wheels are all located on uneven road sections.The instantaneous state of the rear wheel entering the uneven road section is The status of the front and rear wheels in this process is The coefficients obtained by Equations ( 17)- (19) and Equation (12) are Then, the vertical displacements at bus points A and B are Equation ( 21) is still solved by MATLAB and Figure 6 is the net displacement.At this time, the maximum net displacement difference is 3.05 cm at the front wheel and 3.1 cm at the rear wheel.Their displacement difference is not obvious in Figure 6a.This may be due to the fact that the uneven road surface is larger than the body span in a single cycle.The vertical positions of the front wheel and the rear wheel are basically in synchronization.The net vertical acceleration of the front and rear wheels is more than 0.1 m/s 2 , the net peak acceleration is 0.14 m/s 2 at the rear wheels and 0.1 m/s 2 at the front wheels.The acceleration is still higher at the rear wheels than at the front wheels.This result can still be used as an argument for the rear row being bumpier, although the difference is small or not salient.
Equation ( 21) is still solved by MATLAB and Figure 6 is the net displacement.At this time, the maximum net displacement difference is 3.05 cm at the front wheel and 3.1 cm at the rear wheel.Their displacement difference is not obvious in Figure 6a.This may be due to the fact that the uneven road surface is larger than the body span in a single cycle.The vertical positions of the front wheel and the rear wheel are basically in synchronization.The net vertical acceleration of the front and rear wheels is more than 0.1 m/s 2 , the net peak acceleration is 0.14 m/s 2 at the rear wheels and 0.1 m/s 2 at the front wheels.The acceleration is still higher at the rear wheels than at the front wheels.This result can still be used as an argument for the rear row being bumpier, although the difference is small or not salient.

Road Condition III
The front wheel drives into the flat section, but the rear wheel is still in the uneven section.At this time, the bus position is 2ns

Road Condition III
The front wheel drives into the flat section, but the rear wheel is still in the uneven section.At this time, the bus position is 2ns ≤ x ≤ 2ns + (b + c) (n ∈ N*).The differential equation of motion in the road condition III is still Equation (11).When x = 2ns (n ∈ N*), the motion state of road condition III should satisfy Equation ( 21) when x = 2ns (n ∈ N*) in road condition II.The continuity of the bus displacement satisfies The subscript 3 represents road condition III.When x = 2ns + (b + c) (n ∈ N*), the front wheel has just entered the flat section and the rear wheel is still in the uneven section.The instantaneous state of the front wheel entering the flat road section is The statuses of the front and rear wheels in this process are The coefficients can be determined by Equations ( 22)- (24) and Equation (12) as The vertical displacements at bus points A and B are Figure 7 shows the net displacement and vertical acceleration calculated by Equation (26).The maximum net displacement difference is 5.95 cm at the front wheel and 12.1 cm at the rear wheel.When the front wheel drives on a flat road and the rear wheel is still on a bumpy road, the displacement near the rear wheel position is still large.The net vertical acceleration of the front and rear wheels is more than 3.2 m/s 2 , the peak net acceleration of the rear wheels is 5.42 m/s 2 , and the peak acceleration of the front wheels is 3.28 m/s 2 .The rear wheel acceleration is generally higher than the front wheels.This means that the rear wheels are still bumpier at the road condition II.

Mechanical Responses of Bus in Driving under the Coupling Effect of Speed and Uneven Road Surface
With the increasing number of private cars, urban traffic becomes more and more crowded, and the speed of buses is getting slower and slower.Barabino et al. [43] conducted an experimental study on the bus driving conditions in the city of Cagliari, located in Sardinia (Italy).The original data of the bus driving are listed in Table 2.The speed of buses in rural regions is much faster than that in urban areas, so the four groups of buses selected in this section run at the speed of 4 m/s, 6 m/s, 8 m/s, and 10 m/s, respectively.Their mechanical responses when passing through the uneven road surface will be discussed below.
Table 2. Location-at-time raw data sample gathered by the app Torque on Tue Jul 2018 [43].Note: alat means the acceleration component in latitude and along means the acceleration component in longitude.

Road Condition I
The sketch of road conditions for the bus is shown in Figure 3.The front wheel drives into the uneven road section, but the rear wheel is still in the flat road section.The bus travel distance is 0 ≤ x ≤ (b + c).The net displacement of the front and rear wheels of the bus at this point is presented in Figure 8.When the speed is relatively high, the maximum net displacement be reached quicker.This means that the net displacement of the bus is smaller when the bus speed is slower.Thus, a greater bus speed will generate bigger net displacement.Therefore, a bus running at low speed has smaller displacement and higher passenger comfort.Figure 9 shows the acceleration change of the bus at different speeds.

Mechanical Responses of Bus in Driving under the Coupling Effect of Speed and Uneven Road Surface
With the increasing number of private cars, urban traffic becomes more and more crowded, and the speed of buses is getting slower and slower.Barabino et al. [43] conducted an experimental study on the bus driving conditions in the city of Cagliari, located in Sardinia (Italy).The original data of the bus driving are listed in Table 2.The speed of buses in rural regions is much faster than that in urban areas, so the four groups of buses selected in this section run at the speed of 4 m/s, 6 m/s, 8 m/s, and 10 m/s, respectively.Their mechanical responses when passing through the uneven road surface will be discussed below.

Road Condition I
The sketch of road conditions for the bus is shown in Figure 3.The front wheel drives into the uneven road section, but the rear wheel is still in the flat road section.The bus travel distance is 0 ≤ x ≤ (b + c).The net displacement of the front and rear wheels of the bus at this point is presented in Figure 8.When the speed is relatively high, the maximum net displacement be reached quicker.This means that the net displacement of the bus is smaller when the bus speed is slower.Thus, a greater bus speed will generate bigger net displacement.Therefore, a bus running at low speed has smaller displacement and higher passenger comfort.Figure 9 shows the acceleration change of the bus at different speeds.The acceleration also shows that when the bus speed is low, it takes longer to reach the maximum net acceleration of the vehicle, which means that the acceleration changes less at low speed, and the passenger comfort is higher.Table 3 summarizes the maximum net displacement and acceleration at the front and rear wheels of the bus at different speeds.The net displacement basically does not change with speed and is mainly induced by road incentives.When the speed is low, both  Table 3 summarizes the maximum net displacement and acceleration at the front and rear wheels of the bus at different speeds.The net displacement basically does not change with speed and is mainly induced by road incentives.When the speed is low, both Table 3 summarizes the maximum net displacement and acceleration at the front and rear wheels of the bus at different speeds.The net displacement basically does not change with speed and is mainly induced by road incentives.When the speed is low, both maximum net displacement and acceleration are still small, although the gap is small.However, it can be clearly observed from Figure 9 that the faster the speed is, the faster the bus reaches its maximum acceleration.This indicates that the higher the bus speed is, the longer the higher acceleration is maintained.maximum net displacement and acceleration are still small, although the gap is small.However, it can be clearly observed from Figure 9 that the faster the speed is, the faster the bus reaches its maximum acceleration.This indicates that the higher the bus speed is, the longer the higher acceleration is maintained.

Road Condition II
Both the front and rear wheels drive into an uneven road section.At this time, the bus position is (b + c) ≤ x ≤ 2ns (n ∈ N*).N* is a positive integer.The net displacement and the net acceleration at different speeds are presented in Figures 10 and 11, respectively.(   Figure 10 shows the displacement of the bus in the middle section.It clearly shows that there is only one wave peak in (a), one and a half wave peaks in (b), and two wave peaks in (c), but the results are not obvious, while there are two wave peaks in (d).This means that faster speed will generate a higher frequency of large displacement.Figure 10c has obvious fluctuations, which are different from the other three figures.This mainly reflects the large displacement in a short time.This may be because the bus resonates with the ground when v = 8 m/s, resulting in the continuous large displacement of the whole bus.
The change in acceleration still shows that the comfort of the bus will be further reduced when the speed is fast.The change in v = 8 m/s is the most dramatic.At this time, the resonance situation does occur when the bus passes through the section (b + c) ≤ x ≤ 2ns (n ∈ N*), the comfort of passengers is the worst, and the overall displacement of the bus is the largest.Table 4 records the maximum net displacement and acceleration corresponding to different speeds during the whole bus driving in road condition II.The overall acceleration changes uniformly and the displacement changes slightly when the bus passes through this section.When the bus speed is 8 m/s, both net displacement and net acceleration reach a maximum value of a sharp increase.This indicates that the resonance phenomenon between the bus and the road occurs at this time, and the passenger comfort is the lowest.

Road Condition III
The front wheel drives into the flat section, but the rear wheel is still in the uneven section.At this time, the bus position is 2ns Figure 10 shows the displacement of the bus in the middle section.It clearly shows that there is only one wave peak in (a), one and a half wave peaks in (b), and two wave peaks in (c), but the results are not obvious, while there are two wave peaks in (d).This means that faster speed will generate a higher frequency of large displacement.Figure 10c has obvious fluctuations, which are different from the other three figures.This mainly reflects the large displacement in a short time.This may be because the bus resonates with the ground when v = 8 m/s, resulting in the continuous large displacement of the whole bus.
The change in acceleration still shows that the comfort of the bus will be further reduced when the speed is fast.The change in v = 8 m/s is the most dramatic.At this time, the resonance situation does occur when the bus passes through the section (b + c) ≤ x ≤ 2ns (n ∈ N*), the comfort of passengers is the worst, and the overall displacement of the bus is the largest.Table 4 records the maximum net displacement and acceleration corresponding to different speeds during the whole bus driving in road condition II.The overall acceleration changes uniformly and the displacement changes slightly when the bus passes through this section.When the bus speed is 8 m/s, both net displacement and net acceleration reach a maximum value of a sharp increase.This indicates that the resonance phenomenon between the bus and the road occurs at this time, and the passenger comfort is the lowest.

Road Condition III
The front wheel drives into the flat section, but the rear wheel is still in the uneven section.At this time, the bus position is 2ns ≤ x ≤ 2ns + (b + c) (n ∈ N*).The net displacement and net acceleration at different speeds are presented in Figures 12 and 13, respectively.The calculation results in Figure 12 show that the bus has a small net displacement at the beginning when the bus is climbing at a speed of v = 6 m/s.This happens when the net displacements of the front and rear wheels are basically identical, and the spring is near the origin when the second section of the road is approaching the end.and net acceleration at different speeds are presented in Figures 12 and 13, respectively.The calculation results in Figure 12 show that the bus has a small net displacement at the beginning when the bus is climbing at a speed of v = 6 m/s.This happens when the net displacements of the front and rear wheels are basically identical, and the spring is near the origin when the second section of the road is approaching the end.Figures 12c,d and 13 show that the springs of the front and rear wheels of the bus are in a tensile state before climbing, so there will be a very small displacement for a short period of time, then the spring will be compressed rapidly under the impact, and then repeated tensile compression will occur.Therefore, this phenomenon still means that

Figures 12c,d and 13
show that the springs of the front and rear wheels of the bus are in a tensile state before climbing, so there will be a very small displacement for a short period of time, then the spring will be compressed rapidly under the impact, and then repeated tensile compression will occur.Therefore, this phenomenon still means that when the bus is climbing a slope, slowing down is still the wisest choice, which will slow down the large instantaneous impact of the spring.A lower speed means that the passenger's comfort is improved.The acceleration in Figure 13 shows that at different speeds, the acceleration undergoes a process of increasing from large to small, and then increasing.However, faster speeds will result in more wave crests and a more intense frequency of acceleration changes.
Table 5 records the maximum net displacement and acceleration corresponding to different speeds during the whole bus driving in road condition III.The maximum net displacement and acceleration of the bus at different speeds are basically the same when the bus is climbing.The reason is that the center of gravity of the bus is 2/3 away from the front wheel, and the excitation generated during climbing only has a large impact on the front wheel.The rear wheel should be kept in the motion state of road condition II at first, and then changed.The change is more drastic for greater speeds (see Figures 12c,d and 13).

Mechanical Responses of Bus in Driving under the Coupling of Spring Stiffness and Uneven Pavement
For the completeness of the theoretical analysis, this section will discuss the mechanical response of the bus under three different coefficient ratios.The coefficient ratio of these three groups will take 1 as the boundary, and one group will be selected before and after each.The bus speed is given as a constant 6 m/s, and the grouping of stiffness coefficients is listed in Table 6.The low ratio means that the stiffness coefficient is larger at the front wheel than at the rear wheel.The medium ratio means both the front and rear wheels have the same stiffness coefficient.The high ratio means that the stiffness coefficient is smaller at the front wheel than at the rear wheel.

Road Condition I
The front wheel drives into the uneven road section, but the rear wheel is still in the flat road section, and the bus travel distance is 0 ≤ x ≤ (b + c).In the case of road condition I, the displacement and acceleration changes of the front and rear wheels with different stiffness coefficient ratios are shown in Figure 14.The displacement changes little under different stiffness coefficient ratios.When the stiffness ratio is medium, the net displacement of the front and rear wheels is the minimum.The acceleration in Figure 14 obviously shows that the envelope of the black line can include all the curves.This means that the bus with low spring stiffness will cause a large acceleration change when the bus goes downhill.
the displacement and acceleration changes of the front and rear wheels with different stiffness coefficient ratios are shown in Figure 14.The displacement changes little under different stiffness coefficient ratios.When the stiffness ratio is medium, the net displacement of the front and rear wheels is the minimum.The acceleration in Figure 14 obviously shows that the envelope of the black line can include all the curves.This means that the bus with low spring stiffness will cause a large acceleration change when the bus goes downhill.

Road Condition II
Both the front and rear wheels drive into an uneven road section, the bus position is (b + c) ≤ x ≤ 2ns (n ∈ N*).N* is a positive integer.In the case of road condition II, the displacement and acceleration of the front and rear wheels with different spring stiffness coefficient ratios are shown in Figure 15.

Road Condition II
Both the front and rear wheels drive into an uneven road section, the bus position is (b + c) ≤ x ≤ 2ns (n ∈ N*).N* is a positive integer.In the case of road condition II, the displacement and acceleration of the front and rear wheels with different spring stiffness coefficient ratios are shown in Figure 15.It is observed that the displacement is not obvious under different stiffness ratios.This may be that the slope is too small compared with the bus, and the front and rear wheels of the bus basically go up and down the slope synchronously.In the front wheel acceleration analysis in Figure 15, the low stiffness ratio means a large acceleration change, It is observed that the displacement is not obvious under different stiffness ratios.This may be that the slope is too small compared with the bus, and the front and rear wheels of the bus basically go up and down the slope synchronously.In the front wheel acceleration analysis in Figure 15, the low stiffness ratio means a large acceleration change, while the high stiffness ratio means a small change range.On the rear wheel acceleration, the change of the rear wheel of the bus is relatively small when the stiffness ratio is high.In other cases, the acceleration change is large.

Road Condition III
The front wheel drives into the flat section, but the rear wheel is still in the uneven section.At this time, the bus position is 2ns ≤ x ≤ 2ns + (b + c) (n ∈ N*).In the case of road condition III, the changes in displacement and acceleration at the front and rear wheels with different stiffness coefficient ratios are shown in Figure 16.The displacement changes are similar regardless of the stiffness coefficient ratio.The envelope of the black line can include all the curves about acceleration analysis in Figure 16.This implies that the curve with lower stiffness will lead to greater acceleration change when the bus goes uphill.

Conclusions
In this paper, from the perspective of rigid body dynamics, a motion differential equation suitable for bus driving was established by building a bumpy excitation road surface, and the mechanical response of the entire process on the road was analyzed.Secondly, based on the speed of the bus, the displacement and acceleration changes at different speeds were studied.Finally, the stiffness ratio of the front and rear wheels of the bus was analyzed.From these studies, the following conclusions can be drawn: (1) When a bus travels on uneven roads, the rear seats are more bumpy than the front seats.The reason for this phenomenon is that the center of gravity of the bus is closer to the rear wheels, the rear wheels bear more force and are more prone to displacement.This indicates that bus design should pay special attention to the damping of the bus backseat and control its displacement.Passengers who are "old, weak, disabled, pregnant", especially elderly people with osteoporosis, should avoid the seats at the rear of the bus as much as possible.
(2) The bus should not drive too fast on uneven roads.For the simplified bus model and excitation pavement discussed in this paper, faster speed means higher displacement

Conclusions
In this paper, from the perspective of rigid body dynamics, a motion differential equation suitable for bus driving was established by building a bumpy excitation road surface, and the mechanical response of the entire process on the road was analyzed.Secondly, based on the speed of the bus, the displacement and acceleration changes at different speeds were studied.Finally, the stiffness ratio of the front and rear wheels of the bus was analyzed.From these studies, the following conclusions can be drawn: (1) When a bus travels on uneven roads, the rear seats are more bumpy than the front seats.The reason for this phenomenon is that the center of gravity of the bus is closer to the rear wheels, the rear wheels bear more force and are more prone to displacement.This indicates that bus design should pay special attention to the damping of the bus backseat and control its displacement.Passengers who are "old, weak, disabled, pregnant", especially elderly people with osteoporosis, should avoid the seats at the rear of the bus as much as possible.
(2) The bus should not drive too fast on uneven roads.For the simplified bus model and excitation pavement discussed in this paper, faster speed means higher displacement and acceleration variation, which is not conducive to passenger comfort.Regardless of the speed, similar displacements and accelerations will ultimately be achieved, but the speed is higher, and the duration of maximum displacement and acceleration is longer.
(3) The spring stiffness ratio of the bus indicates that the bus is more stable when the rear wheels are high in spring stiffness.A bus model with low stiffness, whether climbing or downhill, will produce significant acceleration changes.During the manufacturing process of cars, it is better to ensure that the rear wheels have a higher stiffness coefficient, while appropriately reducing the stiffness of the front wheels.
The results of this paper have certain significance.Engineers can use the differential equations and analytical solutions in this paper to roughly analyze the mechanical response of a moving bus by establishing an actual excitation road surface.These results can provide some guidance for the design and driving of buses.However, the impact of damping has been ignored in this paper.Damping can be taken into account in future bus models to more accurately analyze the mechanical response of buses during driving.

Figure 1 .
Figure 1.Force model of a bus in general driving conditions.

Figure 1 .
Figure 1.Force model of a bus in general driving conditions.

Figure 2 .
Figure 2. Relationship between ground support forces on front and rear wheels

Figure 2 .
Figure 2. Relationship between ground support forces on front and rear wheels and acceleration.

Figure 3 .
Figure 3.A simplified mechanical analysis of uneven pavement.

Figure 4 .
Figure 4. Sketch of the force analysis of the upper part of the bus.

22 Figure 3 .
Figure 3.A simplified mechanical analysis of uneven pavement.

Figure 4 .
Figure 4. Sketch of the force analysis of the upper part of the bus.

Figure 4 .
Figure 4. Sketch of the force analysis of the upper part of the bus.

Figure 9 .
Figure 9. Net vertical acceleration of front and rear wheels at different speeds 0 ≤ x ≤ (b + c).

Figure 9 .
Figure 9. Net vertical acceleration of front and rear wheels at different speeds 0 ≤ x ≤ (b + c).

Figure 9 .
Figure 9. Net vertical acceleration of front and rear wheels at different speeds 0 ≤ x ≤ (b + c).

Figure 11 .
Figure 11.Net vertical acceleration of front and rear wheels at different speeds (b + c) ≤ x ≤ 2ns (n ∈ N*).

Figure 11 .
Figure 11.Net vertical acceleration of front and rear wheels at different speeds (b + c) ≤ x ≤ 2ns (n ∈ N*).

Figure 14 .
Figure 14.Net displacement and vertical acceleration 0 ≤ x ≤ (b + c) (The black line is the envelope of scattered points at low stiffness ratio; the red line is the envelope of scattered points at medium stiffness ratio; the blue line is the envelope of scattered points at high stiffness ratio).

Figure 14 .
Figure 14.Net displacement and vertical acceleration 0 ≤ x ≤ (b + c) (The black line is the envelope of scattered points at low stiffness ratio; the red line is the envelope of scattered points at medium stiffness ratio; the blue line is the envelope of scattered points at high stiffness ratio).

22 Figure 15 .
Figure 15.Net displacement and vertical acceleration (b + c) ≤ x ≤ 2ns (n ∈ N*) (The black line is the envelope of scattered points at low stiffness ratio; the red line is the envelope of scattered points at medium stiffness ratio; the blue line is the envelope of scattered points at high stiffness ratio).

Figure 15 .
Figure 15.Net displacement and vertical acceleration (b + c) ≤ x ≤ 2ns (n ∈ N*) (The black line is the envelope of scattered points at low stiffness ratio; the red line is the envelope of scattered points at medium stiffness ratio; the blue line is the envelope of scattered points at high stiffness ratio).

Figure 16 .
Figure 16.Net displacement and vertical acceleration 2ns ≤ x ≤ 2ns + (b + c) (n ∈ N*) (The black line is the envelope of scattered points at low stiffness ratio; the red line is the envelope of scattered points at medium stiffness ratio; the blue line is the envelope of scattered points at high stiffness ratio).

Figure 16 .
Figure 16.Net displacement and vertical acceleration 2ns ≤ x ≤ 2ns + (b + c) (n ∈ N*) (The black line is the envelope of scattered points at low stiffness ratio; the red line is the envelope of scattered points at medium stiffness ratio; the blue line is the envelope of scattered points at high stiffness ratio).

Table 2 .
[43]tion-at-time raw data sample gathered by the app Torque on Tue Jul 2018[43].Note: a lat means the acceleration component in latitude and a long means the acceleration component in longitude.

Table 3 .
Maximum net displacement and acceleration of front and rear wheels at different speeds.Both the front and rear wheels drive into an uneven road section.At this time, the bus position is (b + c) ≤ x ≤ 2ns (n ∈ N*).N* is a positive integer.The net displacement and the net acceleration at different speeds are presented in Figures10 and 11, respectively.

Table 3 .
Maximum net displacement and acceleration of front and rear wheels at different speeds.

Table 4 .
Maximum net displacement and acceleration of front and rear wheels at different speeds.

Table 4 .
Maximum net displacement and acceleration of front and rear wheels at different speeds.

Table 5 .
Maximum net displacement and acceleration of front and rear wheels at different speeds.

Table 6 .
The values of stiffness coefficient grouping.