The use of interface elements is widespread in most numerical methods to separate the structural elements from soil elements. It allows slipping and gapping or overlapping between the two media in contact with highly distinct deformability characteristics, such as piles and retaining walls. However, many existing analytical methods still assume no pile-soil slipping before failure in axially loaded piles. This paper attempts to compare the slipping and non-slipping analytical models by introducing an interface with a virtual thickness. Results show that both have different load-displacement characteristics when comparing to the measured data. The non-slipping model would result in a stiffer load-displacement curve if the same input stiffness parameters applied to both models. Different procedures to calibration the input stiffness parameters in slipping and non-slipping models are recommended.