The harmonic balance finite element method proposed by us solves the nonlinear time-periodic solution expressed as the sum of fundamental and harmonic components. This paper describes a new procedure for calculation which emphasizes the field computation in the harmonic domain. As a result, the improvement of reducing memory and calculation time is achieved.