Numerical prediction of the earthquake response of classical columns using the distinct element method