Analysis of Driving Dynamics Considering Driving Resistances in On-Road Driving

Internal combustion engine emissions are a serious worldwide problem. To combat this, emission regulations have become stricter with the goal of reducing the proportion of transportation emissions in global air pollution. In addition, the European Commission passed the real driving emissions–light-duty vehicles (RDE-LDV) regulation that evaluates vehicle emissions by driving on real roads. The RDE test is significantly dependent on driving conditions such as traffic or drivers. Thus, the RDE regulation has the means to evaluate driving dynamics such as the vehicle speed per acceleration (v·apos) and the relative positive acceleration (RPA) to determine whether the driving during these tests is normal or abnormal. However, this is not an appropriate way to assess the driving dynamics because the v·apos and the RPA do not represent engine load, which is directly related to exhaust emissions. Therefore, in the present study, new driving dynamic variables are proposed. These variables use engine acceleration calculated from wheel force instead of the acceleration calculated from the vehicle speed, so they are proportional to the engine load. In addition, a variable of driving dynamics during braking is calculated using the negative wheel force. This variable can be used to improve the accuracy of the emission assessment by analyzing the


Introduction
Greenhouse gas and air pollutant emissions are two of the most important issues currently discussed in the automotive market. Particulate matter and nitrogen oxide are mainly emitted from road transport sources and are key reasons for establishing strict regulatory policies on exhaust and greenhouse gas worldwide [1][2][3]. These regulatory policies have encouraged the development of emission-reduction technologies for automobiles. However, they could not prevent the use of a special software that could detect emission test mode and change its driving pattern to meet emission levels [4]. As a result, there is a significant difference between the amount of exhaust gas that is emitted during certification tests and that is emitted during on-road driving [5][6][7]. Therefore, the European Union has developed and applied the real driving emission-light-duty vehicles (RDE-LDV) test method that measures and evaluates exhaust while driving on actual roads.
The RDE-LDV uses the portable emissions measurement system (PEMS). This system is installed in a vehicle to measure the exhaust gas emitted as it moves on a road. Since this method evaluates vehicle emissions while driving on actual roads under random traffic and/or weather conditions, it has the advantage of measuring emissions under realistic driving situations. In addition, it does not allow automobile manufacturers to target specific driving conditions by setting up controlled driving strategies. However, because of these random conditions, all the emission data from on-road driving are not considered valid. The RDE tests must fulfill standard test procedures and driving conditions for their results   Figure 1 shows the driving route and altitude information utilized in this study. According to RDE standards [18], the driving route must include urban, rural, and motorway roads sequentially, and the distance ratio and velocity of each section must be within a specified range to be accepted as a valid test. The test vehicle entered the urban area on the right side of the map and drove approximately 27 km from the starting point. Afterward, it drove on rural roads, moving from right to left on the map, for approximately 35 km. Finally, the vehicle drove the motorway section, moving from the left to the top right of the map, for approximately 30 km. The driving route configuration met the distance ratio criteria. The start and end points had similar altitudes [19], with a gentle uphill slope in the rural section and a gentle downhill slope in the motorway section.  Figure 1 shows the driving route and altitude information utilized in this study. According to RDE standards [18], the driving route must include urban, rural, and motorway roads sequentially, and the distance ratio and velocity of each section must be within a specified range to be accepted as a valid test. The test vehicle entered the urban area on the right side of the map and drove approximately 27 km from the starting point. Afterward, it drove on rural roads, moving from right to left on the map, for approximately 35 km. Finally, the vehicle drove the motorway section, moving from the left to the top right of the map, for approximately 30 km. The driving route configuration met the distance ratio criteria. The start and end points had similar altitudes [19], with a gentle uphill slope in the rural section and a gentle downhill slope in the motorway section.

Analysis Methods
After the driving test is conducted according to specified test procedure and driving conditions, the results must be verified per RDE standards. Data completeness and normality are verified through the moving average window (MAW) and the power binning method. If the test results are determined as valid, carbon dioxide (CO2) amounts measured during the real road test and measured under the worldwide harmonized light vehicles test procedure (WLTP) mode are compared to correct exhaust emissions.
The present study focused on data analysis, the step before final exhaust emission correction. In on-road driving data analysis, the MAW method configures consecutive average windows of a certain distance at one-second intervals intervals from the starting

Analysis Methods
After the driving test is conducted according to specified test procedure and driving conditions, the results must be verified per RDE standards. Data completeness and normality are verified through the moving average window (MAW) and the power binning method. If the test results are determined as valid, carbon dioxide (CO 2 ) amounts measured during the real road test and measured under the worldwide harmonized light vehicles test procedure (WLTP) mode are compared to correct exhaust emissions.
The present study focused on data analysis, the step before final exhaust emission correction. In on-road driving data analysis, the MAW method configures consecutive average windows of a certain distance at one-second intervals intervals from the starting point to the end point. It also verifies the completeness and normality of the driving data based on the average velocity and exhaust data in each window. The distance of the window is determined so that the CO 2 emission in the window is the same as half of the total CO 2 emission by the test vehicle in the WLTP mode. Vehicle velocity data were obtained from the GPS and the OBD, and they were judged to be reliable since they were nearly the same [17]. However, since the GPS signals were occasionally disconnected in some areas, this study used the velocity data acquired from the OBD for the analysis. Since the acceleration resolution calculated from the velocity data was not high enough (a res > 0.01 m/s 2 ), the velocity data was smoothed with a T4253H filter (as suggested in the RDE regulation). The altitude data were obtained from the GPS at 1-s intervals. The interval of altitude data was then converted into 1 m of vehicle movement to calculate the slope of the road. The converted altitude data were not high in resolution, so they were corrected through a two-step smoothing process. The v·a pos (m 2 /s 3 or W/kg) and RPA (m/s 2 ) were used as parameters for analyzing the driving dynamics. In v·a pos , v refers to the vehicle velocity (m/s) and a pos refers to the vehicle acceleration (m/s 2 ) that exceeds 0.1 m/s 2 . The unit is (W/kg), referring to the power per unit mass (specific power). RPA was calculated separately for the urban, rural, and motorway sections with the following equations: where j: 1 to M k i: 1 to N k k: urban, rural, motorway ∆t: time step (s) (=1 s) d i,k : distance covered in time step i considering the urban, rural and motorway shares (m) M k : number of samples for urban, rural and motorway shares with a pos N k : total number of samples for urban, rural and motorway shares

Results and Discussions
The present study analyzed the on-road driving data with the MAW method to check the completeness and normality of the data. The driving dynamics of the vehicle were analyzed by calculating the v·a pos and RPA according to the method suggested by RDE standards. Afterward, various driving resistances, such as air drag or rolling resistance, were calculated. A new vehicle speed per acceleration (v·a eng ) and RPA were calculated using these driving resistances. These new driving dynamics variables were then compared with the existing results. Finally, one more new parameter (v·a neg ) that can evaluate the driving dynamics during the braking process was suggested. This was done by identifying when the engine does not generate power (e.g., inertial driving or braking).

Data Analysis According to the Real Driving Emission Standard
The present study verified the completeness and normality of the real road test results according to the MAW method. The MAW method uses the average velocity in each window to identify urban, rural, or motorway roads. The window with an average velocity of less than 45 km/h is defined as urban. The window with an average velocity of 45 km/h or more, but less than 80 km/h, is defined as rural. The window with an average velocity of 80 km/h or more, but less than 145 km/h, is defined as motorway. In Figure 2, the CO 2 emissions per unit distance is plotted on the CO 2 characteristics curve. The dot data on the graph are the CO 2 emissions per unit distance (g/km), calculated by dividing the CO 2 emissions in each window by the moving distance. The black solid line is the CO 2 characteristic curve based on CO 2 emissions in WLTP mode, and the two pairs of gray dotted lines refer to the primary and secondary allowable limits. The data within the primary limits are defined as normal testing conditions. Data equal to or higher than the primary limit and below the secondary limit are defined as severe testing conditions, and data equal to or lower than the primary limit and above the secondary limit are Energies 2021, 14, 3408 5 of 16 defined as soft testing conditions. In other words, a higher CO 2 emission per unit distance is considered a more severe condition. If a datum point falls outside of the secondary allowable limits it is considered invalid in the data correction process. dotted lines refer to the primary and secondary allowable limits. The data within the primary limits are defined as normal testing conditions. Data equal to or higher than the primary limit and below the secondary limit are defined as severe testing conditions, and data equal to or lower than the primary limit and above the secondary limit are defined as soft testing conditions. In other words, a higher CO2 emission per unit distance is considered a more severe condition. If a datum point falls outside of the secondary allowable limits it is considered invalid in the data correction process. As is shown if Figure 2, all the urban data were in normal conditions. However, a significant portion of rural data in the 70-80 km/h section were in the severe testing condition or outside of the allowable limit. Although most data were in the normal condition in the motorway section, some fell in severe testing conditions. Table 3 summarizes the quantitative analysis of the data. In the MAW method, the driving distance ratios in the urban, rural, and motorway sections must be configured to be about 34%, 33%, and 33%, respectively, and none can be less than 15%. As is shown in Table 3, the driving distance ratios in the urban, rural, and motorway sections were 34%, 43%, and 23%, respectively. This met the completeness requirement since all were 15% or more. The normality criterion states that at least 50% of the data should be in the normal condition in each window. The rural section showed the lowest normal data rate and was at 54%. This demonstrated that the test data met the normality requirement. Therefore, the on-road driving data fulfilled all validity criteria required by the MAW method. The severity index parameter checks the deviation of data from the CO2 characteristic curve and indicates a severe driving condition when the value is high. As is shown in Figure 2, most datum points fell above the black solid line, indicating a positive severity in all driving windows. The rural section had the highest record of 27.33% as most datum points fell in the severe driving condition area. The validity of the driving dynamics was also checked. There are two representative driving dynamics in RDE standard: v·apos and RPA. Figure 3a,b shows the graphs that confirm the validity of v·apos and RPA, respectively. The urban, rural, and motorway sections were based on spontaneous velocity instead of average velocity. The urban section was defined as 60 km/h or less. The rural section was defined as between 60 km/h and 90 As is shown if Figure 2, all the urban data were in normal conditions. However, a significant portion of rural data in the 70-80 km/h section were in the severe testing condition or outside of the allowable limit. Although most data were in the normal condition in the motorway section, some fell in severe testing conditions. Table 3 summarizes the quantitative analysis of the data. In the MAW method, the driving distance ratios in the urban, rural, and motorway sections must be configured to be about 34%, 33%, and 33%, respectively, and none can be less than 15%. As is shown in Table 3, the driving distance ratios in the urban, rural, and motorway sections were 34%, 43%, and 23%, respectively. This met the completeness requirement since all were 15% or more. The normality criterion states that at least 50% of the data should be in the normal condition in each window. The rural section showed the lowest normal data rate and was at 54%. This demonstrated that the test data met the normality requirement. Therefore, the on-road driving data fulfilled all validity criteria required by the MAW method. The severity index parameter checks the deviation of data from the CO 2 characteristic curve and indicates a severe driving condition when the value is high. As is shown in Figure 2, most datum points fell above the black solid line, indicating a positive severity in all driving windows. The rural section had the highest record of 27.33% as most datum points fell in the severe driving condition area. The validity of the driving dynamics was also checked. There are two representative driving dynamics in RDE standard: v·a pos and RPA. Figure 3a,b shows the graphs that confirm the validity of v·a pos and RPA, respectively. The urban, rural, and motorway sections were based on spontaneous velocity instead of average velocity. The urban section was defined as 60 km/h or less. The rural section was defined as between 60 km/h and 90 km/h, and the motorway section was defined as more than 90 km/h. Figure  urban section was longer than in the other sections and showed more frequent acceleration and deceleration due to traffic signals. Therefore, more data exceeded 0.1 m/s 2 (=a pos ) and were displayed on the urban section than on the other sections [20]. The driving duration was relatively short in the motorway section, and the speed change was not large when compared with the other sections. Therefore, only a small number of data were displayed. Since v·a pos _(95) and v·a pos _(mean) contain the multiplied velocity term, they increase as velocity increases. The test vehicle drove below the upper limit during the whole test procedure, and it was judged that there was no excessively severe driving from the aspect of the driving dynamics.
km/h, and the motorway section was defined as more than 90 km/h. Figure 3a defines the v·apos corresponding to the top 95% in each section (urban, rural, and motorway) as v·apos_(95). It defines the average v·apos in each window as v·apos_(mean). Each small gray dot represents individual v⋅apos values, and the solid black line is the highest line that v·apos_(95) should not exceed for this test to be recognized as valid. The driving time in the urban section was longer than in the other sections and showed more frequent acceleration and deceleration due to traffic signals. Therefore, more data exceeded 0.1 m/s 2 (=apos) and were displayed on the urban section than on the other sections [20]. The driving duration was relatively short in the motorway section, and the speed change was not large when compared with the other sections. Therefore, only a small number of data were displayed. Since v·apos_(95) and v·apos_(mean) contain the multiplied velocity term, they increase as velocity increases. The test vehicle drove below the upper limit during the whole test procedure, and it was judged that there was no excessively severe driving from the aspect of the driving dynamics. The lower limit of acceleration was set for the RPA in Figure 3b to check validity of driving dynamics. The RPA was obtained by summing the v·apos calculated in each driving section and dividing it by the moving distance in the section. Due to this, the highest value was shown in the urban section since it had the largest data set even though the v⋅ apos_(mean) was small. Since the data in the rural and motorway sections were not significantly different, the RPA was determined according to the v·apos_(mean). The actual road driving test conducted in this study confirmed that the PRA was above the lower limit in all windows. Table 4 summarizes the quantitative analysis results for the data. According to RDE standard, the valid moving distance ratio in the urban section is 29-44% and 23-43% in the rural and motorway sections. All the ratios were in the valid range. Since the range of the total valid driving period is 90-120 min, 107 min of it was valid. As is shown in Figure  3a, the amount of data in each section were the highest at 1259 in the urban section and the lowest at 351 in the motorway section. The analysis confirmed that this test had satisfied all the validity criteria required by the RDE standard.  The lower limit of acceleration was set for the RPA in Figure 3b to check validity of driving dynamics. The RPA was obtained by summing the v·a pos calculated in each driving section and dividing it by the moving distance in the section. Due to this, the highest value was shown in the urban section since it had the largest data set even though the v·a pos _(mean) was small. Since the data in the rural and motorway sections were not significantly different, the RPA was determined according to the v·a pos _(mean). The actual road driving test conducted in this study confirmed that the PRA was above the lower limit in all windows. Table 4 summarizes the quantitative analysis results for the data. According to RDE standard, the valid moving distance ratio in the urban section is 29-44% and 23-43% in the rural and motorway sections. All the ratios were in the valid range. Since the range of the total valid driving period is 90-120 min, 107 min of it was valid. As is shown in Figure 3a, the amount of data in each section were the highest at 1259 in the urban section and the lowest at 351 in the motorway section. The analysis confirmed that this test had satisfied all the validity criteria required by the RDE standard.

Driving Resistances
Ruling out excessive exhaust emission test cases (due to excessive power) and minimal exhaust emission test cases (due to unusual cruising speeds), both of which deviate from typical driving patterns, are the reasons behind verifying the validity of driving dynamics in on-road driving tests. Therefore, driving dynamics must reflect engine load, which is directly related to exhaust emissions. However, since the v·a pos calculated by the conventional method does not reflect the deceleration and acceleration due to external forces, it is difficult to say that it represents engine load. For example, air resistance at a higher velocity is greater than that at a lower velocity, even under the same cruising driving conditions with zero acceleration, and it requires a higher engine load. The load that the engine faces even differs when driving uphill and downhill at the same speed. Therefore, the driving dynamics calculation must consider these driving resistances to reflect the engine load. Figure 4 depicts standard driving resistances applied to a moving vehicle, and the resistances have the following correlation.
where F w : wheel force (N) R d : drag resistance (N) R r : rolling resistance (N) R g : gradient resistance (N) R a : acceleration resistance (N) This section calculates the driving resistances of the test vehicle and compares them to examine how these resistances (other than acceleration) affect engine load.

Driving Resistances
Ruling out excessive exhaust emission test cases (due to excessive power) and minimal exhaust emission test cases (due to unusual cruising speeds), both of which deviate from typical driving patterns, are the reasons behind verifying the validity of driving dynamics in on-road driving tests. Therefore, driving dynamics must reflect engine load, which is directly related to exhaust emissions. However, since the v·apos calculated by the conventional method does not reflect the deceleration and acceleration due to external forces, it is difficult to say that it represents engine load. For example, air resistance at a higher velocity is greater than that at a lower velocity, even under the same cruising driving conditions with zero acceleration, and it requires a higher engine load. The load that the engine faces even differs when driving uphill and downhill at the same speed. Therefore, the driving dynamics calculation must consider these driving resistances to reflect the engine load. Figure 4 depicts standard driving resistances applied to a moving vehicle, and the resistances have the following correlation.
where Fw: wheel force (N) Rd: drag resistance (N) Rr: rolling resistance (N) Rg: gradient resistance (N) Ra: acceleration resistance (N) This section calculates the driving resistances of the test vehicle and compares them to examine how these resistances (other than acceleration) affect engine load.

Air Resistance (Ra)
The air resistance applied to a vehicle is calculated as follows: where Cd: drag coefficient A: full projected area (m 2 ) ρ: air density (kg/m 3 ) v: vehicle speed (m/s)

Air Resistance (R a )
The air resistance applied to a vehicle is calculated as follows: where C d : drag coefficient A: full projected area (m 2 ) ρ: air density (kg/m 3 ) v: vehicle speed (m/s) The drag coefficient (C d ) of the test vehicle used in this study was 0.272, and the full projected area (A) was assumed to be 80% of the product of the vehicle's overall width and height. The PEMS measured the air density (ρ) and velocity (v). Figure 5a shows the air resistance calculated using Equation (3) over a wide vehicle velocity range. Since the air resistance is proportional to the square of the vehicle velocity, it rapidly increased as the velocity increased. Air resistance was approximately 283 N at 100 km/h. The drag coefficient (Cd) of the test vehicle used in this study was 0.272, and the full projected area (A) was assumed to be 80% of the product of the vehicle's overall width and height. The PEMS measured the air density (ρ) and velocity (v). Figure 5a shows the air resistance calculated using Equation (3) over a wide vehicle velocity range. Since the air resistance is proportional to the square of the vehicle velocity, it rapidly increased as the velocity increased. Air resistance was approximately 283 N at 100 km/h.

Rolling Resistance (Rr)
It is difficult to find an accurate formula to calculate rolling resistance because tire size, tire type, and even road conditions affect it. It is known that rolling resistance generally increases as the vehicle weight increases, as the tire inflation pressure decreases, and as the velocity increases [21,22]. Since the vehicle weight and inflation pressure were fixed in this study, the parameter that affected the rolling resistance was the vehicle velocity, and it was calculated using the equation below [23].

Rolling Resistance (R r )
It is difficult to find an accurate formula to calculate rolling resistance because tire size, tire type, and even road conditions affect it. It is known that rolling resistance generally increases as the vehicle weight increases, as the tire inflation pressure decreases, and as the velocity increases [21,22]. Since the vehicle weight and inflation pressure were fixed in this study, the parameter that affected the rolling resistance was the vehicle velocity, and it was calculated using the equation below [23].
where C sr = −514.7(P tire /100) 3  The equation uses the road gradient because wheel load varies according to the road gradient. The road gradient was calculated using the altitude data (see Figure 1) obtained with the GPS. Figure 5a shows the rolling resistance calculated with Equation (4) over a wide vehicle velocity range. Although the rolling resistance was much higher than the Energies 2021, 14, 3408 9 of 16 air resistance in low-velocity sections (approximately 180 N), it increased slowly as the velocity increased. It was surpassed by the air resistance at the velocity of approximately 88 km/h. The rolling resistance was roughly 233 N (about 83% of the air resistance) at the vehicle velocity of 100 km/h. Due to this tendency, the rolling resistance variation throughout the test process was very small compared to the other resistances. For example, the air resistance increased by about 283 N when the velocity increased from 0 km/h to 100 km/h, while the rolling resistance only increased by approximately 66 N. The changes in gradient resistance and acceleration resistance were much larger. Therefore, the wheel force in Equation (2) was mostly affected by the gradient resistance and the acceleration resistance, while the impact of the rolling resistance was relatively small. As mentioned above, it is difficult to find an accurate formula to calculate the rolling resistance, so the calculation accuracy can be poor. However, the error from calculation inaccuracy is also relatively small because of its small impact.

Gradient Resistance (R g )
The gradient resistance is the resistance caused by the vehicle weight. It depends on the slope of the road. It becomes a resistance on an uphill slope and an acceleration force on a downhill slope, as the following equation shows.
where W: vehicle weight (N) θ: road gradient (deg) Figure 5b shows the road gradient converted from the GPS altitude data. It also shows the gradient resistance calculated using Equation (5). The maximum road gradient in the driving route was about 3.8 • (6.6%), and its gradient resistance was about 1160 N. The maximum road gradient on a downhill was about −3.2 • (−5.6%), and its acceleration force was about 980 N. The gradient resistance on a typical road gradient of 1 • (1.7%) was approximately 306 N, higher than air resistance or rolling resistance at 100 km/h [24].

Acceleration Resistance (R a )
Acceleration resistance is a force needed for vehicle acceleration (rather than a resistance) and is calculated with the following equation: where m: vehicle mass (kg) a: vehicle acceleration (m/s 2 ) The graph in Figure 5c depicts the measured acceleration and the acceleration resistance. The maximum acceleration that occurred during the driving test was about 2.0 m/s 2 , and its acceleration resistance was approximately 3600 N. This was much larger than the other resistances. Deceleration by braking was controlled within −3 m/s 2 , except for one sudden braking event at 73 min in the graph.
During the entire driving process, the variation ranges in air resistance, rolling resistance, gradient resistance, and acceleration resistance were ±233 N, ±57 N, ±1078 N, and ±1825 N, respectively. Acceleration resistance showed the largest variation range, even though only positive acceleration was considered. Rolling resistance showed the smallest range. This indicates the impact of the wheel force on the calculation. While the rolling resistance does not significantly affect the result even when using a constant value, the acceleration resistance or gradient resistance must be calculated. The air resistance also requires an accurate calculation.

Wheel Force (F w )
Conventional methods calculated v·a pos and RPA using the acceleration (a) of the acceleration resistance. However, since this study aimed to evaluate the driving dynamics based on the load applied to the engine, a new characteristic acceleration was calculated by dividing the wheel force (calculated in Equation (2)) according to the relationship of force = mass times acceleration in Equation (6). Here, the wheel force was the force delivered from the engine through the transmission. Therefore, this new acceleration was defined as the engine acceleration (a eng ) and the v·a eng calculated with it can be considered the vehicle's specific power. Figure 6 shows the v·a eng and RPA eng . Since v·a eng considered the air resistance, rolling resistance, and gradient resistance (in addition to the acceleration resistance), the values were generally larger, as is shown in Figure 6a when compared to Figure 3a. In particular, since the air resistance sharply increased as velocity increased, the v·a eng _(95) increased to near the upper limit in the motorway section.

Data Analysis Considering Driving Resistances
±1825 N, respectively. Acceleration resistance showed the largest variation range, even though only positive acceleration was considered. Rolling resistance showed the smallest range. This indicates the impact of the wheel force on the calculation. While the rolling resistance does not significantly affect the result even when using a constant value, the acceleration resistance or gradient resistance must be calculated. The air resistance also requires an accurate calculation.

Wheel Force (Fw)
Conventional methods calculated v·apos and RPA using the acceleration (a) of the acceleration resistance. However, since this study aimed to evaluate the driving dynamics based on the load applied to the engine, a new characteristic acceleration was calculated by dividing the wheel force (calculated in Equation (2)) according to the relationship of force = mass times acceleration in Equation (6). Here, the wheel force was the force delivered from the engine through the transmission. Therefore, this new acceleration was defined as the engine acceleration (aeng) and the v·aeng calculated with it can be considered the vehicle's specific power. Figure 6 shows the v·aeng and RPAeng. Since v·aeng considered the air resistance, rolling resistance, and gradient resistance (in addition to the acceleration resistance), the values were generally larger, as is shown in Figure 6a when compared to Figure 3a. In particular, since the air resistance sharply increased as velocity increased, the v⋅aeng_(95) increased to near the upper limit in the motorway section. This tendency is clearly demonstrated in Figure 7, which shows the v⋅aeng_(mean) for each driving resistance. As was explained above, the v·a by air resistance sharply increased as the vehicle moved from the urban section to the motorway. However, the v·a by acceleration resistance in the urban section was not significantly different than in the motorway section even though it had a multiplied velocity term. This indicates that vehicle acceleration in the urban section was larger than that in the motorway section. This was also different to the v·apos in Figure 3a that used vehicle acceleration, too. In Figure  3a, the v·apos_(mean) increased as the vehicle speed increased. The reason for the difference between the v·a by acceleration resistance and the v·apos is that while the v·apos excluded vehicle acceleration of 0.1 m/s 2 or less, the v·a included even negative acceleration if the aeng was larger than 0.1 m/s 2 . This can also explain the significantly higher data frequency in Figure 6a than in Figure 3a. For the v·a by rolling resistance, the acceleration by This tendency is clearly demonstrated in Figure 7, which shows the v·a eng _(mean) for each driving resistance. As was explained above, the v·a by air resistance sharply increased as the vehicle moved from the urban section to the motorway. However, the v·a by acceleration resistance in the urban section was not significantly different than in the motorway section even though it had a multiplied velocity term. This indicates that vehicle acceleration in the urban section was larger than that in the motorway section. This was also different to the v·a pos in Figure 3a that used vehicle acceleration, too. In Figure 3a, the v·a pos _(mean) increased as the vehicle speed increased. The reason for the difference between the v·a by acceleration resistance and the v·a pos is that while the v·a pos excluded vehicle acceleration of 0.1 m/s 2 or less, the v·a included even negative acceleration if the a eng was larger than 0.1 m/s 2 . This can also explain the significantly higher data frequency in Figure 6a than in Figure 3a. For the v·a by rolling resistance, the acceleration by rolling resistance did not change much as explained above. However, the v·a by rolling resistance increased as velocity increased since the velocity was multiplied. The v·a by gradient resistance was a characteristic of the driving route. As is shown in Figure 1, the rural section had positive v·a values by gradient resistance because it was a gradual uphill section. However, the v·a values were negative in the motorway section because it was a gradual downhill section. rolling resistance did not change much as explained above. However, the v·a by rolling resistance increased as velocity increased since the velocity was multiplied. The v·a by gradient resistance was a characteristic of the driving route. As is shown in Figure 1, the rural section had positive v·a values by gradient resistance because it was a gradual uphill section. However, the v·a values were negative in the motorway section because it was a gradual downhill section. Similar to the v·aeng, the RPAeng showed a tendency to increase when compared to the conventional calculation method. The RPA in Figure 3b was high in the urban section, and it decreased as velocity increased since the number of valid datum points decreased. However, the PRAeng in Figure 6b did not decrease but even increased with velocity increase because it had more valid datum points. Since the v·aeng and the RPAeng are new parameters defined in this study, it was not feasible to apply existing upper and lower limits to them. It is necessary to establish new allowable limits suitable for the new parameters.

Data Analysis Considering Driving Resistances
It was checked how the newly defined v·aeng reflected the actual engine load characteristics (compared to the existing v·apos) by comparing its relationship with CO2 emissions. Since CO2 emissions are proportional to fuel consumption, and higher fuel consumption indicates higher power, checking the relationship between v·a and CO2 is practically the same as checking the relationship between v·a and engine power. Although fuel consumption is not exactly proportional to engine power since the ignition timing and/or injection strategy are dependent on driving conditions, a linear relationship of some degree was expected. Figure 8a shows the relationship between the v·apos and CO2 emissions. As the graph shows, the v·apos increased as the CO2 emission increased and had a linearity level of 0.5244 based on the determination coefficient (R2). As explained in Section 3.2.4, the v·apos that considered only acceleration resistance can show some level of linearity since it had the most impact among the driving resistances on the calculation of wheel force. However, the data became more distributed as CO2 emissions increased, and it became more difficult to identify a correlation. Figure 8b shows the relationship between the v·aeng and CO2 emissions. Compared to Figure 8a, Figure 8b displays a clear linear relationship without data distribution even as CO2 emissions increased. Since the determination coefficient of 0.7815 was closer to 1 than the coefficient of 0.5244 from the conventional calculation, it was determined that the v·aeng reflects engine load characteristics better than the v·apos. Similar to the v·a eng , the RPA eng showed a tendency to increase when compared to the conventional calculation method. The RPA in Figure 3b was high in the urban section, and it decreased as velocity increased since the number of valid datum points decreased. However, the PRA eng in Figure 6b did not decrease but even increased with velocity increase because it had more valid datum points. Since the v·a eng and the RPA eng are new parameters defined in this study, it was not feasible to apply existing upper and lower limits to them. It is necessary to establish new allowable limits suitable for the new parameters.
It was checked how the newly defined v·a eng reflected the actual engine load characteristics (compared to the existing v·a pos ) by comparing its relationship with CO 2 emissions. Since CO 2 emissions are proportional to fuel consumption, and higher fuel consumption indicates higher power, checking the relationship between v·a and CO 2 is practically the same as checking the relationship between v·a and engine power. Although fuel consumption is not exactly proportional to engine power since the ignition timing and/or injection strategy are dependent on driving conditions, a linear relationship of some degree was expected. Figure 8a shows the relationship between the v·a pos and CO 2 emissions. As the graph shows, the v·a pos increased as the CO 2 emission increased and had a linearity level of 0.5244 based on the determination coefficient (R 2 ). As explained in Section 3.2.4, the v·a pos that considered only acceleration resistance can show some level of linearity since it had the most impact among the driving resistances on the calculation of wheel force. However, the data became more distributed as CO 2 emissions increased, and it became more difficult to identify a correlation. Figure 8b shows the relationship between the v·a eng and CO 2 emissions. Compared to Figure 8a, Figure 8b displays a clear linear relationship without data distribution even as CO 2 emissions increased. Since the determination coefficient of 0.7815 was closer to 1 than the coefficient of 0.5244 from the conventional calculation, it was determined that the v·a eng reflects engine load characteristics better than the v·a pos .

Driving Dynamics during Braking
So far, the process has focused on analyzing the driving dynamics while the driver steps on the accelerator. All driving resistances were considered to calculate wheel force, and only positive wheel forces were used in the driving dynamics analysis. Driving dynamics calculated through these procedures were used to determine whether the on-road driving test was under severe or soft testing conditions. Thus, it was possible to rule out these unusual test cases that could lead to wrong exhaust emission results. However, in reality, braking also affects exhaust emissions. For example, severe braking dissipates significant kinetic energy, so a large amount of fuel is required to compensate for it. Therefore, it is necessary to analyze the braking dynamics, in addition to the driving dynamics, to accurately evaluate exhaust emissions.

Driving Dynamics during Braking
So far, the process has focused on analyzing the driving dynamics while the driver steps on the accelerator. All driving resistances were considered to calculate wheel force, and only positive wheel forces were used in the driving dynamics analysis. Driving dynamics calculated through these procedures were used to determine whether the on-road driving test was under severe or soft testing conditions. Thus, it was possible to rule out these unusual test cases that could lead to wrong exhaust emission results. However, in reality, braking also affects exhaust emissions. For example, severe braking dissipates significant kinetic energy, so a large amount of fuel is required to compensate for it. Therefore, it is necessary to analyze the braking dynamics, in addition to the driving dynamics, to accurately evaluate exhaust emissions.
It was possible to analyze the driving dynamics during braking by using negative wheel force out of the wheel force calculated above. Since the characteristic acceleration (calculated by dividing the negative wheel force by the vehicle mass) is negative, this study defined the absolute value of this acceleration as the negative acceleration (aneg). The graph in Figure 9 demonstrates the v·aneg. As is shown in the graph, the number of v⋅aneg data was large in the urban section because both acceleration and deceleration frequently occur in urban driving. On the other hand, the number of datum points was relatively small in the rural and motorway sections because the brake was rarely used. In addition, most data in these sections remained below 10 W/kg, which means that the aneg values were very small. It was possible to analyze the driving dynamics during braking by using negative wheel force out of the wheel force calculated above. Since the characteristic acceleration (calculated by dividing the negative wheel force by the vehicle mass) is negative, this study defined the absolute value of this acceleration as the negative acceleration (a neg ). The graph in Figure 9 demonstrates the v·a neg . As is shown in the graph, the number of v·a neg data was large in the urban section because both acceleration and deceleration frequently occur in urban driving. On the other hand, the number of datum points was relatively small in the rural and motorway sections because the brake was rarely used. In addition, most data in these sections remained below 10 W/kg, which means that the a neg values were very small.

Driving Dynamics during Braking
So far, the process has focused on analyzing the driving dynamics while the driver steps on the accelerator. All driving resistances were considered to calculate wheel force, and only positive wheel forces were used in the driving dynamics analysis. Driving dynamics calculated through these procedures were used to determine whether the on-road driving test was under severe or soft testing conditions. Thus, it was possible to rule out these unusual test cases that could lead to wrong exhaust emission results. However, in reality, braking also affects exhaust emissions. For example, severe braking dissipates significant kinetic energy, so a large amount of fuel is required to compensate for it. Therefore, it is necessary to analyze the braking dynamics, in addition to the driving dynamics, to accurately evaluate exhaust emissions.
It was possible to analyze the driving dynamics during braking by using negative wheel force out of the wheel force calculated above. Since the characteristic acceleration (calculated by dividing the negative wheel force by the vehicle mass) is negative, this study defined the absolute value of this acceleration as the negative acceleration (aneg). The graph in Figure 9 demonstrates the v·aneg. As is shown in the graph, the number of v⋅aneg data was large in the urban section because both acceleration and deceleration frequently occur in urban driving. On the other hand, the number of datum points was relatively small in the rural and motorway sections because the brake was rarely used. In addition, most data in these sections remained below 10 W/kg, which means that the aneg values were very small. The overall scale of the v·a under braking conditions was not much different from the v·a eng or the v·a pos scales. The graph in Figure 10 compares the scales of the v·a eng and the v·a neg for each section. The hollow marks in the graph represent the data of the on-road driving test case of this study. The graph also shows the v·a pos for reference. The filled marks represent the data of the additional nine on-road driving test cases (for a total of 10 test cases) of various routes in addition to the route shown in Figure 1 (for a total of three driving routes). All test cases used the same vehicle. Figure 10b demonstrates that while the v·a eng _(mean) and the v·a neg _(mean) values were almost the same in the urban section, the v·a neg _(mean) values were much smaller than the v·a eng _(mean) values in the motorway section. This is due to the fact that the driver was forced to continuously press on the accelerator pedal to maintain the velocity on the motorway since the air and rolling resistances were large. In contrast, there were few cases where the driver was required to step on the brake pedal since the traffic flow was smooth on the motorway. Figure 10a shows the 95th quartile data and indicates that the v·a eng _(95) and the v·a neg _(95) values were similar to each other on the motorway. This is a kind of error caused by the lack of valid v·a neg data in the motorway section. The amount of data from the motorway section shown in Figure 6a was large, but the amount in Figure 9 was relatively small. In addition, most of the v·a neg data were distributed below 10 W/kg. However, some data had a high value of 25 W/kg or greater, since the motorway was a high-velocity section and the v·a neg values increased sharply even when the a neg values only increased slightly. Accordingly, the 95th quartile data were overestimated by the small number of large values, although the average value was low. Due to this error, the v·a eng _(95) and the v·a neg _(95) values were almost the same. For the same reason, the deviation between the tests of the v·a neg _(95) in the rural and motorway sections (the filled marks on the graph) was also very large. In the case of a eng _(95), the data of all tests were concentrated within the range of ±5 W/kg, regardless of the driving route. However, the datum points of the v·a neg _(95) were distributed over a wide range of up to ±18 W/kg. driving test case of this study. The graph also shows the v·apos for reference. The filled marks represent the data of the additional nine on-road driving test cases (for a total of 10 test cases) of various routes in addition to the route shown in Figure 1 (for a total of three driving routes). All test cases used the same vehicle. Figure 10b demonstrates that while the v·aeng_(mean) and the v·aneg_(mean) values were almost the same in the urban section, the v·aneg_(mean) values were much smaller than the v·aeng_(mean) values in the motorway section. This is due to the fact that the driver was forced to continuously press on the accelerator pedal to maintain the velocity on the motorway since the air and rolling resistances were large. In contrast, there were few cases where the driver was required to step on the brake pedal since the traffic flow was smooth on the motorway. Figure 10a shows the 95th quartile data and indicates that the v·aeng_(95) and the v·aneg_(95) values were similar to each other on the motorway. This is a kind of error caused by the lack of valid v·aneg data in the motorway section. The amount of data from the motorway section shown in Figure 6a was large, but the amount in Figure 9 was relatively small. In addition, most of the v⋅aneg data were distributed below 10 W/kg. However, some data had a high value of 25 W/kg or greater, since the motorway was a high-velocity section and the v·aneg values increased sharply even when the aneg values only increased slightly. Accordingly, the 95th quartile data were overestimated by the small number of large values, although the average value was low. Due to this error, the v·aeng_(95) and the v·aneg_(95) values were almost the same. For the same reason, the deviation between the tests of the v·aneg_(95) in the rural and motorway sections (the filled marks on the graph) was also very large. In the case of aeng_(95), the data of all tests were concentrated within the range of ±5 W/kg, regardless of the driving route. However, the datum points of the v·aneg_(95) were distributed over a wide range of up to ±18 W/kg. These deviations in braking dynamics between tests are also presented in the average values. In the urban section, drivers must use a brake to decelerate. In rural or motorway driving, on the other hand, some drivers may prefer to slow down naturally through inertial driving without using brakes. Because of these differences in braking patterns, braking dynamics showed the larger deviation than driving dynamics. Although the current RDE standard evaluates driving validity based only on acceleration, not controlling for These deviations in braking dynamics between tests are also presented in the average values. In the urban section, drivers must use a brake to decelerate. In rural or motorway driving, on the other hand, some drivers may prefer to slow down naturally through inertial driving without using brakes. Because of these differences in braking patterns, braking dynamics showed the larger deviation than driving dynamics. Although the current RDE standard evaluates driving validity based only on acceleration, not controlling for the braking dynamics can cause inaccurate evaluation of exhaust emissions because of the braking pattern differences.
The total energy (v·a_(total)) consumed by driving and braking in each section of the on-road road driving test was calculated by adding all data in Figures 6a and 9. The unit of an individual datum point was (W/kg), but since each datum point was acquired at 1-s intervals, the unit of their sum can be considered as (kJ/kg). The graph in Figure 11 shows the total energy consumed for driving and braking in each section, as well as the ratio of the braking energy to the driving energy. In Figure 11a, v·a eng _(total) is the total energy consumed for driving in each section. It was 8.6 kJ/kg in the urban section, 10.4 kJ/kg in the rural section, and 9.8 kJ/kg in the motorway section. The large deviation between the tests was due to driving route differences. The total energy required for each driving route was similar. For example, the two lowest v·a eng _(total) values in the urban and rural sections were from the same driving route. On the other hand, v·a neg _(total) was the total energy lost by braking. It showed a tendency to gradually decrease as the vehicle velocity increased; 4.8 kJ/kg in the urban section, 1.9 kJ/kg in the rural section, and 1.1 kJ/kg in the motorway section. This was due to the small number of valid datum points in the motorway section. A driver was less likely to brake on a motorway than in an urban section, as confirmed in Figure 9. Therefore, the amount of energy lost by braking was less on the motorway. tions were from the same driving route. On the other hand, v·aneg_(total) was the total energy lost by braking. It showed a tendency to gradually decrease as the vehicle velocity increased; 4.8 kJ/kg in the urban section, 1.9 kJ/kg in the rural section, and 1.1 kJ/kg in the motorway section. This was due to the small number of valid datum points in the motorway section. A driver was less likely to brake on a motorway than in an urban section, as confirmed in Figure 9. Therefore, the amount of energy lost by braking was less on the motorway.
As is shown in Figure 11b, the ratios of the energy lost by braking to the energy used for driving were 54%, 18%, and 7% in the urban, rural, and motorway sections, respectively. The deviations between the 10 test cases were ±5% or less. For the urban driving, the deviation between the tests from 49% to 57% is very small. However, for the motorway driving, the data distribution from 2% to 12% is relatively large compared to the values. The acceptable range of deviation between tests is not known yet. In the future work, the relation between the exhaust emission and braking dynamics will be analyzed and the acceptable range of the deviation will also be clearly defined.

Conclusions
The present study conducted on-road road driving tests with a passenger gasoline vehicle according to RDE standards, and driving data were acquired with the PEMS. The data analysis process applied a new method of analyzing driving dynamics that considered all driving resistances to improve the conventional driving dynamics analysis method that utilizes acceleration only. This study also proposed a method of analyzing driving dynamics during the braking process. The conclusions drawn from this study are presented below. As is shown in Figure 11b, the ratios of the energy lost by braking to the energy used for driving were 54%, 18%, and 7% in the urban, rural, and motorway sections, respectively. The deviations between the 10 test cases were ±5% or less. For the urban driving, the deviation between the tests from 49% to 57% is very small. However, for the motorway driving, the data distribution from 2% to 12% is relatively large compared to the values. The acceptable range of deviation between tests is not known yet. In the future work, the relation between the exhaust emission and braking dynamics will be analyzed and the acceptable range of the deviation will also be clearly defined.

Conclusions
The present study conducted on-road road driving tests with a passenger gasoline vehicle according to RDE standards, and driving data were acquired with the PEMS. The data analysis process applied a new method of analyzing driving dynamics that considered all driving resistances to improve the conventional driving dynamics analysis method that utilizes acceleration only. This study also proposed a method of analyzing driving dynamics during the braking process. The conclusions drawn from this study are presented below.

•
The newly calculated v·a eng values are generally larger than the existing v·a pos values because they reflect air, rolling, and gradient resistances in addition to acceleration resistance. For the same reason, the RPA values are also larger overall when compared to the conventional method. Therefore, it is necessary to establish new allowable limit criteria to evaluate the validity of the new driving dynamics.

•
Comparing the correlation between the existing v·a pos and the newly calculated v·a eng with CO 2 emissions, the v·a eng values show a clearer linear relationship, indicating that v·a eng reflects the engine load characteristics better than v·a pos does. This signifies that v·a eng is a more suitable criterion for determining the validity of driving data.

•
The deviation of braking dynamics between tests is particularly large in the motorway due to the differences in braking patterns. Although the current RDE standard evaluates the driving validity based on positive acceleration only, it is necessary to establish a regulation on braking patterns since the braking pattern also can affect the exhaust emissions.

•
The ratio of the total energy lost by braking to the total energy consumed by driving in urban, rural, and motorway areas of each driving section was about 54%, 18%, and 7%, respectively. The deviation between tests was not larger than ±5%. However, the deviation was relatively large compared to the value in the motorway. Therefore, it is necessary to calculate the acceptable energy ratio range by analyzing the relation between exhaust emissions and braking dynamics. Funding: This research received no external funding.