Evaluation of the Crosswind Effect on the Horizontal Curve Safety on Rural Highways Using Vehicle Dynamic Simulation

of


Introduction
Strong winds and storms are some of the natural disasters that cause a lot of damage in the world every year, including in Iran. Winds are classified into different categories based on the speed and impact they have on the sea surface. According to studies by the World Meteorological Organization, winds of more than 15 meters per second are known as storms. In the 1830s and 1840s, public transportation systems were used in most of Europe's major cities, and horse-drawn vehicles became commonplace in urban traffic. In the late nineteenth century, steam engines replaced electric motors, and eventually, internal combustion engines entered the market [1]. Figures 1 and 2 show an example of these changes in pressure. e aerodynamics of road vehicles is quite different from the aerodynamics of aircraft. ree-dimensional airflow, unstable vibrations, and road surface all affect this process [3]. As the fluid passes through the surface, the velocity of the fluid also has properties such as those shown in Figure 1. erefore, a vehicle passing under constant wind conditions has a constant speed ratio according to Figure 2 [2].
How the side wind affects the vehicle and how to apply it to identify existing weaknesses and errors will be investigated by studying books and articles written about calculating the radius of horizontal arcs in combination with longitudinal slopes, and then, the behavior of vehicles in the face of arcs combined with longitudinal slope and wind is investigated using vehicle dynamic simulator software [5]. Cars, trucks, and buses are the most widely used vehicles on the roads. Because of the importance of buses due to the large number of people in them and the high traffic on the roads, this issue was studied for this particular vehicle. However, research on passenger cars and trucks also raises our awareness.

Literature Review. Attention has recently been paid to the aerodynamic loads on vehicles in inverted motion
conditions. An important feature of car roads is that they are usually located in large areas and as a result are exposed to wind [6]. e lateral wind model is based on the narrowbody theory, which examines the character of horizontal lateral wind velocity, in a 1973 study by Hochu and Amal Man.
is model shows the lateral force and torque of inverted motion (Figure 3).
is study showed how the lateral wind speed characteristic affects the size of the offroad and that relying solely on wind tunnel measurement when analysing sensitivity to strong winds will not be sufficient in Figure 4 [7]. In 1991, Ohno and Kohri investigated the effects of lateral air barriers on traction and lift forces. ey showed that traction and elevation can be reduced by zero lateral inversion through the lateral air barrier, but when the inverted angle is increased, elevation will also increase [8]. However, the rise can be reduced even in conditions of zero inverted movements by optimizing the shape of the lateral air barrier as shown in Figure 5 [9]. e models were also equipped with pressure transmitters that made it possible to analyse different types of pressure while creating movements during lateral winds [10]. ey concluded that the transient reaction of the body with sharp parts can be detected by the formation and destruction of the separation bubble, leading to the achievement of the maximum amount of inverted torque. In addition, the creation of radii at sharp edges leads to an increase in constant inverse torque due to increased local pressures [11]. Ryan and Dominy [12] and Amundsen and Ranes [13] obtained the aerodynamic force of the passenger car model by overturning in the wind heading of about 30 [12,13]. Fang [14] adjusted a vehicle dynamic simulation  and model of high-speed aerodynamic characteristics and analysed the speed and the effects of lateral winds on vehicle displacement stability [14]. Maruyama and Yamazaki in 2006 examined the risk of an accident based on vehicle weight, tire load, aerodynamic force, and driving behavior, including strong mental behavior and other factors that are very difficult to determine [15]. Wu [16] studied the safety of highways affected by Typhoon climate and developed a car slip model. In addition, he calculated the speed of the car's safety stimulus under different degrees of Typhoon climate but did not produce a quantitative wind response on paper [16]. Rodriguez and his colleagues (2014) used a stimulus simulator to measure the effect of opposite wind on driver behavior and vehicle lateral displacement and used a calculation model based on cross-cultural strength, driver behavior, and vehicle lateral displacement [17]. Zhou et al. [18] and Ueckermann et al. [19] studied the effect of opposite winds on design indices of alignment with wind tunnel testing but did not consider the effect of aerodynamic torque generated by wind and the unequal coefficient on the lateral stability of torques [18,19].

e Geometry of the Route and Type of Vehicle.
e performance of vehicles in different geometric road conditions and at different speeds was evaluated to examine their stability against overturning and slipping. e design speeds of 40, 60, 80, 100, and 120 km/h, longitudinal slopes of −8%, −4%, 0, 4%, and 8%, and maximum design speed of 8% are considered. Besides, six side wind speeds were selected based on the World Meteorological Organization's classification. erefore, the lateral winds of 0, 16, 34, 55, 80, and 111 km/h were selected. e minimum radius of the arcs is also used according to the AASHETO suggestion. Furthermore, the transition area with length Lr is considered to provide safe entry and exit of vehicles to the horizontal arcs [20]. Table 1 summarizes this explanation. In this study, a bus was selected [21]. e dimensional and mechanical specifications of vehicles are given in Table 2.

e Coefficient of Friction along the Path in Different
Modes. Figure 6 shows the coefficient of lateral friction for the combination of v � 80 km/h and g � −8% along the path. In Figure 6, the coefficient of lateral friction in different parts of the path, including the initial straight path, circular arc, and straight end path, can be observed. e coefficient of friction is constant in the straight path until the start of the round application (beginning of Lr, station 41). During the round up to station 80 or the starting point of the arc, it increases to 0.05 due to the increase in lateral force in the positive direction. At the beginning of the arc, the lateral force changes direction due to the centrifugal force, so the lateral friction is reduced and this trend continues until the  V C V C∞ Figure 4: Lateral force coefficient and inverted torque due to strong wind [7]. Shock and Vibration end of Lr (station 90) and reaches −0.18. en, it remains almost constant along the arc and at the endpoint of the arc (approximately station 257), and the friction jumps at once due to the disappearance of the centrifugal force and increases to the value of 0.055. en, l r ′ has a decreasing trend by the end of it (station 295) due to the reduction in the existing lateral friction coefficient and reaches 0.01. As can be observed, the trend of changes in the coefficient of friction in the case of lateral wind speed of 111 km/h is wider than other wind speeds. erefore, as the lateral wind speed increases, the amount of friction increases due to the increase in lateral force. Figure 7 shows the coefficient of lateral friction available for the combination of w � 34 km/h and v � 80 km/h along the route.

Side-airdam
As shown, the trend of changes in the coefficient of friction in the case of a slope of 8% is wider than that of other slopes. erefore, the slope of −8% has the highest coefficient of friction and then the slope is −4% and the slopes of 0 and 4% are almost the same and the slope of 8% has the lowest coefficient of friction. e higher the coefficient of friction, the greater the road safety. Of course, it should be noted that Figure 7, despite the constant lateral wind, shows w � 34 km/h. Figure 8 shows the coefficient of lateral friction available for the combination of w � 34 km/h and g � −8% along the path. As can be seen, the trend of changes in the coefficient of friction at a speed of 40 km/h has a higher coefficient of friction than that of other speeds, and as the speed increases, the coefficient of friction decreases, resulting in reduced safety. Figure 9 shows the maximum values of the coefficient of friction on a slope of −4% for different wind speeds. According to Figure 9, the maximum value of the coefficient of friction at a speed of 40 km/h is for the time when there is no lateral wind, but at other speeds, the maximum value of the coefficient of friction is for the state where the lateral wind is its maximum value. For a side wind of 111 km/h, there are three data because the bus overturns at speeds of 100 and 120 km/h, and there are three data at a wind speed of 80 km/h because the car overturns when the speed is 120 km/h [22]. Table 3 shows the maximum values of the coefficient of friction, the percentage difference between them, the ground state at different speeds, and the wind speed of 16 km/h for the longitudinal slope state of −4%.

Analysing the Statistical Data of Friction Coefficient.
Since the value of the Kolmogorov-Smirnov is equal to 0.939, the assumption that the residues are normal is confirmed. e second method can be observed from Figure 10, which confirms the assumption that the residues are normal. e third method is to use Figure 10 to achieve the normality of the residues. erefore, the assumption that the residues are normal is confirmed. According to Figure 10, the variance of the residues is constant and its hypothesis is confirmed.
is hypothesis is confirmed using the Durbin-Watson test value. e statistical values of this test are between 0 and 4, in which the values less than 1 will indicate positive correlations. e Durbin-Watson value of this model is 1.753, so the assumption that the residues are uncorrelated is confirmed. As mentioned, this issue is referred to as the alignment problem in the statistical literature. e values of Vif and tolerance are used to investigate the presence or absence of alignment. us, if Vif <10 and tolerance> 0.1, there is no alignment problem and the assumption of explanatory variables is confirmed. According to the values of Vif and tolerance, the assumption of explanatory variables is confirmed and the average value of residues is zero. Finally, the proposed model of maximum lateral friction in the horizontal arc for the bus is shown as follows: Figure 11 shows the lateral acceleration for the combination of v � 80 km/h and g � −8% along the path. As shown, in the case where there is no lateral wind speed, the lateral acceleration is at its maximum and the lowest lateral acceleration is for the situation where the lateral wind is at its   erefore, as the lateral wind decreases, the lateral acceleration also increases, which reduces road safety. Figure 12 shows the lateral acceleration for the combination of w � 34 km/h and v � 80 km/h along the route.

Lateral Acceleration along the Route in Various Modes.
As shown, the −8% slope has the highest lateral acceleration and the 8% slope has the lowest lateral acceleration.
erefore, under constant lateral wind of 34 km/h, the slopes of −8%, −4%, 0%, 4%, and 8% have the highest lateral acceleration and the lowest safety, respectively. Figure 13 shows the lateral acceleration for the combination of w � 34 km/h and g � −8% along the path. As seen, the speed of 40 km/h has the highest lateral acceleration and the reason is because of its low radius, which also has the lowest safety compared with others. erefore, the higher the speed, the lower the lateral acceleration and the greater the safety. Figure 14 shows the maximum values of lateral acceleration on a slope of −4% for different speeds.
Depending on the shape of the maximum lateral acceleration at speeds of 40 and 60 km/h, the lower the lateral wind, the more they are approximately equal at higher speeds. Table 4 shows the maximum lateral acceleration values, the percentage difference between them, the ground state at different speeds, and the wind speed of 16 km/h for the longitudinal slope state of −4%.

Analysing the Statistical Data of Lateral Acceleration.
First, considering that the value of the Kolmogorov-Smirnov is equal to 0.354, the assumption that the residues are normal is confirmed. Second, Figure 15 shows that the assumption that the residues are normal is confirmed. e third method is to use Figure 15 to achieve the normality of the residues. erefore, the assumption that the residues are normal is confirmed. According to Figure 15, it can be seen that the variance of the residues is constant and its hypothesis is confirmed.
is hypothesis is confirmed using the Durbin-Watson test value. e statistical values of this test are between 0 and 4, in which the values less than 1 will indicate positive correlations. e value of the Durbin-Watson value of this model is equal to 1,534, so the assumption that the residues are uncorrelated is confirmed. As mentioned, in the statistical literature, this issue is referred to as the alignment problem. e    Figure 16 shows the overturning angle for the combination of v � 80 km/h and g � −8% along the path. As observed, in the case of lateral wind speed of 111 km/h, the overturning angle is its maximum value and the lowest overturning angle is for the case where there is no lateral wind. erefore, as the lateral wind increases, the overturning angle also increases, which reduces road safety. Figure 17 shows the reverse angle for the combination of w � 34 km/h and v � 80 km/h along the path.
As shown, the slopes of 0%, 4%, and −4% are almost the same, the highest overturning angle has 8% slope and −8% slope having the lowest overturning angle. Figure 18 shows the overturning angle for the combination of w � 34 km/h and g � −8% along the path. As shown, the speed of 120 km/ h has the highest overturning angle due to the high speed when entering the arc. erefore, the higher the speed, the greater the overturning angle, resulting in lower safety. Figure 19 shows the maximum values of the overturning angle at a slope of −4% for different speeds.
According to Figure 19, at equal speeds, the higher the wind speed, the higher the maximum overturning angle. Table 5 shows the maximum values of the overturning angle, the percentage difference between them, the ground state at different speeds, and the wind speed of 16 km/h for the longitudinal slope state of −4%.

Analysing the Statistical Data of Overturning Angle.
First, considering that the value of the Kolmogorov-Smirnov is equal to 0.442, the assumption that the residues are normal is confirmed. e second method can be seen from Figure 20, which confirms the assumption that the residues are normal. ird, using Figure 20 can contribute to achieve the normality of the residues. erefore, the assumption that the residues are normal is confirmed. According to Figure 20, the variance of the residues is constant and its hypothesis is confirmed.
is hypothesis is confirmed using the Durbin-Watson test value. e statistical values of this test are between 0 and 4, in which the values less than 1 will indicate positive correlations.
e Durbin-Watson value of this model is 2.048, so the assumption that the residues are uncorrelated is confirmed. As mentioned, in the statistical literature, this issue is referred to as the alignment problem. e values of Vif are used to study the presence or absence of alignment.
us, if Vif <10 and tolerance> 0.1, there is no alignment problem and the assumption of uncorrelated variables is confirmed. According to the values of Vif, the assumption of explanatory variables is confirmed and the average value of residues is zero. Finally, the proposed model of maximum overturning angle in the horizontal arc for the bus is shown as follows: roll � 0.00005406w 2 − 0.01 g 2 + 1.526 logs + 1.313.
(3) Figure 21 shows the slip angle for the combination of v � 80 km/h and g � −8% along the path. In the case of lateral wind speed of 111 km/h, the slip angle is its maximum and the lowest slip angle is for the case where there is no lateral  wind. erefore, as the lateral wind increases, the slip angle also increases, which reduces road safety. Figure 22 demonstrates the slip angle for the combination of w � 34 km/h and v � 80 km/h along the path.

Slip Angle along the Path in Different Modes.
As can be seen, −8% slope has the highest slip angle and the 8% slope has the lowest slip angle.
erefore, under constant lateral wind of 34 km/h, the slopes of −8%, −4%, 0%, 4%, and 8%, respectively, have the highest slip angle and the lowest safety. Figure 23 shows the slip angle for the combination of w � 34 km/h and g � −8% along the path. As can be seen in the bus, the speed of 120 km/h has the highest angle of slip, followed by the speeds of 60, 80, 40, and 100, respectively.    Figure 24 shows the maximum slip angle values at a slope of −4% for different velocities. According to Figure 24 at equal speeds, the higher the wind speed, the higher the maximum slip angle. Table 6 shows the maximum values of the slip angle, the percentage difference between them, the ground state at different speeds, and the wind speed of 16 km/h for the longitudinal slope state of −4%.

Analysing the Statistical Data of Slip Angle.
First, considering that the Kolmogorov-Smirnov value is equal to 0.2, the assumption that the residues are normal is confirmed. Second, Figure 25 shows that the assumption that the residues are normal is confirmed. ird, using Figure 25 contributes to achieve the normality of the residues. erefore, the assumption that the residues are   normal is confirmed. According to Figure 25, the variance of the residues is constant and its hypothesis is confirmed.
is hypothesis is confirmed using the Durbin-Watson test value. e statistical values of this test are between 0 and 4, in which the values less than 1 will indicate positive correlations.
e Durbin-Watson value of this model is 1.509, so the assumption that the residues are uncorrelated is confirmed. As mentioned, in the statistical literature, this issue is referred to as the alignment problem. To investigate the presence or absence of alignment, the values of Vif and tolerance are used. us, if Vif <10 and tolerance> 0.1, there is no alignment problem and the assumption of uncorrelated variables is confirmed. According to the values of Vif, the assumption of explanatory variables is confirmed and the average value of residues is zero. Finally, the proposed model of maximum slip angle in the horizontal arc for the bus is shown as follows: slip � −0.004w + 0.0000001 g 3 − 82.273 inverse s + 1.024. (4)

Conclusion
In this research, the issue and the subject of the research were discussed and the necessity of conducting research was explained. en, the questions and hypotheses were stated along with the main research objectives to clarify the details of the work. One of the highlights of this chapter is the effects of lateral winds on vehicle aerodynamics. In the following, the literature review has been discussed and an attempt has been made to collect and present a complete collection of researchers' studies on this subject. e studies of Mr Bowman, Takada, and Ohri, as well as Heyat and P, can be introduced as the most comprehensive studies in the field of lateral wind. en, it gives a complete description of the method used for the research and explains the details of the study process. Accordingly, if the vehicle speed and longitudinal slope are constant, the existing lateral friction will increase as the lateral wind increases. For every 10 km/h increase in lateral wind, the lateral friction for each vehicle, SUV, bus, and two-axle truck increases by 3%, 13%, and 27%, respectively. en, if the longitudinal slope and lateral wind intensity are constant, the existing lateral friction coefficient decreases as the vehicle speed increases. Per each 10 km/h increase in vehicle speed, the lateral friction for each vehicle, SUV, bus, and two-axle truck is reduced by 9%, 14%, and 17%, respectively.
If the lateral wind intensity and vehicle speed are constant, as the longitudinal slope shifts from uphill to downhill, the existing lateral friction coefficient increases; i.e., the longitudinal slope is 8% with the lowest available lateral friction coefficient and the longitudinal slope is −8%. e maximum available lateral friction coefficient is available. For every one percent increase in the longitudinal slope to the downside, the lateral friction for each vehicle, SUV, bus, and two-axle truck increases by 4, 5, and 6%, respectively. Accordingly, if the vehicle speed and longitudinal slope are constant, the lateral acceleration increases as the lateral wind increases. For every 10 km/h increase in lateral wind, the lateral acceleration for each vehicle, SUV, bus, and two-axle truck increases by 2, 3, and 5%, respectively. Accordingly, if the vehicle speed and longitudinal slope are constant, the overturning angle will increase as the lateral wind increases. For every 10 km/h increase in lateral wind, the overturning angle for each vehicle, SUV, bus, and two-axle truck increases by 1.5%, 8%, and 12%, respectively. Furthermore, if the longitudinal slope and lateral wind intensity are constant, the overturning angle will increase as the vehicle speed increases. For every 10 km/h increase in vehicle speed, the overturning angle for each vehicle, SUV, bus, and two-axle truck is reduced by 2, 3, and 5%, respectively. Accordingly, if the vehicle speed and longitudinal slope are constant, the deviation rate increases as the lateral wind increases. For every 10 km/h increase in lateral wind, the deviation rate for each vehicle, SUV, bus, and two-axle truck increases by 1.5, 2.5, and 5%, respectively. en, if the longitudinal slope and lateral wind intensity are constant, the deviation rate will decrease as the vehicle speed increases.

Data Availability
Requests for access to these data should be made to the corresponding author (email address: aliabdi@eng.ikiu.ac.ir)