The theoretical model of axial ultrasonic vibration grinding force is built on the basis of a mathematical model of cutting deforming force deduced from the assumptions of thickness of the undeformed debris under Rayleigh distribution and a mathematical model of friction based on the theoretical analysis of relative sliding velocity of abrasive and workpiece. Then, the coefficients of the ultrasonic vibration grinding force model are calculated through analysis of nonlinear regression of the theoretical model by using MATLAB, and the law of influence of grinding depth, workpiece speed, frequency and amplitude of the mill on the grinding force is summarized after applying the model to analyze the ultrasonic grinding force. The result of the above-mentioned law shows that the grinding force decreases as frequency and amplitude increase, while increases as grinding depth and workpiece speed increase; the maximum relative error of prediction and experimental values of the normal grinding force is 11.47% and its average relative error is 5.41%; the maximum relative error of the tangential grinding force is 10.14% and its average relative error is 4.29%. The result of employing regression equation to predict ultrasonic grinding force approximates to the experimental data, therefore the accuracy and reliability of the model is verified.