NONLINEAR NORMAL MODES OF COUPLED SELF-EXCITED OSCILLATORS Jerzy Warminski 1 1 Department of Applied Mechanics, Lublin University of Technology, Lublin, Poland, j.warminski@pollub.pl Abstract: The main purpose of this paper is to present nonlinear modes formulation which can be applied to a class of coupled self-excited nonlinear oscillatory systems under selected resonance conditions. A many degree of freedom system is transformed into the problem with reduced number of degrees by application of nonlinear normal modes. The modes are formulated with respect to the system s vibration amplitude on the basis of free vibrations of a nonlinear conservative system and then applied to systems with force o parametric excitations. The analysis shows effectiveness of the method for three examples: classical viscous damped oscillator and for self-excited models with parametric and external excitations. An important difference between linear and nonlinear mode formulation is presented for selected resonance conditions. Keywords: Vibrations, Nonlinear Normal Modes, Self- Excitation. 1. INTRODUCTION Interaction between different vibration types may lead to very interesting and untypical phenomena. Due to nonlinear damping and nonlinear elasticity characteristic the system can exhibit untypical and intuitively unexpected phenomena. If the system is additionally excited by periodic force or/and excited parametrically then additional interactions near the resonance neighbourhoods occur. Even for the simplest one degree of freedom model many new dynamical phenomena appear [1], []. Interest of this paper, in particular, is to determine Nonlinear Normal Modes of a coupled many degree of freedom system composed of nonlinear self-excited oscillators. The formulation should allow for reduction of the problem without any loss of dynamics of the original system. Due to nonlinear damping, producing self-excitation, the vibrations modes are strongly, both, displacement (amplitude and velocity dependent. Therefore a transformation by linear normal modes does not guarantee a proper decoupling of the model. Moreover, if the self-excited oscillators are driven by parametrically or by external forces, their response can be periodic, quasi-periodic or chaotic. For the one degree of freedom system the regions of mentioned vibration types can be found by analytical approximate methods or directly by numerical methods. The problem arises if the model has many degrees of freedom. For nonlinear systems the Linear Normal Modes ( can be used by making assumption that the model is weakly nonlinear. However, for strong nonlinearities that approach leads to results which are far from real dynamics. For driven parametrically and self-excited systems two possible vibrations modes may appear: (a resonant modes, inside the parametric resonance region where single frequency response takes place, (b non-resonant modes, manifested by quasi-periodic response. Analysis of the resonant modes has been presented in the papers [3], [5], [6] while the nonlinear modes of coupled self-excited oscillators (an autonomous model without internal resonance condition have been investigated in [4]. Despite strong nonlinear damping, the system s motion has been decoupled successfully into two independent modal self-excited oscillators by the centre manifold reduction [8]. The present paper is a review of the former analysis. The nonlinear normal modes ( of coupled self-excited oscillators under resonance conditions are formulated for selected resonance regions. The differences between linear and nonlinear normal modes formulation for selected models and parameters are presented.. MODEL OF THE VIBRATING SYSTEM Let us consider the model of coupled nonlinear oscillators with dynamics represented by a set of differential equations mẍ+kx+f s (X+f d (X, Ẋ+f p (X, Ωτ =f e (ωτ (1 where X and τ are dimensionless coordinate and time respectively, m, k are mass and stiffness matrixes related to the linear part, f d (X, Ẋ, f p (X, Ωτ, f e (ωτ represent nonlinear damping functions, parametric and external excitations. Equation (1 is written in a general form which include nonlinear terms dependent on displacement and velocity as well as time dependent terms. All nonlinear functions couple the system s vibration modes. A linear decoupling procedure Serra Negra, SP - ISSN 178-3667 1
can be applied only if the system is weakly nonlinear or linear with small damping. It would be rather difficult to propose a method enables decoupling of the system and solving the problem in general. Therefore, Eq.(1 is considered only for selected cases. The case of autonomous nonlinear selfexcited system has been studied in [4]. In the present paper we demonstrate a method which can be applied for nonlinear non-autonomous systems near resonance regions. The first considered non-autonomous systems is a classical nonlinear Duffing oscillator with weak viscous damping and then the same technique is applied for self-excited oscillators forced parametrically, and next parametrically and externally. Differential equation (1 can be transformed in a classical way by linear coordinate transformation X = uy ( where u is a linear modal vector and Y vector of normal coordinates. Substituting ( into (1 and next multiplying both equation sides by u T we get M Ÿ + KY = ε [ F(Y, Ẏ+ Q cos ωτ] (3 where: M = diag[m j ]=u T mu, K = diag[m j ω 0j ]= u T ku, and ω 0j, means j natural frequency of a linear model. A small parameter ε is used formally for collecting all nonlinear terms, damping in functions F(Y, Ẏ, and excitation in Q. If the parameter ε is equal to zero, then we get a set of uncoupled linear equations which can be solved independently. However, in our case the system is coupled and the nonlinear normal modes have to be formulated bearing in mind system s nonlinear features. 3. RESONANT NONLINEAR NORMAL MODES FORMULATION In a classical approach we can assume that vibration modes are close to the linear model and then we can solve only one from a set of equations (3, independently finding Y k = Y k (t and putting Y j = 0, j = 1,...,n, j = k. The methodology is based on the assumption that near the resonance zones the response corresponds to one "almost" linear vibrating mode. However, received in this way one mode response may not be correct if nonlinear terms play an essential role. The approach can be extended by making assumption that vibration modes of a forced system are close to the modes of free vibrations of a nonlinear conservative system. If we are able to find nonlinear normal modes of the nonlinear conservative system then on this basis we can find resonance response for "relatively" strong nonlinear models. To formulate nonlinear vibration modes we separate nonlinear elasticity terms in equation (1. mẍ + kx + ε f s (X =0 (4 where f s (X represent nonlinear conservative forces. It is worth to note that in general it is difficult to find a solution of the system (4, therefore the formal small parameter is still kept. We seek the approximate solutions in the form X = au(a cos ω 0 τ (5 where a and ω 0 are amplitude and frequency of free vibrations and u(a is an amplitude dependent eigenvector. For the convenience its coefficients are normalised with respect to the first coefficient, which is assumed to be equal one, u = col[1,u (a,...,u n (a]. Substituting the solution (5 into (4, expanding all nonlinear terms in Fourier series and then grouping them in proper harmonics we can find approximate solutions by the harmonic balance method. In the first order approximation we reject all terms of harmonics of order higher than one. This simplification leads to a set of nonlinear algebraic equations that relates modes coefficients with amplitude, u j = f j (a, j =1,,...,n. Having such formulated amplitude dependent modal coefficients we can transform coordinates (Eq.( X = u(ay (6 Substituting (6 to (4 and multiplying by u T (a we can notice that for a periodic function Y j = acosω 0j τ, the system is reduced to a set of independent nonlinear differential equations M j (a Ÿj + M j (aω 0j(a Y j =0,j=1,,...,n (7 where M j = m 1 u 1j + m u j +... + m nu nj are modal masses, and ω 0j is natural frequency of the nonlinear system found from the first balanced harmonic. The details of formulation of Resonant Nonlinear Modes has been published in papers [5], [6]. 4. RESONANT LINEAR AND NONLINEAR NOR- MAL MODES OF VISCOUS DAMPED NONLIN- EAR MODEL As an example we consider a two degree of freedom system presented in Fig.1 excited by external periodic force. The system may have relatively strong nonlinear elasticity and weak linear viscous damping. x 1 m 1 x f (x 1 1 f (x 1 1,x f (x m c 1 c 1 c qsinωτ Figure 1 Model of coupled nonlinear oscillators with viscous damping and external excitation Serra Negra, SP - ISSN 178-3667
Differential equations of motion are written as m 1 Ẍ 1 + c 1 Ẋ 1 + δ 1 X 1 + γ 1 X 1 3 + c 1 (Ẋ1 Ẋ δ 1 (X 1 X +γ 1 (X 1 X 3 = q cos ωτ m Ẍ + c Ẋ + δ X + γ X 3 c 1 (Ẋ1 Ẋ + δ 1 (X 1 X γ 1 (X 1 X 3 =0 (8 Taking into account transformation (, the differential equation of motion (3 take form M 1 Ÿ 1 + M 1 ω01y 1 = ε (u 11 F1 + u 1 F (9 M Ÿ + M ω0y = ε (u 1 F1 + u F where Note that the frequencies are functions of vibration amplitude and the mode coefficients are normalised in such a way that u 1j =1. Applying nonlinear normal modes the differential equation of motion take form M 1 (aÿ1 + C 1 (aẏ1 + M 1 (aω 01(aY 1 = Q 1 cos ωτ M (aÿ + C (aẏ1 + M (aω 0(aY = Q cos ωτ (1 where modal masses M 1 (a, M (a are defined by Eq(7, and C 1 = c 1 u 11 + c u 1 + c 1 (u 11 u 1 C = c 1 u 1 + c u + c 1 (u 1 u (13 are modal damping coefficients, while and M 1 = m 1 u 11 + m u 1, M = m 1 u 1 + m u F 1 = F s1 + F s1 + F d1 + F d1 + F e1 F = F s F s1 + F d F d1 + F e F s1 = γ 1 (u 11 Y 1 + u 1 Y 3 F s = γ (u 1 Y 1 + u Y 3 F s1 = γ 1 [(u 11 + u 1 Y 1 (u 1 + u Y ] 3 F d1 = α 1 (u 11 Ẏ 1 + u 1 Ẏ F d = α (u 11 Ẏ 1 + u 1 Ẏ [ F d1 = α 1 (u 11 + u 1 Ẏ1 (u 1 + u Ẏ F e1 = u 11 q 1 cosωτ, Fe = u 1 q 1 cosωτ, ] 3 Q 1 = q 1 u 11,Q = q 1 u 1 (14 denote modal amplitudes of excitation. We may notice the difference in a structure of the equation (1 and equation (9. Damping in Eq.(1 is the term which has not been included in the nonlinear modes formulation, therefore should be assumed as small. The solution of the first and the second of Eq.(1 is obtained immediately in form Y 1 = a 1 cos (ωτ φ 1 (15 Y = a cos (ωτ φ where amplitude and phase are determined from equations known from theory of linear oscillators Coefficient γ 1, γ, γ 1 correspond to nonlinear Duffing type stiffness of the first and the second oscillator, γ 1 to nonlinear stiffness of the coupling spring. Viscous damping is denoted by parameters α 1, α, and α 1, respectively. In the classical approach we can assume that vibration modes are close to the linear and then we can solve the first of the equations (9 independently putting Y =0. In the considered example one mode solution has been found by Multiple Scale of time method (see numerical results below. However, received in this way one mode response leads to incorrect results if nonlinear terms are strong. According to the method described in previous chapter the nonlinear modes are formulated on the basis of a nonlinear conservative system. The first balanced harmonic leads to the amplitude and the vibration mode coefficient relationship 3 4 a [ m u ( γ1 u 1 3 + γ 1 (u 1 u 3 + m 1 u 1 ( γ1 (u 1 u 3 γ u 3 ] +(u 1 u (m 1 u 1 + m u δ 1 + u 1 u (m δ 1 m 1 δ =0 (10 and to the natural frequencies of the nonlinear conservative system [ ω0j = 1 m 1 u 1j +m δ u 1 u 1j + δ 1 (u 1j u j + δ u j j + ( 3 4 a γ 1 u 4 1j + γ ] 1(u 1j u j 4 + γ u 4 j (11 a 1 = a = Q 1 M 1 (ω 01 ω +( C1 M 1 ω Q M (ω 0 ω +( C M ω, tan φ 1 =, tan φ = C 1 ω M 1(ω 01 ω C ω M (ω 0 ω (16 We have to remember that Eq.(16 are to be solved together with constrain equation (10 and the frequency equation (11 as functions of amplitude. Exemplary analysis is performed for the following numerical data m 1 =1,m =,δ 1 =1,δ =1,δ 1 =0.3, γ 1 =0.1,γ =0.1,γ 1 =0 (17 Linear vibration modes are represented by vectors, mode I: u 10 = col[1,.3770], mode II: u 10 = col[1, 0.1035]. Index zero points out the linear coefficients. Natural linear frequencies are: ω 01 =0.7661, ω 0 =1.1675. The resonance curves determined by linear ( and nonlinear ( normal modes are presented in Fig.. As can be seen the resonance region found by various modes formulation leads to essentially different results. Near the first natural frequency ω 01 the resonance zone and the curve determined by are much smaller then those received by classical one mode approach. Opposite situation takes place near the second natural frequency ω 0, where linear approach ( reduces real resonance region and its amplitudes. Serra Negra, SP - ISSN 178-3667 3
a 1 a.5 1.5 1 0.5 0 5 4 3 1 (a 0.4 0.6 0.8 1 1. ω (b paper can be used properly under the condition that the response of the system is periodic. The periodic response of the model (18 is observed near the resonance regions. In these regions, the self-excitation frequency is quenched by excitation. Therefore, in spite of existence of the nonlinear self-excited terms motion of the system is periodic. Over those regions the response is quasi-periodic and then another approach is desired. The nonlinear modes are determined near resonance regions making the same assumption as before that a mode shape of a nonlinear conservative system is close to that of a full complex structure. Applying the procedure described in chapters 3 and 4 we received reduced model of the considered problem M j Ÿ j + M j ω 0j Y j C αj Ẏ j + C βj Ẏ 3 j C µj Y j cos ϑτ = Q j cos ωτ (19 Modal mass M j, natural frequency ω 0j and external excitation Q j remain with the same definitions as in the former chapter. New parameters occurring in (19 correspond to self- and parametric excitation and they are defined as C αj = α 1 u 1j + α u j, C βj = β 1 u 4 1j + β u 4 j - coefficients of nonlinear modal damping, C µj =(u 1j u j - modal amplitude of parametric excitation and frequency ϑ. The new terms complicate the solutions comparing with the former example. To get solution an approximate method has to be used. Therefore, equation (19 is written in form with a formal small parameter engaged in the procedure 0 0.8 1. 1.6 Figure Resonance curves for linear ( and nonlinear ( normal modes, (a resonance near ω 01, (b resonance near ω 0 5. NONLINEAR NORMAL MODES OF FORCED SELF-EXCITED OSCILLATORS In this chapter we apply the formulated for a more complex system composed of coupled oscillators presented in Fig.(1, however damping of the oscillators is assumed to be nonlinear and producing self-excitation, moreover apart from external excitation, parametric coupling is assumed due to periodically changing stiffness of the coupling spring. Differential equations of motion for such a system take form m 1 Ẍ 1 + f d1 (Ẋ1 + δ 1 X 1 + γ 1 X 3 1 +(δ 1 µ cos ϑτ (X 1 X =qcos ωτ m Ẍ + f d (Ẋ + δ X + γ X 3 (δ 1 µ cos ϑτ (X 1 X =0 (18 Let us assume that damping is strongly nonlinear represented by Rayleigh s terms, f d1 = α 1 Ẋ 1 + βẋ3 1, f d = α Ẋ +βẋ3. Resonance normal modes formulated in this ( M j Ÿ j + M j ω0j Y j = ε Cαj Ẏ j C βj Ẏj 3 + C µj Y j cos ϑτ + C (0 qj cos ωτ where C αj = ε C αj, C βj = ε C βj, C µj = ε C µj, C qj = ε C qj. We assume small and positive value of ε. According to multiple scale method we seek a periodic solution in series of a small parameter Y j (τ,ε=y j0 (T 0,T 1 +εy j1 (T 0,T 1 +... (1 near the resonance region ϑ = ω 0j + εσ j ( where σ j is a detuning parameter near the first or the second resonance region. Substituting solution (1 into differential equations of motion (0 we get set of recurrent equations in successive perturbation orders D0Y j0 + ϑ Y j0 =0 [ D0Y j1 + ϑ Y j1 = σ j Y j0 D 0 D 1 Y j0 + C αj D 0 D 1 Y j0 C βj (D 0 D 1 Y 10 3 + C µj Y j0 cos ϑτ + C ] qj cos ωτ /M j (3 Further analysis is curried out for two variants, (a a case of parametric and self-excited system, Cqj =0, and (b a case of self-excited system with parametric and external excitation, Cqj > 0 and ω = ϑ. A secular terms elimination Serra Negra, SP - ISSN 178-3667 4
condition leads to the modulation equations M j ϑȧ j = 1 C αj ϑa j 3 8 C βj ϑ 3 a 3 j 1 4 C µj a sin φ j 1 C qj sin φ M j ϑa j φ j = 1 σ jm j a j 1 4 C µj a j cos φ j 1 C qj cos φ (4 In a steady state amplitude ȧ j =0and phase φ j =0. Thus we can determine resonance curve of the system Cqj a jcµj Λ j +C qj Cqj +4a j C µjλ j Λ j =0 (5 where ( λ 1j = a j Cαj ω 3 4 a j C β j ω 3 (, λ j = M j ω 0j ω, Λ j = ( Cqj a j C ( µjλ j +4a j C µj Cµj λ 1j +C qj λ j Equation (5 has to be solved together with constrain equation (10. The approximate solution has form Y j = a [ j cos (ωt + φ+ ε a j 16M j C µj cos (3ωt + φ+ a j C ] βj ω sin (3ωt +3φ (6 If the coefficient C qj is taken as equal to zero, then we get resonance caused by parametric excitation, otherwise both, parametric and externally forced vibration occur. It is worth to note that the formulated nonlinear normal modes are valid for periodic system s response. Modulation equations (4 allows to determine solution also for non stationary regimes which is in contrast to assumption made in nonlinear mode formulation. However, we may expect that for slowly varying amplitude and phase the formulation should remain valid. Numerical example is presented for data m 1 =1,m =,α 1 =0.01, β 1 =0.05, α =0.01, β =0.05,δ 1 =0.3, µ =0.3, q=0, or, µ =0., q=0.1 (7 The modal coefficients take values for first mode; M 1 = 1.3004, C µ1 = 0.379, C α1 =0.0665, C β1 =1.646 for second mode; M = 1.0885, C µ = 0.930, C α =0.01044, C β =0.0501 0.75 0.10 Μ 0.15 0.0 0.80 ϑ 0.5 0.30 0.8 0.6 0.4 a 0. 0.0 0.85 Figure 3 Amplitude a versus amplitude µ and frequency ϑ of parametric excitation received by linear normal modes approach ( near the first natural frequency ω 01 amplitude 0.8 0.6 0.4 0..40 u1.35.30.5.0.15 (a 0.76 0.78 0.80 0.8 0.84 (b 0.0 0.1 0. 0.3 0.4 0.5 0.6 0.7 amplitude Figure 4 Parametric resonance curves for linear ( and nonlinear ( normal modes (a, mode coefficient (b, resonance near ω 01 The resonance curves calculated for the assumed data near the first natural frequency ω 01 are presented in Fig.3. The surface shows resonance region versus parametric excitation amplitude µ and frequency ϑ calculated by the linear normal mode transformation. For comparison a cross-section for µ = 0.3 is taken and presented in Fig.4(a. Amplitudes obtained by nonlinear modes ( are essentially smaller than those got by classical linear modes (. Also the resonance zones fit much better to direct numerical simulations [6]. The mode coefficient u 1 versus amplitude is plotted in Fig.4(b. For a linear system this coefficient is independent of amplitude (the horizontal line. Periodic resonances looks totally different if external force additionally excites the system. An exemplary bifurcation diagram amplitude of a versus external excitation q obtained for a one DOF model is presented in Fig.5. It is clearly visible that for selected parameters there are possible five nontrivial solutions [1], []. For amplitude q =0we get two nontrivial and one trivial solutions (starting points in vertical axis and then increasing q we get five branches of nontrivial solutions which transit into three and eventually to one solution. For the considered reduced order model the resonance curve should exhibit five nontrivial solutions in the resonance zones. Numerical ϑ Serra Negra, SP - ISSN 178-3667 5
(a Figure 5 Bifurcation diagram vibration amplitude versus amplitude of excited force q and fixed excitation frequency ϑ = 1.07; α =0.01, β =0.05, γ =0.1, µ =0. (b simulations shows that in fact for selected set of parameters there are multiple solutions in the resonance zones [1]. Taking into account this observation the set of nonlinear equations (5 and (10 has to be solved together. The solution has been received by solving numerically nonlinear algebraic equations for various starting points and then a continuation method used to get resonance branches. The resonance curve determined by s near the second natural frequency are presented in Fig.6(a. The additional solutions are represented by a loop in the resonance curve and there are five possible nontrivial solutions in this regions. Stability analysis shows that only two upper are stable. This result is in a good agreement with numerical simulations (see [1] for details. Linear normal modes transformation leads to a resonance curve presented in Fig.6(b. We see an important difference between both results. The loop is much bigger and is located in a different place, on the left side of the resonance. Moreover, only three nontrivial solutions occur, instead of expected five. It is clearly seen that the proposed s essentially improve quality of the reduced model. 6. CONCLUSION The paper presents a formulation of resonant nonlinear normal modes which allows to reduce dynamics of a many degrees of freedom system into independent decoupled nonlinear oscillators. The reduction is possible if the response of the system is periodic. Therefore, the method can be applied mainly for the resonance region. The formulated nonlinear modes differs from that presented e.g. in [7], [9]. The modes are not time dependent but related with the system s parameters and vibration amplitude. The mode coefficients are determined as functions of amplitude which leads to an additional constrain equation. In this way the resonance curve has been corrected with regards to classical linear transformation. The modes can be successfully applied for a system with weak damping or a system in which damping is eliminated by system s dynamics [10], [11]. The effectiveness of Figure 6 Resonance curves around ω 0; nonlinear modes formulation (a, linear normal modes (b the method is tested for nonlinear oscillator with weak viscous damping and forced by a harmonic force, and then extended for self-excited oscillators driven parametrically and externally. The presented examples show important differences in resonant response for linear and nonlinear modes approach. The results received by formulated nonlinear normal modes fit very well to direct numerical simulations results. ACKNOWLEDGMENTS Financial support of Structural Funds in the Operational Programme - Innovative Economy (IE OP financed from the European Regional Development Fund Project No POIG.0101.0-00-015/08 is gratefully acknowledged. The author would like to thank Prof. Jose Manual Balthazar for his invitation and making possible the attendance in the conference. Serra Negra, SP - ISSN 178-3667 6
REFERENCES [1] K. Szabelski, J. Warminski, The parametric self excited non-linear system vibrations analysis with the inertial excitation, Int. Journal of Non-Linear Mechamics 30,, pp.179-189, 1995. [] K. Szabelski, J.Warminski, The self-excited system vibrations with the parametric and external excitations, Journal of Sound and Vibration, 187,(4, pp. 595-607, 1995. [3] J.Warminski, Nonlinear Vibration Modes of Coupled Parametrically Self-Excited Oscillators, in Nonlinear Modes of Vibrating Systems, EUROMECH 457, Eds. S. Bellizzi, B.Cochelin, C.-H. Lamarque, Frejus, France, Press of the ENTPE, Lyon, pp. 135-138, 004. [4] J. Warminski, Nonlinear normal modes of coupled self-excited oscillators in regular and chaotic vibration regimes, Journal of Theoretical and Applied Mechanics, No. 3, vol. 46, pp. 693-714, 008. [5] W. Szemplinska-Stupnicka, One-modal Analysis of Resonant Vibrations of Non-linear Systems, vol., (in Polish, IPPT PAN, 1977. [6] J. Warminski, Nonlinear Vibration Modes of Parametrically Coupled Self-Excited Oscillators, Machine Dynamics Problems, 3, 4, pp. 107-118, 008. [7] D. Jiang, C. Pierre, S. W. Shaw, Nonlinear normal modes for vibratory systems under harmonic excitation, Journal of Sound and Vibrations, 88, pp.791-81, 005. [8] S. W. Shaw, C. Pierre, Normal Modes for Non-Linear Vibratory Systems, Journal of Sound and Vibrations, 164(1, pp. 85-14, 1993. [9] A. F. Vakakis, Non-Linear Normal Modes (s and Their Applications in Vibration Theory: an Overview, Mechanical Systems and Signal Processing, 11 (1, pp. 3-, 1997. [10] R. Rand, A. Barcilon, T. Morrison, Parametric Resonance of Hopf Bifurcation, Nonlinear Dynamics, 39, pp. 411-41, 005. [11] M. Pandey, R. Rand, A. Zehnder, Frequency Locking in a Forced Mathieu-van der Pol-Duffing System, Nonlinear Dynamics, 54, pp. 3-1, 008. [1] J. Warminski, Nonlinear normal modes of a selfexcited system driven by parametric and external excitations, Nonlinear Dynamics, DOI 10.1007/s11071-010-9679-5, 010, in print. Serra Negra, SP - ISSN 178-3667 7