Numerical method for the analysis of ground improved by columnar inclusions