This paper proposes an efficient reliability-based optimization method for designing a superconducting magnetic energy system in presence of uncertainty. To evaluate the probability of failure of constraints, samplingbased reliability analysis method is employed, where Monte Carlo simulation is incorporated into dynamic Kriging models. Its main feature is to drastically reduce the numbers of iterative designs and computer simulations during the optimization process without sacrificing the accuracy of reliability analysis. Through comparison with existing methods, the validity of the proposed method is examined with the TEAM Workshop Problem 22.