Estimating the deformation of soil around the pile contributes to reliable design of structures under pullout force. This work presents the results of a series of small-scale physical modelling tests designed to investigate the uplift resistance of piles with diameter of 5 cm and slenderness ratios of 1, 2, 3 and 4 in loose sand. Close photogrammetric technique and particle image velocimetry(PIV) were employed to observe the failure patterns due to uplift force on piles. The results show that the shear zones curve slightly outward near the ground surface. After peak resistance, the shear strain concentrates into a pair of narrow shear bands,then a flow around mechanism is formed accompanied by a reduction in the uplift resistance. The results from the laboratory tests were verified by analytical method proposed by Chattopadhyay and PLAXIS 2D and 3D finite element method software. It is found that the depth and width of the failure surface increase with the increment of the slenderness ratio. A good agreement is observed among the measured bearing capacity and obtained failure surface of the models and the results of numerical modelling. Finally, the maximum deformation of loose and dense sand respectively with densities of 25% and 75% were compared in the stage of fully removing pile. The results shows that the deformation of the soil is related to its density, therefore it depends on its dilatancy. Estimating the deformation of soil around the pile contributes to reliable design of structures under pullout force. This work presents the results of a series of small-scale physical modeling tests designed to investigate the uplift resistance of piles with diameter of 5 cm and slenderness ratios of 1, 2, 3 and 4 in loose sand. Close photogrammetric technique and particle image velocimetry (PIV) were employed to observe the failure patterns due to uplift force on piles. The results show that the shear zones curve slightly outward near the ground surface. After peak resistance, the shear strain concentrates into a pair of narrow shear bands, then a flow around mechanism is formed by a reduction in the uplift resistance. The results from the laboratory tests were verified by analytical method proposed by Chattopadhyay and PLAXIS 2D and 3D finite element method software. It is found that the depth and width of the failure surface increase with the increment of the slenderness A good agreement is observed among the measured bearing capacity and obtained failure surface of the models and the results of numerical modeling. Finally, the maximum deformation of loose and dense sand respectively with densities of 25% and 75% were compared in the stage ofly removing pile. The results shows that the deformation of the soil is related to its density, therefore it depends on its dilatancy.