Significant progress has been made in understanding the behavior and mechanism of sand reliquefaction through physical modelling experiments. However, questions remain unanswered in comprehending the mechanism and factors influencing the reliquefaction potential for a wide range of seismic and ground conditions. In this study, numerical modelling has been performed using a commercially available finite element package assuming plane strain condition. Sand reliquefaction behavior was modelled using a nonlinear UBC3D Plaxis liquefaction model (UBC3D-PLM) which is an effective stress-based elastic–plastic model. This model primarily uses Mohr–Coulomb yield condition in 3D principal stress space. A sand model with dimension 1050 × 600 × 570 mm (length × width × depth) was used for the analysis having a relative density of 50%. Whereas the dimensions are the same as the dimensions of a tank used for physical tests. The model was constructed assuming homogeneous uniform sand deposit. The model domain was discretized using 10-node tetrahedral elements. Size of the element in the finite element mesh is decided based on the largest frequency considered in the analysis. Viscous dampers with Neumann type of boundary condition were used in this modelling as it offers the advantage to nullify the reflected stresses generated during the shaking. To verify the computation technique, a finite element model of the shaking table model was simulated with similar boundary conditions. The results obtained were in good agreement with those obtained by experiments. In this study, three different repeated shaking patterns namely, incremental, uniform, and decremental were employed in the numerical analysis imparting sinusoidal waveform. The applied acceleration amplitude was in the range of (0.1 to 0.4 g) with 2 Hz shaking frequency and 60 s shaking duration. Repeated shaking was applied only after the dissipating the excess pore pressure generated during the previous shaking event. The validation of the numerical approach was achieved by comparing the maximum pore pressure ratio values obtained during the repeated shaking events with the experiments. The results obtained from numerical approach were in good agreement with the available experimental results from 1 g shaking table tests with similar seismic conditions. Among all the shaking patterns and independent shakings, incremental pattern reported lower pore pressure ratio values indicating the beneficial effects of seismic preshaking and the inherent sand resistance offered by the sand specimen against reliquefaction.