Strojniški vestnik - Journal of Mechanical Engineering 60(2014)11, 716-724 © 2014 Journal of Mechanical Engineering. All rights reserved. D0l:10.5545/sv-jme.2013.1588 Original Scientific Paper Received for review: 2013-12-04 Received revised form: 2014-05-09 Accepted for publication: 2014-06-16 Mechanics and Dynamics of Helical Milling Operations Changyi Liu1* - Gui Wang2 - Matthew S. Dargusch2 1 Nanjing University of Aeronautics & Astronautics, College of Mechanical and Electrical Engineering, China 2 The University of Queensland, School of Mechanical and Mining Engineering, Australia This paper presents analytical models dealing with time domain cutting forces and the time and frequency domain dynamics of helical milling operations. The cutting forces both on the side-cutting edges and on the end-cutting edges along the helical feed path are modelled by considering the tangential and the axial motions of the tool. The dual periodicity, which is caused by the spindle rotation, as well as the period of the helical feed of the cutting tool has been included. The models of the machining process dynamics and of the chatter stability problems are decomposed into two parts: the critical axial depth of cut, and the critical radial depth of cut, both of which are solved. The dynamics model allows the selection of processing parameters, including axial cutting depth, radial cutting depth, and spindle rotation velocity in the chatter-free zone and the prediction of chatter vibrations during helical milling. Experiments have also been performed for real machining observations in order to verify the chatter-free zone obtained from the analytical model. Keywords: cutting, process modelling and planning, helical milling, cutting dynamics, chatter 0 INTRODUCTION Helical milling has been applied to generate or enlarge boreholes by means of a milling tool being fed on a helical path into the work piece. It is an advanced hole-making technology, used mainly for typically difficult-to-cut materials. Helical milling is shown to be able to machine H7 quality holes with a surface finish of 0.3 ^m Ra in hardened AISI D2 tool steel, and to be able to enhance the tool life [1]. Helical milling was utilized in the hole-making in carbon fibre-reinforced plastic, and the cutting force coefficients were identified and the mechanistic modelling technique was used to predict cutting forces [2]. The method was employed on CFRP-titanium layer compounds, to model the non-deformed chip geometries for the process, and to explain the effect of the helical feed on the cutting forces [3]. Accuracy and precision are important considerations for machining products along with efficiency and cost. The investigation of machining stability and chatter vibration is significant for process planning. Milling can be one of the most complicated cutting process. A cutting forces prediction model considering the tool's geometrical features and the tool-work piece interaction is the foundation for the chatter modelling. Most research analyses the cutting forces on discrete milling tools elements, and then integrates them to result in the cutting forces on cutting tools. The cutter is divided into a number of discrete disk elements along its axis, and each elements' surface is considered as the plane of cut [4] and [5]. The factors that influence milling forces during the multi-axis milling process, such as the varying lead and tilt angles, the cutting tool deflections and form errors, have been included [6] and [7]. During the helical milling processing, both side-cutting edges and end-cutting edges involved simultaneously is an important feature. The cutting forces of the end-milling process using a flat end milling tool considering both the side edge and the side edge cuttings simultaneously were modelled [8]. For helical milling operations, as well as plunge milling operations, axial feed is another feature that should not be neglected. By considering the rigid body motion of the cutter and three translational and torsional vibrations of the structure, models to predict cutting forces, vibrations, and the chatter stability in the frequency domain for plunge milling process have been presented [9] and [10]. The finite element method (FEM) and finite difference method have been applied as a novel approach to simulate the dynamics of the machining system. A coupled finite element (FE) model was created to simulate the dynamic properties at the tool end point of the whole mechanical system, including the tool, spindle, and machine tool frame [11]. FEM has also been used to model the dynamics and vibration. This approach simulates the process, simultaneously including both vibration and the chip formation [12]. With regard to machining dynamic effects, machining geometric defects have been modelled, and the dynamic displacements due to clamping and machining forces have been defined using FEM [13]. The resonance frequencies of the axial vibration mode of the milling tool were predicted using FEM model of the vibration milling tool and demonstrated that high-frequency vibrations superimposed onto the tool lead to stabilization of the 716 *Corr. Author's Address: College of Mech. and Electrical Engineering, Nanjing University of Aeronautics & Astronautics, 210016, Nanjing, China, liuchangyi@nuaa.edu.cn Strojniski vestnik - Journal of Mechanical Engineering 60(2014)11, 716-724 milling process with superior surface finishes [14]. These FEM models reflected the physical processes in the machining system accurately and visually, but such approaches are time-consuming and require intensive computing resources. Machining chatter was detected using on-process methods, such as acquiring and processing cutting forces, vibration signals and acoustics signals, or by using off-process methods [15]. A new method for the detection of chatter in the end-milling operation based on the wavelet transform has been suggested, which provides various ways to determine chatter characteristics real-time or post process [16]. Sound pressure, machining force and tool displacements are measured during the process to evaluate the stability of high-speed milling, and tool displacements are used as input data to simulate the chip thickness variation during the process [17]. Helical milling has been presented as an alternative enabling technology for drilling operations [1] to [3]. Recently, a few analyses of the cutting forces of the helical milling process have emerged [18] and [19]. Although the mechanisms and dynamics of special operations with axial feed characteristics, such as plunge milling and drilling, or with the tangential feed properties, such as multi-axis milling and circular milling, were investigated, the dynamics and the chatter stability model for helical milling or milling with axial feed have not been found. The characteristics of the helical milling operation are discussed in this paper, including modelling the cutting force influenced by the axial feed; to build the dynamics model in the four degrees of freedom (DOF) cutting tool - machine tool vibration system, and to simulate the chatter limitation for both the axial cutting depth and the radial cutting depth. This research based on the analytical cutting forces model, combining the interactions both on the side-cutting edges and on the end-cutting edges of helical milling operations, focusing on the machining process dynamics model and chatter stability problems, which will dealt with into two parts: the critical axial depth of cut, and the critical radial depth of cut. Both the chatter stability limitations of the critical axial depth of cut and the critical radial depth of cut during helical milling process will be solved. 1 CUTTING FORCE MODEL FOR HELICAL MILLING The motion curve of a certain point on the cutting edge of the milling tool performing helical milling is composed of the helical motion of the tool axis and the circular motion of the point relative to the axis. The diameter of the bore and the endmill diameter Dm, the rotating angular velocity Qh of the helix feed, the axial feed speed fva and the tangential feed speed fvt, the tangential feed rate per tooth fzt, the axial feed rate per tooth fza, the pitch of the helix curve of the reference frame P, and the spindle rotational velocity Q, are depicted in Fig. 1. To define the movement of the cutting tool, and the cutting force, an X Y, Z global coordinate system (GCS) is attached the work piece, and an x, y, z local coordinate system (LCS) is attached to the cutter. While modelling the cutting forces of the helical milling process, there have been two considerations that might also affect the dynamics of the helical milling process. One is the cutting force Dm Fig. 1. Schematics of helical milling Mechanics and Dynamics of Helical Milling Operations 717 Strojniski vestnik - Journal of Mechanical Engineering 60(2014)11, 716-724 fluctuating periodically as a result of the cutting tool circumferential feed; the other is the cutting force component due to the cutting tool axial feed. The axial feed force mostly occurs at the end-cutting edge of the milling tools. The detailed interpretation and the development of the relative equations in this section refer to the authors' previously published literature [18]. The cutting forces loaded on the tool are a combination of side-edge cutting forces and end-edge cutting forces; these two parts are modelled and then summarized: F = F' + F (1) where, F' and F* are cutting force components on the side-cutting edge and on the end-cutting edge, respectively. The cutting force components on the jth end-cutting edge to the LCS coordinate can be represented Fx A - cos^j 0 Fyj 0 sin^. F / 0 0 t; 0 0 0 0 0 0 1 0 0 1 [®][ * * ]f f ,(2) J ■ = FZj TZj \f - + F ' ] F - + F ' yj • F . + F . T* y. (3) cos Qht sin Qht 0 0 - sin Qht cos Qht 0 0 0 0 10 0 0 0 1 Finally, the superposition of all the cutting forces on the Nm cutting teeth is the cutting forces on the cutting tool. The established time domain analytical cutting forces model for helical milling operations reflects both the cutting process on the side-cutting edge and on the end-cutting edge, incorporating the influence of the spindle motion and the tool helical feed. 2 DYNAMICS OF HELICAL MILLING Helical milling is a typical interrupted cutting operation with a tangential feed and an axial feed. The machine tool-spindle-milling tool system was considered to be flexible in four degrees of measurement (DOM), i.e. in the global X, Y, Z and C directions, as shown in Fig. 2. This flexibility could lead to chatter instability if the axial depth of cut, the radial depth of cut and the spindle rotation speed were not selected properly. where K*vc, K*ve, K*nc, K are the cutting force constants on the end-cutting edges along the circumferential and normal directions respectively, [ * * ]= K K've K*c K*e K*c K*e dr cosQ - sin0 0 sin0 cos0 0 0 0 0 r cosd -r sin0 0 A = Nf 2n B = NJz, CQSfflj 2n d = argtan - A B is the relative rotational angle of the cutting tooth j, r is the distance from centre point to a certain point on the end-cutting edge. The side-cutting edge force components are added to the end-cutting force components on the jth tooth, then converted to GCS coordinates. Fig. 2. Four DOM dynamics system of helical milling The proposed cutting force model was used to build the dynamics model of the helical milling system and to predict the limitation of the chatter stability. Based on the regenerative chatter mechanism, the time domain dynamics relationship between the cutting force and dynamic displacement can be expressed as: D D 2 e 718 Liu, C. - Wang, G. - Dargusch, M.S. Strojniski vestnik - Journal of Mechanical Engineering 60(2014)11, 716-724 [M^^ncH^H^On^)} (4) where R is the dynamic displacement or amplitude vector with four components X, Y, Z and C. From the proposed cutting forces model, the cutting forces have a functional relationship with ap, ae, fzt, fat, Q, Qh etc. The objective of chatter prediction during the helical milling operations should be to find the proper relationship between ap, ae, and Q in order to avoid chatter vibration. To find the chatter lobe including the ap, ae, Q directly is difficult. Therefore, a certain level of idealization and assumption has to be introduced. According to the cutting forces model that resulted from cutting forces on the side-cutting edges and cutting forces on the end-cutting edges, the dynamics model and chatter limitation problem was decomposed into two parts: the X and Y direction, and Z and C direction. 2.1 The Chatter Limitation of Axial Depth of Cut The X and Y components of the cutting forces on the cutting tools mainly result from the interaction of the side-cutting edges, because the magnitude of the interaction on the end-cutting edges is smaller, and their sum is almost zero. Hence, the axial depth of cut is affected by the X and Y components of the cutting forces. Regenerative chatter is the result of variation in the dynamic chip thickness as depicted in Fig. 3. The variation of dynamic chip thickness can be described as: Ah (vj ) = K - h (c )} = [G {G>c )] {F (a>c)}, (9) where [g (o)] is the transfer function, or frequency response function. From Eq. (5), a frequency- Y \F' . , : > — a g. E* • ¿—I T P°J I Fy I j=1 2 -K te sin 29j - K rc (1 - cos 29j) -Kte (1 + cos 29) - sin 29 f^ j _ 1 ; ^fAxl Ktc (l - cos 29)- Krcsm2<9 Ktc sin 2^. - Krc (l + cos 2^ )WJ 2 tc P L ^[AyJ' 29j - (l + cos 2^j) , f1. V, e[®„ . ] , ^ , ^ where "jo. V,eK.®.]" K" = KK- Mechanics and Dynamics of Helical Milling Operations 719 Strojniski vestnik - Journal of Mechanical Engineering 60(2014)11, 716-724 dependent dynamic chip thickness vector as in Fig. 3 shown, fAX (©„ AY iA* (»< )}= {ZZfri* »"i (°< »= = (1 -e^){R(«, )}=( -)»G(*, )] (f(»„)}.(10) Substitute Eq. (10) to Eq. (9), we obtain: 1 2 K'tcap (1 - e» )• {(» )} 'I AXY (»c )][G (» )] {F (»c )}. (11) The characteristic function of Eq. (11) is: [I]-2K'tcap (1 -e-1^)) (a>e)][G(g>c)] = 0. (12) The response function has a nontrivial solution only if its determinant is zero, therefore, det|[I] +A[Axy (coc )][G(coc )]| = 0 , and the nontrivial solution of Eq. (11) exists. Where A is the eigenvalue of Eq. (12), and the vibration system has the critical chatter frequency. Experimental modal analysis was performed to extract the modal parameters of the helical milling tool-machine tool system. The frequency response function (FRF) [G (o)] of the milling tool attached to a Mikron 710 machining centre was identified through hammer impact modal tests, as described in the literature [20]. In order to describe the directional dynamic helical milling force coefficients [A^ , the zero order approximation of its Fourier series expansion was applied, and a single frequency solution derived. The predicted stability lobe identifying the critical axial depth of cut (aplim) and spindle speed to construct the lobes can be determined, as depicted in Fig. 4. 2.2 The Chatter Limitation of Radial Depth of Cut The Z and C components of cutting force (FZ, T) on the cutting tools mainly result from the interaction of the end-cutting edges, because the end-cutting edges engage predominantly in the axial feed direction and create continuous cutting. Hence, the radial depth of cut is affected by the Z and C components of cutting forces. Considering the influence of the dynamic chip thickness on the axial vibration and torsional vibration, as presented in the literature, the regenerative chatter in the Z and C direction is the result of the variation in the dynamic chip thickness with period r0 as depicted in Fig. 5. The variation of dynamic chip thickness can be described as: h (p) = hz -Ahv= (AZ (cp;) - fza Ap; / O p )/cos0.(13) a) \ 1 i -c ^ a ] ^ "A -J" 0p b) Fig. 5. Dynamic chip thickness variation in a) Z and b) C direction Neglecting the edge effect on the end-cutting edges, and substituting dynamic chip thickness into the end-cutting edge force model results in the following relation: dFa a dT* Frequency [rev/min] Fig. 4. The stability lobes prediction of axial depth of cut AZ 1 . A„,} * - * * tan0 + k- -((0 + k- )fj®p r - k-r tsad -(r - k-r tan0))/®^ (14) 720 Liu, C. - Wang, G. - Dargusch, M.S. Strojniski vestnik - Journal of Mechanical Engineering 60(2014)11, 716-724 where, K*nc = k*K*c. Integrating Eq. (15), the dynamic cutting forces on an end-cutting edge can be described as, A ln r + B + k*r -(A ln| r + B\ + k*r )) © ; -2 r2- k* A (r - B ln |r + B| ) -( r2- k* A (r - B ln |r + B| )))© t AZ AÇj (15) D Sum up all the flutes and dynamic cutting force components on the cutting tool (FZ, T) as follows: Nm \ F J =1 I J K\c [A9(t)] AZ Aç (16) where the directional dynamic helical milling force coefficients are: Nm j=1 i1 -fj*f A ln + k ^ma,, - t a„2 - k' Aa + k' AB ln (17) The Z and C components of the displacement are: w»- &! f (•)}■ K - $ {R iffl)i- {t ;"[{ • f i"»- f • fi")}- £;")}-{("»-{R» i"»- -(1 - ){R i")} - e" )[g (")] {f (")}. (18) Substitute Eq. (18) to Eq. (16), for the chatter frequency rnc, then: {F (c )} = K \c (1 - ec)[ Az ,(ce )]• •[G (coc )] {F (COc )}. (19) The characteristic function of Eq. (19) is: [I]-Kvc (1 - e-1^ )[Am(coc )][G (COc )] = 0. (20) Theoretically, similar to the aforementioned resolution approach to the chatter limitation of ap, when det [l] + A[AZ0(coc)][G(coc)] = 0, the non- trivial solution of Eq. (19) exists, A is the eigenvalue of Eq. (20). However, for the directional dynamic helical milling force coefficients [AZ4 (tthe situation is more complicated, and the variable ae is not linear. In the range from 0 to the end-edge length (Dm/2, half the tool diameter), [AZ4 (t)] can be simplified to: [ m< )]=: ' Ka, Dm 1 2 j=1 \ — a, - 2 a, 1 -f»/*p}. (21) If the influence of the fza was neglected, the procedure to obtain the approximate solution would be identical with that in literature [9]. It is then necessary to solve the eigenvalue and obtain the predicted stability lobe that identifies the critical radial depth of cut (aelim) and the spindle speed, as shown in Fig. 6. 500 1000 1500 2000 2500 3000 3500 4000 Frequency [rpm] Fig. 6. The s tability lobes prediction of radial depth of cut 3 EXPERIMENTAL METHODS The chatter-free machining experiments of helical milling process were performed under dry conditions. The radial cutting depth of ae, cutting velocity of vc, Mechanics and Dynamics of Helical Milling Operations 721 D 2 D 2 D 2 Strojniski vestnik - Journal of Mechanical Engineering 60(2014)11, 716-724 axial feed rate of fza, and tangential feed rate of fzt are independent processing parameters, on which spindle frequency n and axial cutting depth ap depend. The experiments were designed to orthogonal trails with three variable factors (ae, vc, fza) where each factor had three levels. The range of n, ap and ae based on the solved stability lobe diagrams shown in Figs. 4 and 6, along with the recommendations of the milling tool manual and the consideration that the processing parameters would be suitable for the finishing operation. All the processing parameters were depicted in Table 1. The rough wrought Ti6Al4V alloy blank has the geometry of a cuboid (160 mm width, 160 mm length, and 20 mm height) with a hole in a diameter of 60 mm in the centre as depicted in Fig. 1. M.A. Ford 20 mm 5-flute flat end cylindrical carbide milling tools (17878703A) were employed. The occurrence of chatter could be detected by measuring the amplitude of the cutting forces, as well as the machining accuracy (surface roughness, roundness, and cylindricity). Once chatter occurs, the cutting force amplitude increases drastically, the surface roughness escalates, and the geometric accuracy of the hole is lost. Experiments were performed on a five-axis Mikron UCP-710 machining centre. A Kistler 9265B piezo-electric dynamometer, a Mahr S3P 2D surface profile tester, and a Hexagon 3D coordinate measuring machine were setup to measure the cutting forces, machined surface roughness, roundness and cylindricity, respectively. Table 1. Helical milling experiment parameters Test No. Vc [m/min] n [rpm] ae [mm] ap [mm] fza [mm/tooth] fzt [mm/tooth] 1 100 1591 0.25 1.5 0.48 0.060 2 100 1591 0.5 1.5 0.60 0.075 3 100 1591 0.75 1.5 0.72 0.090 4 120 1910 0.25 1.5 0.72 0.090 5 120 1910 0.5 1.5 0.48 0.060 6 120 1910 0.75 1.5 0.60 0.075 7 135 2149 0.25 1.5 0.60 0.075 8 135 2149 0.5 1.5 0.72 0.090 9 135 2149 0.75 1.5 0.48 0.060 4 EXPERIMENTAL RESULTS AND DISCUSSION The inspection of the operation accuracy, including the surface roughness of the inner surface of the hole, the roundness of the hole cross section, and the cylindricity of the hole were measured, and shown in Table 2. Table 2. Helical milling experiment results Test No. Surface roughness M Roundness M Cylindricity M 1 < 1.6 3.6, 4.7, 5.4 7.01 2 < 1.6 4.8, 5.2, 6.3 8.66 3 < 1.6 6.2, 6.9, 7.5 8.95 4 < 1.6 3.8, 5.1, 7.3 9.72 5 < 1.6 4.5, 5.4, 7.9 11.2 6 < 1.6 4.3, 5.6, 8.1 11.9 7 < 1.6 8.3, 9.5, 10.3 12.3 8 < 1.6 8.0, 8.0, 10.9 15.2 9 < 1.6 8.4, 8.7, 12.0 15.3 During the experimental process, chatter vibration between the cutter and work piece did not occur. In order to assess the chatter-free conditions during the experiment, the cutting forces signal was analysed to confirm that there was no abnormal forces impact and fluctuation. The experimental and simulated cutting force results of helical milling with typical cutting conditions, cutting speed vc 100 m/min, axial feed rate fza 0.5 mm/tooth, tangential feed rate fzt 0.1 mm/tooth, were compared in Fig. 7. The amplitude of the experimental cutting forces in comparison to the simulated cutting forces, which were typical stable force signals in Fig. 7a and b, has a margin of error within 10%. No impact was detected. The amplitude spectrum of the experimental result shown in Fig. 7c properly matches the simulated cutting forces shown in Fig. 7d. The unmatched spectrum in the experimental cutting forces might be due to the serrate chip formation and other environmental influences. From the observation of the experiment processes and the measurement of the machining accuracy and the cutting forces, it could be concluded that the given helical milling processing parameters can deliver the needed machining accuracy, and the operation can be chatter free within the stability limitation predicted by the presented analysis and model. 5 CONCLUSION Based on the modelling of cutting forces, the machining dynamics and the chatter stability of the helical milling process have been modelled in both the time and frequency domains. The cutting forces model predicts both the interaction of the side-cutting 722 Liu, C. - Wang, G. - Dargusch, M.S. Strojniski vestnik - Journal of Mechanical Engineering 60(2014)11, 716-724 X direction E^lOOr éi 75L 8 50 S 25 fe 0„L 100 150 Time [s] Z direction 100 150 Time [s] c) d) Fig. 7. Experiment and simulation of the helical milling cutting forces results; a) experiment; b) simulation; c) amplitude spectrum by Fourier transform of experiment; and d) amplitude spectrum by Fourier transform of simulation a) b) edges with the work piece and the end-cutting edges with the work piece along the helical feed path. The helical milling dynamics model represents the special tangential feed and axial feed. This model includes the dual periodicity which is caused by spindle rotation, as well as the period of the helical feed of the cutting tool. The frequency domain solutions of the critical axial cutting depth and critical radial cutting depth with spindle velocity have also been identified. The 4 DOM vibration system has been decomposed to two 2 DOM systems and the critical axial cutting depth, and critical radial cutting depth have been solved separately. On the given machining condition of the helix milling operation, the experimental results indicate that the simulation from the presented models effectively forecasts the cutting forces and chatterfree process parameters. Chatter did not occur during machining, and the helical milling operation satisfied the required surface roughness, roundness and cylindricity. 6 ACKNOWLEDGEMENTS This work was supported by the Nanjing University of Aeronautics & Astronautics Fundamental Research Funds [grant number NS2013048]. 7 REFERENCES [1] Iyer, R., Koshy, P., Ng, E. (2007). Helical milling: An enabling technology for hard machining precision holes in AISI D2 tool steel. International Journal of Machine Tools & Manufacture, vol. 47, no. 2, p. 205-210, DOI: 10.1016/j.ijmachtools.2006.04.006. [2] Wang, H., Qin, X., Li, H., Ren, C. (2013). Analysis of cutting forces in helical milling of carbon fiber-reinforced plastics. Proceedings of the Institution of Mechanical Engineers, Part B: Journal of Engineering Manufacture, vol. 227, no. 1, p. 62-74, D0I:10.1177/0954405412464328. [3] Denkena, B., Boehnke, D., Dege, J.H. (2008). Helical milling of CFRP-titanium layer compounds. CIRP Mechanics and Dynamics of Helical Milling Operations 723 Strojniski vestnik - Journal of Mechanical Engineering 60(2014)11, 716-724 Journal of Manufacturing Science and Technology, vol. 1, no. 2, p. 64-69, D0L10.1016/j.cirpj.2008.09.009. [4] Larue, A., Altintas, Y. (2005). Simulation of flank milling processes. International Journal of Machine Tools and Manufacture, vol. 45, no. 4-5, p. 549-559, D0L10.1016/j.ijmachtools.2004.08.020. [5] Altintas, Y., Engin, S. (2001). Generalized Modeling of Mechanics and Dynamics of Milling Cutters. CIRP Annals - Manufacturing Technology, vol. 50, no. 1, p. 25-30, D0I:10.1016/s0007-8506(07)62063-0. [6] Ozturk, E., Budak, E. (2007). Modeling of 5-axis milling processes. Machining Science and Technology: An International Journal, vol. 11, no. 3, p. 287 - 311, D0I:10.1080/10910340701554808. [7] Lazoglu, I., Liang, S.Y. (2000). Modeling of ballend milling forces with cutter axis inclination. Journal of Manufacturing Science and Engineering-Transactions of the ASME, vol. 122, no. 1, p. 3-11, D0I:10.1115/1.538885. [8] Dang, J.-W., Zhang, W.-H., Yang, Y., Wan, M. (2010). Cutting force modeling for flat end milling including bottom edge cutting effect. International Journal of Machine Tools and Manufacture, vol. 50, no. 11, p. 986-997, D0I:10.1016/j.ijmachtools.2010.07.004. [9] Ko, J.H., Altintas, Y. (2007). Dynamics and stability of plunge milling operations. Journal of Manufacturing Science and Engineering-Transactions of the ASME, vol. 129, no. 1, p. 32-40, D0I:10.1115/1.2383070. [10] Altintas, Y., Ko, J.H. (2006). Chatter stability of plunge milling. CIRP Annals-Manufacturing Technology, vol. 55, no. 1, p. 361-364, D0I:10.1016/S0007-8506(07)60435-1. [11] Kolar, P., Sulitka, M., Janota, M. (2011). Simulation of dynamic properties of a spindle and tool system coupled with a machine tool frame. International Journal of Advanced Manufacturing Technology, vol. 54, no. 1-4, p. 11-20, D0I:10.1007/s00170-010-2917-7. [12] Mahnama, M., Movahhedy, M.R. (2010). Prediction of machining chatter based on FEM simulation of chip formation under dynamic conditions. International Journal of Machine Tools & Manufacture, vol. 50, no. 7, p. 611-620, D0I:10.1016/j.ijmachtools.2010.03.009. [13] Chaari, R., Abdennadher, M., Louati, J., Haddar, M. (2011). Modelling of the 3D machining geometric defects accounting for workpiece vibratory behaviour. International Journal of Simulation Modelling, vol. 10, no. 2, p. 66-77, D01:10.2507/ijsimm10(2)2.173. [14] Ostasevicius, V., Gaidys, R., Dauksevicius, R., Mikuckyte, S. (2013). Study of vibration milling for improving surface finish of difficult-to-cut materials. Strojniski vestnik-Journal of Mechanical Engineering, vol. 59, no. 6, p. 351-357, D0I:10.5545/ sv-jme.2012.856. [15] Bajic, D., Celent, L., Jozic, S. (2012). Modeling of the influence of cutting parameters on the surface roughness, tool wear and the cutting force in face milling in off-line process control. Strojniški vestnik -Journal of Mechanical Engineering, vol. 58, no. 11, p. 673-682, D01:10.5545/sv-jme.2012.456. [16] Yoon, M.C., Chin, D.H. (2005). Cutting force monitoring in the endmilling operation for chatter detection. Proceedings of the Institution of Mechanical Engineers Part B-Journal of Engineering Manufacture, vol. 219, no. 6, p. 455-465, D0I:10.1243/095440505x32292. [17] Polli, M.L., Weingaertner, W.L., Schroeter, R.B., Gomes, J.d.O. (2012). Analysis of high-speed milling dynamic stability through sound pressure, machining force and tool displacement measurements. Proceedings of the Institution of Mechanical Engineers Part B-Journal of Engineering Manufacture, vol. 226, no. A11, p. 1774-1783, D0I:10.1177/0954405412460128. [18] Liu, C., Wang, G., Dargusch, M. (2012). Modelling, simulation and experimental investigation of cutting forces during helical milling operations. The International Journal of Advanced Manufacturing Technology, vol. 63, no. 9-12, p. 839-850, D0I:10.1007/ s00170-012-3951-4. [19] Wang, H., Qin, X., Ren, C., Wang, Q. (2012). Prediction of cutting forces in helical milling process. The International Journal of Advanced Manufacturing Technology, vol. 58, no. 9-12, p. 849-859, D0I:10.1007/ s00170-011-3435-y. [20] Damir, A., Ng, E.G., Elbestawi, M. (2011). Force prediction and stability analysis of plunge milling of systems with rigid and flexible workpiece. International Journal of Advanced Manufacturing Technology, vol. 54, no. 9-12, p. 853-877, D0I:10.1007/s00170-010-2982-y. 724 Liu, C. - Wang, G. - Dargusch, M.S.