The probabilistic approach was employed to consider the uncertainties in prediction of spudcan penetration in clay. The Latin Hypercube Sampling method was used to generate 300 sample sets of input parameters based on the laboratory and field data. Each set of input parameters was used in conducting Coupled Eulerian Lagrangian finite element analysis. The probabilistic density function of the penetration depth under the maximum applied preload was obtained. The comparison of the results to the measured data showed that the introduced method can produce reasonable prediction of spudcan penetration considering the uncertainties involved in this problem.