Skip welcome & menu and move to editor
Welcome to JS Bin
File
Add library
Share
Login or Register
Blog
Help
Load cached copy from
  • Run with JS
  • Save
  • Share
  • …
 
<!DOCTYPE html>
<meta charset="utf-8">
<head>
<meta name="description" content="[DieselBooster (Arduino) simulation with basic interpolation turned on]">
  <style>
    .xxxline {
      fill: none;
      stroke: steelblue;
      stroke-width: 2px;
    }
​
    .xxxline2 {
      fill: none;
      stroke: red;
      stroke-width: 2px;
    }
​
    pre {
      margin: 0;
      display: flex;
    }
​
    #chart0 svg {
      height: 300px;
    }
  </style>
​
  <link rel="stylesheet" type="text/css" href="https://cdnjs.cloudflare.com/ajax/libs/nvd3/1.8.5/nv.d3.min.css"/>
​
  <!--script src="https://d3js.org/d3.v4.min.js"></script-->
  <script src="https://cdnjs.cloudflare.com/ajax/libs/d3/3.5.2/d3.min.js" charset="utf-8"></script>
  <script src="https://cdnjs.cloudflare.com/ajax/libs/nvd3/1.8.5/nv.d3.js"></script>
  <script src="https://cdnjs.cloudflare.com/ajax/libs/underscore.js/1.8.3/underscore-min.js"></script>
​
  <script>
    /** helpers **/
    var chartNo = 0;
    function chartData(data, xLabel, yLabel, height = 300, multiAxis = false, yLabel2 = '') {
      var chartElement = document.createElement("div");
      var id = 'chart' + chartNo++;
      chartElement.setAttribute('id', id);
      chartElement.setAttribute('style', 'height: ' + height + 'px;');
      document.getElementById('body').appendChild(chartElement);
      //chartElement.appendChild(document.createElement("svg"));
      _.each(data, function (item, index) {
        item.yAxis = index + 1;
        item.type = 'line';
      });
​
      nv.addGraph(function () {
        var origin = multiAxis ? nv.models.multiChart() : nv.models.lineChart()
        var chart = origin
                .margin({left: 70, right: 70})  //Adjust chart margins to give the x-axis some breathing room.
                .useInteractiveGuideline(true)  //We want nice looking tooltips and a guideline!
        //.showLegend(true)       //Show the legend, allowing users to turn on/off line series.
        //.showYAxis(true)        //Show the y-axis
        //.showXAxis(true)        //Show the x-axis
            ;
​
        chart.xAxis     //Chart x-axis settings
            .axisLabel(xLabel)
            .tickFormat(d3.format('.02f'));
​
        if (!chart.yAxis1) {
          chart.yAxis1 = chart.yAxis;
        }
        chart.yAxis1     //Chart y-axis settings
            .axisLabel(yLabel)
            .tickFormat(d3.format('.02f'));
​
        if (multiAxis) {
          chart.yAxis2
              .axisLabel(yLabel2)
              .tickFormat(d3.format(',.02f'));
        }
        /* Done setting the chart up? Time to render it!*/
        //var data = sinAndCos();   //You need data...
​
        d3.select('#' + id).append('svg')    //Select the <svg> element you want to render the chart in.
            .datum(data)         //Populate the <svg> element with chart data...
            .transition().duration(350)  //how fast do you want the lines to transition?
            .call(chart);          //Finally, render the chart!
​
        //Update the chart when window resizes.
        nv.utils.windowResize(function () {
          chart.update()
        });
        return chart;
      });
    }
​
    function print(text) {
      //document.getElementById('printout').innerHTML += text;
      var pre = document.createElement("pre");
      pre.innerHTML = text;
      document.getElementById('body').appendChild(pre);
    }
    function println() {
      document.getElementById('body').appendChild(document.createElement("BR"));
    }
    /** end helpers **/
​
​
    /** car **/
    function originalSensorReferenceOnly() {
      // at min pressure voltage is 0.5, at max pressure voltage is 5V; at idle (~4400psi) is ~1V; 5800psi=1.39V
      // max measured rail pressure: 15200 psi = 1048 bar (104.8MPa) -> max sensor pressure 1500 bar?
      // see http://rb-aa.bosch.com/boaasocs/index.jsp;jsessionid=A20738712FCFB9E51CA919DD7D2F9E91.sundoro2?ccat_id=275&prod_id=516
      // C1 = 0.8 * p / PN (or 0.12 * p / PN)
      // Ua = (C1 + C) * Us
      const C0 = 0.1;
      const C1F = 0.8;
      const US = 16; //Volts, supply voltage max (before pullup?)
      const UV = 5;// Volts, supply voltage
​
      function pressureToVoltage(pressureMPa) {
        var C1 = (C1F) * pressureMPa / PN;
        var UA = (C1 + C0) * US;
        return UA / 5;
      }
    }
​
    const PN = 180; // nominal max sensor pressure (MPa)
    const BAR_TO_PSI = 14.5038;
    const MPA_TO_BAR = 10;
​
    function pressureToVoltage(pressureMPa) {
      // correct voltage range is 0.5 to 4.5V which maps to pressure 0 to PN
      return 0.5 + (pressureMPa * 4 / PN );
    }
​
    function ecuVoltageToPressureMPa(voltage) {
      // the inverse function, used by ECU to command pressure
      return (voltage - 0.5) * PN / 4;
    }
    /** end car **/
​
    /** tuning module **/
    VERSION = 'db2';
    CURVE_POINTS = 15;
    ADC_RATIO = 5 / 1024;
    OFFSET_MIDPOINT = 128;
​
    function ConfigStringToObject(str) {
      str = str.substr(1); // trim 'X'
​
      if (str.substr(0, 3) === VERSION) {
        str = str.substr(3); // m200M750G94F128C100,100,100,100,100,100,100,100,100,100,100,100,100,100,100
​
        let posMin = str.indexOf('m');
        let posMax = str.indexOf('M');
        let posGain = str.indexOf('G');
        let posOffset = str.indexOf('F');
        let posCurve = str.indexOf('C');
​
        let minRange = parseInt(str.substring(posMin + 1, posMax));
        let maxRange = parseInt(str.substring(posMax + 1, posGain));
        let globalGain = parseInt(str.substring(posGain + 1, posOffset));
        let globalOffset = parseInt(str.substring(posOffset + 1, posCurve));
        let curvePointSplit = str.substring(posCurve + 1).split(',');
        let curvePoints = [];
        for (let i = 0; i < curvePointSplit.length; i++) { // TODO: use some map() stuff
          curvePoints.push(parseInt(curvePointSplit[i]));
        }
​
        //return new TuningModuleSettingsStorage(VERSION, minRange, maxRange, curvePoints, globalGain, globalOffset);
        obj =
        {
          CONFIG_VERSION: VERSION,
          minRange: minRange,
          maxRange: maxRange,
          curvePoints: curvePoints,
          globalGain: globalGain,
          globalOffset: globalOffset
        };
​
        return obj;
​
      }
    }
​
    TuningSettings = {
      CONFIG_VERSION: 'db2\0',
      CURVES: 15,
      minRange: 1024,
      maxRange: 1024,
      curvePoints: [100, 100, 100, 100, 100, 100, 100, 100, 100, 100, 100, 100, 100, 100, 100],
      globalGain: 100,
      globalOffset: OFFSET_MIDPOINT
    };
​
    function TuningVoltageToVoltage(voltage, returnObject = false) {
      let settings = TuningSettings;
      let inputValue = voltage / ADC_RATIO;
​
      /** old method without interpolation **/
      function getCurvePoint(inputValue) {
        let result = -1;
        if ((inputValue > settings.minRange) && (inputValue < settings.maxRange)) {
          let binSize = Math.floor((settings.maxRange - settings.minRange) / CURVE_POINTS);
          if (binSize > 0) {
            let point = Math.floor((inputValue - settings.minRange - binSize / 2) / binSize);
            if (point < CURVE_POINTS) {
              result = point;
            }
          }
        }
        return result;
      }
​
      function getCurveGainInterpolated(inputValue) {
        let result = 100;
        if ((inputValue > settings.minRange) && (inputValue < settings.maxRange)) {
          let binSize = Math.floor((settings.maxRange - settings.minRange) / CURVE_POINTS);
          if (binSize > 0) {
            let pointFloat = (inputValue - settings.minRange - binSize / 2) / binSize;
            let point = Math.floor(pointFloat);
            let nearestPoint = point + 1;
​
            let fraction = Math.abs(point - pointFloat);
​
            let nearestGain = (nearestPoint < CURVE_POINTS && nearestPoint >= 0 ) ? settings.curvePoints[nearestPoint] : 100;
            let currentGain = (point < CURVE_POINTS && point >= 0) ? settings.curvePoints[point] : 100;
​
            let gain = (currentGain * (1 - fraction) + nearestGain * fraction);
​
            if (point < CURVE_POINTS) {
              //result = settings.curvePoints[point];
              result = gain;
            }
          }
        }
        return result;
      }
​
      //let curvePoint = getCurvePoint(inputValue);
      let curveGain = getCurveGainInterpolated(inputValue);
      //lastCurvePoint = curvePoint;
      let result = inputValue; // do all calculations in float
​
      // do not alter value if no valid curve point was found
      //if (curvePoint >= 0) {
      if (curveGain > 0) {
        //result = settings.curvePoints[curvePoint] * inputValue / 100;
        result = curveGain * inputValue / 100;
​
​
        if (result == 0) { // sanity check
          console.warn('sanity check');
          result = inputValue;
        }
      }
      result = result * settings.globalGain / 100 + (settings.globalOffset - OFFSET_MIDPOINT);
      if (!returnObject) {
        //return Math.floor(result);
        return Math.floor(result) * ADC_RATIO;
      } else {
        return {adcIn: inputValue, adcOut: result, vin: voltage, vout: result * ADC_RATIO, curveGain: curveGain};
      }
    }
    /** end tuning module **/
​
​
    function doStuff() {
      print("<PRE>");
​
      //TuningSettings = (ConfigStringToObject('sXdb2m180M900G94F128C80,78,80,82,84,86,88,90,91,92,93,94,95,97,99'.substr(1))); //remove trailing 's';
      TuningSettings = (ConfigStringToObject('sXdb2m180M900G100F128C80,78,78,80,82,83,84,85,87,88,90,91,92,95,98'.substr(1))); //remove trailing 's';
​
      print("Sensor simulation\n");
      let data1 = [];
      for (i = 0; i <= PN; i += 10) {
        print(i + " MPa\t" + (i * MPA_TO_BAR * BAR_TO_PSI).toFixed(2) + " psi\t" + pressureToVoltage(i).toFixed(2) + " V\n");
        data1.push({x: i, y: pressureToVoltage(i)});
      }
      //print("\n");
      chartData([{values: data1, key: 'Sensor output'}], 'Pressure (MPa)', 'Volts');
​
      let data2 = [];
      let data2_2 = [];
      print("Tuning module response simulation\n");
      //print("Vin\tPressure\tVtuned\n");
      for (i = 0.8; i < 4.5; i += 0.025) {
        var Vtuned = TuningVoltageToVoltage(i);
        //print(i.toFixed(2) + ' V\t' + ecuVoltageToPressureMPa(i).toFixed(2) + 'MPa\t' + Vtuned.toFixed(2) + 'V');
        data2.push({x: i, y: Vtuned});
        data2_2.push({x: i, y: ecuVoltageToPressureMPa(i)});
        //print("\n");
      }
      //println();
      //plotData(data2, 300, 300);
      chartData(data = [{values: data2, key: 'Tuning module response'}, {
        values: data2_2,
        key: 'ecuVoltageToPressureMPa',
        yAxis: 2
      }], xLabel = 'V', yLabel = 'V', height = 300, multiAxis = true, yLabel2 = 'MPa');
​
      let data3 = [];
      let data3_1 = [];
      let data3_2 = [];
      let data4 = [];
      let data5 = [];
      let tuningDataOut = [];
      print("ECU control loop with tuning\n");
      print("Desired\tVoltage\tVtuned\tNew pressure\tNew V\tp Increase\t\n");
      for (desiredMPa = 10; desiredMPa < PN; desiredMPa += 0.2) {
        var vin = pressureToVoltage(desiredMPa);
        var vout = TuningVoltageToVoltage(vin);
        var correctionFromECU = vin / vout;
        var correctedPressure = desiredMPa * correctionFromECU;
        var newVoltage = pressureToVoltage(correctedPressure);
        var percentageIncrease = ((correctedPressure - desiredMPa) / desiredMPa * 100);
        //print(desiredMPa + ' MPa\t' + vin.toFixed(2) + ' V\t' + vout.toFixed(2) + ' V\t' + correctedPressure.toFixed(1) + ' MPa\t' + newVoltage.toFixed(2) + ' V\t' + percentageIncrease.toFixed(2) + '%\t');
        if (correctedPressure > PN || newVoltage < 0.5 || newVoltage > 4.5) {
          print(desiredMPa + ' MPa\t' + vin.toFixed(2) + ' V\t' + vout.toFixed(2) + ' V\t' + correctedPressure.toFixed(1) + ' MPa\t' + newVoltage.toFixed(2) + ' V\t' + percentageIncrease.toFixed(2) + '%\t');
          print('OVERLOAD');
          print("\n");
        }
        //print("\n");
​
        data3.push({x: desiredMPa, y: vin});
        data3_1.push({x: desiredMPa, y: vout});
        data3_2.push({x: desiredMPa, y: newVoltage});
        data4.push({x: desiredMPa, y: correctedPressure, y2: correctedPressure});
        data5.push({x: desiredMPa, y: percentageIncrease});
        var tuningOutputObj = TuningVoltageToVoltage(vout, true);
        tuningDataOut.push({x: desiredMPa, y: tuningOutputObj.curveGain});
      }
​
      chartData([
        {values: data4, key: 'correctedPressure'},
        //{values: data3_2, key: 'newVoltage'},
        //{values: data3_1, key: 'vout'},
        {values: data3, key: 'vin'},
        {values: data3_1, key: 'vout'},
        {values: data3_2, key: 'newVoltage'},
      ], 'MPa', 'MPa', 500, true, 'V');
​
      chartData([
        //{values: data3, key: 'vin'},
        {values: data3_2, key: 'newVoltage'},
        {values: tuningDataOut, key: 'curveGain'},
      ], 'MPa', 'V', 500, true, '%');
​
      chartData([{values: data5, key: 'Pressure increase'}], 'MPa', '%', 300);
​
      //print("</PRE>");
    }
  </script>
​
</head>
<body onLoad="doStuff();" id='body'>
<!--pre id='printout'></pre-->
</body>
​
Output 300px

You can jump to the latest bin by adding /latest to your URL

Dismiss x
public
Bin info
anonymouspro
0viewers