Appendix 7a Verification of overall trip dynamics
1.INTRODUCTION
This Appendix describes the calculation procedures to verify the overall trip dynamics, to determine the overall excess or absence of dynamics during urban, rural and motorway driving.
2.SYMBOLS, PARAMETERS AND UNITS
RPARelative Positive Acceleration
—
difference
—
larger
—
larger or equal
—
per cent
—
smaller
—
smaller or equal
—
acceleration [m/s2]
—
acceleration in time step i [m/s2]
—
positive acceleration greater than 0,1 m/s2 [m/s2]
—
positive acceleration greater than 0,1 m/s2 in time step i considering the urban, rural and motorway shares [m/s2]
—
acceleration resolution [m/s2]
—
distance covered in time step i [m]
—
distance covered in time step i considering the urban, rural and motorway shares [m]
—
discrete time step
—
discrete time step of positive acceleration datasets
—
refers to the respective category (t=total, u=urban, r=rural, m=motorway)
—
number of samples for urban, rural and motorway shares with positive acceleration greater than 0,1 m/s2
—
total number of samples for the urban, rural and motorway shares and the complete trip
—
relative positive acceleration for urban, rural and motorway shares [m/s2 or kWs/(kg*km)]
—
duration of the urban, rural and motorway shares and the complete trip [s]
—
compound data smoother
—
vehicle speed [km/h]
—
actual vehicle speed in time step i [km/h]
—
actual vehicle speed in time step i considering the urban, rural and motorway shares [km/h]
—
actual vehicle speed per positive acceleration greater than 0,1 m/s2 in time step j considering the urban, rural and motorway shares [m2/s3 or W/kg].
—
95th percentile of the product of vehicle speed per positive acceleration greater than 0,1 m/s2 for urban, rural and motorway shares [m2/s3 or W/kg]
3.TRIP INDICATORS
3.1. Calculations
3.1.1. Data pre-processing
Dynamic parameters like acceleration, or RPA shall be determined with a speed signal of an accuracy of 0,1 % for all speed values above 3 km/h and a sampling frequency of 1 Hz. This accuracy requirement is generally fulfilled by signals obtained from a wheel (rotational) speed sensor.
The speed trace shall be checked for faulty or implausible sections. The vehicle speed trace of such sections is characterised by steps, jumps, terraced speed traces or missing values. Short faulty sections shall be corrected, for example by data interpolation or benchmarking against a secondary speed signal. Alternatively, short trips containing faulty sections could be excluded from the subsequent data analysis. In a second step the acceleration values shall be calculated and ranked in ascending order, as to determine the acceleration resolution .
If , the vehicle speed measurement is sufficiently accurate.
If , data smoothing by using a T4253H Hanning filter shall be performed
The T4235 Hanning filter performs the following calculations: The smoother starts with a running median of 4, which is centred by a running median of 2. The filter then re-smoothes these values by applying a running median of 5, a running median of 3, and hanning (running weighted averages). Residuals are computed by subtracting the smoothed series from the original series. This whole process is then repeated on the computed residuals. Finally, the smoothed final speed values are computed by summing up the smoothed values obtained the first time through the process with the computed residuals.
The correct speed trace builds the basis for further calculations and binning as described in paragraph 8.1.2.
3.1.2. Calculation of distance, acceleration and
The following calculations shall be performed over the whole time based speed trace (1 Hz resolution) from second 1 to second tt (last second).
The distance increment per data sample shall be calculated as follows:
where:
is the distance covered in time step i [m]
is the actual vehicle speed in time step i [km/h]
is the total number of samples
The acceleration shall be calculated as follows:
where:
The product of vehicle speed per acceleration shall be calculated as follows:
where:
3.1.3. Binning of the results
After the calculation of ai and , the values vi , di , ai and shall be ranked in ascending order of the vehicle speed.
All datasets with belong to the ‘urban’ speed bin, all datasets with belong to the ‘rural’ speed bin and all datasets with belong to the ‘motorway’ speed bin.
The number of datasets with acceleration values shall be bigger or equal to 150 in each speed bin.
For each speed bin the average vehicle speed shall be calculated as follows:
where:
is the total number of samples of the urban, rural, and motorway shares.
3.1.4. Calculation of per speed bin
The 95th percentile of the values shall be calculated as follows:
The values in each speed bin shall be ranked in ascending order for all datasets with and the total number of these samples Mk shall be determined.
Percentile values are then assigned to the values with as follows:
The lowest value gets the percentile 1/Mk , the second lowest 2/Mk , the third lowest 3/Mk and the highest value
is the value, with If cannot be met, shall be calculated by linear interpolation between consecutive samples j and j+1 with and .
The relative positive acceleration per speed bin shall be calculated as follows:
where:
is the relative positive acceleration for urban, rural and motorway shares in [m/s2 or kWs/(kg*km)]
is a time difference equal to 1 second
is the sample number for urban, rural and motorway shares with positive acceleration
is the total sample number for urban, rural and motorway shares
4.VERIFICATION OF TRIP VALIDITY
4.1.1. Verification of per speed bin (with v in [km/h])
and
is fulfilled, the trip is invalid.
If and is fulfilled, the trip is invalid.