Introduction to Mechanical Vibrations

Similar documents
Energy Density / Energy Flux / Total Energy in 1D. Key Mathematics: density, flux, and the continuity equation.

Ideal Amplifier/Attenuator. Memoryless. where k is some real constant. Integrator. System with memory

BE.430 Tutorial: Linear Operator Theory and Eigenfunction Expansion

Problems and Solutions for Section 3.2 (3.15 through 3.25)

In this section we will study periodic signals in terms of their frequency f t is said to be periodic if (4.1)

Electrical Engineering Department Network Lab.

State-Space Model. In general, the dynamic equations of a lumped-parameter continuous system may be represented by

Optimization of Rotating Machines Vibrations Limits by the Spring - Mass System Analysis

Section 8 Convolution and Deconvolution

EXISTENCE THEORY OF RANDOM DIFFERENTIAL EQUATIONS D. S. Palimkar

CHAPTER 2 TORSIONAL VIBRATIONS

BEST LINEAR FORECASTS VS. BEST POSSIBLE FORECASTS

On Acoustic Radiation by Rotating Dipole Sources in Frequency Domain

ODEs II, Supplement to Lectures 6 & 7: The Jordan Normal Form: Solving Autonomous, Homogeneous Linear Systems. April 2, 2003

PLASTIC BUCKLING OF SSSS THIN RECTANGULAR PLATES SUBJECTED TO UNIAXIAL COMPRESSION USING TAYLOR-MACLAURIN SHAPE FUNCTION

Comparison between Fourier and Corrected Fourier Series Methods

Walter D. Pilkey. Dmitry V. Balandin. Nikolai N. Bolotnik. Sergey V. Purtsezov. Introduction. Nomenclature

1 Notes on Little s Law (l = λw)

ES 330 Electronics II Homework 03 (Fall 2017 Due Wednesday, September 20, 2017)

Big O Notation for Time Complexity of Algorithms

E will be denoted by n

λiv Av = 0 or ( λi Av ) = 0. In order for a vector v to be an eigenvector, it must be in the kernel of λi

Fresnel Dragging Explained

Solutions Manual 4.1. nonlinear. 4.2 The Fourier Series is: and the fundamental frequency is ω 2π

Dynamic Model and Analysis of Nonlinear Vibration Characteristic of a Curve-Face Gear Drive Cai, Z. Lin, C. Zhiqin Cai Chao Lin*

Introduction to Mobile Robotics Mapping with Known Poses

Bernstein Direct Method for Solving. Variational Problems

ECE-314 Fall 2012 Review Questions

[ m] x = 0.25cos 20 t sin 20 t m

Calculus BC 2015 Scoring Guidelines

th m m m m central moment : E[( X X) ] ( X X) ( x X) f ( x)

Chemistry 1B, Fall 2016 Topics 21-22

ME 501A Seminar in Engineering Analysis Page 1

Harmonic excitation (damped)

Vibration 2-1 MENG331

An Exact Solution for the Free Vibration Analysis of Timoshenko Beams

C(p, ) 13 N. Nuclear reactions generate energy create new isotopes and elements. Notation for stellar rates: p 12

N! AND THE GAMMA FUNCTION

Economics 8723 Macroeconomic Theory Problem Set 2 Professor Sanjay Chugh Spring 2017

Math 6710, Fall 2016 Final Exam Solutions

Clock Skew and Signal Representation

STK4080/9080 Survival and event history analysis

Applications of force vibration. Rotating unbalance Base excitation Vibration measurement devices

Notes 03 largely plagiarized by %khc

Using Linnik's Identity to Approximate the Prime Counting Function with the Logarithmic Integral

Spectral Simulation of Turbulence. and Tracking of Small Particles

Sliding Mode Control for Robust Stabilization of Uncertain Input-Delay Systems

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

Chapter 3 Moments of a Distribution

CLOSED FORM EVALUATION OF RESTRICTED SUMS CONTAINING SQUARES OF FIBONOMIAL COEFFICIENTS

Lecture 9: Polynomial Approximations

A quadratic convergence method for the management equilibrium model

B. Maddah INDE 504 Simulation 09/02/17

1. Solve by the method of undetermined coefficients and by the method of variation of parameters. (4)

12 Getting Started With Fourier Analysis

Optimum design of complementary transient experiments for estimating thermal properties

Dynamic h-index: the Hirsch index in function of time

Calculus Limits. Limit of a function.. 1. One-Sided Limits...1. Infinite limits 2. Vertical Asymptotes...3. Calculating Limits Using the Limit Laws.

SHOCK AND VIBRATION RESPONSE SPECTRA COURSE Unit 21 Base Excitation Shock: Classical Pulse

Paper 3A3 The Equations of Fluid Flow and Their Numerical Solution Handout 1

Moment Generating Function

Actuarial Society of India

Extremal graph theory II: K t and K t,t

WILD Estimating Abundance for Closed Populations with Mark-Recapture Methods

The Recovery of Rockets.

EXAMPLE SHEET B (Partial Differential Equations) SPRING 2014

SUMMATION OF INFINITE SERIES REVISITED

The Eigen Function of Linear Systems

5.74 Introductory Quantum Mechanics II

ME 3210 Mechatronics II Laboratory Lab 6: Second-Order Dynamic Response

CHARACTERIZATIONS OF THE NON-UNIFORM IN TIME ISS PROPERTY AND APPLICATIONS

Application examples and potential for explicit finite elements in the analysis of friction brakes

Three Point Bending Analysis of a Mobile Phone Using LS-DYNA Explicit Integration Method

Practice Problems - Week #4 Higher-Order DEs, Applications Solutions

Lecture 23 Damped Motion

10.3 Autocorrelation Function of Ergodic RP 10.4 Power Spectral Density of Ergodic RP 10.5 Normal RP (Gaussian RP)

CHAPTER 2. Problem 2.1. Given: m k = k 1. Determine the weight of the table sec (b)

Chemical Engineering 374

MALAYSIAN JOURNAL OF MATHEMATICAL SCIENCES. Boundary Value Problem for the Higher Order Equation with Fractional Derivative

Boundary-to-Displacement Asymptotic Gains for Wave Systems With Kelvin-Voigt Damping

Let s express the absorption of radiation by dipoles as a dipole correlation function.

ENGINEERING MECHANICS

Chapter 2. Foundations of quantum mechanics

INVESTMENT PROJECT EFFICIENCY EVALUATION

Transverse Vibrations of Elastic Thin Beam Resting on Variable Elastic Foundations and Subjected to Traveling Distributed Forces.

An interesting result about subset sums. Nitu Kitchloo. Lior Pachter. November 27, Abstract

Mechanical Vibrations

If boundary values are necessary, they are called mixed initial-boundary value problems. Again, the simplest prototypes of these IV problems are:

The Importance of Ordering the Number of Lattice Points Inside a Rational Polyhedron Using Generating Functions

Stationarity and Error Correction

David Randall. ( )e ikx. k = u x,t. u( x,t)e ikx dx L. x L /2. Recall that the proof of (1) and (2) involves use of the orthogonality condition.

DETERMINATION OF NATURAL FREQUENCY AND DAMPING RATIO

14.02 Principles of Macroeconomics Fall 2005

Solutions to selected problems from the midterm exam Math 222 Winter 2015

Introduction to Numerical Analysis. In this lesson you will be taken through a pair of techniques that will be used to solve the equations of.

Basic Magnetism Thorsten Glaser, University of Münster. Basic Magnetism. 1. Paramagnetism and Diamagnetism (macroscopic)

CSE 5311 Notes 12: Matrices

Robust Dynamic Output Feedback Second-Order Sliding Mode Controller for Uncertain Systems

Pure Math 30: Explained!

Online Supplement to Reactive Tabu Search in a Team-Learning Problem

Transcription:

CHAPTER 1 Iroducio o Mechaical Vibraios Vibraio is he oio of a paricle or a body or syse of coeced bodies displaced fro a posiio of equilibriu. Mos vibraios are udesirable i achies ad srucures because hey produce icreased sresses, eergy losses, cause added wear, icrease bearig loads, iduce faigue, creae passeger discofor i vehicles, ad absorb eergy fro he syse. Roaig achie pars eed careful balacig i order o preve daage fro vibraios. Vibraio occurs whe a syse is displaced fro a posiio of sable equilibriu. The syse eds o reur o his equilibriu posiio uder he acio of resorig forces (such as he elasic forces, as for a ass aached o a sprig, or graviaioal forces, as for a siple pedulu). The syse keeps ovig back ad forh across is posiio of equilibriu. A syse is a cobiaio of elees ieded o ac ogeher o accoplish a objecive. For exaple, a auoobile is a syse whose elees are he wheels, suspesio, car body, ad so forh. A saic elee is oe whose oupu a ay give ie depeds oly o he ipu a ha ie while a dyaic elee is oe whose prese oupu depeds o pas ipus. I he sae way we also speak of saic ad dyaic syses. A saic syse coais all elees while a dyaic syse coais a leas oe dyaic elee. A physical syse udergoig a ie-varyig ierchage or dissipaio of eergy aog or wihi is eleeary sorage or dissipaive devices is said o be i a dyaic sae. All of he elees i geeral are called passive, i.e., hey are icapable of geeraig e eergy. A dyaic syse coposed of a fiie uber of sorage elees is said o be luped or discree, while a syse coaiig elees, which are dese i physical space, is called coiuous. The aalyical descripio of he dyaics of he discree case is a se of ordiary differeial equaios, while for he coiuous case i is a se of parial differeial equaios. The aalyical foraio of a dyaic syse depeds upo he kieaic or geoeric cosrais ad he physical laws goverig he behaviour of he syse. 1.1 CLASSIFICATION OF VIBRATIONS Vibraios ca be classified io hree caegories: free, forced, ad self-excied. Free vibraio of a syse is vibraio ha occurs i he absece of exeral force. A exeral force ha acs o he syse causes forced vibraios. I his case, he exciig force coiuously supplies eergy o he syse. Forced vibraios ay be eiher deeriisic or rado (see Fig. 1.1). Selfexcied vibraios are periodic ad deeriisic oscillaios. Uder cerai codiios, he 1

2 SOLVING VIBRATION ANALYSIS PROBLEMS USING MATLAB equilibriu sae i such a vibraio syse becoes usable, ad ay disurbace causes he perurbaios o grow uil soe effec liis ay furher growh. I coras o forced vibraios, he exciig force is idepede of he vibraios ad ca sill persis eve whe he syse is preveed fro vibraig. x x = x() Fig. 1.1(a) A deeriisic (periodic) exciaio. x Fig. 1.1(b) Rado exciaio. 1.2 ELEMENTARY PARTS OF VIBRATING SYSTEMS I geeral, a vibraig syse cosiss of a sprig (a eas for sorig poeial eergy), a ass or ieria (a eas for sorig kieic eergy), ad a daper (a eas by which eergy is gradually los) as show i Fig. 1.2. A udaped vibraig syse ivolves he rasfer of is poeial eergy o kieic eergy ad kieic eergy o poeial eergy, aleraively. I a daped vibraig syse, soe eergy is dissipaed i each cycle of vibraio ad should be replaced by a exeral source if a seady sae of vibraio is o be aiaied.

INTRODUCTION TO MECHANICAL VIBRATIONS 3 Sprig k Daper c Saic equilibriu posiio 0 Mass Exciaio force F() Displacee x Fig. 1.2 Eleeary pars of vibraig syses. 1.3 PERIODIC MOTION Whe he oio is repeaed i equal iervals of ie, i is kow as periodic oio. Siple haroic oio is he siples for of periodic oio. If x() represes he displacee of a ass i a vibraory syse, he oio ca be expressed by he equaio x = A cos ω = A cos 2π τ where A is he apliude of oscillaio easured fro he equilibriu posiio of he ass. The repeiio ie τ is called he period of he oscillaio, ad is reciprocal, f = 1, is called he τ frequecy. Ay periodic oio saisfies he relaioship x () = x ( + τ) 2π Tha is Period τ = s/cycle ω Frequecy f = 1 τ = ω cycles/s, or Hz 2π ω is called he circular frequecy easured i rad/sec. The velociy ad acceleraio of a haroic displacee are also haroic of he sae frequecy, bu lead he displacee by π/2 ad π radias, respecively. Whe he acceleraio X && of a paricle wih reciliear oio is always proporioal o is displacee fro a fixed poi o he pah ad is direced owards he fixed poi, he paricle is said o have siple haroic oio. The oio of ay vibraig syses i geeral is o haroic. I ay cases he vibraios are periodic as i he ipac force geeraed by a forgig haer. If x() is a periodic fucio wih period τ, is Fourier series represeaio is give by x() = a 0 2 + (a cos ω + b si ω) = 1 where ω = 2π/τ is he fudaeal frequecy ad a 0, a 1, a 2,, b 1, b 2, are cosa coefficies, which are give by: zτ a 0 = 2 τ 0 x() d

4 SOLVING VIBRATION ANALYSIS PROBLEMS USING MATLAB zτ a = 2 x() cos ω d τ 0 b = 2 zτ x() si ω d τ 0 The expoeial for of x() is give by: x() = iω ce = The Fourier coefficies c ca be deeried, usig c = 1 z τ (x) e iω d τ 0 The haroic fucios a cos ω or b si ω are kow as he haroics of order of he periodic fucio x(). The haroic of order has a period τ/. These haroics ca be ploed as verical lies i a diagra of apliude (a ad b ) versus frequecy (ω) ad is called frequecy specru. 1.4 DISCRETE AND CONTINUOUS SYSTEMS Mos of he echaical ad srucural syses ca be described usig a fiie uber of degrees of freedo. However, here are soe syses, especially hose iclude coiuous elasic ebers, have a ifiie uber of degree of freedo. Mos echaical ad srucural syses have elasic (deforable) elees or copoes as ebers ad hece have a ifiie uber of degrees of freedo. Syses which have a fiie uber of degrees of freedo are kow as discree or luped paraeer syses, ad hose syses wih a ifiie uber of degrees of freedo are called coiuous or disribued syses. 1.5 VIBRATION ANALYSIS The oupus of a vibraig syse, i geeral, deped upo he iiial codiios, ad exeral exciaios. The vibraio aalysis of a physical syse ay be suarised by he four seps: 1. Maheaical Modellig of a Physical Syse 2. Forulaio of Goverig Equaios 3. Maheaical Soluio of he Goverig Equaios 1. Maheaical odellig of a physical syse The purpose of he aheaical odellig is o deerie he exisece ad aure of he syse, is feaures ad aspecs, ad he physical elees or copoes ivolved i he physical syse. Necessary assupios are ade o siplify he odellig. Iplici assupios are used ha iclude: (a) A physical syse ca be reaed as a coiuous piece of aer (b) Newo s laws of oio ca be applied by assuig ha he earh is a ieral frae (c) Igore or eglec he relaivisic effecs All copoes or elees of he physical syse are liear. The resulig aheaical odel ay be liear or o-liear, depedig o he give physical syse. Geerally speakig, all physical syses exhibi o-liear behaviour. Accurae aheaical odel-

INTRODUCTION TO MECHANICAL VIBRATIONS 5 lig of ay physical syse will lead o o-liear differeial equaios goverig he behaviour of he syse. Ofe, hese o-liear differeial equaios have eiher o soluio or difficul o fid a soluio. Assupios are ade o liearise a syse, which peris quick soluios for pracical purposes. The advaages of liear odels are he followig: (1) heir respose is proporioal o ipu (2) superposiio is applicable (3) hey closely approxiae he behaviour of ay dyaic syses (4) heir respose characerisics ca be obaied fro he for of syse equaios wihou a deailed soluio (5) a closed-for soluio is ofe possible (6) uerical aalysis echiques are well developed, ad (7) hey serve as a basis for udersadig ore coplex o-liear syse behaviours. I should, however, be oed ha i os o-liear probles i is o possible o obai closed-for aalyic soluios for he equaios of oio. Therefore, a copuer siulaio is ofe used for he respose aalysis. Whe aalysig he resuls obaied fro he aheaical odel, oe should realise ha he aheaical odel is oly a approxiaio o he rue or real physical syse ad herefore he acual behaviour of he syse ay be differe. 2. Forulaio of goverig equaios Oce he aheaical odel is developed, we ca apply he basic laws of aure ad he priciples of dyaics ad obai he differeial equaios ha gover he behaviour of he syse. A basic law of aure is a physical law ha is applicable o all physical syses irrespecive of he aerial fro which he syse is cosruced. Differe aerials behave differely uder differe operaig codiios. Cosiuive equaios provide iforaio abou he aerials of which a syse is ade. Applicaio of geoeric cosrais such as he kieaic relaioship bewee displacee, velociy, ad acceleraio is ofe ecessary o coplee he aheaical odellig of he physical syse. The applicaio of geoeric cosrais is ecessary i order o forulae he required boudary ad/or iiial codiios. The resulig aheaical odel ay be liear or o-liear, depedig upo he behaviour of he elees or copoes of he dyaic syse. 3. Maheaical soluio of he goverig equaios The aheaical odellig of a physical vibraig syse resuls i he forulaio of he goverig equaios of oio. Maheaical odellig of ypical syses leads o a syse of differeial equaios of oio. The goverig equaios of oio of a syse are solved o fid he respose of he syse. There are ay echiques available for fidig he soluio, aely, he sadard ehods for he soluio of ordiary differeial equaios, Laplace rasforaio ehods, arix ehods, ad uerical ehods. I geeral, exac aalyical soluios are available for ay liear dyaic syses, bu for oly a few oliear syses. Of course, exac aalyical soluios are always preferable o uerical or approxiae soluios. 4. Physical ierpreaio of he resuls The soluio of he goverig equaios of oio for he physical syse geerally gives he perforace. To verify he validiy of he odel, he prediced perforace is copared wih he experieal resuls. The odel ay have o be refied or a ew odel is developed ad a ew predicio copared wih he experieal resuls. Physical ierprea-

6 SOLVING VIBRATION ANALYSIS PROBLEMS USING MATLAB io of he resuls is a ipora ad fial sep i he aalysis procedure. I soe siuaios, his ay ivolve (a) drawig geeral ifereces fro he aheaical soluio, (b) develope of desig curves, (c) arrive a a siple ariheic o arrive a a coclusio (for a ypical or specific proble), ad (d) recoedaios regardig he sigificace of he resuls ad ay chages (if ay) required or desirable i he syse ivolved. 1.5.1 COMPONENTS OF VIBRATING SYSTEMS (a) Siffess elees Soe ies i requires fidig ou he equivale sprig siffess values whe a coiuous syse is aached o a discree syse or whe here are a uber of sprig elees i he syse. Siffess of coiuous elasic elees such as rods, beas, ad shafs, which produce resorig elasic forces, is obaied fro deflecio cosideraios. The siffess coefficie of he rod (Fig. 1.3) is give by k = EA l The cailever bea (Fig.1.4) siffess is k = 3 EI l3 The orsioal siffess of he shaf (Fig.1.5) is K = GJ l E,A, l EA k= l u F Fig.1.3 Logiudial vibraio of rods. E,I, l F v k= 3EI l 3 Fig.1.4 Trasverse vibraio of cailever beas.

INTRODUCTION TO MECHANICAL VIBRATIONS 7 G,J, l T k= GJ l Fig. 1.5 Torsioal syse. Whe here are several sprigs arraged i parallel as show i Fig. 1.6, he equivale sprig cosa is give by algebraic su of he siffess of idividual sprigs. Maheaically, k eq = i = 1 k i k 1 k 2 k Fig. 1.6 Sprigs i parallel. Whe he sprigs are arraged i series as show i Fig. 1.7, he sae force is developed i each sprig ad is equal o he force acig o he ass. k 1 k 2 k 3 k Fig. 1.7 Sprigs i series. The equivale siffess k eq is give by: 1 1/k eq = 1 i = 1 k i Hece, whe elasic elees are i series, he reciprocal of he equivale elasic cosa is equal o he reciprocals of he elasic cosas of he elees i he origial syse. (b) Mass or ieria elees The ass or ieria elee is assued o be a rigid body. Oce he aheaical odel of he physical vibraig syse is developed, he ass or ieria elees of he syse ca be easily ideified.

8 SOLVING VIBRATION ANALYSIS PROBLEMS USING MATLAB (c) Dapig elees I real echaical syses, here is always eergy dissipaio i oe for or aoher. The process of eergy dissipaio is referred o i he sudy of vibraio as dapig. A daper is cosidered o have eiher ass or elasiciy. The hree ai fors of dapig are viscous dapig, Coulob or dry-fricio dapig, ad hyseresis dapig. The os coo ype of eergy-dissipaig elee used i vibraios sudy is he viscous daper, which is also referred o as a dashpo. I viscous dapig, he dapig force is proporioal o he velociy of he body. Coulob or dry-fricio dapig occurs whe slidig coac ha exiss bewee surfaces i coac are dry or have isufficie lubricaio. I his case, he dapig force is cosa i agiude bu opposie i direcio o ha of he oio. I dry-fricio dapig eergy is dissipaed as hea. Solid aerials are o perfecly elasic ad whe hey are defored, eergy is absorbed ad dissipaed by he aerial. The effec is due o he ieral fricio due o he relaive oio bewee he ieral plaes of he aerial durig he deforaio process. Such aerials are kow as visco-elasic solids ad he ype of dapig which hey exhibi is called as srucural or hysereic dapig, or aerial or solid dapig. I ay pracical applicaios, several dashpos are used i cobiaio. I is quie possible o replace hese cobiaios of dashpos by a sigle dashpo of a equivale dapig coefficie so ha he behaviour of he syse wih he equivale dashpo is cosidered ideical o he behaviour of he acual syse. 1.6 FREE VIBRATION OF SINGLE DEGREE OF FREEDOM SYSTEMS The os basic echaical syse is he sigle-degree-of-freedo syse, which is characerized by he fac ha is oio is described by a sigle variable or coordiaes. Such a odel is ofe used as a approxiaio for a geerally ore coplex syse. Exciaios ca be broadly divided io wo ypes, iiial exciaios ad exerally applied forces. The behavior of a syse characerized by he oio caused by hese exciaios is called as he syse respose. The oio is geerally described by displacees. 1.6.1 FREE VIBRATION OF AN UNDAMPED TRANSLATIONAL SYSTEM The siples odel of a vibraig echaical syse cosiss of a sigle ass elee which is coeced o a rigid suppor hrough a liearly elasic assless sprig as show i Fig. 1.8. The ass is cosraied o ove oly i he verical direcio. The oio of he syse is described by a sigle coordiae x() ad hece i has oe degree of freedo (DOF). k L Fig. 1.8 Sprig ass syse.

INTRODUCTION TO MECHANICAL VIBRATIONS 9 The equaio of oio for he free vibraio of a udaped sigle degree of freedo syse ca be rewrie as &&x() + kx () = 0 Dividig hrough by, he equaio ca be wrie i he for &&x() + ω 2 x () = 0 i which ω = codiios k / is a real cosa. The soluio of his equaio is obaied fro he iiial x(0) = x 0, &x(0) = v 0 where x 0 ad v 0 are he iiial displacee ad iiial velociy, respecively. The geeral soluio ca be wrie as x() = A 1 e i ω + A e i ω 2 i which A 1 ad A 2 are cosas of iegraio, boh coplex quaiies. I ca be fially siplified as: x() = X e + e = X cos (ω 2 φ) so ha ow he cosas of iegraio are X ad φ. This equaio represes haroic oscillaio, for which reaso such a syse is called a haroic oscillaor. There are hree quaiies defiig he respose, he apliude X, he phase agle φ ad he frequecy ω, he firs wo depedig o exeral facors, aely, he iiial exciaios, ad he hird depedig o ieral facors, aely, he syse paraeers. O he oher had, for a give syse, he frequecy of he respose is a characerisic of he syse ha says always he sae, idepedely of he iiial exciaios. For his reaso, ω is called he aural frequecy of he haroic oscillaor. The cosas X ad φ are obaied fro he iiial codiios of he syse as follows: X = x i ( ω φ) i ( ω φ) F v + H G I ω K J 0 2 0 v ad φ = a 1 x L NM 0 ω 0 2 O QP The ie period τ, is defied as he ie ecessary for he syse o coplee oe vibraio cycle, or as he ie bewee wo cosecuive peaks. I is relaed o he aural frequecy by τ = 2 π = 2π ω k Noe ha he aural frequecy ca also be defied as he reciprocal of he period, or f = 1 1 k = τ 2π i which case i has uis of cycles per secod (cps), where oe cycle per secod is kow as oe Herz (Hz).

10 SOLVING VIBRATION ANALYSIS PROBLEMS USING MATLAB 1.6.2 FREE VIBRATION OF AN UNDAMPED TORSIONAL SYSTEM A ass aached o he ed of he shaf is a siple orsioal syse (Fig. 1.9). The ass of he shaf is cosidered o be sall i copariso o he ass of he disk ad is herefore egleced. k l I G Fig. 1.9 Torsioal syse. The orque ha produces he wis M is give by M = GJ l where J = he polar ass oe of ieria of he shaf J = diaeer d I K G = shear odulus of he aerial of he shaf. l = legh of he shaf. The orsioal sprig cosa k is defied as k = T GJ θ = l The equaio of oio of he syse ca be wrie as: I G && θ + k θ = 0 F HG πd 4 32 The aural circular frequecy of such a orsioal syse is ω = The geeral soluio of equaio of oio is give by &θ 0 θ() = θ 0 cos ω + si ω ω 1.6.3 ENERGY METHOD for a circular shaf of F HG k I I GKJ 1/ 2 Free vibraio of syses ivolves he cyclic ierchage of kieic ad poeial eergy. I udaped free vibraig syses, o eergy is dissipaed or reoved fro he syse. The kieic eergy T is sored i he ass by virue of is velociy ad he poeial eergy U is sored i he for of srai eergy i elasic deforaio. Sice he oal eergy i he syse

INTRODUCTION TO MECHANICAL VIBRATIONS 11 is cosa, he priciple of coservaio of echaical eergy applies. Sice he echaical eergy is coserved, he su of he kieic eergy ad poeial eergy is cosa ad is rae of chage is zero. This priciple ca be expressed as T + U = cosa d or (T + U) = 0 d where T ad U deoe he kieic ad poeial eergy, respecively. The priciple of coservaio of eergy ca be resaed by T 1 + U 1 = T 2 + U 2 where he subscrips 1 ad 2 deoe wo differe isaces of ie whe he ass is passig hrough is saic equilibriu posiio ad selec U 1 = 0 as referece for he poeial eergy. Subscrip 2 idicaes he ie correspodig o he axiu displacee of he ass a his posiio, we have he T 2 = 0 ad T 1 + 0 = 0 + U 2 If he syse is udergoig haroic oio, he T 1 ad U 2 deoe he axiu values of T ad U, respecively ad herefore las equaio becoes T ax = U ax I is quie useful i calculaig he aural frequecy direcly. 1.6.4 STABILITY OF UNDAMPED LINEAR SYSTEMS The ass/ieria ad siffess paraeers have a affec o he sabiliy of a udaped sigle degree of freedo vibraory syse. The ass ad siffess coefficies eer io he characerisic equaio which defies he respose of he syse. Hece, ay chages i hese coefficie will lead o chages i he syse behavior or respose. I his secio, he effecs of he syse ieria ad siffess paraeers o he sabiliy of he oio of a udaped sigle degree of freedo syse are exaied. I ca be show ha by a proper selecio of he ieria ad siffess coefficies, he isabiliy of he oio of he syse ca be avoided. A sable syse is oe which execues bouded oscillaios abou he equilibriu posiio. 1.6.5 FREE VIBRATION WITH VISCOUS DAMPING Viscous dapig force is proporioal o he velociy &x of he ass ad acig i he direcio opposie o he velociy of he ass ad ca be expressed as F = c &x where c is he dapig cosa or coefficie of viscous dapig. The differeial equaio of oio for free vibraio of a daped sprig-ass syse (Fig. 1.10) is wrie as: x&& + c & x + k x =0

12 SOLVING VIBRATION ANALYSIS PROBLEMS USING MATLAB k c k( + x) cẋ x g (a) (b) Fig. 1.10 Daped sprig-ass syse. By assuig x() = Ce s as he soluio, he auxiliary equaio obaied is c s s k + + = which has he roos 2 0 s 1, 2 = c c 2 2 2 ± F H G I K J The soluio akes oe of hree fors, depedig o wheher he quaiy (c/2) 2 k/ is zero, posiive, or egaive. If his quaiy is zero, c = 2ω This resuls i repeaed roos s 1 = s 2 = c/2, ad he soluio is x() = (A + B)e (c/2) As he case i which repeaed roos occur has special sigificace, we shall refer o he correspodig value of he dapig cosa as he criical dapig cosa, deoed by C c = 2ω. The roos ca be wrie as: e j s 1, 2 = ζ ± ζ2 1 ω where ω = (k/) 1/2 is he circular frequecy of he correspodig udaped syse, ad ζ = c c C = 2 ω c is kow as he dapig facor. If ζ < 1, he roos are boh iagiary ad he soluio for he oio is x() = Xe ζω k si ( ω + φ) where ω d = 1 ζ2 ω is called he daped circular frequecy which is always less ha ω, ad φ is he phase agle of he daped oscillaios. The geeral for of he oio is show i Fig. 1.11. For oio of his ype, he syse is said o be uderdaped. d

INTRODUCTION TO MECHANICAL VIBRATIONS 13 x() Xe <1 Fig. 1.11 The geeral for of oio. If ζ = 1, he dapig cosa is equal o he criical dapig cosa, ad he syse is said o be criically daped. The displacee is give by x() = (A + B)e ω The soluio is he produc of a liear fucio of ie ad a decayig expoeial. Depedig o he values of A ad B, ay fors of oio are possible, bu each for is characerized by apliude which decays wihou oscillaios, such as is show i Fig. 1.12. x() =1 Fig. 1.12 Apliude decayig wihou oscillaios. I his case ζ > 1, ad he syse is said o be overdaped. The soluio is give by: x() = Ce ( ζ + ζ2 1) ω ( ζ ζ2 1) ω + Ce 1 2 The oio will be o-oscillaory ad will be siilar o ha show i Fig. 1.13. x() >1 Fig. 1.13 No-oscillaory oio. 1.6.6 LOGARITHMIC DECREMENT The logarihic decree represes he rae a which he apliude of a free daped vibraio decreases. I is defied as he aural logarih of he raio of ay wo successive apliudes.