A new reliability calculation method is proposed based on design sensitivity analysis by the finite element method for nonlinear performance constraints in the optimal design of electromagnetic devices. In the proposed method, the reliability of a given design is calculated by using the Monte Carlo simulation (MCS) method after approximating a constraint function to a linear one in the confidence interval with the help of its sensitivity information. The validity and numerical efficiency of the proposed sensitivity-assisted MCS method are investigated by comparing its numerical results with those obtained by using the conventional MCS method and the first-order reliability method for analytic functions and the TEAM Workshop Problem 22.