A full-discretization method for prediction of milling stability

This paper presents a full-discretization method based on the direct integration scheme for prediction of milling stability. The fundamental mathematical model of the dynamic milling process considering the regenerative effect is expressed as a linear time periodic system with a single discrete time delay, and the response of the system is calculated via the direct integration scheme with the help of discretizing the time period. Then, the Duhamel term of the response is solved using the full-discretization method. In each small time interval, the involved system state, time-periodic and time delay items are simultaneously approximated by means of linear interpolation. After obtaining the discrete map of the state transition on one time interval, a closed form expression for the transition matrix of the system is constructed. The milling stability is then predicted based on Floquet theory. The effectiveness of the algorithm is demonstrated by using the benchmark examples for one and two degrees of freedom milling models. It is shown that the proposed method has high computational efficiency without loss of any numerical precision. The code of the algorithm is also attached in the appendix.