Hello I am modeling a vehicle that has a seperately excited DC motor as the power plant utilizing the constant-torque and constant-power regions for traction. The acceleration is governed by the following equation $$ M \frac{dv}{dt} = F_t -(\beta_1+\beta_2v^{2}+\beta_3v) $$
The problem with my model is that the $\beta_1$ component, the rolling resistance, is a constant value and causes the velocity to go negative for the first 0.7 seconds when the model is executed . Any suggestion on how to fix this issue?