The liquefaction of gravelly soil has been demonstrated in many earthquakes, but the constitutive model to predict its liquefaction behavior is limited. The dilatancy equation plays an important role in establishing elastoplastic models of geotechnical materials. Recognizing the unique dilatancy characteristics of gravelly soil, this study introduces an improved nonlinear dilatancy equation into the generalized plasticity model. Consequently, a unified elastoplastic constitutive model is proposed, aimed at accurately predicting both the monotonic and cyclic behaviors of saturated gravelly soil under diverse drainage conditions. A series of drained and undrained triaxial tests of gravelly soil under monotonic and cyclic loadings are carried out to validate the model performance. In addition, the model performance is also examined through monotonic and liquefaction tests of gravelly soil obtained from prior publications. Furthermore, the proposed model is implemented into the nonlinear finite element program GEODYNA and applied to the numerical simulation of the Kobe Port Island liquefaction to better validate the capabilities of the proposed model. From the comprehensive analysis of triaxial tests and numerical simulations, the proposed model performs well in predicting mechanical behavior of gravelly soil and can provide a powerful tool for seismic liquefaction analysis of gravelly soil sites.