Shell Eco-Marathon Vehicle Simulation — V9

GPS track integration, zone-aware driver strategy, lap-based race simulation across the full IMS Infield Road Course.


Overview

V9 moves the simulation from a fixed-time run to a competition-accurate race model. Four laps of the IMS Infield Road Course are simulated under Shell Eco-Marathon Article 226 rules: 4 laps × 3832.5 m = 15,330 m total, 35-minute maximum, 7.30 m/s minimum average speed.

The main additions are a GPS-derived track profile at 5-meter resolution (767 breakpoints), a grade force term in the vehicle dynamics equation, and a zone-aware driver strategy that picks between burn-coast, cruise-hold, and corner-hold based on where the vehicle is on the lap. Output is now built around lap-level metrics — per-lap time, fuel, and km/L reported alongside Article 226 pass/fail validation.


V9 Changelog

NEW Lap-based simulation replacing fixed simulation time — runs until x_Vehicle ≥ 15,330 m (4 full laps). Per-lap breakdowns reported consistent with how Shell scoring is actually applied.
NEW Article 226 pass/fail validation — total race time checked against the 2100 s (35 min) limit. attempt_valid flag and time margin reported in command window and plotted as title annotations on Fig 1.
NEW IMS track profile — GPS + elevation data integrated — track arrays at 5 m resolution derived from two official Shell CSV sources. Curvature computed with a 15 m sliding window; 19 turn complexes identified producing a three-zone classification at 767 breakpoints.
NEW Grade force added to vehicle dynamicsF_grade = m × g × grade included in the force balance at every simulation step. Track elevation range is 3.8 m with max practical grade of ~0.9%.
NEW Zone-aware three-mode driver strategy — Zone 1 (long straights, 34.6%): burn-coast hysteresis. Zone 2 (short straights, 47.0%): cruise-hold at 7.45 m/s. Zone 3 (turn complexes, 18.4%): corner-hold at per-turn speed limit.
UPDATE Speed band shifted up from V8v_target raised to 7.45 m/s to guarantee minimum average speed ≥ 7.30 m/s (Article 226). Estimated lap time 497.7 s × 4 = 1991 s, giving a 109 s margin against the 2100 s limit.
NEW Competition context comparison — Shell Score benchmarked against the 2025 SEM Americas ICE Prototype field in command window output.
PENDING BSFC map — Engine efficiency fixed at 20% pending dyno runs on the assembled vehicle.
PENDING Confirmed mass and gear ratiom = 100 kg and G = 10 pending final drivetrain assembly measurements.

Track Profile — IMS Infield Road Course

Total lap distance3,832.5 m
Turn complexes detected19 (curvature threshold 0.5 deg/m)
Zone 1 — long straight (burn-coast)1,325 m (34.6%) — includes main 1325 m straight
Zone 2 — short straight (cruise-hold)1,800 m (47.0%)
Zone 3 — turn zones (corner-hold)705 m (18.4%)
Elevation range3.8 m total | net change ~0 m (closed loop)
Max practical grade~0.9% (brief) — grade force term included in dynamics
Tightest turnR = 19.3 m at x ≈ 3132–3176 m → v_corner = 7.54 m/s
Zone 1 — Long Straight (burn-coast)
Zone 2 — Short Straight (cruise-hold)
Zone 3 — Turn Complex (corner-hold)
Track grade and elevation profile Fig 5 — Grade % and relative elevation vs position around the lap. Pink shading marks turn zones.
Track zone map overhead view GPS track overhead with zone color coding. White = burn-coast, grey = cruise-hold, pink = turn complexes.

Zone-Aware Driver Strategy

Zone Condition Strategy Speed Reference
1Long straight (>200 m)Burn-Coast (Stateflow hysteresis)v_low=7.5 m/s / v_high=8.0 m/s
2Short straight (<200 m)Cruise-Holdv_target = 7.45 m/s
3Turn / corner zoneCorner-Holdtrack_v_limit (per-turn, 7.54–8.50 m/s)
Zone-aware driver strategy Simulink block Driver strategy block: track zone lookup feeds the mode selector, which routes to Stateflow burn-coast or cruise/corner hold.
Throttle and RPM vs time Fig 4 — Throttle state (burn/coast, top) and engine RPM with idle/max limits (bottom). Lap markers shown as dashed lines.

Race-Level Results

Baseline: v_low=7.2, v_high=7.7, v_target=7.45 m/s, G=10, m=100 kg, Cd=0.30, 4 laps × 3832.5 m.

15,330 m
Total Distance (4 laps)
2016.2 s
Total Race Time (33.6 min)
✓ +83.8 s
Article 226 Valid — Time Margin
7.604 m/s
Avg Speed (min 7.300 m/s)
18.38 g
Fuel Consumed
621.3 km/L
Shell Score
Engine efficiency note: Score of 621.3 km/L computed with fixed 20% thermal efficiency. The BSFC sensitivity sweep shows this is the single highest-leverage unknown in the model — see Parameter Sweeps below. Engine RPM at v_target = 7.45 m/s with G = 10 is ~2800 RPM, within the valid dyno map range (2708–8708 RPM).

Simulation Figures

Fig 1 — Vehicle Speed vs Time (Full Race)

Vehicle speed vs time full race Full race speed trace across all four laps. v_high_on (red), v_low_on (blue), v_target (green), v_min_avg (magenta). Lap markers dashed. Title annotations show avg speed, Shell Score, pass/fail status, and time margin.

Fig 2 — Speed vs Distance, All 4 Laps Overlaid

All 4 laps overlaid speed vs distance All four laps on a single distance axis with zone shading in the background. Consistent lap profiles confirm steady-state is reached by lap 1.

Fig 3 — Fuel & Rolling Shell Score

Fuel consumption and rolling Shell Score vs distance Cumulative fuel (g, top) and rolling Shell Score (km/L, bottom) vs total distance. The stepped fuel profile confirms no energy leakage during coast phases. km/L converges after lap 2.

Fig 4 — Throttle Command & Engine RPM

Throttle and RPM vs time Throttle command (top) and engine RPM with idle/max limits marked (bottom).
Powertrain subsystem block diagram Powertrain subsystem: torque map lookup, gear ratio G, drivetrain efficiency, wheel force output.

Fig 5 — Track Grade & Elevation Profile

Track grade and elevation profile Grade % (left axis) and relative elevation in meters (right axis) vs position. Pink shading marks turn zones.
Grade force vs position around lap F_grade = m × g × grade vs lap position. Shows peak resistive and assistive grade force values around the circuit.

Track Animations

Animation 1 — Single Run Track Playback

Single run track playback, V9 baseline — 4× speed

Animation 2 — Ghost Comparison: Cd=0.12 vs Cd=0.35

TOP-class aerodynamics (Cd=0.12) vs current estimated drag (Cd=0.35) running simultaneously on the IMS lap. Lower drag means longer coast phases, fewer burn cycles, and a real position gap that opens over the lap.

▶ TOP-class   Cd = 0.12
Fuel burned3.656 g
Shell Score777.9 km/L
Lap time~497 s
Gap at finish+2.3 m ahead
▶ Current   Cd = 0.35
Fuel burned5.247 g
Shell Score541.6 km/L
Lap time~505 s
Gap at finish
Reducing Cd from 0.35 → 0.12 saves 1.591 g of fuel per lap and gains 236.3 km/L (+43.6%). Over 4 competition laps: 6.4 g saved.
Ghost comparison: TOP-class Cd=0.12 (yellow) vs current Cd=0.35 (cyan) — 8× speed

Competition Context — 2025 SEM Americas Field

RankTeamScore (km/L)Notes
BYU (2024 reference)995Fully covered wheels, professional build
1Mater Dei HS7182025 Americas champion
2IFRS Erechim629
3Michigan Tech323
4Cedarville University315
5Penn St Behrend177
6Schurr HS92
This vehicle (V9 simulation)621.3 km/Lm=100 kg, Cd=0.30, G=10, effic=0.20 (fixed)

Parameter Sweep Results — SweepV9_track

Five independent sweeps on the full IMS lap-distance simulation. All scored in km/L. Baseline: 621.3 km/L.

Sweep 1 — Burn-Coast Speed Band

Optimum: v_low = 7.0 m/s, v_high = 7.5 m/s628.0 km/L | t_total = 2039.5 s (margin +61 s)
Speed band sweep heat mapSweep 1 — 2D heat map: v_low (y-axis) vs v_high (x-axis), color mapped to km/L.
Gear ratio vs Shell ScoreSweep 2 — Gear ratio has minimal efficiency impact across G=7–14. Set it to keep RPM in range and move on.

Sweep 3 — Vehicle Mass

Best: m = 80 kg → 697.9 km/L | Worst: m = 180 kg → 432.7 km/L | Sensitivity: 61.2%
Mass vs Shell ScoreSweep 3 — Monotonically decreasing. Each 10 kg costs roughly 25–30 km/L.
Cd vs Shell ScoreSweep 4 — Near-linear improvement as drag drops. Cd=0.15 reaches Mater Dei territory.

Sweep 4 — Drag Coefficient

Best: Cd = 0.05 → 1003.3 km/L | Worst: Cd = 0.55 → 449.9 km/L | Sensitivity: 123.0% — highest of all five sweeps

Sweep 5 — Engine Efficiency (BSFC Sensitivity)

Engine efficiency swept from 8% to 34%. The relationship is perfectly linear — Shell score scales directly with efficiency. Every 2 percentage points is worth approximately 62 km/L.

Key implication: At Cd=0.30 and m=100 kg, matching Mater Dei's 718 km/L requires roughly 23% engine efficiency — achievable with a well-tuned carburetor and ignition timing. Matching BYU's 995 km/L requires ~32% efficiency, Cd near 0.10, and mass well below 100 kg.
Engine efficiency vs Shell ScoreSweep 5 — Perfectly linear. Score scales proportionally with efficiency. Baseline at 20% = 621 km/L. Each 2% gain is worth ~62 km/L.
η_engkm/LContext
8%248.5
16%497.1
20%621.3← baseline
22%683.5Mater Dei range
24%745.6Mater Dei+
32%994.1BYU 2024 range

Design Sensitivity Ranking (V9)

Drag Coeff (Cd) 123.0%
Engine Effic. 124.7%
Mass 61.2%
Speed Band 1.1%
Gear Ratio 0.8%

Key takeaways:

Command Window Output

=== Shell Eco Run Results (V9) ===
Total distance:    15330.0 m
Total time:        2016.2 s  (33.6 min)
Time margin:       +83.8 s
Attempt valid:     YES ✓
Avg speed:         7.604 m/s
Shell Score:       621.3 km/L

--- 2025 Field Comparison ---
  vs Mater Dei (1)   718.0 km/L  → 86.5%
  vs IFRS Erechim    629.0 km/L  → 98.8%
  vs Michigan Tech   323.0 km/L  → 192.4%

--- Per-Lap Breakdown ---
Lap  Time(s)  LapTime  AvgSpd  Fuel(g)  km/L
1    505.8    505.8    7.578   4.911    581.4
2    1009.3   503.5    7.611   4.506    633.6
3    1512.8   503.5    7.612   4.486    636.5
4    2016.2   503.4    7.613   4.478    637.5

F_Engine check:    CLEAN — no leakage during coast
Coast shape:       CURVED (exponential) — physics correct

Lap 1 fuel (4.91 g) is higher than laps 2–4 (~4.48–4.51 g) because the vehicle accelerates from rest at t=0. Laps 2–4 are within 0.6% of each other — steady-state cycling is established by lap 2. The max uphill grade of 1.66% produces a grade force (16.28 N) that is 1.8× flat rolling resistance, making it the dominant terrain feature on the IMS infield circuit and the main driver of the lap 1 vs lap 2+ fuel difference.


Model Architecture (V9)

Full system Simulink overview V9 Full system overview: DriverStrat (zone-aware), Powertrain, VehicleDynamics with F_grade, EnergyModeling, and lap counter.
Powertrain subsystem V9 Powertrain subsystem: torque map lookup, gear ratio, drivetrain efficiency, wheel force output.
Zone-aware driver strategy V9 Zone-aware driver strategy: track zone lookup feeds the mode selector, which routes to Stateflow burn-coast or cruise/corner hold.

MATLAB & Simulink Scripts

View Source Code on GitHub


Future Work (V10)