This paper presents a systematic optimization process for designing an electromagnetic vibration energy harvester using FEA(finite element analysis) to improve computational accuracy and efficiency. A static FEA is used in the optimization process where trend analysis in a short period of time is rather important than precise computation, while a dynamic FEA is used in the verification step for the final result where precise computation is more important. An electromechanical transduction factor can be calculated efficiently by using an approach to use the radial component of magnetic flux density directly instead of an approach to compute the flux density gradient. The proposed optimization process was verified through a case study where simulation and experiment results were compared.