The objective of this study is to perform finite element analyses using the anisotropic hardening constitutive model on the basis of the total stress concept. An anisotropic hardening constitutive model had been developed in a companion paper, and was then formulated by implicit stress integration and consistent tangent moduli. A nonlinear finite element analysis program was coded including the algorithm, and as a result, the nonlinear solution was accurately calculated and converged to be asymptotically quadratic. In the analysis of a test embankment it was found that the proposed model could predict the displacement of soils more reasonably than the analysis with von Mises type model. In addition the proposed model could predict accurately the actual behavior through the reanalysis of the problem by a reasonable evaluation of the strength parameter.