To obtain accurate optimal design results in electric machines, the finite element analysis (FEA) technique should be used; however, it is time-consuming. In addition, when the design of experiments (DOE) is conducted in the optimal design process, mechanical design, analysis, and post process must be performed for each design point, which requires a significant amount of design cost and time. This study proposes an automated DOE procedure through linkage between an FEA program and optimal design program to perform DOE easily and accurately. Parametric modeling was developed for the FEA model for automation, the files required for automation were generated using the macro function, and the interface between the FEA and optimal design program was established. Shape optimization was performed on permanent magnet synchronous motors (PMSMs) for small electric vehicles to maximize torque while maintaining efficiency, torque ripple, and total harmonic distortion of the back EMF using the built-in automation program. Fifty FEAs were performed for the experimental points selected by optimal Latin hypercube design and their results were analyzed by screening. Eleven metamodels were created for each output variable using the DOE results and root mean squared error tests were conducted to evaluate the predictive performance of the metamodels. The optimization design based on metamodels was conducted using the hybrid metaheuristic algorithm to determine the global optimum. The optimum design results showed that the average torque was improved by 2.5% in comparison to the initial model, while satisfying all constraints. Finally, the optimal design results were verified by FEA. Consequently, it was found that the proposed optimal design method can be useful for improving the performance of PMSM as well as reducing design cost and time.