The object of this study is the process of end milling, taking into account the discontinuity of the process, simultaneous cutting with several flutes arranged in a spiral, tool runout, and feedback in the elastic machining system, in particular, for the depth of cutting. The subject of the study is the cutting force and identification of its empirical model. During identification, the cutting force coefficient is automatically determined when matching the theoretical and experimental oscillograms of the cutting force component. The reported results related to forecasting the cutting force at end milling are based on a mechanistic approach and involve the process modeling method for forecasting. The simulation uses an algorithm for representing the interaction of the cutter flutes workpiece engagement, based on the scan of the cutter according to the rotation angle coordinate. The algorithm makes it possible to identify empirical coefficients and exponents of the cutting force model based on experimental oscillograms of cutting force components. The built model is implemented in an application program and owing to the representation of the machining system in the form of a closed structural diagram, it allows predicting the elastic displacement, which will determine the actual cutting depth. The developed program under an interactive mode using digital files of experimental cutting force components makes it possible to perform model identification and predict cutting force components with an error of 4.6 %. The adequacy of the algorithms was confirmed by measuring the profile of the machined surface in the places where the cutting mode changed with the feed stopped. The developed simulation algorithm makes it possible to take into account the simultaneous cutting by several flutes arranged in a spiral, the runout of the tool, and the feedback in the elastic machining system, in particular, the depth of cutting