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·apos 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·aeng) 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·aneg) 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).
3.1. 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 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.
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 3a defines the v·a
pos corresponding to the top 95% in each section (urban, rural, and motorway) as v·a
pos_(95). It defines the average v·a
pos in each window as v·a
pos_(mean). Each small gray dot represents individual v⋅a
pos values, and the solid black line is the highest line that v·a
pos_(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 (=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.
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.
3.2. 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
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.
3.2.1. Air Resistance (Ra)
The air resistance applied to a vehicle is calculated as follows:
where
Cd: drag coefficient
A: full projected area (m2)
ρ: air density (kg/m3)
v: vehicle speed (m/s)
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.
3.2.2. 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].
where
Csr = −514.7(Ptire/100)3 + 53.72(Ptire/100)2 − 1.877(Ptire/100) + 0.03051
Cdr = −793.1(Ptire/100)3 + 83.98(Ptire/100)2 − 2.977(Ptire/100) + 0.03759
v: vehicle speed (km/h)
W: vehicle weight (N)
θ: road gradient (deg)
Ptire: inflation pressure (bar)
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 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.
3.2.3. Gradient Resistance (Rg)
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].
3.2.4. Acceleration Resistance (Ra)
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/s2)
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.
3.2.5. 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.
3.3. Data Analysis Considering Driving Resistances
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.
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.
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.
3.4. 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 (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.
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.
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
Figure 6a and
Figure 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.
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.