The Screw Driving Sounding (SDS) is a recently developed in-situ testing method for soil/site characterisation. To better understand the mechanism of testing, it is essential to create a simulation model for the SDS test to assist in better understanding the response and in improving the testing process, not to mention the development of empirical correlations to estimate geotechnical parameters for use in design. Complex problems involving large deformations are usually difficult to solve with the classical finite element (FE) method because large deformations can lead to large mesh distortions and contact problems. In this paper, a computational model using the finite element method is developed to simulate the drilling process involved in SDS tests; the Coupled Eulerian-Lagrangian (CEL) approach is used to deal with large deformation problems. The SDS drilling process in a sandy deposit with defined stress states (i.e., relative density and effective confining pressure) is simulated, and the SDS-derived parameters are monitored. Based on the simulation results, a chart was established to correlate the measured SDS parameter with the internal friction angle for different vertical effective overburden stresses. The derived chart is validated with the results of laboratory tests performed on samples taken from several sites in Christchurch, NZ, adjacent to the locations of the SDS tests. The results show that the CEL FE framework can model complex physical processes encountered during the SDS drilling. Moreover, the developed chart can be used to estimate the friction angle of the sandy soil based on the SDS-measured torque at a given depth.