On the dynamics of ball end milling: modeling of cutting forces and stability analysis

This paper presents a dynamic force model and a stability analysis for ball end milling. The concept of the equivalent orthogonal cutting conditions, applied to modeling of the mechanics of ball end milling, is extended to include the dynamics of cutting forces. The tool is divided into very thin slices and the cutting force applied to each slice is calculated and summed for all the teeth engaged. To calculate the instantaneous chip thickness of each tooth slice, the method of regenerative chip load calculation which accounts for the effects of both the surface undulations and the instantaneous deflection is used. To include the effect of the interference of the flank face of the tool with the finished surface of the work, the plowing force is also considered in the developed model. Experimental cutting forces are obtained using a five-axis milling machine with a rotary dynamometer. The developed dynamic model is capable of generating force and torque patterns with very good agreement with the experimental data. Stability of the ball end milling in the semi-finishing operation of die cavities is also studied in this paper. The tangential and radial forces predicted by the method of equivalent orthogonal condition are fitted by the equations Ft=Kt(Z)bhav and Fr=Kr(Z)Ft, where b is the depth of cut and hav is the average chip thickness along the cutting edge and Z is the tool axis coordinate. The polynomial functions Kt(Z) and Kr(Z) are the cutting force constants. The interdependency of the axial and radial depths of cut in ball end milling results in an iterative solution of the characteristic equation for the critical width of cut and spindle speed. In addition, due to different cutting characteristics of the cutting edge at different heights of the ball nose, stability lobes are represented by surfaces. Comparison of the time domain simulation for the shoulder removal process in die cavity machining with the analytical predictions shows that the proposed method is capable of accurate prediction of the stability lobes.