The bottom interaction stress distribution of circular foundation on rock subgrade is generally nonlinear. Different from in soil, it is more complex and have no common conclusion. But it is important for design of this kind foundation. According to the circular foundation’s axial symmetry, using the homogeneous Lame equation in the circular cylindrical coordinates, and the Boussinesq-Galerkin’s general solution, selecting the Love’s strain function as a trial function, the analytic expressions of the bottom stress distribution of the circular foundation on rock subgrade is obtained by the trial Love’s strain function satisfying the boundary conditions. The results of field model tests and numerical simulation analysis are discussed here. It is provided about the bottom interaction stress distribution of circular foundation that is based on rock subgrade by compared the theoretical analysis with the test and numerical calculation. This study conclusion is valuable for design and practice of this foundation engineering.