ENGINEERING FOR RURAL DEVELOPMENT Jelgava,

Similar documents
APPROXIMATE LIMIT CYCLES FOR THE RAYLEIGH MODEL

Ordinary differential equations

NUMERICAL STUDY OF COEXISTING ATTRACTORS FOR THE HÉNON MAP

Applied Mathematics, 2011, 2, doi: /am Published Online December 2011 (

QUADRATURE is an old-fashioned word that refers to

The Active Universe. 1 Active Motion

Chaos in drive systems

Linear and Non-linear Feedback Control Strategies for a 4D Hyperchaotic System

THE INTERVAL LATTICE BOLTZMANN METHOD FOR TRANSIENT HEAT TRANSFER IN A SILICON THIN FILM

DYNAMICAL ANALYSIS OF AN NON-ŠIL NIKOV CHAOTIC SYSTEM

DISTRIBUTION OF SUB AND SUPER HARMONIC SOLUTION OF MATHIEU EQUATION WITHIN STABLE ZONES

The Application of Duffing Oscillator in Weak Signal Detection

Remark on boundary value problems arising in Ginzburg-Landau theory

Solutions of Klein - Gordan equations, using Finite Fourier Sine Transform

u( t) + K 2 ( ) = 1 t > 0 Analyzing Damped Oscillations Problem (Meador, example 2-18, pp 44-48): Determine the equation of the following graph.

ORDER REDUCTION USING POLE CLUSTERING AND FACTOR DIVISION METHOD

MAC-solutions of the nonexistent solutions of mathematical physics

dx dt dy = G(t, x, y), dt where the functions are defined on I Ω, and are locally Lipschitz w.r.t. variable (x, y) Ω.

221A Lecture Notes WKB Method

Jim Lambers MAT 169 Fall Semester Lecture 4 Notes

Physics 116C Solution of inhomogeneous ordinary differential equations using Green s functions

Application of Exp-Function Method to. a Huxley Equation with Variable Coefficient *

Bend Forms of Circular Saws and Evaluation of their Mechanical Properties

Analytical Based Truncation Principle of Higher- Order Solution for a x 1/3 Force Nonlinear Oscillator

Linear Systems with Constant Coefficients

Reversible magnetization processes in scalar Preisachtype models of hysteresis

SOUND INTENSITY PROBE CALIBRATOR FOR FIELD USE: CALCULATING THE SOUND FIELD IN THE CALIBRATOR USING BOUNDARY ELEMENT MODELLING

State space systems analysis (continued) Stability. A. Definitions A system is said to be Asymptotically Stable (AS) when it satisfies

Lecture 20: Numerical Integration III

Polynomials and Division Theory

HOMOCLINIC ORBITS IN A POROUS PELLET

Time delay in a basic model of the immune response

IDENTIFICATION AND MODIFICATION OF FRAME STRUCTURE

221B Lecture Notes WKB Method

COEXISTENCE OF POINT, PERIODIC AND STRANGE ATTRACTORS

PERFORMANCE ANALYSIS OF HARMONICALLY FORCED NONLINEAR SYSTEMS

Travelling Profile Solutions For Nonlinear Degenerate Parabolic Equation And Contour Enhancement In Image Processing

1.9 C 2 inner variations

Special Relativity solved examples using an Electrical Analog Circuit

Fredholm Integral Equations of the First Kind Solved by Using the Homotopy Perturbation Method

Nonlinear Resonances Analysis of a RLC Series Circuit Modeled by a Modified Van der Pol Oscillator

A Model of two mutually interacting Species with Mortality Rate for the Second Species

New Expansion and Infinite Series

13.4 Work done by Constant Forces

ADAPTIVE CONTROL AND SYNCHRONIZATION OF RÖSSLER PROTOTYPE-4 SYSTEM

1B40 Practical Skills

BIFURCATIONS IN ONE-DIMENSIONAL DISCRETE SYSTEMS

Review of Gaussian Quadrature method

A REVIEW OF CALCULUS CONCEPTS FOR JDEP 384H. Thomas Shores Department of Mathematics University of Nebraska Spring 2007

ADAPTIVE CONTROL AND SYNCHRONIZATION OF LIU S FOUR-WING CHAOTIC SYSTEM WITH CUBIC NONLINEARITY

NUMERICAL INTEGRATION. The inverse process to differentiation in calculus is integration. Mathematically, integration is represented by.

Diverse modes of eco-evolutionary dynamics in communities of antibiotic-producing microorganisms

Harmonic Mean Derivative - Based Closed Newton Cotes Quadrature

International Journal of Solids and Structures

Jack Simons, Henry Eyring Scientist and Professor Chemistry Department University of Utah

On the Uncertainty of Sensors Based on Magnetic Effects. E. Hristoforou, E. Kayafas, A. Ktena, DM Kepaptsoglou

Solving the (3+1)-dimensional potential YTSF equation with Exp-function method

STABILITY AND BIFURCATION ANALYSIS OF A PIPE CONVEYING PULSATING FLUID WITH COMBINATION PARAMETRIC AND INTERNAL RESONANCES

Chaotic Dynamics of an Extended Duffing Oscillator Under Periodic Excitation

ADAPTIVE CONTROL AND SYNCHRONIZATION OF A HIGHLY CHAOTIC ATTRACTOR

Fig. 1. Open-Loop and Closed-Loop Systems with Plant Variations

ODE: Existence and Uniqueness of a Solution

LINEAR ALGEBRA APPLIED

7.1 Integral as Net Change and 7.2 Areas in the Plane Calculus

A discrete chaotic multimodel based on 2D Hénon maps

INTRODUCTION. The three general approaches to the solution of kinetics problems are:

CHAPTER 4a. ROOTS OF EQUATIONS

arxiv:math/ v1 [math.ds] 23 Aug 2006

Lecture Solution of a System of Linear Equation

Section 6.1 INTRO to LAPLACE TRANSFORMS

Properties of Integrals, Indefinite Integrals. Goals: Definition of the Definite Integral Integral Calculations using Antiderivatives

63. Representation of functions as power series Consider a power series. ( 1) n x 2n for all 1 < x < 1

Studies on Nuclear Fuel Rod Thermal Performance

Simulated climate vegetation interaction in semi-arid regions affected by plant diversity

Chapter 4 Contravariance, Covariance, and Spacetime Diagrams

A Critical Path Problem Using Intuitionistic. Trapezoidal Fuzzy Number

ENGI 9420 Lecture Notes 7 - Fourier Series Page 7.01

A Brief Note on Quasi Static Thermal Stresses In A Thin Rectangular Plate With Internal Heat Generation

SOLUTION OF QUADRATIC NONLINEAR PROBLEMS WITH MULTIPLE SCALES LINDSTEDT-POINCARE METHOD. Mehmet Pakdemirli and Gözde Sarı

Model Reduction of Finite State Machines by Contraction

Chapter 3 Solving Nonlinear Equations

Improper Integrals, and Differential Equations

When a force f(t) is applied to a mass in a system, we recall that Newton s law says that. f(t) = ma = m d dt v,

Global organization of spiral structures in bi-parameter space of dissipative systems with Shilnikov saddle-foci

Note 12. Introduction to Digital Control Systems

Designing finite automata II

How to simulate Turing machines by invertible one-dimensional cellular automata

On the Decomposition Method for System of Linear Fredholm Integral Equations of the Second Kind

OPEN NEWTON - COTES QUADRATURE WITH MIDPOINT DERIVATIVE FOR INTEGRATION OF ALGEBRAIC FUNCTIONS

Conservation Law. Chapter Goal. 5.2 Theory

On the application of explicit spatial filtering to the variables or fluxes of linear equations

Industrial Electrical Engineering and Automation

ANALYSIS OF MECHANICAL PROPERTIES OF COMPOSITE SANDWICH PANELS WITH FILLERS

STEP FUNCTIONS, DELTA FUNCTIONS, AND THE VARIATION OF PARAMETERS FORMULA. 0 if t < 0, 1 if t > 0.

Math 1B, lecture 4: Error bounds for numerical methods

Bifurcation structures in maps of Hénon type

Impact of the tribological characteristics on the dynamics of the ultrasonic piezoelectric motor

Review of Calculus, cont d

Pressure Wave Analysis of a Cylindrical Drum

Method of Localisation and Controlled Ejection of Swarms of Likely Charged Particles

Transcription:

RESERCH IN LOL DYNMICS O DUIN-TYPE OSCILLTOR Y METHOD O COMPLETE IURCTION NLYSIS Ris Smirnov, Jurij Ivnov, Vldimir Nikishin Rig Technicl University, Ltvi rj@df.rtu.lv, ijm@df.rtu.lv, nikis@df.rtu.lv strct. Numerous reserches of the Duffing-type oscilltor hve reveled ig vriety of dynmics of this nonliner system. Nonliner effects, which result from chnges of initil conditions nd prmeters, cn e mnifested in vrious, sometimes seemingly unexpected, wys. The tsk of studying nonliner effects hs importnt prcticl vlue for engineering dynmicl systems working y the Duffing-type oscilltor. or the reserch of nonliner effects we usully use the method of hrmonic lnce, nlog simultion or the method of direct scnning. Such techniques re little limited regrding the ility to predict ll dynmicl chrcteristics, including ifurctions tht led to occurrence of unstle, suhrmonic nd chotic solutions. The given work is devoted to the study of complicted dynmics of the Duffing-type oscilltor nd the resons of its occurrence t the chnge of the system prmeters on the sis of the system pproch y the method of complete ifurction nlysis. The results re received y direct numericl simultion with the use of softwre NLO nd SPRIN. ifurction digrms t the reserch of dynmics of the Duffing-type oscilltor re constructed y mens of the prmeter continution pproch of stedy nd unstle periodic solution. The exmple of the Duffing-type oscilltor illustrtes high efficiency of the method of complete ifurction nlysis for providing n explntion of the resons of the irth of regulr nd chotic solutions, their coexistence nd influence on trnsient processes. t the systemtic pproch to certin prmeters, t which the complete ifurction nlysis specifies the preconditions of the irth of unexpected nonliner effects, in the investigted system n dditionl nlysis of vrious periodic nd chotic ttrctors, forecsting of rre ttrctors nd rre phenomen, studying of scenrios of their irth is crried out. Keywords: nonliner effect, dynmicl system, ifurction digrm, rre ttrctor, chos. Introduction n importnt tsk for ll nonliner dynmicl systems is the study of the nonliner effects tht llow predicting nd explining complicted ehviour of such systems. Modern reserches of nonliner phenomen in the nture nd engineering re relted to the study of typicl structures in the ehviour of nonliner dynmicl systems. This work is devoted to the reserch of regulr nd chotic dynmics of the Duffing-type oscilltor nd systemtic reveling of nonliner effects. The work illustrtes the dvntge of the method of complete ifurction groups for the reserch of nonliner effects in comprison to the pplied methods: hrmonic lnce nd direct scnning. Dynmicl model nd ifurction nlysis for regulr nd chotic dynmics reserch The im of the given work is to reserch in the dynmics of the Duffing-type oscilltor (), study regulr nd chotic solutions, nonliner effects [-]. The uthors pply up-to-dte numericl nd nlyticl methods of the nlysis of nonliner dynmicl systems. irstly, these re the points mpping method nd the prmeter continution method. It is offered to pply systemtic pproch to the reserch of periodic nd chotic forced oscilltions in nonliner systems: serch of periodic regimes, contour mpping, constructing of regime ttrction domins, construction of ifurction digrms. In this pper, periodic nd non-periodic stedy or trnsition oscilltions of the influence of the frequency of externl hrmonic forcing re oserved in computer-simulted solutions of the secondorder nonliner ordinry differentil eqution of the Duffing-type oscilltor (): ( t) & x& x& x x h cos, () where x displcement; nd coefficients of the Duffing-type restoring force; dissiption coefficient of the liner dissiptive force; (t) hcos(t) hrmonic driving force with the period Т π ; h mplitude; w frequency. DOI: 0.66/ERDev07.6.N

Let us show the ppliction of the method of hrmonic lnce [] for the reserch of the Duffingtype oscilltor. t clcultion y the method of hrmonic lnce on the first hrmonics the solution of eqution () is set in the form of (). Prmeters nd re unknown vlues. Solution () is plced in eqution () nd the mplitudes re equted t the sine nd cosine of the first hrmonics. fter the trnsformtions we receive system of two nonliner lgeric equtions (). The system of equtions () is solved numericlly t zero initil pproch in softwre MTL. System () is reduced to one cuic eqution () for the mplitude of the sinus components of the solution. We find the root of cuic eqution () nd determine the mplitude of the cosine component of the solution sed on expression (5). Cuic eqution () is recorded in the cnonicl form (6). or the cnonicl form (6) the discriminnt of the eqution is (7). t positive discriminnt (7) cuic eqution () hs one rel root; t negtive discriminnt the cuic eqution will hve three rel roots (ig..). The system descried y eqution () possesses one or three possile regimes ccording to the quntity of rel roots of cuic eqution (). In cse of coexistence of three regimes of the sic group there is, so-clled in the electricl engineering, the trigger or hysteresis effect when there is discontinuous chnge of the regime t the chnge (increse nd reduction) of the prmeters (see ig. ). 6 ( h) ( t) cos( t) x sin. () ) ( ) h ) ( ) ) 0. () h ( ) h 0 ) 6 s 6 z. () h. (5) r s z 0. (6) 8 r p q D, p 7 ) h ( h) ( ). h, ( ( ) ), s r, q r r s z 7. (7) When clculting with the help of the hrmonic lnce method on the first nd third hrmonics the solution of eqution () is set in the form of (8). Prmeters,,, re unknown vlues. Solution (8) is plced in eqution () nd the mplitudes re equted t the sine nd cosine of the first nd third hrmonics. fter the trnsformtions we receive system of four nonliner lgeric equtions (). The system of equtions () is solved numericlly t zero initil pproch in softwre MTL. ( t) cos( t) sin( t) cos( t) x sin (8) 0

0 0 0 h () ig.. ifurction digrms of the Duffing-type oscilltor of sic periodic regimes t the chnge of the frequency of externl hrmonic forcing pplying the hrmonic lnce method on the first hrmonics (left fig.) nd pplying the method of complete ifurction nlysis of the prmeter continution pproch of stle (drk gry) nd unstle (light gry) periodic solutions with the sic ifurction group Т (right fig.). System prmeters:,, 0.08, h 5, w vr The method of hrmonic lnce (HM) is the pproched method of definition of periodic regimes in nonliner oscilltory systems, sed on the theory tht t presence of nonlinerity the estlished oscilltions re frequently close to the hrmonic. The HM is lso suitle for the study of oscilltions of utonomous systems. The HM in mny cses cn give successful results in terms of qulity nd often in quntity not only t smll nonlinerities, ut lso for essentilly nonliner systems. High efficiency HM in the reserch of the sic regime cn e explined y filtering properties of mny nonliner systems, which llows pproximting the solution y the sum of hrmonics. or the reserch of su-hrmonic nd chotic solutions the pplicility of the method is strongly limited y the complexity of trnsformtions (see ig. nd ig. ). or the disply of regulr nd chotic sttionry regimes (ttrctors) in ifurction digrms the scnning method ( method of direct scnning) is usully used [5]. This method llows defining quntittive vlues of phse coordintes of ttrctors tht correspond to the results of the integrtion fter the completion of the trnsition process t the movement on the prmeter of one or severl initil conditions elonging to the re of ttrction of the sttionry regime (see ig. ). To receive complete ifurction imge we need to crry out whole set of dditionl reserches t fixed vlues of prmeters in order to define the quntity of ttrctors, which crries occsionl rther thn systemtic chrcter. The ide of the RTU professor M. Zkrzhevsky of continution of stedy nd unstle solutions of one ifurction group efore the construction of complete ifurction imge llows crrying out complete ifurction nlysis nd systemticlly find nd explin nonliner effects. Under the supervision of M. Zkrzhevsky the Spring nd NLO softwre hve een developed for the reserch of dynmicl systems [6-8]. T w x

The method of complete ifurction groups nd the complex of lgorithms developed on its sis nd the softwre llow finding qulittively new previously unfmilir regulr nd chotic regimes, rre ttrctors (R) nd rre phenomen, new ifurction groups for typicl nd widely pplied clssicl nonliner dynmicl models nd studying the interction of vrious ifurction groups [8-]. The completeness of ifurction group is defined, if we find for it ll stedy nd unstle regimes nd lso connected with this ifurction group in spce of prmeters of dissiptive nonliner dynmicl systems. t vrition of prmeters in nonliner oscilltory systems qulittive chnges in ehviour of the system tke plce, which leds to the loss of stility of the periodic solution, the cscde ifurctions of period douling, to the formtion of res with infinite numer of unstle regimes UPI (Unstle Periodic Infinitium) nd to the irth of chotic ttrctors [0]. In nonliner dynmicl systems the irth of chotic oscilltions is connected with the sic ifurction group Т, suhrmonic ifurction groups nt nd their coexistence []. E T E T E E E 5 ig.. ifurction digrm of the Duffing-type oscilltor pplying the method of complete ifurction nlysis of the prmeter continution pproch of stedy nd unstle periodic solution with the sic ifurction group Т nd suhrmonic ifurction group T (left fig.). System prmeters:,, 0.08, h 5, w vr. Nonliner multistility phenomenon of the coexistence of ttrctors: t w 5 the right figure displys ll the domins of ttrction for stle periodic regimes E, E, E 5 nd coordintes for unstle periodic regimes E, E (see ig. ) E E E 5 E E E E E 5 ig.. Phse portrits nd phse trjectories of stle periodic regimes (ttrctors E, E, E 5 ) nd phse portrits of unstle periodic regimes E, E from the ifurction groups Т nd T (see ig. ). Initil conditions: E (x 5.0, v 6.758); E (x.7577, v 0.56); E (x -.6, v 0.5085); E (x -.75, v 5.86); E 5 (x -0.08568, v 0.07). The nonliner driven Duffing-type oscilltor prmeters:,, 0.08, h 5, w 5

or disply of the chotic trnsient process, regulr nd chotic ttrctors in ifurction digrms the scnning method is used to determine the phse spce of the quntittive vlue. ig.. ifurction digrm (dependence of sttionry regulr periodic nd chotic solutions displcement x nd velocity v) of the Duffing-type oscilltor pplying the scnning method ( method of direct scnning) on the excittion force frequency w. System prmeters: -,, 0., h 7, w vr T T T 5T twins UPI 5T T twins R P5 T twins T R P5 5T twins UPI 5T ig. 5. Complete ifurction digrm of stle (drk gry) nd unstle (light gry) periodic solutions (dependence displcement x nd velocity v) nd chotic motions or ttrctors (UPI regions) of 7 ifurction groups (T, three T nd three 5T) of the driven Duffing-type oscilltor pplying the method of complete ifurction nlysis of the prmeter w (frequency of externl hrmonic forcing) continution. New nonliner phenomen hve een discovered: rre ttrctors (R) (see ig. 6). System prmeters: -,, 0., h 7, w vr w.5 w.666 w.0 ig. 6. Chotic ttrctors (w.5 nd w.0) from initil conditions (x -, v 0) during 50000 periods on the Poincré plne nd the phse portrit of suhrmonic solution - ttrctor P 5 of the Duffing-type oscilltor. t the vlue of the prmeter w.666, we cn oserve trnsient chos (chotic oscilltions) during 70 periods from initil conditions (x -, v 0). Then, periodic ttrctor P 5 (rre ttrctor R) is formed in the system. It is displyed in the middle imge. The nonliner driven Duffing-type oscilltor prmeters: -,, 0., h 7, w vr Efficient lgorithms of softwre Spring llow conducing reserches of regulr nd chotic solutions of ll types ccording to the method of complete ifurction groups. The exmples in ig. 5 of the ifurction nlysis for the Duffing-type oscilltor y mens of softwre Spring illustrte the dvntge nd systemticity of the reserch of dynmicl systems pplying the method of complete ifurction groups.

Results nd discussion This work shows the dvntge of the method of complete ifurction groups in comprison to usully pplied methods of hrmonic lnce nd direct scnning, which do not llow complete predicting of mny nonliner effects, including ifurctions tht led to the occurrence of unstle, su-hrmonic nd chotic solutions. With the exmple of the Duffing-type oscilltor the uthors hve considered nd illustrted new possiilities in the reserch of the irth nd existence of nonliner effects, for exmple, such s multistility, chotic oscilltions nd rre ttrctors. These new possiilities hve ppered owing to the method of complete ifurction groups, which is the sis of the serch of ll sic periodic solutions nd their continution in the investigted prmeter, studying the topology of their ifurctions. Thus, ifurction digrms disply ll ifurction groups, which llow studying the irth of nonliner effects, nlyzing interction of ifurction groups t the chnge of prmeters. typicl element of ifurction group with chotic ttrctor is the res in spce of prmeters with infinite numer of unstle periodic regimes resulting from the cscde of ifurctions of period douling. The systemtic pproch to crrying out ifurction nlysis on the sis of the method of complete ifurction groups llows explining the origin of rre ttrctors - nonliner effects tht re hrd to e dignosed nd which exist in smll intervl of chngele prmeter. Conclusions. The exmple of the Duffing-type oscilltor hs een used to show the dvntge of crrying out ifurction nlysis with the ppliction of the method of complete ifurction groups.. The nlysis of coexistence of vrious ifurction groups hs enled the systemtic study of mechnisms of the irth of nonliner effects in the investigted dynmicl system.. It hs een shown tht the method of complete ifurction groups llows defining wy to chotic oscilltions nd rre ttrctors.. The results of the reserch of ttrctors, the irth of which is closely connected with unstle solutions, cn e used to forecst nd void ctstrophic situtions in electric mechnicl nd other nonliner dynmicl systems. References. Duffing., Erzwungene Schwingungen ei vernderlicher Eigenfrequenz und ihre technische edeutung,. Vieweg, runschweig, /, 8.. Kovcic I., rennn M.J. The Duffing Eqution: Nonliner Oscilltors nd their ehviour. John Wiley nd Sons, Chichester, 0, 0 p.. Ued Y. The Rod to Chos - II, eril Press Inc. - Snt Cruz, 00.. Virtions in Engineering: Reference ook. In 6 volumes. M.: Mchine Engineering, V.. Oscilltions of Nonliner Mechnicl Systems/Edited y I.I. lehmn. 7, p.5 (in Russin). 5. Lkshmnn M., Rjsekr S. Nonliner Dynmics. Integrility, Chos nd Ptterns. Springer- Verlg, erlin, Heidelerg, 00. 6. Zkrzhevsky M., Ivnov Yu., rolov V. NLO: Universl Softwre for lol nlysis of Nonliner Dynmics nd Chos. In Proceeding of the nd ENOC, Prgue 6. v.. pp.6-6. 7. Zkrzhevsky M. New concepts of nonliner dynmics: complete ifurction groups, protuernces, unstle periodic infinitiums nd rre ttrctors. Journl of Viroengineering, Volume 0, Issue, 008, pp. -. 8. Schukin I., Zkrzhevsky M., Ivnov Yu., Kugelevich V., Mlgin V., rolov V. ppliction of softwre SPRIN nd method of complete ifurction groups for the ifurction nlysis of nonliner dynmicl system. Journl of Viroengineering, Volume 0, Issue, 008, pp. 50-58.. Smirnov R.S., Zkrzhevsky M.V., Schukin I.T., Yevstignejev V.Yu. The Influence of Nonliner Dissiption on the irth of Rre ttrctors in Nonliner Dynmics. Interntionl Symposium R on Rre ttrctors nd Rre Phenomen in Nonliner Dynmics, Rīg, 0, pp. -. 0. Smirnov R. The trnsition to chos nd criteri of chotic ehvior in dynmicl systems with liner nd nonliner dmping: DYVIS-0 - the Dynmics of Viroimpct (Strongly Nonliner) Systems, Russin cdemy of Sciences, Moscow- Klin, 0, pp. 6-67.. Smirnov R., Zkrzhevsky M., Schukin I. lol nlysis of the Nonliner Duffing-vn der Pol Type Eqution y ifurction Theory nd Complete ifurction roups Method. Viroengineering Procedi, 0, Vol.6, Iss.,.-.lpp. ISSN 5-05.