Based on the slip-line field theory, a two-dimensional slip failure mechanism with mesh-like rigid block system was constructed to analyze the ultimate bearing capacity problems of rough foundation within the framework of the upper bound limit analysis theorem. In the velocity discontinuities in transition area, the velocity changes in radial and tangent directions are allowed. The objective functions of the stability problems of geotechnical structures are obtained by equating the work rate of external force to internal dissipation along the velocity discontinuities, and then the objective functions are transformed as an upper-bound mathematic optimization model. The upper bound solutions for the objective functions are obtained by use of the nonlinear sequential quadratic programming and interior point method. From the numerical results and comparative analysis, it can be seen that the method presented in this work gives better calculation results than existing upper bound methods and can be used to establish the more accurate plastic collapse load for the ultimate bearing capacity of rough foundation.