A systematic approach of the energy method is proposed for analysis of axisymmetric deep drawing in which the total deforming region is divided into five sections by the geometric characteristic. The corresponding solution is found through optimization of the total energy dissipation with respect to some parameters assumed in the kinematically admissible velocity field defined over each region. The sheet blank is divided into three-or five-layers to consider the bending effect. For the evaluation of frictional energy, it is assumed that the blank holding force acts on the outer rim of the flange and that the contact pressure acting on punch shoulder or die shoulder has uniform distributions, respectively. The computed results by the present method are compared with the experiment and the computed results by the elastic-plastic finite element method for the distribution of thickness strain and the relation between the punch stroke and punch load. The results for the case of multi-layers show better agreements than for the case of a single layer in load vs. stroke relation and strain distribution. It is thus shown that the multi-layer technique can be effectively employed in analyzing axisymmetric deep drawing in connection with the energy method.