Abstract: To adress the new β titanium alloy Ti13Nb13Zr as a good medical material for bio-implants but prone to chattering vibration during machining due to its high strength and toughness, this paper comprehensively considers regeneration effect, modal coupling effect, process damping effect and other factors in the milling process, simplifies the milling process into a two-degree-of-freedom milling model, establishes the system dynamics equations, and adopts the second-order full-discrete method for the solving. By comparing the simulation results, the validity of the model was proved, and the influence of various factors on the milling stability was explored. Finally, the orthogonal experiment of three factors and four degrees of freedom was designed, and the milling mode, tool wear band bandwidth, tool immersion ratio and other factors were ranked by stability influence, and it was concluded that the optimal milling scheme of Ti13Nb13Zr under the premise of ensuring stability was using up milling, tool immersion ratio of 0.25, and tool wear band bandwidth of 35 μm.