Numerical studies on bearing capacity problems of layered clay are performed for smooth and rough strip footings. The finite element method and finite difference method (FLAC) are used for computations of the bearing capacity, entire load-displacement curve and the failure mechanism. The presented results show that it is possible to analyze the bearing capacity of layered clay and to give a progressive failure mechanism clearly. To obtain high quality solutions, it is necessary to review the results on control parameters(e.g., yield function, number of calculation) and compare the results by two numerical methods.