This paper presents the pullout behavior of suction pile using finite difference method; and the commercial software, FLAC3D, was employed for the numerical analyses. The ultimate pullout capacity of suction pile was predicted using conventional equations, and the results were compared with the results from numerical analyses with varying pile diameter, pile length, and the undrained shear strength of clays. Based on the results from 24 analyses, it was found that the failure pattern depends not only on the drainage condition of suction pile, but also on the pile dimensions and the material properties of surrounding soils. The developed side shear (DSS) along the internal surface of the suction pile was collected from numerical analyses, which was used to classify the failure type between sliding failure and tensile failure. Regardless of the external DSS, the high internal DSS tends to result in sliding failure in the numerical analyses, which conforms well to the estimation from conventional equations.