Considering the initial shear stiffness depth effect of the pile-soil interface and keeping the uniformity along the contact surface thickness direction of the shear stress and shear strain,the paper put forward the zero-thickness zone contact surface constitutive model based on the statistical damage constitutive model.According to the shear experiment data ,we obtained the simulation curves fitting parameters of the improved constitutive mod-el under different normal stresses .Calculation results indicated that the zero-thickness zone contact surface con-stitutive model developed was correct...