Mises stress method and interfacial energy method were used to predict the rafting type of nickel based single crystal superalloy in different orientations at 950℃. The rafting prediction results of above two methods were consistent. The rafting type is N-type in[001] orientation

P-type in[011] orientation and there is no rafting in[111] orientation. Based on the rafting type prediction results and crystal slip theory

a creep constitutive model of nickel based single crystal superalloy considering Kachanov-Robotnov (K-R) equation of damage evolution was introduced. This creep constitutive model was incorporated as a UMAT user defined subroutine in Abaqus and applied to description of the creep behaviors of CMSX-4 single crystal superalloy in[001]

[011] and[111] orientations at 950℃

180~450MPa. This model could predict the rafting and slip type of nickel based single crystal superalloy and couple well with the creep deformation physical mechanism. Therefore

this model could simulate the second and third stages of creep curve of nickel based single crystal superalloy

as this has been verified by the experimental data.