Skip to content
Go back

FCHEV Longitudinal Dynamics Model: From Mechanics to Electrics

|

Written at somewhere on the Earth

In our previous discussions on Dynamic Programming and Fuel Cell Structure, we explored the control algorithms and the power source itself. However, to execute these algorithms effectively, we require the most critical input variable: PdemandP_{demand} (Power Demand).

The Energy Management System (EMS) does not inherently “know” if the car is climbing a steep hill or facing a strong headwind; it simply asks: “How many kilowatts of electricity does the electric machine (E-machine) need right now?”

Therefore, this article will construct a fundamental Longitudinal Dynamics Model for a vehicle to convert Mechanical Power (such as velocity, acceleration, and road grade) into Electrical Power (PelecP_{elec}).

Step 1: Calculating Resistive Forces

When a vehicle moves, it must overcome four fundamental resistive forces. The sum of these forces is known as the total tractive force at the wheels (FtotalF_{total}).

Ftotal=Froll+Faero+Fgravity+FaccF_{total} = F_{roll} + F_{aero} + F_{gravity} + F_{acc}

1.1. Rolling Resistance Force

This is the friction between the tires and the road surface.

Froll=mgCrcos(α)F_{roll} = m \cdot g \cdot C_r \cdot \cos(\alpha)

1.2. Aerodynamic Drag Force

This is the air resistance acting against the vehicle. It increases with the square of the velocity (the faster the car moves, the greater the drag).

Faero=12ρACdv2F_{aero} = \frac{1}{2} \cdot \rho \cdot A \cdot C_d \cdot v^2

1.3. Road Gravity Force (Gradient Resistance)

This is the component of gravity pulling the vehicle back when climbing a slope (or assisting when descending).

Fgravity=mgsin(α)F_{gravity} = m \cdot g \cdot \sin(\alpha)

1.4. Acceleration Force (Inertial Force)

The force required to change the vehicle’s speed, derived from Newton’s Second Law.

Facc=ma=mdvdtF_{acc} = m \cdot a = m \cdot \frac{dv}{dt}

Nomenclature Table:

SymbolDescriptionStandard Unit
mmVehicle Masskgkg
ggGravitational Acceleration (9.819.81)m/s2m/s^2
CrC_rRolling Resistance Coefficient (typ. 0.010.0150.01 - 0.015)-
α\alphaRoad Grade Angleradrad
ρ\rhoAir Density (1.2251.225)kg/m3kg/m^3
AAFrontal Aream2m^2
CdC_dAerodynamic Drag Coefficient-
vvVehicle Velocitym/sm/s
aaAccelerationm/s2m/s^2

Step 2: From Force to Torque & Mechanical Power

Once we have the total resistive force FtotalF_{total}, we can calculate the required Torque at the wheels (TwheelT_{wheel}) and the Mechanical Power Demand (PdemandP_{demand}).

2.1. Power Demand at Wheels

This represents the instantaneous mechanical power required at the tire-road interface to maintain the vehicle’s current trajectory. It is worth noting that if the vehicle is stationary (v=0v=0), the power demand is zero, regardless of the magnitude of resistive forces (e.g., holding position on a steep gradient).

Pdemand=FtotalvP_{demand} = F_{total} \cdot v

2.2. Wheel Torque

This is the rotational moment required at the wheel axle to overcome the linear resistive force FtotalF_{total}. Based on the principles of mechanics, the wheel acts as a lever arm with a length equal to the tire radius (rtirer_{tire}). This relationship highlights the impact of wheel geometry. For a constant resistive force, a larger tire radius (rtirer_{tire}) necessitates a higher torque output from the powertrain to propel the vehicle.

Twheel=FtotalrtireT_{wheel} = F_{total} \cdot r_{tire}

(rtirer_{tire}: Tire radius in meters)

Step 3: Drivetrain Transmission

The Electric Machine (Motor) is not connected directly to the wheels because the rotational speed and torque of the motor are significantly higher and lower, respectively, than those required at the wheels. Therefore, the motor transmits power through a reduction gear (Transmission/Differential) to multiply torque.

3.1. Torque at the Motor Shaft

After calculating the required torque at the wheels (TwheelT_{wheel}), we must convert it back to the motor shaft. Why?

Because there is always a Transmission/Gearbox between the wheels and the motor. As mentioned, electric motors typically spin very fast but produce relatively low torque. The gearbox amplifies this torque to enable the vehicle to move. Thus, we calculate TmotorT_{motor} to determine how much effort the motor must exert.

Tmotor=TwheelNgηgT_{motor} = \frac{T_{wheel}}{N_g \cdot \eta_g}

Where:

3.2. Motor Angular Velocity

The efficiency of an electric motor (ηmotor\eta_{motor}) varies greatly depending on its rotational speed (RPM). To look up the Efficiency Map in Step 4, we must know exactly how fast the motor is spinning.

First, convert the vehicle’s linear velocity (vv) to the wheel’s angular velocity (ωwheel\omega_{wheel}):

ωwheel=vrtire\omega_{wheel} = \frac{v}{r_{tire}}

Then, convert this to the motor shaft via the gear ratio:

ωmotor=ωwheelNg\omega_{motor} = \omega_{wheel} \cdot N_g

(Standard Unit: rad/srad/s. To convert to RPM, multiply by 60/2π60/2\pi).

3.3. Mechanical Power at the Motor Shaft

This is the actual physical power at the motor’s rotating shaft. It is also the power that the drivetrain must transmit from the motor to the wheels. It is the product of torque (TT) and angular velocity (ω\omega).

Pmech=TmotorωmotorP_{mech} = T_{motor} \cdot \omega_{motor}

Note: If we ignore transmission losses (ηg=1\eta_g = 1), then PmechP_{mech} would be exactly equal to PdemandP_{demand} at the wheels (Law of Conservation of Energy: Power remains constant through gears; only Torque and Speed change inversely).

Step 4: Electric Motor Model

This is the final step for the Energy Management System to determine how much electrical power the vehicle requires. We must convert Mechanical Power (PmechP_{mech}) into Electrical Power (PelecP_{elec}) that the Battery or Fuel Cell must supply. However, this process incurs losses due to motor efficiency (ηmotor\eta_{motor}).

Motor efficiency is not constant; it is a function of Torque and Speed: ηmotor=f(Tmotor,ωmotor)\eta_{motor} = f(T_{motor}, \omega_{motor}).

Calculation Logic for PelecP_{elec}

Case 1: Motoring Mode (Pmech>0P_{mech} > 0) The vehicle consumes electricity to move.

Pelec=Pmechηmotor(T,ω)P_{elec} = \frac{P_{mech}}{\eta_{motor}(T, \omega)}

Case 2: Regenerative Braking Mode (Pmech<0P_{mech} < 0) The vehicle generates electricity to recharge the battery.

Pelec=Pmechηmotor(T,ω)P_{elec} = P_{mech} \cdot \eta_{motor}(T, \omega)

Case Study: Toyota Mirai Calculation 🚗

Let’s calculate the electrical power required for a Toyota Mirai (Gen 2) cruising on a highway.

Assumed Specifications:

Scenario: Constant velocity of 100km/h100 \, km/h (27.78m/s27.78 \, m/s) on a flat road (α=0,a=0\alpha = 0, a = 0).

1. Calculate Resistive Forces

2. Calculate Mechanical Power at Wheels

Pdemand=515.1527.78=14,310W14.31kWP_{demand} = 515.15 \cdot 27.78 = 14,310 \, W \approx 14.31 \, kW

3. Calculate Power at Motor Shaft (Including Gear Losses)

Pmech=Pdemandηg=14.310.9714.75kWP_{mech} = \frac{P_{demand}}{\eta_g} = \frac{14.31}{0.97} \approx 14.75 \, kW

4. Calculate Electrical Power Demand (PloadP_{load})

This is the final figure fed into the EMS.

Pelec=Pmechηm=14.750.9216.03kWP_{elec} = \frac{P_{mech}}{\eta_m} = \frac{14.75}{0.92} \approx 16.03 \, kW

Conclusion: To maintain a speed of 100 km/h, the powertrain system (Fuel Cell + Battery) must continuously supply 16.03 kW. The EMS controller will then decide how to split this 16.03 kW between Hydrogen and the Battery.

I hope this post clarifies the “origin” of the PloadP_{load} or PdemandP_{demand} variable often seen in simulation models.


👨🏻‍💻🏀

Previous Post
Lithium-ion Batteries in Hybrid Vehicles: Modeling and Technical Characteristics
Next Post
Formulating Dynamic Programming for EMS: State-Space and Cost Function