Analytical modeling of a thinwalled cylindrical workpiece during the turning process. Stability analysis of a cutting process
Artem Gerasimenko^{1} , Mikhail Guskov^{2} , Alexander Gouskov^{3} , Philippe Lorong^{4} , Alexander Shokhin^{5}
^{1, 3}Bauman Moscow State Technical University, 2nd Baumanskaya, 5, Moscow 105005, Russia
^{1, 2, 4}PIMM Laboratory, Arts et Metiers ParisTech, 151 bd de l’Hopital, Paris 75013, France
^{3, 5}IMASH RAN, 4 Maly Kharitonyevsky Pereulok, Moscow 101990, Russia
^{1}Corresponding author
Journal of Vibroengineering, Vol. 19, Issue 8, 2017, p. 58255841.
https://doi.org/10.21595/jve.2017.18061
Received 28 November 2016; received in revised form 6 February 2017; accepted 21 February 2017; published 31 December 2017
Copyright © 2017 JVE International Ltd. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
JVE Conferences
The purpose of this work is to develop a mathematical model of the dynamics of turning a thinwalled cylindrical shell. This model uses a finite number of degrees of freedom and takes into account the variability of dynamic compliance. It is then possible to obtain estimates of the boundaries of stability of the continuous cutting process. The model is based on the theory of shells with application of Galerkin’s method in conjunction with the expansion of the displacement field in beam and trigonometric functions. On the basis of the developed model, an algorithm designed for constructing boundaries of stability of turning of the thinwalled cylindrical parts is presented and compared to experimental results. A strategy to define matter removal sequences is proposed.
Keywords: cylindrical shell, chatter, turning, selfoscillations.
1. Introduction
Investigating the stability of the cutting process has been the subject of studies for many years now. The first works appeared at the beginning of the 50s and explained the basic mechanism of occurrence of selfoscillations in the machining process “over the trace on the previously cut surface”. Currently, this effect is considered the main cause of chatter vibrations. This regenerative mechanism of excitation of oscillations is widely used in the method of analysis of the stability of the stationary cutting process. This method is based on the analysis of the characteristic equation with the further construction of the stability diagram [1, 2]. Most of the works devoted to the study of the stability of the process of turning use as an object a cylindrical full workpiece [3, 4] or beams [5, 6]. However, in [7], the author conducted an experiment of turning thinwalled cylindrical part and showed the presence of areas that were vibrating during the process. More recently a stability study of turning aluminum tubes, finite element modeled, for a variety of processing conditions was carried out [8, 9]. A time domain simulation of the turning process is presented in [10] where the influence of damping on emergence of unstable cutting process was shown. All these works show the influence of damping on the occurrence of vibrations during cutting. Another feature of the turning of thinwalled cylindrical parts is a strong change in the dynamic stiffness and in the natural frequencies of the workpiece during the pass.
In this paper, the turning of a thinwalled cylindrical shell using an analytical approach to model the dynamics is studied. This analytical approach is based on the KirhgoffLove shell theory. The assumption of the middle surface inextensibility is used [11]. The system of differential equations of dynamics of the shell, expressed in terms of the axial, circumferential and radial components of displacement, was conducted by using variational principles. The research of solutions of this system is performed by Fourier expansion of the components of displacement along the circumferential direction and beam functions in axial direction [1214] and uses the Galerkin’s method [1517]. The obtained analytical model allows estimating the boundaries of stability zones for the continuous cutting process. It takes into account the continuous changing of dynamics due to matter removal and motion of the tool over a flexible part. An experimental validation is presented and a strategy to optimize the turning process is proposed.
In Section 2, we present the analytical model that contains a finite number of degrees of freedom. In Section 3 we present the way we use to estimate the stability of the process and in Section 4 we describe the design of the experiment uses for validation. The experiment was chosen in order to have an increase in compliance of workpiece while the tool progresses until chatter occurs. A presentation of the obtained results and a discussion are given in Section 5. Finally, a way to optimize to sequence of tool passes is proposed in Section 6. The conclusions take place in Section 7.
2. Analytical simulation of the turning process
A cylindrical thinwalled shell (Fig. 1) is considered. The shell is in rotation around its axis at a constant angular velocity $\mathrm{\Omega}$. The shell is loaded with a moving force, the cutting force, which has, with respect to the shell, an angular velocity –$\mathrm{\Omega}$. The left edge of the shell is fixed, the right one is free. The geometric parameters of the shell are $L\text{,}\text{}R\text{,}\text{}h$, which account for length, radius and thickness, respectively.
The following coordinates systems (see Fig. 1) are used: $\mathbf{X}\mathbf{Y}\mathbf{Z}$ mobile system of coordinates, linked to the rotating part; $\stackrel{\u02d8}{\mathbf{X}}\stackrel{\u02d8}{\mathbf{Y}}\stackrel{\u02d8}{\mathbf{Z}}$ fixed system of coordinates, connected to the frame of the machine tool.
Fig. 1. Schema of the model
2.1. Basic system of equations of the thinwalled shell
The system of differential equations of dynamics of the shell, expressed in terms of the components $u(x,\phi ,t)$, $v(x,\phi ,t)$, $w(x,\phi ,t)$, on axial, circumferential and radial displacement, can be represented as the follows [18]:
The expression of the differential operators ${L}_{ij}^{\text{'}}$ are:
${L}_{12}^{\text{'}}=\frac{1+\mu}{2R}{\partial}_{x}{\partial}_{\phi},{L}_{23}^{\text{'}}=\frac{1}{{R}^{2}}{\partial}_{\phi}\frac{{h}^{2}}{12{R}^{4}}{\partial}_{\phi}^{3},{L}_{13}^{\text{'}}=\frac{\mu}{R}{\partial}_{x},$
${L}_{21}^{\text{'}}={L}_{12}^{\text{'}},{L}_{31}^{\text{'}}={L}_{13}^{\text{'}},{L}_{32}^{\text{'}}={L}_{23}^{\text{'}},$
where ${\partial}_{*}=\partial /{\partial}^{*}$ denotes the partial derivative with respect to $*$, ${p}_{u}$, ${p}_{v}$, ${p}_{w}$ represent the components of the surface forces acting on the shell in the axial, circumferential and radial directions, respectively.
The rotationinduced load terms are expressed as follows:
As shown in [19], these components have a small impact on the dynamics of the shell, and can be neglected in further calculations.
The components of the surface forces, applied to the shell, modelled by an interaction with the tool, are expressed as follows:
where ${x}_{P}$, ${\phi}_{P}$ are the coordinates of the position of the tool, in the mobile frame, in axial and circumferential direction, respectively (both are time dependent); ${P}_{g}\left(t,\phi ,x\right)$ refers to the components of the external point cutting force, $\delta $ is the Dirac function.
2.2. Model of the cutting force
In this work, a linear cutting law (cutting force model), linearized around the nominal cutting parameters, is used. Therefore, the force, which models the effect of the tool on the shell for each direction, is expressed as follows:
where $\chi $ is the coefficient of cutting stiffness, ${p}_{0g}$ is the cutting force in nominal conditions, ${h}_{w}=\left(w\left(tT,\phi ,x\right)w\left(t,\phi ,x\right)\right)$ represents dynamic perturbation of the cutting depth, due to vibrations, $T$ is the period of rotation and ${e}_{u}$, ${e}_{v}$, ${e}_{w}$ are the direction cosines of the force vector.
To study the stability of the cutting process, the system is considered in a perturbed state and therefore the part of the cutting forces associated with the perturbation, which is used, can be written as follows:
2.3. Discretization of the system of equations
In order to find the solution of the Eq. (1), a Fourier expansion of the components of displacement along the axial and circumferential coordinate is carried out, with the further use of the Galerkin’s method.
The displacement field is written as follows:
where ${\mathrm{\Psi}}_{mk}^{u}\left(x\right)$, ${\mathrm{\Psi}}_{mk}^{v}\left(x\right)$, ${\mathrm{\Psi}}_{mk}^{w}\left(x\right)$ are the unknown axial basic functions, $k\in \{1,\dots ,N\}$ is the number of harmonics in the circumferential direction, $m\in \{1,\dots ,M\}$ is the number of harmonics in the axial direction, ${\alpha}_{km}^{c}\left(t\right)$, ${\alpha}_{km}^{s}\left(t\right)$ are unknown amplitude functions of time.
The axial basic functions are beam functions written as follows [14, 15, 20]:
The parameters ${\lambda}_{m}$, ${C}_{1}$, ${C}_{2}$, ${C}_{3}$, ${C}_{4}$ are determined according to a set of boundary conditions.
The relationship between functions ${\mathrm{\Psi}}^{u}\left(x\right)={\mathrm{\Psi}}^{v}\left(x\right)={\mathrm{\Psi}}^{w}\left(x\right)$ is determined by the state of nonextensible middle surface, namely ${\epsilon}_{\phi}=0$, ${\gamma}_{x\phi}=0$ [11].
According to boundary conditions, at hand clampedfree, the function ${\mathrm{\Psi}}^{w}$ will have the following form [14, 15, 20]:
Note that the ${\lambda}_{m}$, ${\sigma}_{m}$ do not depend on the number of circumferential harmonics $k$.
The final system of equations of motion of the shell is presented in the fixed coordinate system $\stackrel{\u02d8}{u}$, $\stackrel{\u02d8}{v}$, $\stackrel{\u02d8}{w}$, because of simplicity of representation of the point of interaction (point A, Fig. 2) between the tool and the shell for further analysis of stability.
The components of displacement of the shell in the moving frame are $u$, $v$, $w$ and $\stackrel{\u02d8}{u}$, $\stackrel{\u02d8}{v}$, $\stackrel{\u02d8}{w}$ are the components of this displacement in the fixed frame. The angle ${\psi}_{P}$ gives the position of the tool in the fixed frame. This angle is not time dependent.
Fig. 2. Relation between the different coordinate systems
2.4. Writing systems in matrix form
The final system of differential equations of dynamics of the shell is written in matrix form for the unknown amplitude dimensionless variables ${\alpha}_{km}^{c}\left(\tau \right)$, ${\alpha}_{km}^{s}\left(\tau \right)$ depending on the time and for one particular circumferential harmonic $k$:
The vector of solution ${\underset{\_}{\mathbf{q}}}_{Mk}\left(\tau \right)$ in Eq. (10) consists of a set of vectors ${\underset{\_}{\mathbf{d}}}_{1k}^{T}\left(\tau \right),\dots ,{\underset{\_}{\mathbf{d}}}_{Mk}^{T}\left(\tau \right)$ for each axial harmonic $m\in \{\mathrm{1,2},3,\dots ,M\}$ and one particular circumferential harmonic $k$, and these vectors consist of a set of symmetric ${\alpha}_{mk}^{s}\left(\tau \right)$ and antisymmetric ${\alpha}_{mk}^{c}\left(\tau \right)$ amplitude components:
where ${\alpha}_{mk}^{c}\left(\tau \right)$, ${\alpha}_{mk}^{s}\left(\tau \right)$ are unknown amplitude functions of time.
The vector of load ${\underset{\_}{\mathbf{p}}}_{k}$ (dimension $2M\times 1$) is written:
The matrices ${\underset{\_}{\underset{\_}{\mathbf{M}}}}_{k}$, ${\underset{\_}{\underset{\_}{\mathbf{K}}}}_{k}$, ${\underset{\_}{\underset{\_}{\mathbf{P}}}}_{Mk}^{f}$ are derived via Galerkin method. Their structures and their components can be found in the appendix. All these matrices are dependent on the tool position along the axial coordinate ${x}_{P}$. Changes in the thickness of the shell is taken into account with the help of the integrals Eq. (29). Also, dynamic compliance variation of the cutting process is taken into account by changing the position of the tool Eq. (32).
The matrix ${\underset{\_}{\underset{\_}{\mathbf{C}}}}_{k}$ is determined with help of the Rayleigh damping model:
The constant of Rayleigh $\alpha $, $\beta $ have been determined from experimental results given in Table 2 In the present case we have $\alpha =\text{7.25}$, $\beta =\text{6.71}{\text{0}}^{\text{8}}$ for ${\varsigma}_{1}=\text{0.0007}$, ${\varsigma}_{2}=\text{0.0008}$ and their corresponding frequencies ${f}_{1}=\text{1661}$ Hz, ${f}_{2}=\text{2798}$ Hz.
The system of Eq. (10) describes the turning process of a thinwalled cylindrical shell whose interaction with the tool is defined by a cutting force concentrated on a mobile point.
3. Stability evaluation of cutting process
The stability of the turning process is investigated, using the previous model and the procedure described in the book of Cheng [2].
The solution of Eq. (10) will be sought in the form:
After substituting Eq. (14) into Eq. (10), the following equation is obtained:
where:
At the boundaries of the stability areas, the $\lambda $ is purely imaginary. So, the critical parameters of the cut are found via characteristic Eq. (15):
By substituting $\lambda =i2\pi {f}_{c}$, the characteristic equation is divided in real and imaginary parts. The parameter ${f}_{c}$ represents the oscillation frequency at the limit of the zone of stability, when the characteristic factor $\lambda $ becomes purely imaginary. Subsequently, the oscillation frequency ${f}_{c}$ will be called chatter frequency.
It is possible to combine the eigenvalues of Eq. (17) into a vector ${\underset{\_}{\mathbf{\Lambda}}}_{k}=\left({\mathrm{\Lambda}}_{k}^{1}\dots {\mathrm{\Lambda}}_{k}^{2M}\right)$. Each component of this vector leads, according to Eq. (16), to the resolution of:
where ${\mathrm{\Lambda}}_{k}^{j}$ is the component of the eigenvalues vector, which is determined from the solution of the eigenvalue problem, i.e $\left[{\underset{\_}{\underset{\_}{\mathbf{L}}}}_{k}{\mathrm{\Lambda}}_{k}^{j}{\underset{\_}{\underset{\_}{\mathbf{P}}}}_{f}\right]{\underset{\_}{\mathbf{v}}}_{j}=0\text{,}$$j=1,\dots ,2M.$ The coefficients ${a}_{k}^{j}$ and ${b}_{k}^{j}$ represent the real and imaginary components of ${\mathrm{\Lambda}}_{k}^{j}$.
This leads to:
The Eq. (19) is divided in the real part and imaginary part:
From the resulting system of Eq. (20), the stiffness of cutting is expressed [2]:
and also [22]:
where ${\psi}_{k}^{j}$ is the angle phase.
The number of oscillation per revolution [14] is given by:
where $z$, ${\epsilon}_{k}^{j}/2\pi $ are the integer and the fractional number of incisions on the surface, respectively.
Taking into account Eqs. (22) and (23) leads to:
The solution of Eq. (20) is a pair of parameters $\left(\stackrel{~}{T},\stackrel{~}{\chi}\right)$ for each set of frequency value ${f}_{c}\in \left({f}_{c0},{f}_{cf}\right)$ and position of the tool ${x}_{P}$, where the parameters ${f}_{c0}$, ${f}_{cf}$ are arbitrary constants.
This approach gives place to the following algorithm:
• On the first step of the algorithm, the frequency parameter ${f}_{c}={{f}_{c}}_{0}$ is fixed (the smallest eigenvalue of the eigenvalue problem of conservative system ${\underset{\_}{\underset{\_}{\mathbf{M}}}}_{k}{\underset{\_}{\ddot{\mathbf{q}}}}_{Mk}+{\underset{\_}{\underset{\_}{\mathbf{K}}}}_{k}{\underset{\_}{\mathbf{q}}}_{Mk}=0$).
• In a second step, the eigenvalue problem is solved, namely $\left[{\underset{\_}{\underset{\_}{\mathbf{L}}}}_{k}{\mathrm{\Lambda}}_{k}^{j}{\underset{\_}{\underset{\_}{\mathbf{P}}}}_{f}\right]{\underset{\_}{\mathbf{v}}}_{j}=0$, $j=1,\dots ,2M$, from which the vector of eigenvalues ${\underset{\_}{\mathbf{\Lambda}}}_{k}=\left({\mathrm{\Lambda}}_{k}^{1}\dots {\mathrm{\Lambda}}_{k}^{2M}\right)$ is determined.
• Then, for each eigenvalue ${\mathrm{\Lambda}}_{k}^{j}$, using the Eqs. (21), (22) and (24), the following parameters ${\stackrel{~}{\chi}}_{k}^{j}$, ${\psi}_{k}^{j}$ and ${\epsilon}_{k}^{j}$ are calculated. The final step enables to calculate the period of rotation ${\stackrel{~}{T}}_{z}^{j}$.
• The cycle is repeated for the next value of the frequency parameter ${f}_{c}={{f}_{c}}_{0}+\u2206{{f}_{c}}_{0}$ (where $\u2206{{f}_{c}}_{0}$ is an arbitrary constant).
The obtained data set enables the drawing of a multiDof lobe stability diagram.
4. Experimental setup and measurements
As an illustration of the above designed approach, an experiment has been carried out, namely. The turning of the outer cylindrical surface of a tubular part, as can be seen on the drawing (Fig. 3), was conducted. The experimental set up consists of a part, fixed by a threejaw chuck to the machine tool spindle, and a set of sensors. This experiment has been partly addressed in [21, 22].
In the framework of the experiment, several types of sensors were used. The experimental set up is shown in Fig. 4. For the vibration measurement of the workpiece two types of sensor are used: i) displacement eddy current sensors (3a, 3b, Fig. 4) and ii) accelerometers (2a, 2b, 2c, 2d, Fig. 4). For the measurement of the frequency of rotation of the workpiece an optical sensor was used (1a, Fig. 4).
On Fig. 4, one can also see the eddy current sensor ’3a’ but, due to presence of chips inside the tube, it was not exploitable. The accelerometer ’2c’ was positioned on ’3a’ sensor support to measure the vibrations of this support.
The pass of the tool is repeated several times until occurrence of vibrations of large amplitudes generating significant damage on the machined surface. The generated defects on the machined surface during the last tool pass are shown on Fig. 5. This occurred for a thickness of the tube going from 5.4 mm to 4.4 mm.
Fig. 3. Geometry of the machined tube
Fig. 4. Experimental setup
Fig. 5. Surface of the tube after the last tool pass
An impact test was carried out before and after the machining. The main aim of this test was the definition of modal characteristics of the workpiece before and after the turning process. The respective transfer functions are presented on Fig. 6. The experimental data given in Table 1 and Table 2 are coming from these transfer functions.
Fig. 6. Transfer function before and after the last tool path
Table 1. Experimental eigenfrequencies before the passage of the tool
Freq. experiment (Hz)

1108.4

1937.2

3384.5

Modal Damping

0.0105

0.0076

0.00051

Modal form

2 lobes

3 lobes

Flexion 2

Table 2. Experimental eigenfrequencies after the passage of the tool
Freq. experiment (Hz)

1100.2

1661.1

2798

Modal Damping

0.0077

0.00070

0.00078

Modal form

2 lobes

3 lobes

4 lobes

The signal of the eddy current sensor ’3b’, fixed on the machine frame, is presented on Fig. 7. This sensor is placed in front of the right extremity of the workpiece where the largest displacements are present.
On the right side ($t>\text{45}\text{}\text{s}$) of the graphic (Fig. 7) the observed oscillations correspond to those visible on the workpiece surface on Fig. 5. On the spectrogram computed from this displacement signal (Fig. 8) one can see the evolution of the resonances of the tube close to the evolution of its eigenfrequencies. On the right of the Fig. 8 the result of the hammer test after the last tool pass is given (solid line from Fig. 6).
And, as shown by the analysis of the spectrogram and FRF (Fig. 8), the dominant vibration component is close to a natural frequency of the workpiece.
Fig. 7. Signal recorded using the sensor 3b
Fig. 8. Spectrogram issued from the measurement of the eddy current sensor ’3b’ with the result of hammer test after the final pass
5. Results and discussion
As previously said, the solution of Eq. (20) is the pair of parameters $\left(\stackrel{~}{T},\stackrel{~}{\chi}\right)$ for each frequency value ${f}_{c}$ and depends on the tool position ${x}_{P}$. We build the stability curves for different positions of the tool ${x}_{P}$, with the following notations: cutting stiffness coefficient $\kappa =\chi /{\chi}_{\mathrm{e}\mathrm{x}\mathrm{p}}$ and rotation frequency ${f}_{rot}=1/T$, with: $\chi $ – cutting stiffness of the system, [N/m]; ${\chi}_{\mathrm{e}\mathrm{x}\mathrm{p}}$ – nominal cutting stiffness, obtained from experiment (here, ${\chi}_{\mathrm{e}\mathrm{x}\mathrm{p}}=$569 N/m), for the toolworkpiece set studied, [N/m];
The Fig. 9 shows the threedimensional stability diagram of the continuous cutting process in the threeparameter space $\left({x}_{P},\kappa ,{f}_{rot}\right)$. This graph shows that conditions of stability change, depending of the tool position. The number of different degrees of freedom in the system considered is the number of sets of lobes representing the limits of stability, highlighted by different colors.
Fig. 9. Stability diagram in the parameter space $\left({x}_{P},\kappa ,{f}_{rot}\right)$
The Fig. 10(а) represents a clipping plan of the stability diagram in the twoparameters $\left(\kappa ,{f}_{rot}\right)$ plane, for the position of the tool ${x}_{p}=\text{0.739}L$. It corresponds to the loss of stability of the quasistatic cutting process. For the considered state of the system, and the specified range of technological parameters, three sets of lobes are obtained, each of the sets having a different color. In practice, the border of the lower zone of stability is our only concern because, from the moment the parameters of the system reach the lower stability boundary, the continuous cutting process becomes unstable. The red vertical line represents the rotation frequency of the workpiece, corresponding to the experiment ${f}_{rot}=\text{12,63}$ Hz. On the same line, corresponding to the rotational frequency of interest, there are several different sets of lobes representing a limit of the stability zone. The red vertical line intersects the first set of lobes in point А. To evaluate the stability of the quasistatic cutting process, the relative position of the point А (the intersection point of the rotational frequency of interest with the lowest set of lobes representing the border of the stability area) has to be estimated by comparing it to the limit value of the cutting stiffness coefficient $\kappa =\chi /{\chi}_{\mathrm{e}\mathrm{x}\mathrm{p}}=\text{1}$ (horizontal blue line). If the point A is above the blue line, then the quasistatic cutting process is stable. If it is below the blue line, the process is unstable. Fig. 10(a) shows that the quasistatic cutting process is located on the border of the zone of stability. One of the characteristics of the boundary of the stability zone of the cutting process consists in the oscillation frequency, called chatter frequency. It is at this frequency that the selfexcited oscillations of the system begin, if the process parameters are within the instability area of the system. The chatter frequency may vary depending on the encountered stability zone and on the positions with respect to this zone. On the Fig. 10(b) is shown the variation of the critical stiffness parameter of the system depending on the chatter frequency ${f}_{chatter}$.
To determine the chatter frequency, which corresponds to the oscillation frequency on the border of the zone of stability for the rotational frequency considered, the value of parameter $\kappa $ must be known at the intersection with the boundary of the zone of stability at the relevant rotational frequency (the point А, Fig. 10(а)). For the lowest limit of the stability area, which corresponds to the black lines in Fig. 10(a), the value of the parameter is ${\kappa}_{A}=\text{1}$ (Fig. 10(а)). Then, reporting this value ${\kappa}_{A}=\text{1}$ on the Fig. 10(b), the value of the chatter frequency is obtained: ${f}_{chatter}=\text{1930}$ Hz. This means that at the limit of the zone of instability selfoscillations begin to develop at the frequency ${f}_{chatter}=\text{1930}$ Hz. It is important to note that each set of boundaries of the stability area has its own chatter function, as shown on Fig. 11(a) and Fig. 11(b).
The stability boundaries for the position of the tool ${x}_{P}=\text{0.574}\text{}\text{L}$ in the plane of parameters $\left(\kappa ,{f}_{rot}\right)$ are presented on Fig. 11(a). There are 6 families of these boundaries in total (each assigned with a different color). On Fig. 11(b), for this position of the tool and this field of parameter $\kappa $, the evolution of the stiffness parameter of the system $\kappa $, depending on the chatter frequency, is plotted. The 6 sets of different limits of stability zone are clearly seen, in different colors (Fig. 11(а)), and each set has its own function of chatter frequency (Fig. 11(b)). Thus, each set of lobes has its own variation of the stiffness parameter of cutting and of the chatter frequency during the process. The lowest boundary of stability at given for the position of the tool ${x}_{P}=\text{0.574}L$ is the set of lobes shown in blue (Fig. 11(а)).
Fig. 10. Clipping plane of the stability diagram in the $\left(\kappa ,{f}_{rot}\right)$ and $\left(\kappa ,{f}_{chat}\right)$ planes, for the tool position ${x}_{p}=\text{0.739}L$
a) Variation of the $\kappa $ as a function of the workpiece rotational frequency
b) Variation of the $\kappa $ with respect to the chatter frequency
Fig. 11. Clipping plane of the stability diagram in the plans $\left(\kappa ,{f}_{rot}\right)$ and $\left(\kappa ,{f}_{chat}\right)$ for the tool position ${x}_{P}=\text{0.574}L$
a) Variation of the $\kappa $ as a function of the workpiece rotational frequency
b) Variation of the $\kappa $ with respect to the chatter frequency
The frequency of oscillations on the stability boundary shows some variations, depending on the position of the tool (Fig. 12). The values shown in this diagram correspond to the frequency values, obtained by the intersection of the vertical line of the rotational frequency with the boundaries of stability regions for each of the positions of the tool. It was found that the frequencies at the boundaries of stability domains are very close to the natural frequencies of the shell.
Fig. 12. Variation of the frequency of oscillations on the border of the zone of stability depending on the position of the tool
During the pass of the tool, a chip has been removed, thus leading to modifications in the dynamic characteristics of the system, including the natural frequencies of the shell. The results obtained by modeling the system during the 16th pass of the tool can be compared with the experimental data (Fig. 13).
During the experiment, it was observed that the system responded for some frequencies (Fig. 13(а)). In the range of 06000 Hz, there are 6 frequencies that lead to a response of the system.
The spectrogram during the final pass of the tool is presented on Fig. 13. The frequencies, which lead to a response of the system when the tool passes on the shell, are shown in red and numbered (Fig. 13(а)). The oscillation frequencies, shown in black and numbered on the borders of the stability areas of the continuous cutting are obtained by calculation (Fig. 13(b)). When crossing the border of the stability area of the limit stiffness parameter $\kappa $, one of the oscillation frequencies becomes the chatter frequency, from which begins the process of selfexcited oscillations.
Fig. 13. Comparison between the chatter frequency of the experiment and the frequencies calculated using the analytical model
a)
b)
The chatter frequency that dominates at the end of the pass in the experiment is 1670 Hz. The frequencies numbered 1, 2, 3 and 5 are close to the experimental result, with an error below 5 %, while the frequencies 4 and 6 are assigned with an error of 15 to 17 %. This difference shows that for $M=\text{2}$, which is the number of axial harmonics that we use for the approximation Eq. (7) introduced in the Galerkin method, is much better for modes whose axial form is close to the first harmonic ($m=\text{1}$) than for the modes whose axial shape is close to the second harmonic ($m=\text{2}$). It is possible to improve this result if we take $M>\text{2}$. Nevertheless, the chatter frequency, which is 1670 Hz and corresponds to the loss of stability, is well described by this model.
The inferred chatter frequency variation alongside the instability zone number 153 change (Fig. 15), results in a discontinuous chatter frequency evaluation. This is in qualitative agreement on the observed stairslike pattern in the signal spectrum, as shown on Fig 14.
Fig. 14. Comparison of the chatter frequency
a) Computed ${f}_{c}$ for lowest critical stifness
b) Spectrogram detailed for the frequency that dominates the response
Fig. 15. Critical relative rigidity of the system variation during the pass
Fig. 15 shows a sectional view of the stability diagram in three dimensions (plans – 2, Fig. 9) in the parameter area $\left(\kappa ,{x}_{P}\right)$. On this diagram, one can see how the limit stiffness of the system changes during the pass of the tool and therefore when the nominal cutting process becomes unstable.
Fig. 15 shows only the lower limit of the zone of stability (for the studied rotational frequency), because ultimately it is the only important one for the stability analysis of the cut. As can be seen, the first loss of stability of the cutting process takes place when the tool is at position ${x}_{p}=\text{0.739}L$. The stability diagram for this position of the tool in the parameters area $\left(\kappa ,{f}_{rot}\right)$ is shown in Fig. 10.
6. Optimizing the operations of turning
Our approach allows estimating the boundary conditions of the cutting for each pass of the tool, which lead the turning process to the vibrations. So, using this information, an estimation of the maximum cutting depth for each pass can be made, for which the cutting will take place without chatter. For example, in the experiment presented above, the last pass should have a cutting depth limit below ${A}_{p}^{\mathrm{l}\mathrm{i}\mathrm{m}}=\text{0,3}$ mm. Thus, in order to avoid vibrations during the turning process, the remaining passes and the depth of cut should be determined on the basis of this finding.
Table 3 presents a strategy of the turning without vibration to produce a part having a final thickness of ${h}_{final}=\text{4,4}$ mm.
The cutting depth specified is the maximum allowed for each of the passes, in order to determine the fastest way to get the piece with a specified wall thickness without vibration. These parameters are determined from the last pass to the first one. This means that the maximum cutting depth for each (current) pass depends on the conditions of the following passes.
Based on the approach described above, a diagram of the turning sequence can be worked out (Fig. 16).
Table 3. Technological plan of the turning without vibration
Number of the pass

Thickness before the pass
${h}_{before}$, (mm)

Thickness after the pass
${h}_{after}$, (mm)

Cutting depth limit
${A}_{p}^{\mathrm{l}\mathrm{i}\mathrm{m}}$, (mm)

N5

7,7

6,9

0,81

N4

6,9

6,2

0,68

N3

6,2

5,6

0,57

N2

5,6

5,1

0,47

N1

5,1

4,7

0,38

N

4,7

4,4

0,3

Fig. 16. Limit thickness sequence diagram of the turning
This diagram is drawn along the axes ${h}_{i}$ and ${h}_{i+1}$, which represent the current wall thickness and the future one, respectively. This plot can not only be used for the point values from Table 3, but can also be practical for intermediate values via interpolation.
For example, for a piece with an initial thickness of ${h}_{i}=\text{7,3}$ mm, a piece with a final wall thickness of ${h}_{final}=\text{5,4}$ mm is obtained. The horizontal segment OA of the diagram shows the initial thickness of the piece, this segment goes through the black line at point A, and from the intersection point, draws a vertical segment $AB$, thus leading to the maximum cutting depth, which is ${A}_{p}^{\mathrm{l}\mathrm{i}\mathrm{m}}=\text{0,8}$ mm for the first pass. The segment $CD=\text{0,65}$ mm illustrates the cutting depth limit for the second pass. For the third and final pass, the cutting depth is $EF=\text{0,45}$ mm, from which a piece with a ${h}_{final}=\text{5,4}$ mm wall thickness is obtained, illustrated by the segment $GF$.
7. Conclusions
A mathematical model of the turning was proposed and developed, with a finite number of degrees of freedom for a reliable estimate of the limits of stability of the continuous cutting. Based on the proposed model, an algorithm was developed to define the limits of the stability regions of the turning process of a thinwalled cylindrical part. It is shown that the boundaries of the stability zones of the cutting process occur at oscillation frequencies close to the natural frequencies of the part. The importance of the use of models with several degrees of freedom is enhanced by the fact that, upon progression of the tool, the dynamic characteristics of the part may vary, in particular, the vibration modes and its Eigen frequencies, which can lead to significant changes within the areas of instability. Since the frequency at which the stability of the cutting process will be lost is not known a priori, as well as the form of vibration corresponding to this frequency, several modes must be taken into account. Thus, with the help of the limit thickness sequence diagram, some recommendations on the choice of technological parameters of the cutting were made, for the turning, to avoid the vibrations observed during the cutting process and to get the final result as rapidly as possible.
Acknowledgements
The research was funded by the financial support of Ministry of Education and Science, NIR No. 9.1073.2014K under the design part of the Stateguaranteed order in scientific research area. The work was carried out with the financial support of the RFBR grant (Project No. 1658150001 NCNI_a).
References
 Altintas Y. Manufacturing Automation: Metal Cutting Mechanics. Machine Tool Vibrations, and CNC Design. Cambridge University Press, 2012. [CrossRef]
 Cheng K. Machining Dynamics: Fundamentals, Applications and Practices. Springer, 2008. [CrossRef]
 Urbikain G., López de Lacalle L. N., Campa F. J., Fernández A., Elías A. Stability prediction in straight turning of a flexible workpiece by collocation method. International Journal of Machine Tools and Manufacture, Vol. 54, Issue 55, 2012, p. 7381. [CrossRef]
 Zoltan Dombovari, Barton David A. W., Wilson R. Eddie, GaborStepan On the global dynamics of chatter in the orthogonal cutting model. International Journal of NonLinear Mechanics, Vol. 46, Issue 1, 2011, p. 330338. [CrossRef]
 Chen C. K., Tsao Y. M. A stability analysis of turning a tailstock supported flexible workpiece, International Journal of Machine Tools and Manufacture, Vol. 46, 2006, p. 1825. [CrossRef]
 Chen C. K., Tsao Y. M. A stability analysis of regenerative chatter in turning process without using tailstock. The International Journal of Advanced Manufacturing Technology, Vol. 29, 2006, p. 648654. [CrossRef]
 Arnold Rn Chatter patterns formed on the surface of thin cylindrical tubes during machining. Journal of Mechanical Engineering Science, Vol. 3, 1961, p. 714. [CrossRef]
 Mehdi K., Rigal J., Play D. Dynamic behavior of a thinwalled cylindrical workpiece during the turning process, Part 2: Experimental approach and validation. Journal of Manufacturing Science and Engineering, Vol. 124, Issues 3569, 2002, p. 58011. [CrossRef]
 Mehdi K., Rigal J., Play D. Dynamic behavior of a thinwalled cylindrical workpiece during the turning process, Part 1: Cutting process simulation. Journal of Manufacturing Science and Engineering, Vol. 124, Issue 3, 2002, p. 569580. [CrossRef]
 Lorong P., Larue A., Duarte A. P. Dynamic materials research. Advanced Materials Research, Vol. 223, 2011, p. 591599. [CrossRef]
 Matveyev V. A., Lipatnikov V. I., Alekhin A. V. The Designing of Wave Horoscope. MGTU, 1997, (in Russian). [CrossRef]
 Lam K. Y., Loy C. T. Effects of boundary conditions on frequencies of a multilayered cylindrical shell. Journal of Sound and Vibration, Vol. 188, Issues 37, 1995, p. 363384. [CrossRef]
 Zhang X. M., Liu G. R. Vibration Analysis of thin cylindrical shells using wave propagation approach. Journal of Sound and Vibration, Vol. 239, Issue 3, 2001, p. 397403. [CrossRef]
 Blevins R. D. Formulas for Natural Frequency and Mode Shape. Van Nostrand Reinhold, New York, 1979. [CrossRef]
 Li Hua, Lam ShinYong Rotating shell dynamics. Studies in Applied Mechanics, Vol. 50, 2005, p. 284. [CrossRef]
 Lopatin A. V., Morozov E. V., Shatov A. V. An analytical expression for fundamental frequency of the composite lattice cylindrical shell with clamped edges. Composite Structures, Vol. 141, 2016, p. 232239. [CrossRef]
 Jam J., Zadeh M., Taghavian H., Eftari B. Vibration analysis of gridstiffened circular cylindrical shells with full free edges. Polish Maritime Research, Vol. 18, Issue 4, 2011, p. 2327. [CrossRef]
 Biderman V. L. Mechanics of ThinWalled Structures. Mashinostroyeniye, 1977, (in Russian). [CrossRef]
 Gerasimenko A. A., Gouskov A. M. Determination of the Natural Frequencies of the Rotating Cylindrical Shell. Mashinostroyeniye, 2012, (in Russian). [CrossRef]
 Kolkunov N. V. Basis of Calculation of the Elastic Shells. 2d. Edition, High School, 1972, (in Russian). [CrossRef]
 Gerasimenko A., Guskov M., Lorong P., Duchemin J., Gouskov A. Experimental investigation of chatter dynamics in thinwalled tubular parts turning. Procedia CIRP, 2016. [CrossRef]
 Gerasimenko A., Guskov M., Duchemin J., Lorong P., Gouskov A. Variable compliancerelated aspects of chatter in turning thinwalled tubular parts. Procedia CIRP, Vol. 31, 2015, p. 5863. [CrossRef]
 Altintas Y. Manufacturing Automation. Metal Cutting Mechanics, Machine Tool Vibrations, and CNC Design. Cambridge University Press, Cambridge, 2000. [CrossRef]
 Bolotin V. V. Vibrations in the Techniques. Mashinostroenie, 1978, p. 352, (in Russian). [CrossRef]