We prove an estimate for spherical functions $phi_lambda(a)$ on $mathrm{SL}(3,mathbb{R})$, establishing uniform decay in the spectral parameter $lambda$ when the group parameter $a$ is restricted to a compact subset of the abelian subgroup $mathrm{A}$. In the case of $mathrm{SL}(3,mathbb{R})$, it improves a result by J.J. Duistermaat, J.A.C. Kolk and V.S. Varadarajan by removing the limitation that $a$ should remain regular. As in their work, we estimate the oscillatory integral that appears in the integral formula for spherical functions by the method of stationary phase. However, the major difference is that we investigate the stability of the singularities arising from the linearized phase function by classifying their local normal forms when the parameters $lambda$ and $a$ vary.