HARMONIC BALANCE SOLUTION OF COUPLED NONLINEAR NON-CONSERVATIVE DIFFERENTIAL EQUATION

Similar documents
MCH T 111 Handout Triangle Review Page 1 of 3

A NEW COMBINED BRACKETING METHOD FOR SOLVING NONLINEAR EQUATIONS

A Primer on Continuous-time Economic Dynamics

Section 2.3. Matrix Inverses

TIME-VARYING AND NON-LINEAR DYNAMICAL SYSTEM IDENTIFICATION USING THE HILBERT TRANSFORM

Counting Paths Between Vertices. Isomorphism of Graphs. Isomorphism of Graphs. Isomorphism of Graphs. Isomorphism of Graphs. Isomorphism of Graphs

A STUDY OF Q-CONTIGUOUS FUNCTION RELATIONS

Eigenvectors and Eigenvalues

for all x in [a,b], then the area of the region bounded by the graphs of f and g and the vertical lines x = a and x = b is b [ ( ) ( )] A= f x g x dx

k and v = v 1 j + u 3 i + v 2

Solutions for HW9. Bipartite: put the red vertices in V 1 and the black in V 2. Not bipartite!

PAIR OF LINEAR EQUATIONS IN TWO VARIABLES

Factorising FACTORISING.

DIFFERENCE BETWEEN TWO RIEMANN-STIELTJES INTEGRAL MEANS

CS 491G Combinatorial Optimization Lecture Notes

Solids of Revolution

Edexcel Level 3 Advanced GCE in Mathematics (9MA0) Two-year Scheme of Work

SOME INTEGRAL INEQUALITIES FOR HARMONICALLY CONVEX STOCHASTIC PROCESSES ON THE CO-ORDINATES

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

Numbers and indices. 1.1 Fractions. GCSE C Example 1. Handy hint. Key point

50 AMC Lectures Problem Book 2 (36) Substitution Method

Now we must transform the original model so we can use the new parameters. = S max. Recruits

SOME COPLANAR POINTS IN TETRAHEDRON

A Study on the Properties of Rational Triangles

Tikrit Journal of Pure Science 21 (1) 2016 ISSN:

Particle Physics. Michaelmas Term 2011 Prof Mark Thomson. Handout 3 : Interaction by Particle Exchange and QED. Recap

Electromagnetism Notes, NYU Spring 2018

I 3 2 = I I 4 = 2A

Necessary and sucient conditions for some two. Abstract. Further we show that the necessary conditions for the existence of an OD(44 s 1 s 2 )

Hyers-Ulam stability of Pielou logistic difference equation

AP Calculus BC Chapter 8: Integration Techniques, L Hopital s Rule and Improper Integrals

Co-ordinated s-convex Function in the First Sense with Some Hadamard-Type Inequalities

SECTION A STUDENT MATERIAL. Part 1. What and Why.?

where the box contains a finite number of gates from the given collection. Examples of gates that are commonly used are the following: a b

GENERALIZED OSTROWSKI TYPE INEQUALITIES FOR FUNCTIONS WHOSE LOCAL FRACTIONAL DERIVATIVES ARE GENERALIZED s-convex IN THE SECOND SENSE

MATH 122, Final Exam

Electrical Circuits II (ECE233b)

On a Class of Planar Graphs with Straight-Line Grid Drawings on Linear Area

Linear Algebra Introduction

Lecture 1 - Introduction and Basic Facts about PDEs

Math 211A Homework. Edward Burkard. = tan (2x + z)

ERT 316: REACTION ENGINEERING CHAPTER 3 RATE LAWS & STOICHIOMETRY

Chapter Gauss Quadrature Rule of Integration

Consider a potential problem in the half-space dened by z 0, with Dirichlet boundary conditions on the plane z = 0 (and at innity).

The study of dual integral equations with generalized Legendre functions

Reverse Engineering Gene Networks with Microarray Data

Harmonic Balance Solution of Mulholland Equation

1 PYTHAGORAS THEOREM 1. Given a right angled triangle, the square of the hypotenuse is equal to the sum of the squares of the other two sides.

f a L Most reasonable functions are continuous, as seen in the following theorem:

3D modelling of heating of thermionic cathodes by high-pressure arc plasmas

Synchronization of different 3D chaotic systems by generalized active control

CIT 596 Theory of Computation 1. Graphs and Digraphs

Chapter 6 Techniques of Integration

CSE 332. Sorting. Data Abstractions. CSE 332: Data Abstractions. QuickSort Cutoff 1. Where We Are 2. Bounding The MAXIMUM Problem 4

CS 2204 DIGITAL LOGIC & STATE MACHINE DESIGN SPRING 2014

Calculus Cheat Sheet. Integrals Definitions. where F( x ) is an anti-derivative of f ( x ). Fundamental Theorem of Calculus. dx = f x dx g x dx

Algebra 2 Semester 1 Practice Final

MULTIPLE CHOICE QUESTIONS

Inequalities of Olympiad Caliber. RSME Olympiad Committee BARCELONA TECH

F / x everywhere in some domain containing R. Then, + ). (10.4.1)

Antiderivatives Introduction

Lecture 6: Coding theory

18.06 Problem Set 4 Due Wednesday, Oct. 11, 2006 at 4:00 p.m. in 2-106

LESSON 11: TRIANGLE FORMULAE

Generalized Kronecker Product and Its Application

SEMI-EXCIRCLE OF QUADRILATERAL

Math Lesson 4-5 The Law of Cosines

APPROXIMATION AND ESTIMATION MATHEMATICAL LANGUAGE THE FUNDAMENTAL THEOREM OF ARITHMETIC LAWS OF ALGEBRA ORDER OF OPERATIONS

1 nonlinear.mcd Find solution root to nonlinear algebraic equation f(x)=0. Instructor: Nam Sun Wang

Lecture Notes No. 10

Tutorial Worksheet. 1. Find all solutions to the linear system by following the given steps. x + 2y + 3z = 2 2x + 3y + z = 4.

Sturm-Liouville Theory

VISIBLE AND INFRARED ABSORPTION SPECTRA OF COVERING MATERIALS FOR SOLAR COLLECTORS

( ) { } [ ] { } [ ) { } ( ] { }

EXTENSION OF THE GCD STAR OF DAVID THEOREM TO MORE THAN TWO GCDS CALVIN LONG AND EDWARD KORNTVED

On the Co-Ordinated Convex Functions

Chapter Five - Eigenvalues, Eigenfunctions, and All That

SIMPLE NONLINEAR GRAPHS

How many proofs of the Nebitt s Inequality?

How do we solve these things, especially when they get complicated? How do we know when a system has a solution, and when is it unique?

ANALYSIS AND MODELLING OF RAINFALL EVENTS

NON-DETERMINISTIC FSA

On the Spectra of Bipartite Directed Subgraphs of K 4

Compression of Palindromes and Regularity.

TOPIC: LINEAR ALGEBRA MATRICES

APPENDIX. Precalculus Review D.1. Real Numbers and the Real Number Line

Pressure-sensitivity Effects on Toughness Measurements of Compact Tension Specimens for Strain-hardening Solids

More Properties of the Riemann Integral

MAT 403 NOTES 4. f + f =

University of Sioux Falls. MAT204/205 Calculus I/II

C. C^mpenu, K. Slom, S. Yu upper boun of mn. So our result is tight only for incomplete DF's. For restricte vlues of m n n we present exmples of DF's

1.2. Linear Variable Coefficient Equations. y + b "! = a y + b " Remark: The case b = 0 and a non-constant can be solved with the same idea as above.

CS 573 Automata Theory and Formal Languages

ILLUSTRATING THE EXTENSION OF A SPECIAL PROPERTY OF CUBIC POLYNOMIALS TO NTH DEGREE POLYNOMIALS

Integration. antidifferentiation

POSITIVE IMPLICATIVE AND ASSOCIATIVE FILTERS OF LATTICE IMPLICATION ALGEBRAS

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

The DOACROSS statement

Technische Universität München Winter term 2009/10 I7 Prof. J. Esparza / J. Křetínský / M. Luttenberger 11. Februar Solution

Surds and Indices. Surds and Indices. Curriculum Ready ACMNA: 233,

Transcription:

GNIT J. nglesh Mth. So. ISSN - HRMONIC LNCE SOLUTION OF COUPLED NONLINER NON-CONSERVTIVE DIFFERENTIL EQUTION M. Sifur Rhmn*, M. Mjeur Rhmn M. Sjeur Rhmn n M. Shmsul lm Dertment of Mthemtis Rjshhi University of Engineering n Tehnology RUET, Rjshhi, nglesh *msr_mth_@yhoo.om Reeive.. ete.. STRCT moifie hrmoni lne metho is emloye to etermine the seon roximte solutions to oule nonliner ifferentil eution ner the limit yle. The solution shows goo greement with the numeril solution. Keywors: Hrmoni lne metho, Coule vn er Pol eution, Nonliner osilltor. Introution The system of two eutions freuently rises in nonliner osilltions, nonliner ynmis n mthemtil hysis et. Rn n Holmes [] first roose the system of two oule vn er Pol eutions with liner iffusive ouling. They investigte the roerties of ertin erioi motion of two ientil vn er Pol osilltors with wek nonliner ouling. Reently, Neem et. l. [] hs foun roximte first integrl for system of two oule vn er Pol osilltors with liner iffusive ouling. There re mny nlytil rohes for roximting erioi solutions of the nonliner systems. The most wiely use methos re the erturtion methos, in whih the solution is exne in ower series of smll rmeter. The LP metho [], KM metho [-] n multi-time exnsion metho [-] re imortnt mong them. Usully, lower orer e.g., first or seon roximte solution is etermine y the erturtion methos ue to voi lgeri omlexities. To tkle similr nonliner rolems, there re more imortnt roximtion tehniues. One of them is the itertion tehniue see [-]. The hrmoni lne H metho [-] is nother tehniue for etermining erioi solutions of nonliner ifferentil eutions y using the trunte Fourier series. Sine the erivtion of higher roximtion is omlite, the first n seon roximte solutions re usully lulte. The vntge of H metho is tht the solution gives esire result though nonlinerities eome signifint. The im of this rtile is to fin seon roximtion to oule vn er Pol s eution following moifie hrmoni lne tehniue.

Rhmn et l.. The solution metho Let us onsier nonliner oule ifferentil eution u u ε f U, u ε u u, ε << where n x x u, U, x x α, β. erioi solution of E. is hosen in the form x os x os os os sin sin os os sin L sin L where, n re onstnts. In generl the unknown funtions, j, j, j, j, j,, L re etermine together with n the initil hse, n. Now sustituting E. into E. n exning the funtion f U, u in Fourier series, we otin os os sin L ε[ F,,,, L os F,,,, Los L G,, f g,,,,,,,,,, ε os ε ε sin ε ε os ε ε sin ε n Los Lsin G Lsin g, f,,,,,,, os L sin L os L,,, sin L, Lsin L Lsin L Lsin L os os sin L ε[ F,,,, L os F,,,, Los L ]

Hrmoni lne Solution G,, f g,,,,,,,,, ε sin ε, ε os ε ε sin ε Los Lsin G Lsin g ε os ε, f,,,,, os L,,, sin L,,, Lsin L Lsin L os L sin L Lsin L y omring the oeffiients of eul hrmoni from E. n E., we otin εf ε, ε F ε ε G ε, ε G ε n ε f ε, ε f ε ε g ε, ε g ε Utilizing the first eution of E., we eliminte tkes the following form G ε ε F ] from ll the rest. Thus E. ε ε F εf ε }/ { { ε ε F εg ε }/ gin utilizing the first eution of E., we eliminte E. tkes the following form g ε f ε { ε f ε ε f ε }/ { from ll the rest. Thus ε f ε ε g ε }/ We use new rmeter ε, << n ε, << with ε Ο n solve the thir-, fourth-et. eutions of E. n E. in owers of n resetively s

Rhmn et l. j j n j j j, j, j, j, j, j, j, j, j, j, j, j, LL LL LL LL j,,l j,,l Now sustituting the vlues of,, L,, L n,, L,, Lfrom E. n E. resetively into the first eution of E. n E. resetively, we etermine n, n then n. To etermine stey stte solution we strt from x n x. Thus we otin sin sin os L sin sin os L where n. Finlly, sustituting the vlues of,,, L n,,, Linto seon,, eution of E., E. n E., we solve them for,, n. Then sustituting the vlues of,, n into the eution, x x os os os os sin sin L L We otin the vlue of n, whih reresents the initil vlue of x n x for the stey-stte solution.. Exmle Consier oule vn er Pol eution u u ε U u ε u u, ε << x x where u, U x x n α, β Let us onsier erioi solution of the form

Hrmoni lne Solution sin os sin os os sin os sin os os x x Sustituting E. into E., we get sin os os [ sin os ε sin sin sin os os ] os HOH HOH ] os os os sin sin sin [ ε n [ sin os sin os os ε sin sin sin os

Rhmn et l. os ] os HOH HOH ] os os os sin sin sin [ ε where HOH stns for the higher orer hrmonis. Comring the oeffiients of eul hrmonis, we otin ε ε ε ε ε ε ε ε ε ε n ε ε ε ε ε ε

Hrmoni lne Solution ε ε ε ε ε ε ε ε ε With the hel of the first eution of E. n E., we eliminte n resetively from the seon, thir, fifth n sixth eution of E. n E., we otin ε ε ε ε ε ε ε ε

Rhmn et l. ε ε n ε ε ε ε ε ε ε ε ε ε ε ε ε ε ε To etermine stey stte solution we strt from x n x. Thus we otin

Hrmoni lne Solution sin sin sin sin os os sin sin os os Herein fourteen unknown untities,,,,,,,,,,,, n will e lulte from fourteen nonliner eutions esrie y E., E., E. n E.. Certinly E. n E. reresents set of nonliner ε ε lgeri eutions with smll rmeter n. ε ε ε oring to [], we shll e le to fin n roximte solution of E. n E. in the form n,,,,,,,,,,,,,,,,,,,,,,,, L L L L L L L L Sustituting E. n E. into E. n E. resetively n euting the oeffiients of n on oth sies, we will hve system of liner eutions of, L,,, L,,, L,,, L,,,,,,,,,,,,,,,,,,,,,,, L, L,,,L. Solving these eutions, we otin,, L L L L

Rhmn et l. n L L L L Finlly, sustituting the vlues of,,, n,,, into first eution of E., E. n E., we solve them for,, n. Then sustituting the vlues of,, n into the eution x os x os os os sin sin os os sin sin We otin the vlue of x n x, whih reresents the initil vlue of x n x for the stey-stte solution.. Results n Disussion In orer to test the ury of n roximte solution, some uthors [,, ] omre nlytil solutions to those otine y numeril tehniues. We hve omre suh n roximte solution of E. to numeril solution for ε. n ε.. First of ll we lot in Fig. n Fig., the seon roximte solution for ε., with initil onitions x., x, x., x, where unknown L

Hrmoni lne Solution oeffiients,,,,,,, ; mlitues, n initil hses, re lulte y the first eutions of E. n E. together with Es., - n then sustituting these vlues of,,,,,,,,,,, in E., we otin x n x. Then orresoning numeril solution hs een omute y Runge-Kutt fourth-orer metho. From the figure it is ler tht the nlytil solution shows goo oiniene with the numeril solution. In Fig., we hve shown the orresoning hse ifferene etween two osilltors. In Fig. n Fig., the seon roximte solution of E. for ε. with the initil onitions x., x, x., x,,,,,,,, where unknown oeffiients re ; mlitues, n initil hses, re lulte y the first eutions of E. n E. together with eutions, - n then sustituting these vlues of,,,,,,,,,,, in E., we otin x n x. Then orresoning numeril solution hs een omute y Runge-Kutt fourth-orer metho. From the figure it is ler tht the nlytil solution shows goo oiniene with the numeril solution. In Fig., we hve shown the orresoning hse ifferene etween two osilltors.... x -. - -. - -. Fig. t... x -. - -. - -. Fig. t Fig. Fig. Fig. shows the hrmoni lne solution of x of E. whih is enote y.. n the orresoning numeril solution is enote y. Here ε.,.,.,.,.,.,.,.,.,.. n.,.,.,.,.,.,. In Fig., we oserve the hrmoni lne solution of x of E. whih is enote y. n the orresoning numeril solution is enote y. Here ε.,.,.,.,.,.,.,

Rhmn et l..,.,.,. n.,.,.,.,.,.,.... x -. - -. - -. Fig. t Fig. Fig. shows the freueny ifferene etween x n. n.. x, when ε.,.,... x -. - -. - -. Fig. t... x -. - -. - -. t Fig. Fig. Fig. In Fig. we oserve the hrmoni lne solution of x of E. whih is enote y.. n the orresoning numeril solution is enote y. Here ε.,.,.,.,.,.,.,.,.,.,. n.,.,.,.,.,.,. In Fig. we oserve the hrmoni lne solution of x of E. whih is enote y.. n the orresoning numeril solution is enote y. Here ε.,.,.,.,.,.,.,.,.,.,. n.,.,.,.,

Hrmoni lne Solution.,.,.... x -. - -. - -. Fig. t Fig. Fig. shows the freueny ifferene etween x n x, when ε.,.,. n.. REFERENCES [] Rn, R. H. n Holmes, P. J., ifurtion of erioi motions in two wekly oule vn er Pol osilltors, Int. J. Nonliner Meh., Vol.. -,. [] Neem, I. n Mhome, F. M., roximte first integrls for system of two oule vn er Pol osilltors with liner iffusive ouling, Mthemtil n omuttionl lition, Vol., No.. -,. [] Mrion, J.., Clssil Dynmis of Prtiles n System Sn Diego, C: Hrourt re Jovnovih,. [] Krylov, N.N. n N.N. ogoliuov, Introution to Nonliner Mehnis, Prineton University Press, New Jersey,. [] ogoliuov, N. N. n Yu.. Mitroolskii, symtoti Methos in the Theory of Nonliner Osilltions, Gorn n reh, New York,. [] Nyfeh,.H., Perturtion Methos, J. Wiley, New York,. [] Nyfeh,. H., Introution to Perturtion Tehniue, John Wylie,. [] Lu He-Xing n Li Wei-Dong n itertion metho for single-egree-of-freeom nonliner ynmil eution vne in Virtion Engineering, Vol.,. [] H. Hu n J. H. Tng, lssil itertion roeure vli for ertin strongly nonliner osilltors, J. Soun Vi., Vol.,. -,. [] Mikens, R.E., Osilltion in Plnr Dynmi Systems, Worl Sientifi, Singore,. [] West, J. C., nlytil Tehniues for Nonliner Control Systems, English Univ. Press, Lonon,. [] Mikens, R. E., Comments on the metho of hrmoni lne, J. Soun Vi., Vol.,. -,.

Rhmn et l. [] Mikens, R. E., generliztion of the metho of hrmoni lne, J. Soun Vi., Vol.,. -,. [] Lim, C. W. n. S. Wu, new nlytil roh to the Duffing-hrmoni osilltor, Physis Letters, Vol.,. -,. [] Lim, C. W. n. S. Wu, urte roximte nlytil solutions to nonliner osillting systems with non-rtionl restoring fore, vne in Virtion Engineering, Vol.,. -,. [] Wu,.S., W.P. Sun n C.W. Lim, n nlytil roximte tehniue for lss of strongly nonliner osilltors Int. J. Nonliner Meh., Vol.,. -,. [] He, J. H., Moifie Linste-Poinre methos for some strongly nonliner osilltions-i: exnsion of onstnt, Int. J. Nonliner Meh., Vol.,. -,. [] Lim, C. W. n S. K. Li, urte higher-orer nlytil roximte solutions to nononservtive nonliner osilltors n lition to vn er Pol me osilltors, Int. Journl of Mehnil Sienes, Vol.,. -,. [] Ymgoue, S.. n T. C. Kofne, On the nlytil roximtion of me osilltions of utonomous single egree of freeom osilltors, Int. J. Nonliner Meh., Vol.. -,. [] M. Shmsul lm, M. E. Hue n M.. Hossin new nlytil tehniue to fin erioi solutions of nonliner systems, Int. J. Nonliner Meh., Vol.. -,.