Conical intersections (CIs) of ethylene have been successfully determined using spin-flip density functional theory (SFDFT) combined with a penalty-constrained optimization method. We present in detail three structures, twisted-pyramidalized, hydrogen-migrated, and ethylidene CIs. In contrast to the linear response time-dependent density functional theory, which predicts a purely twisted geometry without pyramidalization as the S1 global minimum, SFDFT gives a pyramidalized structure. Therefore, this is the first correct optimization of CI points of twisted ethylene by the DFT method. The calculated energies and geometries are in good agreement with those obtained by the multireference configuration interaction (MR-CI) method and the multistate formulation of second-order multireference perturbation theory (MS-CASPT2).