The aim of this paper is to obtain the stability lobes for milling operations with a nonlinear cutting force model. The work is focused on the generation of stability lobes based on a formulation with Chebyshev polynomials and the semidiscretization method, considering a nonlinear cutting force model. Comparisons were conducted between experimental data at 5% radial immersion with aluminum workpiece and predictions based on Chebyshev and semidiscretization. In all cases, the use of nonlinear cutting force model provides better prediction of process stability conditions.