On the robustness of stable turning processes

Size: px
Start display at page:

Download "On the robustness of stable turning processes"

Transcription

1 3 Int. J. Machining and Machinability of Materials, Vol. 4, No. 4, 8 On the robustness of stable turning processes Zoltan Dombovari* Department of Applied Mechanics, Budapest University of Technology and Economics, Budapest, Hungary dombo@mm.bme.hu *Corresponding author R. Eddie Wilson Department of Engineering Mathematics, University of Bristol, England, UK re.wilson@bristol.ac.uk Gabor Stepan Department of Applied Mechanics, Budapest University of Technology and Economics, Budapest, Hungary stepan@mm.bme.hu Abstract: Self-excited non-linear vibrations occurring in the machining processes are investigated in this paper. Our treatment applies analytical techniques to a one Degree of Freedom (DOF but strongly non-linear mechanical model of the turning process. This tool enables us to describe and analyse the highly non-linear dynamics of the appearing periodic motions. Using normal form calculations for the Delay-Differential Equation (DDE model, we prove that the low-amplitude vibrations are unstable all along the stability lobes due to the subcriticality of Hopf bifurcations. This means that self-excited vibrations of the machine tool may occur below the stability boundaries predicted by the linear theory. Consequently, stable stationary cutting may not be robust enough for external perturbations close to the linear stability limits determined during the parameter optimisation of turning processes. Robustness is characterised by the amplitude of unstable oscillations along the stability lobes for non-linear cutting force characteristics having essential inflection points against chip thickness. Keywords: orthogonal cutting; Hopf bifurcation; subcritical; bi-stable. Reference to this paper should be made as follows: Dombovari, Z., Wilson, R.E. and Stepan, G. (8 On the robustness of stable turning processes, Int. J. Machining and Machinability of Materials, Vol. 4, No. 4, pp Biographical notes: Zoltan Dombovari received his MSc from the Budapest University of Technology and Economics Collaborating with Bristol University in 6. Currently, he is a PhD student in the Department of Applied Mechanics at the Budapest University of Technology and Economics. Copyright 8 Inderscience Enterprises Ltd.

2 On the robustness of stable turning processes 3 R. Eddie Wilson was educated at the University of Oxford (Lincoln College where he received a BA in Mathematics. He also received a PhD in Applied Mathematics from the same institution, and spent two years as a Postdoctoral Research Fellow in the Department of Mathematics at the University of Surrey. He joined the Department of Engineering Mathematics at the University of Bristol as a Lecturer in and he is currently in a Reader position in this department. Gabor Stepan received his PhD in 98 from the Technical University of Budapest. He has held visiting positions at the Universities in England, Denmark, the Netherlands and France and a Fulbright Fellowship at CalTech, Pasadena. Currently, he is the Head of the Department of Applied Mechanics at the Budapest University of Technology and Economics and a Member of the Hungarian National Academy of Sciences. Among several other journals, he serves on the Editorial Board of the Phil Trans of the Royal Society. Introduction Nowadays, high speed cutting comes to the front due to high quality and economic requirements. The corresponding machine tools are expensive compared to the usual machines because of their complex structure. The so-called regenerative chatter is not just harmful for these machines, but it also has a negative effect on the machined surface quality. The regenerative effect causes a kind of self-excited high-frequency vibration originated in the chip formation process (Tlusty and Spacek, 954; Tobias, 965. The stability of the regenerative vibrations has been studied extensively during the last few decades (Altintas and Budak, 995; Stepan, 989. The study of the corresponding non-linear vibrations started in the last years only (Insperger et al., 3; Mann et al., 3; Stepan and Kalmar-Nagy, 997; Szalai et al., 4; Zatarain et al., 6. These results led to the final conclusion that unstable vibrations exist in the vicinity of the otherwise stable stationary cutting in the stable parameter regions. These unstable vibrations separate two attractors in the sense of dynamical systems. Apart from the stable stationary cutting, there must exist another attractor, that is, a large amplitude stable oscillation outside the unstable vibration. This large amplitude non-linear vibration is stable in dynamical sense, but these vibrations are of little interest from technology view-point. However, the amplitudes of the unstable vibrations are important: they define that region in the state space also called domain of attraction, or basin where the stable cutting process is not sensitive for perturbations like non-homogeneous workpiece material. The larger the amplitude of the unstable oscillation is, the more robust the cutting process is. In this paper, we prove the existence of the unstable vibrations for a large set of cutting force characteristics. Also, the amplitudes of these oscillations are calculated in the parameter space of the technological parameters like cutting speed and chip width.

3 3 Z. Dombovari, R.E. Wilson and G. Stepan Model construction In order to investigate the fundamental effects of non-linearity during cutting, we have to apply as simple model as possible. Our one Degree of Freedom (DOF orthogonal cutting model (Figure might be far from the real cutting conditions, but it will help to follow analytically the influence of the cutting parameters on the dynamical behaviour. The desired chip thickness h equals the feed per revolution and the chip width w corresponds to the depth of cut. We disregard the geometric non-linearities of the machine tool structure since their influences are usually far smaller than those of cutting force characteristics. Thus, the equation of motion of the chosen model assumes the form: qt ( + κωn qt ( + ωn qt ( = Fx ( t ( m where ω n and κ are the natural angular frequency and the damping ratio of the essential vibration mode described by the general coordinate q that refers to the tool position. We can express these parameters with the modal mass m, stiffness k and the damping factor b. k b ωn =, κ = m mω n F x (t is the essential component of the actual cutting force. Dot refers to time derivative. Figure Orthogonal cutting with regenerative effect. Cutting force characteristics To build up the dynamical model, we need a simple expression of the cutting force characteristics. There are a number of traditional, empirical formulas for the cutting force depending on the chip width and thickness. The most popular and generally applied one is the power-law (Taylor, 97 ν F ( h = Kwh, K, ν x + K and ν are empirical parameters. The exponent ν may vary from /5 (Kalmar-Nagy and Pratt, 999 through 3/4 (Kienzle, 957 to 4/5 (Tlusty and Spacek, 954. The origin of

4 On the robustness of stable turning processes 33 these power-law expressions is in the linear optimisation techniques in the parameter space of the logarithms of the cutting parameters w and h. There is a less frequently used expression of the cutting force, which is essentially a cubic polynomial curve fitted on the experimental data (Shi and Tobias, 984 x ( 3 F ( h = w ρ h+ ρ h + ρ h, ρ ( 3,,3 (e.g. ρ = [N/m ], ρ = [N/m 3 ], ρ 3 = [N/m 4 ]; (Shi and Tobias, 984. Originally, these measurement data were obtained for full immersion milling with a face mill with even number (z = 4 of teeth. For this case, the cutting process nearly corresponds to the orthogonal cutting since the parametric excitation can be averaged in the system. From dynamical view-point, there are essential differences between the presented empirical interpretations of the cutting forces. The power-law in Figure (a has a vertical tangent at the origin, where the tool just touches the surface of the workpiece, and it cannot have inflection point. Also, the bifurcation calculations require Taylor series expansion of the power-law function at the stationary cutting. In turn, the Tobias cubic curve Figure (b has a finite gradient at the origin and an inflection at h ρ inf = (3 3 ρ3 while it has more parameters and so greater class of functions can be identified this way. Figure (a and (b show the curves of the power-law and the Tobias cubic curve cutting force characteristics. Equation of motion As shown in Figure, the so-called regenerative effect arises during the chip separation. Through the relative vibration between the tool and the workpiece, the tool leaves its motion pattern on the workpiece that excites the system after one revolution of the workpiece. The actual chip thickness can be expressed as a function of the present and the delayed motions of the tool ht ( = qt ( τ qt ( + h (4

5 34 Z. Dombovari, R.E. Wilson and G. Stepan where h is the theoretical chip thickness and τ = π/ω is the time delay which is equal to the time period of one revolution of the workpiece of angular velocity Ω. The chip thickness variation has the form χ( t = ht ( h = qt ( τ qt ( The power series of the cutting force with respect to the chip thickness variation is given by F ( h( t = F + k χ( t + k χ ( t + k χ ( t + x 3 x 3 l F + kl =, l l! h x l h Note that in the case of Tobias cubic curve, this is a finite polynomial function of third degree. Let us consider the perturbed system around the steady state solution with the new coordinate x(t defined by q(t = q + x(t. With substitution into (, we obtain q = Fx /( mωn for the equilibrium, which is the static deformation of the tool. Thus, we can obtain the perturbed equation of motion: where 3 xt ( + κωnxt ( + ωnxt ( = kχ ( t + kχ ( t + k3χ ( t m χ( t = x( t τ x( t ( is the chip thickness variation. After the introduction of the dimensionless time we obtain a more general form ( ( + κ ( + ( + ( = ( τ + η ( ( τ ( x t x t w x t w x t x t x t 3 ( ( τ ( + w η x t x t 3 t = ω t n where prime is the derivative with respect to t and k w =, k = w + h + 3 h mω ( ρ ρ ρ 3 n and the coefficients of the square and cubic terms in the non-linear excitation are ρ + 3ρ h ρ η =, η = ρ + ρh + 3ρ3h ρ + ρh + 3ρ3h (5 By abuse of notation, we drop the tilde, thus, the equation of motion is ( 3 x ( t + κx ( t + ( + wxt ( = w xt ( τ + η ( xt ( τ xt ( + η ( xt ( τ xt ( (6 3

6 On the robustness of stable turning processes 35 Repeat this equation of motion for the delayed time t τ to obtain ( ( xt ( xt ( ( xt ( xt ( 3 x ( t τ + κx ( t τ + ( + w x( t τ = w x( t τ + η τ τ + η τ τ 3 (7 and subtract (6 from (7. Then the non-linear equation of motion for the chip thickness variation will be simplified to ( ( ( χ ( t + κχ ( t + ( + w χ( t = w χ( t τ + η χ ( t τ χ ( t + η χ ( t τ χ ( t where much less number of non-linear terms will show up in the algebraic calculations. In first order Delay-Differential Equations (DDE form, we have y ( t = Ly( t + Ry( t τ + g( y(, t y ( t τ (8 where y ( t = col( y( t, y( t = col( χ( t, χ ( t and the linear non-delayed and delayed coefficient matrices and the non-linear term are given by L =,, ( w κ R = w + gy ( ( t, y( t τ wη y t y t ( ( τ ( wη3( y( t τ y( t =.3 Operator formulation In order to investigate the non-linear DDE, we transform the system into the space of continuously differentiable functions by the shift, u ( θ = u( t + θ:[ τ,], θ [ τ,] t In this space, the whole system can be rewritten as an Operational Differential Equation (OpDE ( y = Ay + F y (9 t t t with the linear operator A defined by o yt ( θ, if θ [ τ, A t ( θ = ( Lyt( + Ryt( τ, if θ = ( y and the non-linear operator F defined by F, if θ [ τ,] ( θ = g( yt(, yt( τ, if θ = ( ( ( yt

7 36 Z. Dombovari, R.E. Wilson and G. Stepan It can be shown that (9 corresponds to (8 since the two kinds of derivates are equivalent o d d ut = ut( θ = u ( t + θ, θ [ τ,] dθ dt 3 Linear stability Consider only the linear part y t = A yt of (9 with the exponential trial solution: y ( θ = b( θ, θ [ τ,], λ ( t e λt The trivial solution of (9 is exponentially stable if the characteristic exponents λ k, that is, the eigenvalues of the linear operator A satisfy Re λ k <, k + With the substitution of the general solution into the linear equation, we get a boundary value problem b o ( θ = λb ( θ ( L λi b( + Rb ( τ = With the trial solution b(θ = B exp(λθ, this leads to an eigenvalue problem leading to the characteristic function λτ ( L λi Re D( λ = det + λ κλ w = we λτ The damping factor is fixed for a certain machine tool, so the stability chart will be constructed for the variable parameters w and τ, which are determined by the chip width and the cutting speed. The stability boundaries for these parameters are expressed as a function of the dimensionless vibration frequency ω after the substitution λ = iω :π w( ω = ( ω + 4 ( ω κ ω ω τω ( = jπ+ arctan, j ω κω Then the dimensionless angular velocity Ω of the workpiece is calculated as (3 (4 Ω ( ω = π τω ( (5 The jth parametric function formed from (4 to (5 is called jth lobe. The stability limit consists of the sections of the lobes in the (w, Ω plane beneath, which the cutting process is asymptotically stable (Stepan, 989; Tobias, 965 (see Figure 3.

8 On the robustness of stable turning processes 37 Figure 3 The linear stability chart (a and the vibration frequency at the stability limits (b; κ =. 4 Non-linear investigation At the stability limits, there are two complex conjugate critical characteristic roots. We will choose the dimensionless chip width w as a bifurcation parameter and follow the motion of these critical roots in the complex plane as w increases through the stability limit. 4. Overview of the Hopf bifurcation calculus We briefly present here the procedure of the Hopf bifurcation calculation following the algorithm given in Dombovari (6 and Stepan and Kalmar-Nagy ( Variation of critical eigenvalues According to the Hopf bifurcation theory, we need the derivates of the critical eigenvalues with respect to the bifurcation parameter w, which comes from the implicit differentiation of characteristic function (3: where d λ ( ω γ n( ω γω ( = Re = d w λ= iω w( ω γ d ( ω ( ( γ ω = κ τ ω ω + κ ω + + τ ω ω (6 ( 4 ( 4 ( n ( 8 κτ ( ω ( ω ( 3ω 8κ ( ω ( τ ( ω ω ( ω 6 ω τ ( ω ( ω γ ( ω = 6 κτ ( ωω + 3 κτωω ( ω d ( + + (7

9 38 Z. Dombovari, R.E. Wilson and G. Stepan 4.. Centre manifold At the linear stability limit, the Hopf bifurcation can be studied on a two-dimensional centre manifold embedded in the infinite dimensional phase space. The tangent subspace of the centre manifold at the origin is spanned by the real and imaginary part of the critical eigenvectors of the linear operator A (Figure 4. Figure 4 Tangent subspace and centre manifold at the steady state The critical eigenvectors are calculated from ( A s( θ = i ωs( θ, s( θ = s ( θ + is ( θ R I Substituting the operator A according to (, and solving the corresponding boundary value problem, we obtain the real and the imaginary part of the eigenvectors s cosωθ sinωθ R( θ = and I( θ = ωsinωθ ωcosωθ s (8 Since the base of the tangent space defined by (8 is not orthogonal, to make the projection, we need its reciprocal base spanned by n R and n I, which satisfy the adjoint problem ( n = i n n = nr + ini A ( ϑ ω ( ϑ, ( ϑ ( ϑ ( ϑ (9 where ϑ [, τ] and the adjoint operator A is defined as d A u ( ϑ if ϑ (, τ d ϑ ( u ( ϑ = L u( + R u( τ if ϑ = The orthonormality of the critical eigenvectors are prescribed by the conditions ( n s ( n s, =,, = R R R I where the scalar product is defined by * * ( uv, = u( θ + τ Rv( θd θ + u( v ( τ

10 On the robustness of stable turning processes 39 The solution of the linear boundary value problem (9 provides the critical normed adjoint eigenvectors nr( ϑ = b( ω cos ωϑ b( ω sin ωϑ ϑ [, τ ] n ( ϑ = b ( ωsin ωϑ+ b ( ωcosωϑ I where the coefficient vectors b (ω and b (ω are not presented here. With the help of the new coordinates z ( t = col( z(, t z( t we can decompose the phase space in the following way: y t ( θ = z( t sr ( θ + z ( t s I ( θ + vt ( θ = ( T z( θ + vt ( θ, where the transformation operator is [ ] ( T z( θ = s ( θ s ( θ z( t R and the corresponding projector operator will be T p consequently z (, (, T u = n u n u t R t I t ( t = T p y t I With the new coordinates, we can express the transformed OpDE as ( t ( p z ( t Jb z( t T F z(, t v = + t( θ o t( θ v A v G zv, t ( where ( p ( zv, = ( zv, ( zv, G F T T F t t t and the Jordan block at the critical parameters assumes the form J b ω = ω The centre manifold is bent, so its second degree approximation with the new coordinates z(t is vt ( θ v( θ, z( t = h ( θ z ( t + h ( θ z ( t z ( t + h ( θ z ( t ( 3 where the unknown coefficient functions are calculated from a linear boundary value problem, again. During the lengthy calculations, the restriction of the non-linear operator F ( on the centre manifold requires the following substitution of the non-linear function

11 33 Z. Dombovari, R.E. Wilson and G. Stepan gy ( ( t, y( t τ gz ( ( t, v (, z( t η + η ( a ( z z + ( a z v ( τ, z z v (, z 3 3 3( a ( z z = w( ω a z w ω ω w ω z κωz and η, η 3 come from (5. ( = ( (( + ( + 4. Poincaré-Ljapunov constant The first two scalar equations of ( describe the flow on the centre manifold, where the Hopf bifurcation takes place. Its second and third order terms assume the form ( (, (( (, ( p p T F z t vt T F z t v z t j k a,,,3 jkz( t z( t jk j+ k= = j k b,,,3 jkz( t z( t jk j+ k= The Poincaré-Ljapunov Constant (PLC can be expressed directly with the help of the coefficients from ( (Stepan 989 ( ω = + + 8ω ω + ( 3 a3 + a + b + 3 b3 8 ( a a ( a b b ( b b ( a a b If ( ω >, then the Hopf bifurcation is subcritical, otherwise, it is supercritical. Since the PLC can be calculated along the lobes as a function of the vibration frequency ω, we can follow the criticality of the Hopf bifurcation along the lobes. If we substitute the coefficients a jk, b jk, ( we obtain δn ( ω η + δn( ω η3 ( ω = ( ω (3 ( ωγ ( ω u d d ( ( where ( ( 4 4 ( ω 4 ( ω ( ( ( 4 ( δ ω κω κ τ ω ω κ ω ω n( = 6 48 ( κτωω ( + 48 κτωω ( κω ω + 7ω 3κω ω 4ω (4

12 On the robustness of stable turning processes 33 n d n ( ( ( δ ( ω = 3 u ( ω γ ( ω ω u ( ω = 36κ ω + ω 4ω d (5 and γ n (ω, γ d (ω can be found in (6 and (7. We can determine the first harmonic component of the arising periodic motion analytically. Its amplitude can be expressed as a function of the bifurcation parameter w: γω ( r( ω, w = ( w w( ω ( ω δ ( ω w = 3 δ ω η δ ω η w( ω n n( + n( 3 (6 4.3 The criticality of the Hopf bifurcation It has great technical relevance to what kind of limit cycle exists in the vicinity of the stability boundaries, that is, whether the arisen Hopf bifurcation is subcritical or supercritical. In the subcritical case, an unstable oscillation exists around the stable stationary cutting, which can still lead to chatter for perturbations larger than the amplitude of the unstable oscillation. Consequently the sign of (3 refers to the nature of the stability. Since the chatter frequencies ω (, along the stability boundaries, w(ω, γ d (ω and u d (ω are positive in (3 in accordance with (4, (7 and (5. Thus, the sense of the bifurcation depends only on the nominator of the PLC, more exactly on the signs of δ n (ω defined in (4 and η 3 defined in (5: δ ( ω η + δ ( ω η > n n 3 This way, we get the condition δ ( ω η > η (7 n 3 δn( ω Substituting (5 into the condition (7, the problem leads to a second order polynomial condition with respect to the theoretical chip thickness h : ( ( 3 δ ω n ρ ρ n h h δ ω ρ > δn( ω 3 ρ3 3ρ3 3 δn( ω ρ3 This condition is always satisfied, if it has no real root for h, that is, if the discriminant D satisfies ρ δ ( ω D = ρ + < n 3 ρ 3 3 ρ3 δn( ω Since ρ 3 > for all reasonable cutting force characteristics ( and the last factor is always positive for all ω (,, κ [,, we are left with the simple condition 3ρ ρ > ρ 3

13 33 Z. Dombovari, R.E. Wilson and G. Stepan This condition is equivalent with the condition of the positive gradient at the inflection point (3 since F x ρ ( h w h inf = ρ > 3ρρ3 > ρ 3 ρ3 Thus, if the gradient of the cubic cutting force function at the inflection point is positive, then the Hopf bifurcation is subcritical. Otherwise, the subcriticality is not proven mathematically, but again, having negative derivative at a possible inflection point on the cutting force characteristic is physically unreasonable. 5 Robustness of stable cutting The unstable period one branch determined by the Hopf calculation gives us the good approximation (6 of the amplitude r of the chip thickness variation χ as a function of the chip width values w close enough to the critical value w(ω from below. Figure 5 clearly shows that the otherwise stable cutting process is much more robust at the left side of the lobes than it is at their right side. The parabolas are flat; they have large curvature at the right side of the lobes, so the stationary cutting is more sensitive to the perturbations there. It is safer to tune the parameters to the left side of the lobes, these regions are smaller, though, with respect to the cutting speeds. This conclusion is similar to the one obtained for the power law by Wahi and Chatterjee (5. Figure 5 Structure of periodic orbits; j =, κ =. The subcriticality of the Hopf bifurcation means that the unstable limit cycle separates two coexisting stable motions, the desired stable stationary cutting and an undesired large amplitude chatter that is also stable in the sense of dynamical systems theory. This is why this region is called the region of bi-stability. Since we proved that the Hopf bifurcation is subcritical along the lobes, the bi-stable regions are located below these lobes, covering a large region of the parameter domain of stable stationary cutting.

14 On the robustness of stable turning processes 333 According to (6, the first harmonic component of the limit cycle can be calculated as a projection of the real orbit to the subspace tangent to the centre manifold (Figure 4: z ( t r( ω, wcosωt z( t r( ω, wsinωt t ( θ v This way, the actual chip thickness is given by (4: ht ( = χ( t + h = yt,( + h = z( t + h = r( ω, wcos( ωt + h The tool leaves the surface of the workpiece when the chip thickness becomes zero: ht ( = r( ω, w = h (8 Expressing the bifurcation parameter w from (8, we get the so-called bi-stable limit where the periodic orbits satisfy the switching condition: 3 n( n( 3 wbs ( w( h δ ω η + δ ω η ω = ω δn ( ω Figure 6 shows the bi-stable regions as dark strips in the grey stable cutting domain, just below the white unstable parameter domain. With the appropriate choice of the parameters η, η 3 from (6, the results are shown for both 3/4 law and for the Tobias cubic curve. While the dimensionless linear stability boundary is not affected at all by the applied cutting force characteristic function, the size of the bi-stable region strongly depends on it. Clearly, this bi-stable region is much larger if the cutting force characteristics follow the cubic curve rather than the power-law. Figure 6 Both panels show the stability chart with the bi-stable region in case of the 3/4 powerlaw (a and Tobias cubic curve (b; h =.6 [mm], κ =. 6 Conclusion During the design of the technological parameters, it is much more convenient and common to use the power-law approximation of the real measured cutting force characteristics. From the view-point of the stability prediction of the stationary cutting, the linearisation of this power-law is still appropriate. However, from non-linear

15 334 Z. Dombovari, R.E. Wilson and G. Stepan vibrations view-point, the use of the power-law approximation has great deficiencies. As we showed, the prediction of the bi-stable parameter domain is important if we want to secure stable stationary cutting robust enough for perturbations. However, this domain cannot be predicted correctly, in a conservative way, if we use the power-law approximation. We proved that the Tobias cubic curve approximation of the cutting force characteristics leads to subcritical Hopf bifurcations along all the stability limits in the same way as it was shown for the power-law by Kalmar-Nagy and Pratt (999 and Wahi and Chatterjee (5. Namely, the stable cutting process is more robust against perturbations originated in non-homogeneous workpiece material at the left side of the lobes than it is at their right side. In the meantime we also showed that the bi-stable domain can be substantially larger for more realistic cutting force characteristics involving inflection points than the ones predicted by the power-law. Acknowledgement This research was supported by the Hungarian Nation Sciences Foundation under grant no. OTKA F4738 and OTKA T References Altintas, Y. and Budak, E. (995 Analytical prediction of stability lobes in milling, CIRP, , pp Dombovari, Z. (6 Bifurcation analysis of a cutting process, MSc Thesis, BME, Hungary. Insperger, T., Stepan, G., Bayly, P.V. and Mann, B.P. (3 Multiple chatter frequencies in milling processes, Journal of Sound and Vibration, Vol. 6, pp Kalmar-Nagy, T. and Pratt, R.J. (999 Experimental and analytical investigation of the subcritical instability in metal cutting, ASME, Nevada, USA. Kienzle, O. (957 Spezifische Schnittkräfte bei der Metallbearbeitung (in German, Werkstattstechnik und Maschinenbau, Vol. 47. Mann, B.P., Bayly, P.V., Davies, M.A. and Halley, J.E. (3 Limit cycles, bifurcations, and accuracy of the milling process, Journal of Sound and Vibration, Vol. 77, pp Shi, H.M. and Tobias, S.A. (984 Theory of finite amplitude machine tool instability, International Journal of Machine Tool Design and Research, Vol. 4, pp Stepan, G. (989 Retarded Dynamical Systems, London: Longman. Stepan, G. and Kalmar-Nagy, T. (997 Nonlinear regenerative machine tool vibrations, ASME, California, USA. Szalai, R., Stepan, G. and Hogan, S.J. (4 Global dynamics of low immersion high-speed milling, Chaos, Vol. 4, pp Taylor, F.W. (97 On the art of cutting metals, Transaction of ASME, Vol. 8, pp Tlusty, J. and Spacek, L. (954 Self-excited vibrations on machine tools (in Czech, Nakl CSAV, Prague. Tobias, S.A. (965 Machine Tool Vibrations, London: Blackie. Wahi, P. and Chatterjee, A. (5 Regenerative tool chatter near a Codimension hopf point using multiple scales, Nonlinear Dynamics, Vol. 4, pp Zatarain, M., Muñoa, J., Peigné, G. and Insperger, T. (6 Analysis of the influence of mill helix angle on chatter stability, Annals of the CIRP, Vol. 55, No..

STATE-DEPENDENT, NON-SMOOTH MODEL OF CHATTER VIBRATIONS IN TURNING

STATE-DEPENDENT, NON-SMOOTH MODEL OF CHATTER VIBRATIONS IN TURNING Proceedings of the ASME 215 International Design Engineering Technical Conferences & Computers and Information in Engineering Conference IDETC/CIE 215 August 2-5, 215, Boston, Massachusetts, USA DETC215-46748

More information

SURFACE PROPERTIES OF THE MACHINED WORKPIECE FOR HELICAL MILLS

SURFACE PROPERTIES OF THE MACHINED WORKPIECE FOR HELICAL MILLS Machining Science and Technology, 13:227 245 Copyright 29 Taylor & Francis Group, LLC ISSN: 191-344 print/1532-2483 online DOI: 1.18/191349312167 SURFACE PROPERTIES OF THE MACHINED WORKPIECE FOR HELICAL

More information

On the analysis of the double Hopf bifurcation in machining processes via center manifold reduction

On the analysis of the double Hopf bifurcation in machining processes via center manifold reduction rspa.royalsocietypublishing.org Research Article submitted to journal On the analysis of the double Hopf bifurcation in machining processes via center manifold reduction T. G. Molnar, Z. Dombovari, T.

More information

Available online at ScienceDirect. IFAC-PapersOnLine (2015)

Available online at   ScienceDirect. IFAC-PapersOnLine (2015) Available online at www.sciencedirect.com ScienceDirect IFAC-PapersOnLine 48-12 (2015) 322 327 Regenerative delay, parametric forcing and machine tool chatter: A review Tamas Insperger David Lehotzky Gabor

More information

AN EXAMINATION OF SURFACE LOCATION ERROR AND SURFACE ROUGHNESS FOR PERIOD-2 INSTABILITY IN MILLING

AN EXAMINATION OF SURFACE LOCATION ERROR AND SURFACE ROUGHNESS FOR PERIOD-2 INSTABILITY IN MILLING AN EXAMINATION OF SURFACE LOCATION ERROR AND SURFACE ROUGHNESS FOR PERIOD-2 INSTABILITY IN MILLING Andrew Honeycutt and Tony L. Schmitz Mechanical Engineering and Engineering Science University of North

More information

PERIOD-N BIFURCATIONS IN MILLING: NUMERICAL AND EXPERIMENTAL VERIFICATION

PERIOD-N BIFURCATIONS IN MILLING: NUMERICAL AND EXPERIMENTAL VERIFICATION PERIOD-N BIFURCATIONS IN MILLING: NUMERICAL AND EXPERIMENTAL VERIFICATION Andrew Honeycutt and Tony L. Schmitz Department of Mechanical Engineering and Engineering Science University of North Carolina

More information

Evaluation of active structural vibration control strategies in milling process

Evaluation of active structural vibration control strategies in milling process Evaluation of active structural vibration control strategies in milling process Monnin, J. (a); Wegener, K. (a) a) Institute of Machine Tools and Manufacturing, Zurich, Switzerland Keywords: Mechatronics,

More information

Analytical estimations of limit cycle amplitude for delay-differential equations

Analytical estimations of limit cycle amplitude for delay-differential equations Electronic Journal of Qualitative Theory of Differential Equations 2016, No. 77, 1 10; doi: 10.14232/ejqtde.2016.1.77 http://www.math.u-szeged.hu/ejqtde/ Analytical estimations of limit cycle amplitude

More information

Subcritical Hopf Bifurcation in the Delay Equation Model for Machine Tool Vibrations

Subcritical Hopf Bifurcation in the Delay Equation Model for Machine Tool Vibrations Nonlinear Dynamics 26: 2 42, 2. 2 Kluwer Academic Publishers. Printed in the Netherlands. Subcritical Hopf Bifurcation in the Delay Equation Model for Machine Tool Vibrations TAMÁS KALMÁR-NAGY Department

More information

Approximation of Top Lyapunov Exponent of Stochastic Delayed Turning Model Using Fokker-Planck Approach

Approximation of Top Lyapunov Exponent of Stochastic Delayed Turning Model Using Fokker-Planck Approach Approximation of Top Lyapunov Exponent of Stochastic Delayed Turning Model Using Fokker-Planck Approach Henrik T. Sykora, Walter V. Wedig, Daniel Bachrathy and Gabor Stepan Department of Applied Mechanics,

More information

Analysis of the Chatter Instability in a Nonlinear Model for Drilling

Analysis of the Chatter Instability in a Nonlinear Model for Drilling Sue Ann Campbell Department of Applied Mathematics, University of Waterloo, Waterloo, ON N2L 3G1 Canada e-mail: sacampbell@uwaterloo.ca Emily Stone Department of Mathematical Sciences, The University of

More information

ON STABILITY PREDICTION FOR LOW RADIAL IMMERSION MILLING

ON STABILITY PREDICTION FOR LOW RADIAL IMMERSION MILLING Machine Science and Technology, 9:117 130 Copyright # 2005 Taylor & Francis Inc. ISSN: 1091-0344 print/1532-2483 online DOI: 10.1081/MST-200051378 ON STABILITY PREDICTION FOR LOW RADIAL IMMERSION MILLING

More information

Stability analysis of milling with irregular pitch tools by the implicit subspace iteration method

Stability analysis of milling with irregular pitch tools by the implicit subspace iteration method Int. J. Dynam. Control (214) 2:26 34 DOI 1.17/s4435-13-52-7 Stability analysis of milling with irregular pitch tools by the implicit subspace iteration method Mikel Zatarain Zoltan Dombovari Received:

More information

DETC DYNAMICS OF DRILL BITS WITH CUTTING EDGES OF VARYING PARAMETERS

DETC DYNAMICS OF DRILL BITS WITH CUTTING EDGES OF VARYING PARAMETERS Proceedings of the ASME 2013 International Design Engineering Technical Conferences & Computers and Information in Engineering Conference IDETC/CIE 2013 August 4-August 7, 2013, Portland, OR, USA DETC2013-12594

More information

Balancing of the skateboard with reflex delay

Balancing of the skateboard with reflex delay ENOC 214, July 6-11, 214, Vienna, Austria Balancing of the skateboard with reflex delay Balazs Varszegi, Denes Takacs, Gabor Stepan and S. John Hogan Department of Applied Mechanics, Budapest University

More information

HOPF BIFURCATION CALCULATIONS IN DELAYED SYSTEMS

HOPF BIFURCATION CALCULATIONS IN DELAYED SYSTEMS PERIODICA POLYTECHNICA SER. MECH. ENG. VOL. 48, NO. 2, PP. 89 200 2004 HOPF BIFURCATION CALCULATIONS IN DELAYED SYSTEMS Gábor OROSZ Bristol Centre for Applied Nonlinear Mathematics Department of Engineering

More information

Collocation approximation of the monodromy operator of periodic, linear DDEs

Collocation approximation of the monodromy operator of periodic, linear DDEs p. Collocation approximation of the monodromy operator of periodic, linear DDEs Ed Bueler 1, Victoria Averina 2, and Eric Butcher 3 13 July, 2004 SIAM Annual Meeting 2004, Portland 1=Dept. of Math. Sci.,

More information

UPDATED TOOL PATH MODELLING WITH PERIODIC DELAY FOR CHATTER PREDICTION IN MILLING.

UPDATED TOOL PATH MODELLING WITH PERIODIC DELAY FOR CHATTER PREDICTION IN MILLING. ENOC-2005, Eindhoven, Netherlands, 7-12 August 2005 UPDATED TOOL PATH MODELLING WITH PERIODIC DELAY FOR CHATTER PREDICTION IN MILLING. Ronald Faassen Department of Mechanical Engineering Eindhoven University

More information

Application of Spindle Speed Variation for Chatter Suppression in Turning

Application of Spindle Speed Variation for Chatter Suppression in Turning Application of Spindle Speed Variation for Chatter Suppression in Turning Andreas Otto a,, Günter Radons a a Institute of Physics, Chemnitz University of Technology, 917 Chemnitz, Germany Abstract In the

More information

This is an author-deposited version published in: Eprints ID: 13508

This is an author-deposited version published in:   Eprints ID: 13508 Open Archive Toulouse Archive Ouverte (OATAO) OATAO is an open access repository that collects the work of Toulouse researchers and makes it freely available over the web where possible. This is an author-deposited

More information

2:1 Resonance in the delayed nonlinear Mathieu equation

2:1 Resonance in the delayed nonlinear Mathieu equation Nonlinear Dyn 007) 50:31 35 DOI 10.1007/s11071-006-916-5 ORIGINAL ARTICLE :1 Resonance in the delayed nonlinear Mathieu equation Tina M. Morrison Richard H. Rand Received: 9 March 006 / Accepted: 19 September

More information

Phase Synchronization

Phase Synchronization Phase Synchronization Lecture by: Zhibin Guo Notes by: Xiang Fan May 10, 2016 1 Introduction For any mode or fluctuation, we always have where S(x, t) is phase. If a mode amplitude satisfies ϕ k = ϕ k

More information

Active in-process chatter control E.J.J. Doppenberg*, R.P.H. Faassen**, N. van de Wouw**, J.A.J. Oosterling*, H. Nijmeijer**,

Active in-process chatter control E.J.J. Doppenberg*, R.P.H. Faassen**, N. van de Wouw**, J.A.J. Oosterling*, H. Nijmeijer**, Active in-process chatter control E.J.J. Doppenberg*, R.P.H. Faassen**, N. van de Wouw**, J.A.J. Oosterling*, H. Nijmeijer**, *TNO Science and Industry, P.O. Box 155, 2600 AD Delft, The Netherlands 1 **Eindhoven

More information

Difference Resonances in a controlled van der Pol-Duffing oscillator involving time. delay

Difference Resonances in a controlled van der Pol-Duffing oscillator involving time. delay Difference Resonances in a controlled van der Pol-Duffing oscillator involving time delay This paper was published in the journal Chaos, Solitions & Fractals, vol.4, no., pp.975-98, Oct 9 J.C. Ji, N. Zhang,

More information

Experimental investigation of micro-chaos

Experimental investigation of micro-chaos Experimental investigation of micro-chaos Gergely Gyebrószki Gábor Csernák Csaba Budai Department of Applied Mechanics Budapest University of Technology and Economics Budapest Hungary MTA-BME Research

More information

DETC DELAY DIFFERENTIAL EQUATIONS IN THE DYNAMICS OF GENE COPYING

DETC DELAY DIFFERENTIAL EQUATIONS IN THE DYNAMICS OF GENE COPYING Proceedings of the ASME 2007 International Design Engineering Technical Conferences & Computers and Information in Engineering Conference IDETC/CIE 2007 September 4-7, 2007, Las Vegas, Nevada, USA DETC2007-3424

More information

Additive resonances of a controlled van der Pol-Duffing oscillator

Additive resonances of a controlled van der Pol-Duffing oscillator Additive resonances of a controlled van der Pol-Duffing oscillator This paper has been published in Journal of Sound and Vibration vol. 5 issue - 8 pp.-. J.C. Ji N. Zhang Faculty of Engineering University

More information

Dynamical Systems with Varying Time Delay in Engineering Applications

Dynamical Systems with Varying Time Delay in Engineering Applications Budapest University of Technology and Economics Department of Applied Mechanics Dynamical Systems with Varying Time Delay in Engineering Applications PhD dissertation Tamás G. Molnár Supervisor: Tamás

More information

The Higgins-Selkov oscillator

The Higgins-Selkov oscillator The Higgins-Selkov oscillator May 14, 2014 Here I analyse the long-time behaviour of the Higgins-Selkov oscillator. The system is ẋ = k 0 k 1 xy 2, (1 ẏ = k 1 xy 2 k 2 y. (2 The unknowns x and y, being

More information

Three ways of treating a linear delay differential equation

Three ways of treating a linear delay differential equation Proceedings of the 5th International Conference on Nonlinear Dynamics ND-KhPI2016 September 27-30, 2016, Kharkov, Ukraine Three ways of treating a linear delay differential equation Si Mohamed Sah 1 *,

More information

Research Article The Microphone Feedback Analogy for Chatter in Machining

Research Article The Microphone Feedback Analogy for Chatter in Machining Shock and Vibration Volume 215, Article ID 976819, 5 pages http://dx.doi.org/1.1155/215/976819 Research Article The Microphone Feedback Analogy for Chatter in Machining Tony Schmitz UniversityofNorthCarolinaatCharlotte,Charlotte,NC28223,USA

More information

Local and Global Dynamics of Cutting Processes

Local and Global Dynamics of Cutting Processes Local and Global Dynamics of Cutting Processes Zoltán Dombóvári Supervisor: Prof. Gábor Stépán Department of Applied Mechanics Budapest University of Technology and Economics A thesis submitted for the

More information

Chapter 14 Three Ways of Treating a Linear Delay Differential Equation

Chapter 14 Three Ways of Treating a Linear Delay Differential Equation Chapter 14 Three Ways of Treating a Linear Delay Differential Equation Si Mohamed Sah and Richard H. Rand Abstract This work concerns the occurrence of Hopf bifurcations in delay differential equations

More information

Math Ordinary Differential Equations

Math Ordinary Differential Equations Math 411 - Ordinary Differential Equations Review Notes - 1 1 - Basic Theory A first order ordinary differential equation has the form x = f(t, x) (11) Here x = dx/dt Given an initial data x(t 0 ) = x

More information

Analytic solution of simplified Cardan s shaft model

Analytic solution of simplified Cardan s shaft model Applied and Computational Mechanics 8 (2014) 215 228 Analytic solution of simplified Cardan s shaft model M. Zají cek a,, J. Dupal a a Faculty of Applied Sciences, University of West Bohemia, Univerzitní

More information

Active control of time-delay in cutting vibration

Active control of time-delay in cutting vibration THEORETICAL & APPLIED MECHANICS LETTERS 3, 633 (213) Active control of time-delay in cutting vibration Pingxu Zheng a) and Xinhua Long b) State Key Laboratory of Mechanical System and Vibration, Shanghai

More information

DYNAMICS ASPECT OF CHATTER SUPPRESSION IN MILLING

DYNAMICS ASPECT OF CHATTER SUPPRESSION IN MILLING 11th World Congress on Computational Mechanics (WCCM XI) 5th European Conference on Computational Mechanics (ECCM V) 6th European Conference on Computational Fluid Dynamics (ECFD VI) E. Oñate, J. Oliver

More information

Oscillation Control in Delayed Feedback Systems

Oscillation Control in Delayed Feedback Systems Oscillation Control in Delayed Feedback Systems Fatihcan M. Atay (atay@member.ams.org) Preprint. Final version in Lecture Notes in Control and Information Sciences Vol. 273, pp. 3 6 (22) Abstract. The

More information

A Robust Semi-Analytical Method for Calculating the Response Sensitivity of a Time Delay System

A Robust Semi-Analytical Method for Calculating the Response Sensitivity of a Time Delay System A Robust Semi-Analytical Method for Calculating the Response Sensitivity of a Time Delay System Mohammad H. Kurdi 1 Postdoctoral Associate e-mail: mhkurdi@gmail.com Raphael T. Haftka Distinguished Professor

More information

Dr.Vinod Hosur, Professor, Civil Engg.Dept., Gogte Institute of Technology, Belgaum

Dr.Vinod Hosur, Professor, Civil Engg.Dept., Gogte Institute of Technology, Belgaum STRUCTURAL DYNAMICS Dr.Vinod Hosur, Professor, Civil Engg.Dept., Gogte Institute of Technology, Belgaum Overview of Structural Dynamics Structure Members, joints, strength, stiffness, ductility Structure

More information

MODELING OF REGENERATIVE CHATTER OF A MILLING PROCESS TO DELINEATE STABLE CUTTING REGION FROM UNSTABLE REGION

MODELING OF REGENERATIVE CHATTER OF A MILLING PROCESS TO DELINEATE STABLE CUTTING REGION FROM UNSTABLE REGION International Journal of Mechanical Engineering and Technology (IJMET) Volume 9, Issue, November 208, pp. 748 757, Article ID: IJMET_09 076 Available online at http://www.iaeme.com/ijmet/issues.asp?jtype=ijmet&vtype=9&itype=

More information

STABILITY. Phase portraits and local stability

STABILITY. Phase portraits and local stability MAS271 Methods for differential equations Dr. R. Jain STABILITY Phase portraits and local stability We are interested in system of ordinary differential equations of the form ẋ = f(x, y), ẏ = g(x, y),

More information

Stability and Hopf Bifurcation for a Discrete Disease Spreading Model in Complex Networks

Stability and Hopf Bifurcation for a Discrete Disease Spreading Model in Complex Networks International Journal of Difference Equations ISSN 0973-5321, Volume 4, Number 1, pp. 155 163 (2009) http://campus.mst.edu/ijde Stability and Hopf Bifurcation for a Discrete Disease Spreading Model in

More information

Evaluation of Cutting Forces and Prediction of Chatter Vibrations in Milling

Evaluation of Cutting Forces and Prediction of Chatter Vibrations in Milling H. B. Lacerda and V. T. Lima H. B. Lacerda and V. T. Lima Federal University of Uberlândia School of Mechanical Engineering Av. João Naves de Ávila, Bl. M 38400-90 Uberlândia, MG. Brazil helder@mecanica.ufu.br

More information

This is an author-deposited version published in: Eprints ID: 6397

This is an author-deposited version published in:   Eprints ID: 6397 Open Archive Toulouse Archive Ouverte (OATAO) OATAO is an open access repository that collects the work of Toulouse researchers and makes it freely available over the web where possible. This is an author-deposited

More information

B5.6 Nonlinear Systems

B5.6 Nonlinear Systems B5.6 Nonlinear Systems 4. Bifurcations Alain Goriely 2018 Mathematical Institute, University of Oxford Table of contents 1. Local bifurcations for vector fields 1.1 The problem 1.2 The extended centre

More information

Influence of Torsional Motion on the Axial Vibrations of a Drilling Tool

Influence of Torsional Motion on the Axial Vibrations of a Drilling Tool Sergey A. Voronov e-mail: voronov@rk5.bmstu.ru Alexander M. Gouskov e-mail: gouskov@rk5.bmstu.ru Alexey S. Kvashnin e-mail: a-kvashnin@yandex.ru Bauman Moscow State Technical University, 5, 2-nd Baumanskaya,

More information

Strange dynamics of bilinear oscillator close to grazing

Strange dynamics of bilinear oscillator close to grazing Strange dynamics of bilinear oscillator close to grazing Ekaterina Pavlovskaia, James Ing, Soumitro Banerjee and Marian Wiercigroch Centre for Applied Dynamics Research, School of Engineering, King s College,

More information

2 Discrete growth models, logistic map (Murray, Chapter 2)

2 Discrete growth models, logistic map (Murray, Chapter 2) 2 Discrete growth models, logistic map (Murray, Chapter 2) As argued in Lecture 1 the population of non-overlapping generations can be modelled as a discrete dynamical system. This is an example of an

More information

Nonlinear Stability and Bifurcation of Multi-D.O.F. Chatter System in Grinding Process

Nonlinear Stability and Bifurcation of Multi-D.O.F. Chatter System in Grinding Process Key Engineering Materials Vols. -5 (6) pp. -5 online at http://www.scientific.net (6) Trans Tech Publications Switzerland Online available since 6//5 Nonlinear Stability and Bifurcation of Multi-D.O.F.

More information

Isolated large amplitude periodic motions of towed rigid wheels

Isolated large amplitude periodic motions of towed rigid wheels Nonlinear Dyn 2008) 52: 27 34 DOI 10.1007/s11071-007-9253-y ORIGINAL PAPER Isolated large amplitude periodic motions of towed rigid wheels Dénes Takács Gábor Stépán S. John Hogan Received: 8 August 2006

More information

B5.6 Nonlinear Systems

B5.6 Nonlinear Systems B5.6 Nonlinear Systems 5. Global Bifurcations, Homoclinic chaos, Melnikov s method Alain Goriely 2018 Mathematical Institute, University of Oxford Table of contents 1. Motivation 1.1 The problem 1.2 A

More information

Multi Degrees of Freedom Systems

Multi Degrees of Freedom Systems Multi Degrees of Freedom Systems MDOF s http://intranet.dica.polimi.it/people/boffi-giacomo Dipartimento di Ingegneria Civile Ambientale e Territoriale Politecnico di Milano March 9, 07 Outline, a System

More information

CALCULATION OF NONLINEAR VIBRATIONS OF PIECEWISE-LINEAR SYSTEMS USING THE SHOOTING METHOD

CALCULATION OF NONLINEAR VIBRATIONS OF PIECEWISE-LINEAR SYSTEMS USING THE SHOOTING METHOD Vietnam Journal of Mechanics, VAST, Vol. 34, No. 3 (2012), pp. 157 167 CALCULATION OF NONLINEAR VIBRATIONS OF PIECEWISE-LINEAR SYSTEMS USING THE SHOOTING METHOD Nguyen Van Khang, Hoang Manh Cuong, Nguyen

More information

Identification Methods for Structural Systems. Prof. Dr. Eleni Chatzi Lecture March, 2016

Identification Methods for Structural Systems. Prof. Dr. Eleni Chatzi Lecture March, 2016 Prof. Dr. Eleni Chatzi Lecture 4-09. March, 2016 Fundamentals Overview Multiple DOF Systems State-space Formulation Eigenvalue Analysis The Mode Superposition Method The effect of Damping on Structural

More information

University of Bristol - Explore Bristol Research. Early version, also known as pre-print

University of Bristol - Explore Bristol Research. Early version, also known as pre-print Orosz, G., & Stepan, G. (4). Hopf bifurcation calculations in delayed systems with translational symmetry. https://doi.org/1.17/s33-4- 65-4 Early version, also known as pre-print Link to published version

More information

Stability and bifurcation of a simple neural network with multiple time delays.

Stability and bifurcation of a simple neural network with multiple time delays. Fields Institute Communications Volume, 999 Stability and bifurcation of a simple neural network with multiple time delays. Sue Ann Campbell Department of Applied Mathematics University of Waterloo Waterloo

More information

Stabilizability diagram for turning processes subjected to digital PD control

Stabilizability diagram for turning processes subjected to digital PD control Int. J. Dynam. Control (2014) 2:46 54 DOI 10.1007/s40435-013-0047-4 Stabilizability diagram for turning processes subjected to digital PD control David Lehotzky Janos Turi Tamas Insperger Received: 10

More information

Research Article A Method for Stability Analysis of Periodic Delay Differential Equations with Multiple Time-Periodic Delays

Research Article A Method for Stability Analysis of Periodic Delay Differential Equations with Multiple Time-Periodic Delays Hindawi Mathematical Problems in Engineering Volume 217, Article ID 949142, 8 pages https://doiorg/11155/217/949142 Research Article A Method for Stability Analysis of Periodic Delay Differential Equations

More information

(8.51) ẋ = A(λ)x + F(x, λ), where λ lr, the matrix A(λ) and function F(x, λ) are C k -functions with k 1,

(8.51) ẋ = A(λ)x + F(x, λ), where λ lr, the matrix A(λ) and function F(x, λ) are C k -functions with k 1, 2.8.7. Poincaré-Andronov-Hopf Bifurcation. In the previous section, we have given a rather detailed method for determining the periodic orbits of a two dimensional system which is the perturbation of a

More information

One Dimensional Dynamical Systems

One Dimensional Dynamical Systems 16 CHAPTER 2 One Dimensional Dynamical Systems We begin by analyzing some dynamical systems with one-dimensional phase spaces, and in particular their bifurcations. All equations in this Chapter are scalar

More information

2 Lecture 2: Amplitude equations and Hopf bifurcations

2 Lecture 2: Amplitude equations and Hopf bifurcations Lecture : Amplitude equations and Hopf bifurcations This lecture completes the brief discussion of steady-state bifurcations by discussing vector fields that describe the dynamics near a bifurcation. From

More information

Lecture 3 : Bifurcation Analysis

Lecture 3 : Bifurcation Analysis Lecture 3 : Bifurcation Analysis D. Sumpter & S.C. Nicolis October - December 2008 D. Sumpter & S.C. Nicolis General settings 4 basic bifurcations (as long as there is only one unstable mode!) steady state

More information

Chapter 2 Machining Stability

Chapter 2 Machining Stability Chapter 2 Machining Stability Mark J. Jackson, H. Zhang, and J. Ni Abstract Since the early days of regenerative chatter theory, it has been noticed that the phase difference between the current and the

More information

Advanced Vibrations. Elements of Analytical Dynamics. By: H. Ahmadian Lecture One

Advanced Vibrations. Elements of Analytical Dynamics. By: H. Ahmadian Lecture One Advanced Vibrations Lecture One Elements of Analytical Dynamics By: H. Ahmadian ahmadian@iust.ac.ir Elements of Analytical Dynamics Newton's laws were formulated for a single particle Can be extended to

More information

AA242B: MECHANICAL VIBRATIONS

AA242B: MECHANICAL VIBRATIONS AA242B: MECHANICAL VIBRATIONS 1 / 50 AA242B: MECHANICAL VIBRATIONS Undamped Vibrations of n-dof Systems These slides are based on the recommended textbook: M. Géradin and D. Rixen, Mechanical Vibrations:

More information

Oscillatory Motion. Simple pendulum: linear Hooke s Law restoring force for small angular deviations. small angle approximation. Oscillatory solution

Oscillatory Motion. Simple pendulum: linear Hooke s Law restoring force for small angular deviations. small angle approximation. Oscillatory solution Oscillatory Motion Simple pendulum: linear Hooke s Law restoring force for small angular deviations d 2 θ dt 2 = g l θ small angle approximation θ l Oscillatory solution θ(t) =θ 0 sin(ωt + φ) F with characteristic

More information

COMPUTER AIDED NONLINEAR ANALYSIS OF MACHINE TOOL VIBRATIONS AND A DEVELOPED COMPUTER SOFTWARE

COMPUTER AIDED NONLINEAR ANALYSIS OF MACHINE TOOL VIBRATIONS AND A DEVELOPED COMPUTER SOFTWARE Mathematical and Computational Applications, Vol. 10, No. 3, pp. 377-385, 005. Association for Scientific Research COMPUTER AIDED NONLINEAR ANALYSIS OF MACHINE TOOL VIBRATIONS AND A DEVELOPED COMPUTER

More information

Oscillatory Motion. Simple pendulum: linear Hooke s Law restoring force for small angular deviations. Oscillatory solution

Oscillatory Motion. Simple pendulum: linear Hooke s Law restoring force for small angular deviations. Oscillatory solution Oscillatory Motion Simple pendulum: linear Hooke s Law restoring force for small angular deviations d 2 θ dt 2 = g l θ θ l Oscillatory solution θ(t) =θ 0 sin(ωt + φ) F with characteristic angular frequency

More information

Math 1553, Introduction to Linear Algebra

Math 1553, Introduction to Linear Algebra Learning goals articulate what students are expected to be able to do in a course that can be measured. This course has course-level learning goals that pertain to the entire course, and section-level

More information

On stability prediction for milling

On stability prediction for milling International Journal of Machine Tools & Manufacture 45 (2005) 769 781 www.elsevier.com/locate/ijmactool On stability prediction for milling Janez Gradišek a, *, Martin Kalveram b, Tamás Insperger c, Klaus

More information

Remark on Hopf Bifurcation Theorem

Remark on Hopf Bifurcation Theorem Remark on Hopf Bifurcation Theorem Krasnosel skii A.M., Rachinskii D.I. Institute for Information Transmission Problems Russian Academy of Sciences 19 Bolshoi Karetny lane, 101447 Moscow, Russia E-mails:

More information

ELEC 3035, Lecture 3: Autonomous systems Ivan Markovsky

ELEC 3035, Lecture 3: Autonomous systems Ivan Markovsky ELEC 3035, Lecture 3: Autonomous systems Ivan Markovsky Equilibrium points and linearization Eigenvalue decomposition and modal form State transition matrix and matrix exponential Stability ELEC 3035 (Part

More information

Chatter reduction through active vibration damping

Chatter reduction through active vibration damping Université Libre de Bruxelles F a c u l t y o f A p p l i e d S c i e n c e s Chatter reduction through active vibration damping 8 x 105 7 K cut 6 5 4 With active damping Without active damping 3 2 1 0

More information

CHALMERS, GÖTEBORGS UNIVERSITET. EXAM for DYNAMICAL SYSTEMS. COURSE CODES: TIF 155, FIM770GU, PhD

CHALMERS, GÖTEBORGS UNIVERSITET. EXAM for DYNAMICAL SYSTEMS. COURSE CODES: TIF 155, FIM770GU, PhD CHALMERS, GÖTEBORGS UNIVERSITET EXAM for DYNAMICAL SYSTEMS COURSE CODES: TIF 155, FIM770GU, PhD Time: Place: Teachers: Allowed material: Not allowed: August 22, 2018, at 08 30 12 30 Johanneberg Jan Meibohm,

More information

7 Planar systems of linear ODE

7 Planar systems of linear ODE 7 Planar systems of linear ODE Here I restrict my attention to a very special class of autonomous ODE: linear ODE with constant coefficients This is arguably the only class of ODE for which explicit solution

More information

Autonomous system = system without inputs

Autonomous system = system without inputs Autonomous system = system without inputs State space representation B(A,C) = {y there is x, such that σx = Ax, y = Cx } x is the state, n := dim(x) is the state dimension, y is the output Polynomial representation

More information

5.2.2 Planar Andronov-Hopf bifurcation

5.2.2 Planar Andronov-Hopf bifurcation 138 CHAPTER 5. LOCAL BIFURCATION THEORY 5.. Planar Andronov-Hopf bifurcation What happens if a planar system has an equilibrium x = x 0 at some parameter value α = α 0 with eigenvalues λ 1, = ±iω 0, ω

More information

6. Well-Stirred Reactors III

6. Well-Stirred Reactors III 6. Well-Stirred Reactors III Reactors reaction rate or reaction velocity defined for a closed system of uniform pressure, temperature, and composition situation in a real reactor is usually quite different

More information

Numerical Continuation of Bifurcations - An Introduction, Part I

Numerical Continuation of Bifurcations - An Introduction, Part I Numerical Continuation of Bifurcations - An Introduction, Part I given at the London Dynamical Systems Group Graduate School 2005 Thomas Wagenknecht, Jan Sieber Bristol Centre for Applied Nonlinear Mathematics

More information

Theoretical physics. Deterministic chaos in classical physics. Martin Scholtz

Theoretical physics. Deterministic chaos in classical physics. Martin Scholtz Theoretical physics Deterministic chaos in classical physics Martin Scholtz scholtzzz@gmail.com Fundamental physical theories and role of classical mechanics. Intuitive characteristics of chaos. Newton

More information

Application of numerical continuation to bifurcation analysis of Rijke tube. Kanpur , India Corresponding author:

Application of numerical continuation to bifurcation analysis of Rijke tube. Kanpur , India Corresponding author: n 3 l - Int l Summer School and Workshop on Non-Normal and Nonlinear Effects in Aero- and Thermoacoustics, May 17-21, 2010, Munich n³l Application of numerical continuation to bifurcation analysis of Rijke

More information

Dynamical systems tutorial. Gregor Schöner, INI, RUB

Dynamical systems tutorial. Gregor Schöner, INI, RUB Dynamical systems tutorial Gregor Schöner, INI, RUB Dynamical systems: Tutorial the word dynamics time-varying measures range of a quantity forces causing/accounting for movement => dynamical systems dynamical

More information

STABILITY ANALYSIS OF DAMPED SDOF SYSTEMS WITH TWO TIME DELAYS IN STATE FEEDBACK

STABILITY ANALYSIS OF DAMPED SDOF SYSTEMS WITH TWO TIME DELAYS IN STATE FEEDBACK Journal of Sound and Vibration (1998) 214(2), 213 225 Article No. sv971499 STABILITY ANALYSIS OF DAMPED SDOF SYSTEMS WITH TWO TIME DELAYS IN STATE FEEDBACK H. Y. HU ANDZ. H. WANG Institute of Vibration

More information

Nonsmooth systems: synchronization, sliding and other open problems

Nonsmooth systems: synchronization, sliding and other open problems John Hogan Bristol Centre for Applied Nonlinear Mathematics, University of Bristol, England Nonsmooth systems: synchronization, sliding and other open problems 2 Nonsmooth Systems 3 What is a nonsmooth

More information

Models Involving Interactions between Predator and Prey Populations

Models Involving Interactions between Predator and Prey Populations Models Involving Interactions between Predator and Prey Populations Matthew Mitchell Georgia College and State University December 30, 2015 Abstract Predator-prey models are used to show the intricate

More information

UNCERTAINTY PROPAGATION FOR SELECTED ANALYTICAL MILLING STABILITY LIMIT ANALYSES

UNCERTAINTY PROPAGATION FOR SELECTED ANALYTICAL MILLING STABILITY LIMIT ANALYSES UNCERTAINTY PROPAGATION FOR SELECTED ANALYTICAL MILLING STABILITY LIMIT ANALYSES G. Scott Duncan, Mohammad H. Kurdi, Tony L. Schmitz Department of Mechanical and Aerospace Engineering University of Florida

More information

1859. Forced transverse vibration analysis of a Rayleigh double-beam system with a Pasternak middle layer subjected to compressive axial load

1859. Forced transverse vibration analysis of a Rayleigh double-beam system with a Pasternak middle layer subjected to compressive axial load 1859. Forced transverse vibration analysis of a Rayleigh double-beam system with a Pasternak middle layer subjected to compressive axial load Nader Mohammadi 1, Mehrdad Nasirshoaibi 2 Department of Mechanical

More information

Hardening nonlinearity effects on forced vibration of viscoelastic dynamical systems to external step perturbation field

Hardening nonlinearity effects on forced vibration of viscoelastic dynamical systems to external step perturbation field Int. J. Adv. Appl. Math. and Mech. 3(2) (215) 16 32 (ISSN: 2347-2529) IJAAMM Journal homepage: www.ijaamm.com International Journal of Advances in Applied Mathematics and Mechanics Hardening nonlinearity

More information

Research Article An Efficient Approach for Identifying Stable Lobes with Discretization Method

Research Article An Efficient Approach for Identifying Stable Lobes with Discretization Method Hindawi Publishing Corporation Advances in Mechanical Engineering Volume 213, Article ID 682684, 5 pages http://dx.doi.org/1.1155/213/682684 Research Article An Efficient Approach for Identifying Stable

More information

Dynamical Analysis of a Hydraulic Pressure Relief Valve

Dynamical Analysis of a Hydraulic Pressure Relief Valve Proceedings of the World Congress on Engineering 9 Vol II WCE 9, July -, 9, London, U.K. Dynamical Analysis of a Hydraulic Pressure Relief Valve Gábor Licskó Alan Champneys Csaba Hős Abstract A mathematical

More information

Response of A Hard Duffing Oscillator to Harmonic Excitation An Overview

Response of A Hard Duffing Oscillator to Harmonic Excitation An Overview INDIN INSTITUTE OF TECHNOLOGY, KHRGPUR 710, DECEMBER 8-0, 00 1 Response of Hard Duffing Oscillator to Harmonic Excitation n Overview.K. Mallik Department of Mechanical Engineering Indian Institute of Technology

More information

CHAPTER 6 HOPF-BIFURCATION IN A TWO DIMENSIONAL NONLINEAR DIFFERENTIAL EQUATION

CHAPTER 6 HOPF-BIFURCATION IN A TWO DIMENSIONAL NONLINEAR DIFFERENTIAL EQUATION CHAPTER 6 HOPF-BIFURCATION IN A TWO DIMENSIONAL NONLINEAR DIFFERENTIAL EQUATION [Discussion on this chapter is based on our paper entitled Hopf-Bifurcation Ina Two Dimensional Nonlinear Differential Equation,

More information

Chaotic motion. Phys 420/580 Lecture 10

Chaotic motion. Phys 420/580 Lecture 10 Chaotic motion Phys 420/580 Lecture 10 Finite-difference equations Finite difference equation approximates a differential equation as an iterative map (x n+1,v n+1 )=M[(x n,v n )] Evolution from time t

More information

Periodic Skeletons of Nonlinear Dynamical Systems in the Problems of Global Bifurcation Analysis

Periodic Skeletons of Nonlinear Dynamical Systems in the Problems of Global Bifurcation Analysis Periodic Skeletons of Nonlinear Dynamical Systems in the Problems of Global Bifurcation Analysis M Zakrzhevsky, I Schukin, A Klokov and E Shilvan PERIODIC SKELETONS OF NONLINEAR DYNAMICAL SYSTEMS IN THE

More information

Analysis of a lumped model of neocortex to study epileptiform ac

Analysis of a lumped model of neocortex to study epileptiform ac of a lumped model of neocortex to study epileptiform activity Sid Visser Hil Meijer Stephan van Gils March 21, 2012 What is epilepsy? Pathology Neurological disorder, affecting 1% of world population Characterized

More information

Math 312 Lecture Notes Linear Two-dimensional Systems of Differential Equations

Math 312 Lecture Notes Linear Two-dimensional Systems of Differential Equations Math 2 Lecture Notes Linear Two-dimensional Systems of Differential Equations Warren Weckesser Department of Mathematics Colgate University February 2005 In these notes, we consider the linear system of

More information

Stability of Feedback Solutions for Infinite Horizon Noncooperative Differential Games

Stability of Feedback Solutions for Infinite Horizon Noncooperative Differential Games Stability of Feedback Solutions for Infinite Horizon Noncooperative Differential Games Alberto Bressan ) and Khai T. Nguyen ) *) Department of Mathematics, Penn State University **) Department of Mathematics,

More information

Long-wave Instability in Anisotropic Double-Diffusion

Long-wave Instability in Anisotropic Double-Diffusion Long-wave Instability in Anisotropic Double-Diffusion Jean-Luc Thiffeault Institute for Fusion Studies and Department of Physics University of Texas at Austin and Neil J. Balmforth Department of Theoretical

More information

DYNAMIC BIFURCATION THEORY OF RAYLEIGH-BÉNARD CONVECTION WITH INFINITE PRANDTL NUMBER

DYNAMIC BIFURCATION THEORY OF RAYLEIGH-BÉNARD CONVECTION WITH INFINITE PRANDTL NUMBER DYNAMIC BIFURCATION THEORY OF RAYLEIGH-BÉNARD CONVECTION WITH INFINITE PRANDTL NUMBER JUNGHO PARK Abstract. We study in this paper the bifurcation and stability of the solutions of the Rayleigh-Bénard

More information