An elastoplastic, dynamic, finite-difference method was applied to study the effects of nonlinear seismic soil–pile interaction on the liquefaction potential of marine sand with piles. The developed model was well validated using the centrifuge test. The results showed that acceleration, bending moment, and excess pore water pressure complied well with centrifuge test results. The effect of different affecting parameters on liquefaction potential was investigated using parametric study. Using a sensitivity analysis, the pile embedment parameter was shown to be the most influential parameter. Finally, applying the evolutionary polynomial regression technique, a new model for predicting the liquefaction potential was presented.