We present a consolidation analysis method coupled with an elastoplastic constitutive model for unsaturated soils. The elastoplastic model includes two suction effects related to the increase in effective stresses, yield stress and resistance to plastic deformation due to the increase in suction. The smooth transition from elasticity to plasticity was also achieved by introducing the subloading surface model. The consolidation analysis method employs force equilibrium equations for soil and the Richards' mass conservation equation of pore water as field equations and is formulated by using FEM. Three simulations of typical element test series, a seepage model test and a virtual earth dam during reservoir filling were performed using the consolidation analysis method. The simulation results showed that the consolidation analysis method could be used to analyze geomechanical static problems relating to pore water pressure changes, especially saturation collapse behavior.