This study is aimed at developing a formulation of hyperelastic model for geomaterials and combining it with plasticity. Several existing representative isotropic hyperelastic models with pressure-dependent bulk and shear moduli are reformulated to extend them to finite strains. To be compatible with the framework of multiplicative finite strain elastoplasticity, the stress versus the elastic strain relation, together with the fourth-order elastic tangent moduli tensor, in the description relative to the intermediate configuration is derived for each of the hyperelastic models. Their spatial and material descriptions are also derived. A systematic parametric study with a particular focus on the pressure-dependent property of the elastic moduli is performed to examine and compare the constitutive response of the hyperelastic models under typical cases of simple shear and triaxial compression. Significant differences in the hyperelastic responses depending on the types of model are observed in the analysis. Notably, some models exhibit unexpected unreasonable decrease in stress in the process of triaxial compression. The hyperelastic models are then combined with plasticity. The Cam-clay plasticity with rotational hardening is adopted as a specific prototype model for geomaterials exhibiting induced anisotropy. A formulation of the elastoplastic model based on the multiplicative framework, as well as a stress-point algorithm using fully implicit return-mapping scheme, is developed. Numerical examples are presented to demonstrate the significance of hyperelastic model in the analysis of elastoplasticity. The analysis result reveals that the property of hyperelastic model has significant influences not only in the responses of elastoplastic model but also in the computational aspect to ensure stable analysis.