Smoothing the Brush Model
GM3's longitudinal force, lateral force, and aligning moment are discontinuous at the adhesion–sliding boundary, and the sign function kills gradients entirely. We replace the hard switch with a sigmoid gate, the sign function with $\tanh(k\alpha)$, and absolute values with a softened norm.
\[ G(\sigma, \psi) = \frac{1}{1 + e^{-k\left(\sigma - \frac{1}{\psi}\right)}} \]The gate steepness $k$ acts as a temperature. The resulting worst-case force error decays with $(1/k)^3$, so the smoothed model is the exact model plus a small bounded disturbance — and receding-horizon MPC re-plans from the measured state each step, so the mismatch cannot accumulate.