Deformation of soil skeleton and flow of pore water must be coupled for deformation analyses of saturated soil structures. On the other hand, a rational treatment for contact boundaries connecting each body is needed for analyses of engineering problems since friction and/or sliding are included in almost all technical applications. In this article, in order to strictly consider not only the deformation phenomenon but also the frictional sliding phenomenon, the soil-water coupled finite element program incorporating both the subloading surface and the subloading friction models is developed. Subsequently, simulations of consolidation test by a constant strain rate control are performed. It is revealed by the present program that the frictional sliding phenomenon of the contact boundary influences the deformation behavior as well as flow of pore water phenomenon in the specimen.
