SOLUTION OF DIFFERENTIAL EQUATION FOR THE EULER-BERNOULLI BEAM

Similar documents
APPLICATION OF DIFFERENCE EQUATIONS TO CERTAIN TRIDIAGONAL MATRICES

Numerical Solutions of Fredholm Integral Equations Using Bernstein Polynomials

A GENERAL METHOD FOR SOLVING ORDINARY DIFFERENTIAL EQUATIONS: THE FROBENIUS (OR SERIES) METHOD

Variational Iteration Method for Solving Volterra and Fredholm Integral Equations of the Second Kind

Power Series Solutions to Generalized Abel Integral Equations

sin m a d F m d F m h F dy a dy a D h m h m, D a D a c1cosh c3cos 0

Schrödinger Equation Via Laplace-Beltrami Operator

POWER SERIES SOLUTION OF FIRST ORDER MATRIX DIFFERENTIAL EQUATIONS

SOLUTION OF SYSTEM OF LINEAR EQUATIONS. Lecture 4: (a) Jacobi's method. method (general). (b) Gauss Seidel method.

CALCULATION OF CONVOLUTION PRODUCTS OF PIECEWISE DEFINED FUNCTIONS AND SOME APPLICATIONS

Review of Sections

Certain sufficient conditions on N, p n, q n k summability of orthogonal series

Reduction of Higher Order Linear Ordinary Differential Equations into the Second Order and Integral Evaluation of Exact Solutions

Remarks: (a) The Dirac delta is the function zero on the domain R {0}.

Digital Signal Processing, Fall 2006

Some New Iterative Methods Based on Composite Trapezoidal Rule for Solving Nonlinear Equations

Linear Programming. Preliminaries

Eigenfunction Expansion. For a given function on the internal a x b the eigenfunction expansion of f(x):

A general theory of minimal increments for Hirsch-type indices and applications to the mathematical characterization of Kosmulski-indices

CHAPTER 2: Boundary-Value Problems in Electrostatics: I. Applications of Green s theorem

Chapter 7 Infinite Series

FACULTY OF MATHEMATICAL STUDIES MATHEMATICS FOR PART I ENGINEERING. Lectures

Convergence rates of approximate sums of Riemann integrals

Canonical Form and Separability of PPT States on Multiple Quantum Spaces

0 otherwise. sin( nx)sin( kx) 0 otherwise. cos( nx) sin( kx) dx 0 for all integers n, k.

=> PARALLEL INTERCONNECTION. Basic Properties LTI Systems. The Commutative Property. Convolution. The Commutative Property. The Distributive Property

Graphing Review Part 3: Polynomials

Interpolation. 1. What is interpolation?

Chapter System of Equations

Numerical Solution of Fuzzy Fredholm Integral Equations of the Second Kind using Bernstein Polynomials

METHOD OF FUNDAMENTAL SOLUTIONS FOR HELMHOLTZ EIGENVALUE PROBLEMS IN ELLIPTICAL DOMAINS

National Quali cations AHEXEMPLAR PAPER ONLY

Reversing the Arithmetic mean Geometric mean inequality

Theorem 5.3 (Continued) The Fundamental Theorem of Calculus, Part 2: ab,, then. f x dx F x F b F a. a a. f x dx= f x x

Convergence rates of approximate sums of Riemann integrals

Northwest High School s Algebra 2

Comparison of B-spline Method and Finite Difference Method to Solve BVP of Linear ODEs

Fast Fourier Transform 1) Legendre s Interpolation 2) Vandermonde Matrix 3) Roots of Unity 4) Polynomial Evaluation

Double Sums of Binomial Coefficients

Section 7.3, Systems of Linear Algebraic Equations; Linear Independence, Eigenvalues, Eigenvectors (the variable vector of the system) and

Green s Function Approach to Solve a Nonlinear Second Order Four Point Directional Boundary Value Problem

Prior distributions. July 29, 2002

10.5 Power Series. In this section, we are going to start talking about power series. A power series is a series of the form

BC Calculus Review Sheet

BULLETIN OF MATHEMATICS AND STATISTICS RESEARCH

S(x)along the bar (shear force diagram) by cutting the bar at x and imposing force_y equilibrium.

A GENERALIZATION OF GAUSS THEOREM ON QUADRATIC FORMS

2.1.1 Definition The Z-transform of a sequence x [n] is simply defined as (2.1) X re x k re x k r

Chapter 10: The Z-Transform Adapted from: Lecture notes from MIT, Binghamton University Dr. Hamid R. Rabiee Fall 2013

Taylor Polynomials. The Tangent Line. (a, f (a)) and has the same slope as the curve y = f (x) at that point. It is the best

International Journal of Advancements in Research & Technology, Volume 4, Issue 2, February ISSN

INTEGRATION TECHNIQUES (TRIG, LOG, EXP FUNCTIONS)

[ 20 ] 1. Inequality exists only between two real numbers (not complex numbers). 2. If a be any real number then one and only one of there hold.

National Quali cations SPECIMEN ONLY

Crushed Notes on MATH132: Calculus

Numerical Methods (CENG 2002) CHAPTER -III LINEAR ALGEBRAIC EQUATIONS. In this chapter, we will deal with the case of determining the values of x 1

REVISION SHEET FP1 (AQA) ALGEBRA. E.g., if 2x

INFINITE SERIES. ,... having infinite number of terms is called infinite sequence and its indicated sum, i.e., a 1

Orthogonality, orthogonalization, least squares

Calculus II Homework: The Integral Test and Estimation of Sums Page 1

Chapter 10: The Z-Transform Adapted from: Lecture notes from MIT, Binghamton University Hamid R. Rabiee Arman Sepehr Fall 2010

Options: Calculus. O C.1 PG #2, 3b, 4, 5ace O C.2 PG.24 #1 O D PG.28 #2, 3, 4, 5, 7 O E PG #1, 3, 4, 5 O F PG.

Algebra II, Chapter 7. Homework 12/5/2016. Harding Charter Prep Dr. Michael T. Lewchuk. Section 7.1 nth roots and Rational Exponents

DEPARTMENT OF ELECTRICAL &ELECTRONICS ENGINEERING SIGNALS AND SYSTEMS. Assoc. Prof. Dr. Burak Kelleci. Spring 2018

ALGEBRA. Set of Equations. have no solution 1 b1. Dependent system has infinitely many solutions

Discrete Mathematics I Tutorial 12

MA6351-TRANSFORMS AND PARTIAL DIFFERENTIAL EQUATIONS SUBJECT NOTES. Department of Mathematics FATIMA MICHAEL COLLEGE OF ENGINEERING & TECHNOLOGY

Damped Vibration of a Non-prismatic Beam with a Rotational Spring

Lecture 38 (Trapped Particles) Physics Spring 2018 Douglas Fields

Frequency-domain Characteristics of Discrete-time LTI Systems

EXAM-3 -B2 MATH 261: Elementary Differential Equations MATH 261 FALL 2012 EXAMINATION COVER PAGE Professor Moseley

Indices and Logarithms

MTH 146 Class 16 Notes

In an algebraic expression of the form (1), like terms are terms with the same power of the variables (in this case

EXAM-3 -A4 MATH 261: Elementary Differential Equations MATH 261 FALL 2010 EXAMINATION COVER PAGE Professor Moseley

MATRIX ALGEBRA, Systems Linear Equations

INTEGRAL SOLUTIONS OF THE TERNARY CUBIC EQUATION

Closed Newton-Cotes Integration

General properties of definite integrals

Sequence and Series of Functions

A Constructive Proof of Fundamental Theory for Fuzzy Variable Linear Programming Problems

1. (25 points) Use the limit definition of the definite integral and the sum formulas to compute. [1 x + x2

Chapter 2 Infinite Series Page 1 of 9

Supplemental Handout #1. Orthogonal Functions & Expansions

International Journal of Mathematical Archive-5(1), 2014, Available online through ISSN

Next we encountered the exponent equaled 1, so we take a leap of faith and generalize that for any x (that s not zero),

Section 6.3: Geometric Sequences

( a n ) converges or diverges.

Fourier Series and Applications

Chap8 - Freq 1. Frequency Response

SOME SHARP OSTROWSKI-GRÜSS TYPE INEQUALITIES

Discrete-Time Signals & Systems

Limits and an Introduction to Calculus

y udv uv y v du 7.1 INTEGRATION BY PARTS

Lecture 2. Rational Exponents and Radicals. 36 y. b can be expressed using the. Rational Exponent, thus b. b can be expressed using the

is continuous at x 2 and g(x) 2. Oil spilled from a ruptured tanker spreads in a circle whose area increases at a

M/G /n/0 Erlang queueing system with heterogeneous servers and non-homogeneous customers

Lesson 4 Linear Algebra

Classical Electrodynamics

Important Facts You Need To Know/Review:

Transcription:

Jourl of Applied Mthemtics d Computtiol Mechics () 57-6 SOUION O DIERENIA EQUAION OR HE EUER-ERNOUI EAM Izbel Zmorsk Istitute of Mthemtics Czestochow Uiversit of echolog Częstochow Pold izbel.zmorsk@im.pcz.pl Abstrct. he pper presets the solutio of fourth order differetil equtio with vrious coefficiets occurrig i the vibrtio problem of the Euler-eroulli bem. he cocerig equtio is writte s first order mtri differetil equtio. o solve the equtio the power series method is proposed. Kewords: Euler-eroulli bem power series method mthemticl modellig Itroductio or certi cses of differetil equtios with vrible coefficiets it is possible to determie their ect solutios [-] usig e.g. homotop lsis [] or the Gree's fuctios method []. However i most cses i order to obti solutio it is ecessr to ppl pproimte methods such s fiite differece method [5] the power series method [6] differetil trsformtio method (DM) - improved lor method [5 7 8] or b the use of the grge multiplier formlism [9]. his pper is cotiutio of cosidertio show t [] reltig to the use of mtri d power series methods for solvig ordir differetil equtios. he work presets s emple for proposed method the solutio to the equtio of motio of the o-uiform bem described ccordig to the Euler-eroulli theor b the equtio of the fourth order. ormultio d solutio of the problem sic cocepts of the procedure At the begiig let us recll the procedure schem for fourth order lier differetil equtio: () ( ) ( ) ( ) ( ) f( ) () ()

I. Zmorsk 58 completed b iitil coditios. itroducig fuctios: ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) d ( ) ( ) ( ) ( ) s s s f ssumig tht the fuctios pperig i the equtio () re C clss d c be preseted i the form of power series we c rewrite () s first order mtri differetil equtio []: ( ) ( ) ( ) ( ) () I the bove equtio the followig desigtios were chrged: ( ) ( ) ( ) ( ) ( ) [ ] [ ] ( ) ( ) [ ] [ ] f f ( ) ( ) ( ) ( ) ( ) for d boudr coditio is ( ). A solutio of ihomogeeous equtio () i the power series form c be epressed s sum ( ) ( ) ( ) () where E d re determied from the recursive reltios. he first few vlues of those prmeters re: ( ) ( ) ( ) ()

em s equtio of motio Solutio of differetil equtio for the Euler-eroulli bem 59 Accordig to the Euler-eroulli theor the motio of the bem legth with vrious cross-sectio re A() d momet of ierti I() (ig. ) is described b the prtil differetil equtio: EI u ( ) A( ) u ρ (5) t where u is the fuctio of deflectio ρ is the mss desit d E is the oug s modulus. Equtio (5) is complemeted b pproprite boudr coditios depedig o the method of fiig the eds of the bem. ig.. A sketch of cosidered bem iωt ( ) ( ) Assumig siusoidl rottio of fuctio u t e the equtio of motio c be rewritte s: d d EI d d ( ) A( ) ρω (6) with turl frequec ω. Net let us ssume the cross-sectio re d the momet of ierti i the form of polomils: A α α I ( ) A I( ) (7) where A A() I I() α is proportiolit fctor of the bem s crosssectio d [ ]. After few trsformtios equtio (6) c be epressed s equtio () where: for ( ) ( ) Ω ( ) ( ) ( ) α ( α ). ( ) ( ) 8 ( α ) ( ) (8)

6 I. Zmorsk Prmeter Ω chrcterizes the vibrtio frequec of the bem d is give b formul Ω series: ω ρ A. Note tht fuctios () d () c be preseted s power EI α ( ) ( ) ( ) where (9) Of course we hve to esure the covergece of the series hece: for α > is d for α < is. α α Mtri () occurrig i () is i the form: ( ) () ( ) ( α) ( ) 8( α) Ω d solutio of (6) is give b reltio () for d s i (). I the cse α we ve got vibrtio problem of the Euler-eroulli bem with costt prmeters chrcterizig its phsicl properties. Equtio (5) hs the form: We c ote tht ( ) ( ) ( ) ( ) () Ω () Ω. After coversio ccordig to discussed method mtrices re s follows: O () Ω he solutio of the boudr problem: ( ) ( ) ( ) ( ) ( ) ( ) where: is

Solutio of differetil equtio for the Euler-eroulli bem 6 Ω k Ω Ω Ω Ω k k k k ( Ω ) ( Ω ) k k k k ( Ω ) ( Ω ) ( Ω ) E k... Ω E () he geerl solutio of equtio (6) depeds o four costts which re determied from the boudr coditios. or emple the boudr coditios for ctilever bem re ( ) ( ) ( ) ( ) () he coditios t the bem s boudr gives us d t gives so the iitil coditio to the () is [ ] (5) Coclusios I the preset stud the method reducig the fourth order differetil equtio to the form of the first order mtri equtio ws cosidered o the emple of o-uiform bem s vibrtio equtio. Cosidertios show tht it is reltivel simple method of solvig the boudr problem reduced to form which is es for computer implemettio. If ol the fuctios occurrig i the equtios re the pproprite clss d it s possible to epd them i power series the i the itervl of covergece of these series the preseted method c lso be used to solve more comple problems described b differetil equtios. Refereces [] Elishkoff I. ecquet R. Closed-form solutios for turl frequec for ihomogeeous bems with oe slidig support d the other pied Jourl of Soud d Vibrtio 8() 59-59. [] Che D.-W. Wu J.-S. he ect solutios for the turl frequecies d mode shpes of ouiform bems with multiple sprig - mss sstems Jourl of Soud d Vibrtio 55() 99-. [] Hss H.N. El-wil M.A. A ew techique of usig homotop lsis method for secod order olier differetil equtios Applied Mthemtics d Computtio 9 78-78.

6 I. Zmorsk [] Kukl S. Zmojsk I. Applictio of the Gree s fuctio method i free vibrtio lsis of o-uiform bems Scietific Reserch of the Istitute of Mthemtics d Computer Sciece 5 () 87-9. [5] eh -. Jg M-J. Wg C-C. Alzig the free vibrtios of plte usig fiite differece d differetil trsformtio method Applied Mthemtics d Computtio 6 78 9-5. [6] Qisi M.I. A power series solutio for the o-lier vibrtio of bems Jourl of Soud d Vibrtio 997 99() 587-59. [7] Ozgumus O.O. K M.O. lpwise bedig vibrtio lsis of double tpered rottig Euler-eroulli bem b usig the differetil trsform method Meccic 6 66-67. [8] Mei C. Applictio of differetil trsformtio techique to free vibrtio lsis of cetrifugll stiffeed bem Computers d Structures 8 86 8-8. [9] Cekus D. ree vibrtio of ctilever tpered imosheko bem Scietific Reserch of the Istitute of Mthemtics d Computer Sciece () -7. [] Kukl S. Zmorsk I. Power series solutio of first order mtri differetil equtios Jourl of Applied Mthemtics d Computtiol Mechics () -8.