.. highlightlang:: us .. index:: PEMS .. index:: RDE (Real Driving Emission) .. index:: ISC (In-Service-Conformity) .. index:: y^pems-evaluation .. _in_service_conformity_data_evaluation: Real Driving Emission Data Evaluation ===================================== .. note:: The addin is a beta version. The release version is planned for 2021. The Addin calculates the emission conformity factor from data measured by a PEMS systen during an on-the-road test. A dry/wet correction after ISO 16183:2002(E) and a gas delay correction can be applied to the data. The data from the pems device and the ECU must be time aligned and merged into one file. The following channels are mandatory: .. list-table:: :header-rows: 1 * - Name - Alias - Unit * - Time - time_s - s * - Engine Speed - speed_rpm - 1/min * - Engine Torque - torque_nm - Nm * - Torque Reference - Torque_Reference - Nm * - Torque Percent - Torque_Percent - % * - Torque Friction Loss - Torque_Friction - % * - Ambient Temperature - ambient_temperature_grdC - °C * - Ambient Pressure - ambient_pressure_hpa - hPa * - Exhaust gas flow - ExhaustMassFlow__kgs - kg/s * - Engine coolant temperature - CoolantTemperature - °C * - Fuel Mass Flow - fuel_mass_flow_kg_s - kg/s The following channels are optional: .. list-table:: :header-rows: 1 * - Name - Alias - Unit * - Valid Data Point Flag (valid=1) - ValidPoint - 0/1-Signal * - Vehicle Speed - VehicleSpeed - km/h The Addin uses the method "Moving Averaging Window Principle" (as used by EMROAD) to determine the conformity factors. The emission data is integrated using a moving average window, based on the reference engine work or CO2 mass. The data must be measured with a frequency of at least 1 Hz. If the data was measured with a higher frequency, it will be resampled to 1 Hz. Data in time gaps will be set to zero and ignored by the MAW calculation. Dry/Wet correction ISO 16183:2002(E) ------------------------------------ The ``Dry/wet correction`` in the main dialog box is executed during the data import. Da dialog box will be displayed to select the mandatory channels. Test Data Verification ---------------------- The environment temperature is checked against the lower limit of -7 °C and the upper limit. The data evaluation will start after the coolant temperature has reached 30 °C for the first time or after the coolant temperature is stabilised within +/– 2°C over a period of 5 minutes whichever comes first but no later than 10 minutes after engine start. The valid windows are the windows whose average power exceeds the power threshold of 20 % of the maximum engine power. It is checked, if the percentage of valid windows is equal or greater than 50 %. Gas Delay Times --------------- Installation ------------ * Choose :ref:`tools-add-in-manager`. * Check the "Real Driving Emissions" addin. * Choose OK. .. _usage_isc: Usage ----- * Start UniPlot. * Choose RDE Evaluation=>Configuration to select emission channels: .. image:: S:/uniplot-obj/images/isc-config-dialog-en.* * Choose ISC Evaluation=>Evalution to start a new evalutation: .. image:: S:/uniplot-obj/images/isc-main-dialog-en.* Data Selection -------------- The data can be imported from all supported data formats. The unit converter is used to convert the units used for the RDE calculation. Channel Selection ----------------- The channels mandatory for the calculation can be selected using the following dialog box: .. image:: S:/uniplot-obj/images/isc-channel-dialog-en.* Calculation ----------- All calculated results will be written into the result data file. Its name ist displayed in the Event Window. A double click loads the file into the data browser. The following channels will be added to the data file: .. list-table:: :header-rows: 1 * - Name - Unit - Description * - e.Trip_Work_int - kWh - Total trip work. * - e.Trip_Work_int_a - kWh - Total trip work of valid data points. * - e.AmbientTemp_max_limit - °C - Calculated maximum valid ambient temperature. * - e.MAW_Time - s - Lenght of time of each moving average window. * - e.MAW_Power_percent - % - Percent of each moving average window of the specified max. power. * - e.MAW_start_time - s - Start time of each moving average window. * - e.MAW_end_time - s - End time of each moving average window. * - e.MAW_duration_time - s - Length of each moving average window * - e.MAW_invalid_points - - Number of invalid points in each moving average window. * - e.MAW_XXX - g/kWh or #/kWh - Emission values for each moving average window. * - e.CF_all_XXX - - Conformity factor for each moving average window for all emissions. * - e.CF_XXX_Distribution_y - - Sorted conformity factor for each moving average window for all emissions. * - e.CF_XXX_Distribution_y_unsorted - - Unsorted conformity factor for each moving average window for all emissions. * - e.CF_Distribution_x - % - X-Coordinates for ``e.CF_XXX_Distribution_y`` and ``e.CF_XXX_Distribution_y_unsorted``. * - e.vehicle_speed - km/h - Vehicle velocity * - e.Trip_Cummulativ - km - trip distance cummulativ * - e.Trip_Fuel_Consumption - g - Fuel consumption for the trip * - e.XXX - ppm or # - Copy of specified emission channel. * - e.Trip_flow_XXX - #/s - Particle number * - e.Trip_int_XXX - # - Particle number integral. * - e.Trip_mass_flow_XXX - g/s - Emission mass. * - e.Trip_mass_int_XXX - g - Emission mass integral. * - e.torque - Nm - Torque channel (copy of the specified channel). * - e.Time - s - time channel (copy of the specified channel). * - e.speed - 1/min - engine speed (copy of the specified channel). * - e.exhaust_mass_flow - kg/s - exhaust mass flow (copy of the specified channel). * - e.coolant_t - °C - Coolant temperatur (copy of the specified channel). * - e.ambient_t - °C - Ambient temperatur (copy of the specified channel). * - e.ambient_p - kPa - Ambient Pressure (copy of the specified channel). * - e.valid - - Valid points (0 or 1) (copy of the specified channel). * - e.Power - kW - Calculated power * - e.Power_percent - % - Power in percent (power / power_max * 100). The following attribute values will be added to the data file: (XXX stands for emission channel) .. list-table:: :header-rows: 1 * - Name - Unit - Description * - Name - Einheit - Beschreibung * - e.work_window_limit - kWh - Reference Work, specified in the main dialog box. * - e.coolant_min_limit - °C - Minimum coolant temperature (Default: -7 °C). * - e.power_max - kW - Maximum Power, specified in the main dialog box. * - e.percentile - % - Percentile value. * - e.power_limit_percent - % - Moving Average Percentil for which the emission values are calculated. * - e.ugas_XXX - - UGas values for the different emissions (XXX = NOx, CO2, etc.) * - e.limit_XXX - Unit of XXX - Limits for the given emissions Test cycle (WHTC). (XXX = NOx, CO2, etc.) * - e.Trip_Duration - s - Tripe Time for the complete trip from engine start to stop (engine speed > 100 rpm at the beginning and end of test). * - e.Trip_AV_Speed - 1/min - Average engine speed for the complete trip. * - e.Trip_AV_Torque - Nm - Average torque for the complete trip. * - e.Trip_AV_Power - kW - Average power for the complete trip. * - e.Trip_Work - kWh - Work for the complete trip. * - e.Trip_Work_a - kWh - Work for the complete trip for all valid points (_a: calibration during test is ignored). * - e.AmbientTemp_invalid_points - - Number of invalid points, (T_ambient < Tmin OR T_ambient > Tmin). * - e.AmbientTemp_min_limit - °C - Ambient minimum temperature (T_min), Tmax is calculated: -0.4514 .* (101.3 - rvAP) + (311-273.15)), rvAP is the ambient pressure in kPa. * - e.MAW_Work - kWh - Work of all moving average windows * - e.MAW_Duration - s - Time of all moving average windows. * - e.Trip_MAW_AV_speed - 1/min - Average engine speed for all moving average windows. * - e.Trip_MAW_AV_torque - Nm - Average engine speed for all moving average windows. * - e.Trip_MAW_AV_power - kW - Average power for all moving average windows. * - e.Trip_MAW_AV_coolant_t - °C - Average coolant temperature for all moving average windows. * - e.Trip_MAW_AV_ambient_t - °C - Average ambient temperature for all moving average windows. * - e.Trip_MAW_AV_ambient_p - kPa - Average ambient pressure for all moving average windows. * - e.Trip_MAW_AV_XXX - Unit of XXX - Average emission for all moving average windows. (XXX = NOx, CO2, etc.) * - e.Q90_XXX - g or #/kWh for particle in window - Quantile value for emission XXX at the specified percentil value, default value 90%. (e.percentile) * - e.CF_Q90_XXX - - Conformity factor for emission XXX at the specified percentil value. * - e.MAW_valid_windows - - Number of valid windows. * - e.MAW_all_windows - - Number of windows. * - e.MAW_T70_Windows - - Number of windows where coolant temperatur is above 70 °C. * - e.MAW_valid_windows_percent - % - Number of valid windows in percent (nValidWindows/nWindows * 100). * - e.MAW_valid_Power_Min_percent - % - minimum power of all valid windows. * - e.MAW_valid_Power_Max_percent - % - maximum power of all valid windows. * - e.speed_limit_zero - km/h - vehicle velocity limit for zero velocity (Default: 2 km/h) * - e.speed_limit_urban - km/h - vehicle speed limit for urban driving (Default: 50 km/h) * - e.speed_limit_interurban - km/h - vehicle speed limit for inter urban driving (Default: 75 km/h) * - e.Trip_Zero - % - Percent of trip time with zero vehicle velocity. * - e.Trip_Urban - % - Percent of trip time with urban vehicle velocity. * - e.Trip_InterUrban - % - Percent of trip time with inter urban vehicle velocity. * - e.Trip_Highway - % - Percent of trip time with highway vehicle velocity (v > v_inter_urban). * - e.Trip_Velocity_Mean - km/h - Average vehicle velocity. * - e.Trip_Vehicley_Max - km/h - Maximum vehicle velocity for the complete trip. * - e.Trip_Distance - km - Trip distance. * - e.Trip_Fuel_Mass - g - complete fuel consumption for the trip. * - e.Trip_bsfc - g/kWh - Specific fuel consumption for the trip. * - e.ColdStartDuration - s - Time at the test start until the coolant temperatur is valid. * - e.Trip_AV_XXX - Unit of XXX - Average of all valid data points. * - e.Trip_total_XXX - ppm of # - Total emissions for the trip. * - e.Trip_spez_XXX" - g/kWh or #/kWh - Specific emission for alle valid values (a). * - e.Trip_spez_distance_XXX - g/km or #/km - Distance specific emission for all valid values (a). * - e.Trip_spez_fuel_XXX - g/kg or #/kg - Fuel consumption specific emission for all valid values (a). * - e.Trip_mass_XXX - g - Complete emission mass. (XXX = NOx, CO2, etc.) * - e.channel.XXX - - Emission channel name. * - e.Trip_Duration - s - Tripe Time for the complete trip from engine start to engine stop (engine speed > 100 rpm). Template -------- The template will be filled with the data from the selected datafile. The result file contains all data of the selected file and the result data. The result data and attributes starts with the prefix ``e.``. If global attributes from the NC file are used to be displayed in the IPZ document via placeholder, the period (.) must be replaced by an underscore. .. code-block:: none Attribute: e.CF_Final_NOx --> Placeholder: $e_CF_Final_NOx$ An example file can be found in the addin rs_pems folder. Programming ----------- The source code can be found in ``\\addin\rs_pems``:: def _pems_calculate_test_example1() { obj = _pems_create(); obj.ssNCFile = "s:/projects/pems/pems_example.nc2" if (obj.ssNCFile == "") { return FALSE; } obj.rsWork_Ref_kWh = 30; obj.rsCoolant_Limit_grdC = 70; obj.rsPower_Max_kW = 250; obj.rsPower_Ref_percent = 20; obj.rsPercentile = 90; obj.uGasValues.NOx = 0.001587; obj.limit.Nox = 0.46; obj.oVar["time_s"] = "Time"; obj.oVar["speed_rpm"] = "iENG_SPEED"; obj.oVar["torque_nm"] = ""; obj.oVar["Torque_Reference"] = "sREF_ENG_TORQ"; obj.oVar["Torque_Percent"] = "iPCNT_TORQUE"; obj.oVar["Torque_Friction"] = "iFRICT_TORQUE"; obj.oVar["CoolantTemperature"] = "iCOOL_TEMP"; obj.oVar["NOx"] = "ikNOxzw"; obj.oVar["ExhaustMassFlow__kgs"] = "icMASS_FLOW"; obj.oVar["ambient_temperature_grdC"] = "iSCB_LAT"; obj.oVar["ambient_pressure_hpa"] = "iSCB_LAP"; return _pems_calculate_main(obj); } .. us.history **History** .. list-table:: :header-rows: 1 * - Version - Description * - R2017.0 - New. * - R2020.6 (0.2.5) - The evaluation has been adapted to "COMMISSION REGULATION (EU) 2019/1939 of 7 November 2019" The avarage value windows are now determined in forward direction. The start of the evaluation is adjusted according to the regulation 2019/1939. The CF factor is calculated from the cold phase and hot phase in the ratio 14 to 86. The dialog box has been modified: The windows can now be calculated using the work or the cumulated CO2 mass emission. :sub:`id-1162246`