A robust time integration scheme for snapthrough buckling of shallow arches is proposed. The algorithm is a composite method that consists of three sub-steps. Numerical damping is introduced to the system by employing an algorithm similar to the backward differentiation formulas (BDF) method in the last substep. Optimal algorithmic parameters are established based on stability criteria and minimization of numerical damping. The proposed method is accurate, numerically stable, and efficient as demonstrated through several examples involving loss of stability, large deformation, large displacements and large rotations.