A numerical modelling methodology for the in-situ axial pressure-displacement response of soft soils is developed based on the Discrete Element Method. The resulting models achieve good agreement with the idealized Beirnstein form of the ‘p-z’ equation describing soft soil uniaxial pressure-displacement response. A method of determination of appropriate input parameters for the models from bevameter measurements of soft soils is proposed.