In order to simulate the grinding surface more accurately, a novel modeling method is proposed based on the ubiquitiform theory. Combined with the power spectral density (PSD) analysis of the measured surface, the anisotropic characteristics of the grinding surface are verified. Based on the isotropic fractal Weierstrass–Mandbrot (W-M) function, the expression of the anisotropic fractal surface is derived. Then, the lower bound of scale invariance δmin is introduced into the anisotropic fractal, and an anisotropic W-M function with ubiquitiformal properties is constructed. After that, the influence law of the δmin on the roughness parameters is discussed, and the δmin for modeling the grinding surface is determined to be 10−8 m. When δmin = 10−8 m, the maximum relative errors of Sa, Sq, Ssk, and Sku of the four surfaces are 5.98%, 6.06%, 5.77%, and 4.53%, respectively. In addition, the relative errors of roughness parameters under the fractal method and the ubiquitiformal method are compared. The comparison results show that the relative errors of Sa, Sq, Ssk, and Sku under the ubiquitiformal modeling method are 5.36%, 6.06%, 5.84%, and 4.53%, while the maximum relative errors under the fractal modeling method are 23.21%, 7.03%, 83.10%, and 7.25%. The comparison results verified the accuracy of the modeling method in this paper.