Modeling Nonlinear Evolution of Baryon Acoustic Oscillations: Convergence Regime of $N$-body Simulations and Analytic Models