This paper is concerned with numerical simulation of plastic collapse behavior of square boxes. The finite element limit analysis program was established using shell element. The simulation incorporates with finite element limit analysis based on the upper bound method and the minimization technique. The simulation is performed with various height to side length ratios(h/b) with fixed thickness to side length ratio(t/b). The numerical simulation results are compared with experiments, and is shows good agreement in collapse load and collapse mode.