RAYLEIGH'S METHOD Revision D

Similar documents
Engineering Mechanics Dynamics & Vibrations. Engineering Mechanics Dynamics & Vibrations Plane Motion of a Rigid Body: Equations of Motion

Exercise 8 CRITICAL SPEEDS OF THE ROTATING SHAFT

The Binomial Multi- Section Transformer

BENDING FREQUENCIES OF BEAMS, RODS, AND PIPES Revision S

The Binomial Multi-Section Transformer

Mechanical Vibrations

5.6 Binomial Multi-section Matching Transformer

5.6 Binomial Multi-section Matching Transformer

DYNAMICS VECTOR MECHANICS FOR ENGINEERS: Mechanical Vibrations. Seventh Edition CHAPTER. Ferdinand P. Beer E. Russell Johnston, Jr.

The state space model needs 5 parameters, so it is not as convenient to use in this control study.

EXPERIMENT OF SIMPLE VIBRATION

The Differential Transform Method for Solving Volterra s Population Model

(s)h(s) = K( s + 8 ) = 5 and one finite zero is located at z 1

Application of Homotopy Analysis Method for Solving various types of Problems of Ordinary Differential Equations

EN40: Dynamics and Vibrations. Final Examination Friday May : 2pm-5pm

DETERMINATION OF NATURAL FREQUENCY AND DAMPING RATIO

Problem 1. Problem Engineering Dynamics Problem Set 9--Solution. Find the equation of motion for the system shown with respect to:

Perturbation Theory, Zeeman Effect, Stark Effect

Lecture 11. Solution of Nonlinear Equations - III

EFFECTIVE MODAL MASS & MODAL PARTICIPATION FACTORS Revision I

Project 10.3 Vibrating Beams and Diving Boards

Contents Two Sample t Tests Two Sample t Tests

Studying Interaction of Cotton-Raw Material with Working Bodies of Cotton-Cleaning Machines

PARTIAL DIFFERENTIAL EQUATIONS SEPARATION OF VARIABLES

Discrete population models

Castiel, Supernatural, Season 6, Episode 18

Transfer Function Analysis

Lecture 10: Bounded Linear Operators and Orthogonality in Hilbert Spaces

The Hypergeometric Coupon Collection Problem and its Dual

Complete Solutions to Supplementary Exercises on Infinite Series

School of Mechanical Engineering Purdue University. ME375 Transfer Functions - 1

Chapter 2: Numerical Methods

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

Energy Identities of ADI-FDTD Method with Periodic Structure

Acoustic Field inside a Rigid Cylinder with a Point Source

Application 10.5D Spherical Harmonic Waves

x !1! + 1!2!

fiziks Forum for CSIR-UGC JRF/NET, GATE, IIT-JAM, GRE in PHYSICAL SCIENCES

Chapter 4. Fourier Series

Supplementary Information

Answer Key, Problem Set 1, Written

(5x 7) is. 63(5x 7) 42(5x 7) 50(5x 7) BUSINESS MATHEMATICS (Three hours and a quarter)

A New Type of q-szász-mirakjan Operators

Fourier Series and the Wave Equation

A PROBABILITY PROBLEM

FREE VIBRATIONS OF SIMPLY SUPPORTED BEAMS USING FOURIER SERIES

(4 pts.) (4 pts.) (4 pts.) b) y(x,t) = 1/(ax 2 +b) This function has no time dependence, so cannot be a wave.

Conversion Factors. Dynamic Measurements. Displacement, Velocity, Acceleration Relationships. Acceleration. Pressure. Sinusoids (only for sinusoids)

Lebesgue Constant Minimizing Bivariate Barycentric Rational Interpolation

Mathematics Extension 1

A New Family of Multipoint Iterative Methods for Finding Multiple Roots of Nonlinear Equations

FREE VIBRATION RESPONSE OF A SYSTEM WITH COULOMB DAMPING

2.2 Limits Involving Infinity AP Calculus

A Steady State Heat Conduction Problem in. a Thick Annular Disc Due to Arbitrary. Axisymmetric Heat Flux

MATH 10550, EXAM 3 SOLUTIONS

Bertrand s postulate Chapter 2

NONLOCAL THEORY OF ERINGEN

PHY 181 Test 2 Practice Solutions Spring 2013

19.1 The dictionary problem

After the completion of this section the student

Quiz. Use either the RATIO or ROOT TEST to determine whether the series is convergent or not.

Jacobi symbols. p 1. Note: The Jacobi symbol does not necessarily distinguish between quadratic residues and nonresidues. That is, we could have ( a

Recurrence Relations

and then substitute this into the second equation to get 5(11 4 y) 3y

Implicit Splitting Finite Difference Scheme for Multi-dimensional Wave Simulation

J 10 J W W W W

The driven Rayleigh-van der Pol oscillator

Solution: APPM 1360 Final Spring 2013

Lecture 20 - Wave Propagation Response

Numerical Solution of Nonlinear Integro-Differential Equations with Initial Conditions by Bernstein Operational Matrix of Derivative

Orthogonal Function Solution of Differential Equations

NICK DUFRESNE. 1 1 p(x). To determine some formulas for the generating function of the Schröder numbers, r(x) = a(x) =

Wavelet Transform Theory. Prof. Mark Fowler Department of Electrical Engineering State University of New York at Binghamton

Chapter 10 Partial Differential Equations and Fourier Series

Homework 6: Forced Vibrations Due Friday April 6, 2018

Ma 530 Introduction to Power Series

Flight and Orbital Mechanics. Exams

Bernoulli Polynomials Talks given at LSBU, October and November 2015 Tony Forbes

Most text will write ordinary derivatives using either Leibniz notation 2 3. y + 5y= e and y y. xx tt t

Paper-II Chapter- Damped vibration

Null Spaces, Column. Transformations. Remarks. Definition ( ) ( ) Null A is the set of all x's in R that go to 0 in R 3

Adv. Studies Theor. Phys., Vol. 5, 2011, no. 4, A Numerical Computation of Non-Dimensional. Form of a Mathematical Model of Soil Salinity

CHAPTER 6 RESISTANCE FACTOR FOR THE DESIGN OF COMPOSITE SLABS

Calculus 2 - D. Yuen Final Exam Review (Version 11/22/2017. Please report any possible typos.)

Note that the argument inside the second square root is always positive since R L > Z 0. The series reactance can be found as

3. Supppose the amount of information available on the web is multiplied by 27 every year. How much information will be available a.

Integrals of Functions of Several Variables

Numerical Methods in Fourier Series Applications

Define a Markov chain on {1,..., 6} with transition probability matrix P =

Chapter 5 Problem Solutions


Results of Final Exam

Principle Of Superposition

A Generalization of Ince s Equation

Course Outline. 2. Motion of systems that can be idealized as particles

Finite element analysis of nonlinear structures with Newmark method

Statistics and Data Analysis in MATLAB Kendrick Kay, February 28, Lecture 4: Model fitting

PHYS 1443 Section 003 Lecture #22

Automated Proofs for Some Stirling Number Identities

Transcription:

RAYGH'S METHOD Revisio D B To Irvie Eail: toirvie@aol.co Noveber 5, Itroductio Daic sstes ca be characterized i ters of oe or ore atural frequecies. The atural frequec is the frequec at which the sste would vibrate if it were give a iitial disturbace ad the allowed to vibrate freel. There are a available ethods for deteriig the atural frequec. eaples are Soe. Newto s aw of Motio. Raleigh s Method. Eerg Method. agrage s Equatio Note that the Raleigh, Eerg, ad agrage ethods are closel related. Soe of these ethods directl ield the atural frequec. Others ield a goverig equatio of otio, fro which the atural frequec a be deteried. This tutorial focuses o Raleigh s ethod, which ields the atural frequec. Raleigh's ethod requires a assued displaceet fuctio. The ethod thus reduces the daic sste to a sigle-degree-of-freedo sste. Furtherore, the assued displaceet fuctio itroduces additioal costraits which icrease the stiffess of the sste. Thus, Raleigh's ethod ields a upper liit of the true fudaetal frequec. Defiitio Raleigh's ethod ca be suarized as where a PE a KE = total eerg of the sste () KE = kietic eerg PE = potetial eerg

Note that potetial eerg is also referred to as strai eerg for the case of certai sstes, such as beas. Equatio () ca ol be satisfied if the sste is vibratig at its atural frequec. Refereces. T. Irvie, Natural Frequecies of Bea Bedig Modes, Revisio E, Vibratiodata, 99.

APPENDIX A Pedulu Eaple Cosider a coservative sste. A eaple is the pedulu show i Figure A-. g Figure A-. The kietic eerg becoes zero whe the pedulu reaches its aiu agular displaceet. The kietic eerg reaches its aiu value whe the pedulu passes through =, which is also the "static equilibriu poit." O the other had, the potetial eerg reaches its aiu level as the pedulu reaches its aiu agular displaceet. The potetial eerg reaches its iiu value whe the pedulu is at its static equilibriu poit. For siplicit, the potetial eerg ca be cosidered as zero at the static equilibriu poit. et = pedulu ass = legth = agular displaceet Assue a sall agular displaceet.

The potetial eerg is PE g( cos ) (A-) The kietic eerg is Assue a displaceet equatio of KE ( ) (A-) ( t) si t (A-) a (A-) The velocit equatio is ( t) cos t (A-5) a (A-6) PE g( cos ) a Cosider the epasio (A-7) cos (A-)!! Cosider the aiu potetial eerg. Substitute equatio (A-) ito (A-7). PE a g (A-9)!! PE a g (A-)!!

PE a g (A-) Cosider the aiu kietic eerg. Substitute equatio (A-6) ito (A-). a (A-) KE g (A-) Siplifig, g (A-) The pedulu atural frequec is thus g (A-5) 5

APPENDIX B Catilever Bea with Ed Mass Cosider a ass outed o the ed of a catilever bea, as show i Figure B-. Assue that the ed-ass is uch greater tha the ass of the bea. Figure B-. E is the odulus of elasticit I is the area oet of iertia is the legth g is gravit is the ass is the displaceet The static stiffess at the ed of the bea is k (B-) Equatio (B-) is derived i Referece. The potetial eerg is PE (B-) The kietic eerg is KE (B-) Assue a ed displaceet of Asi t (B-) 6

The correspodig velocit is Acos t (B-5) The aiu displaceet is A. The aiu velocit is A. Thus, KEa A (B-6) PEa A (B-7) Applig Raleigh s ethod, A A (B-) The atural frequec of the ed ass supported b the catilever bea is thus (B-9) (B-) 7

APPENDIX C Catilever Bea with Iteral Distributed Mass Cosider a catilever bea with ass per legth. Assue that the bea has a uifor cross sectio. Deterie the atural frequec., Figure C-. The goverig differetial equatio is t (C-) The boudar coditios at the fied ed = are () = (zero displaceet) (C-) d (zero slope) (C-) The boudar coditios at the free ed = are d (zero bedig oet) (C-) d (zero shear force) (C-5) Propose a quarter cosie wave solutio.

9 o ( ) cos (C-6) d o si (C-7) d o cos (C-) d o si (C-9) The proposed solutio eets all of the boudar coditios epect for the zero shear force at the right ed. The proposed solutio is accepted as a approiate solutio for the deflectio shape, despite oe deficiec. Agai, Raleigh s ethod is used to fid the atural frequec. The total potetial eerg ad the total kietic eerg ust be deteried. The total potetial eerg P i the bea is P d (C-) B substitutio, P o cos (C-) P o cos (C-) P o cos (C-)

P o si (C-) P o (C-5) o 6 P (C-6) The total kietic eerg T is T (C-7) T o cos (C-) T o cos cos (C-9) T o cos cos (C-) T o cos cos (C-) T o cos cos (C-) T o si si (C-) T o (C-)

T o (C-5) Now equate the potetial ad the kietic eerg ters. 6 o o (C-6) 6 (C-7) 6 (C-) 6 / (C-9) f 6 / (C-) f 6 / (C-)

/ f (C-) f 66. (C-)

APPENDIX D Catilever Bea with Iteral Distributed Mass ad Ed Mass Figure D-. The followig is based o a quarter cosie wave solutio. The kietic eerg is o o T (D-) T o (D-) The potetial eerg is o 6 P (D-) Now equate the potetial ad the kietic eerg ters. o o 6 (D-),

6 (D-5) 6 (D-6) / 6 (D-7) / (D-) / f (D-9).6. f (D-)

APPENDIX E Raleigh s Quotiet Raleigh s ethod ca also be applied to ulti-degree-freedo-sstes, as follows. X T K X (E-) X T MX where K is the stiffess atri M is the ass atri X is a assued ode shape with arbitrar scale Equatio (E-) is essetiall a uerical approiatio. It overestiates the true fudaetal frequec. Thus, it should be used i a trial-ad-error aer. Note that the uerator i equatio (E-) is equal to twice the potetial eerg. The deoiator is equal to twice the kietic eerg if first ultiplied b. As a eaple, cosider the sste defied i Figure E- ad Table E-. 5

k k k Figure E-. Table E-. Paraeters Variable k k k. kg. kg Value N/ N/ N/ The hoogeeous equatio of otio is k k k k k k (E-) 6

The ass atri is M kg (E-) The stiffess atri is K N / 5 (E-) The atural frequecies ca be coputed usig the eigevalue ethod. The eigevalues are the roots of the followig equatio. det K M (E-5) Equatio (E-5) ca be solved eactl for sstes with up to four degrees-of-freedo. The first atural frequec is thus. rad /sec (E-6) The et task is to test the Raleigh quotiet ethod. Several cadidate ode shape are evaluated as show i Table E-. Note that X (E-7) 7

Table E-. Raleigh Quotiet Trials X.5.5 (rad/sec).6.6.6 6.5 5.77 Agai, the Raleigh quotiet overestiates the true fudaetal frequec. Thus, the best estiate for the atural frequec after five trials is.6 rad /sec (estiate) (E-) The estiated value is.% higher tha the eact value. The Raleigh quotiet ethod thus gives ver good results for this eaple.