Pure and Practical Resonance in Forced Vibrations

Similar documents
Undetermined Coefficents, Resonance, Applications

An Illustration. The differential equation y y = x + xe x will be solved, verifying that y h = c 1 e x + c 2 e x and y p = x 1 4 xex x2 e x.

Forced Mechanical Vibrations

Forced Damped Vibrations

Math Assignment 5

Differential Equations 2280 Sample Midterm Exam 3 with Solutions Exam Date: 24 April 2015 at 12:50pm

The Harmonic Oscillator

Notes on the Periodically Forced Harmonic Oscillator

Damped Oscillation Solution

Thursday, August 4, 2011

Phase-Amplitude Conversion, Cafe Door, Pet Door, Damping Classifications

Laboratory handouts, ME 340

The student will experimentally determine the parameters to represent the behavior of a damped oscillatory system of one degree of freedom.

4.9 Free Mechanical Vibrations

Math 211. Substitute Lecture. November 20, 2000

M A : Ordinary Differential Equations

Lecture 6: Differential Equations Describing Vibrations

Lecture 14: Forced/free motion; variation of parameters

Solutions for homework 5

Applications of Second-Order Differential Equations

4.2 Homogeneous Linear Equations

MAT 275 Laboratory 6 Forced Equations and Resonance

3.4.1 Distinct Real Roots

Section Mass Spring Systems

Damped & forced oscillators

Laboratory 10 Forced Equations and Resonance

Lab 11 - Free, Damped, and Forced Oscillations

F = ma, F R + F S = mx.

Chapter 2: Linear Constant Coefficient Higher Order Equations

2. Determine whether the following pair of functions are linearly dependent, or linearly independent:

MAT187H1F Lec0101 Burbulla

Mechanical Oscillations

Vibrations and waves: revision. Martin Dove Queen Mary University of London

Second order linear equations

Lecture Notes for Math 251: ODE and PDE. Lecture 16: 3.8 Forced Vibrations Without Damping

Section 3.7: Mechanical and Electrical Vibrations

APPLICATIONS OF SECOND-ORDER DIFFERENTIAL EQUATIONS

dx n a 1(x) dy

Damped harmonic motion

Honors Differential Equations

MAT 275 Laboratory 6 Forced Equations and Resonance

Forced Oscillation and Resonance

Lab 1: Damped, Driven Harmonic Oscillator

Linear Differential Equations. Problems

Differential Equations

Lab 1: damped, driven harmonic oscillator

CHAPTER 4 TEST REVIEW

Vibrations and Waves MP205, Assignment 4 Solutions

Damped Harmonic Oscillator

The Phasor Analysis Method For Harmonically Forced Linear Systems

Double Spring Harmonic Oscillator Lab

MATH 251 Week 6 Not collected, however you are encouraged to approach all problems to prepare for exam

Unforced Oscillations

Oscillations. Simple Harmonic Motion of a Mass on a Spring The equation of motion for a mass m is attached to a spring of constant k is

Math 54. Selected Solutions for Week 10

Work sheet / Things to know. Chapter 3

This work is licensed under a Creative Commons Attribution-Noncommercial-Share Alike 4.0 License.

Dynamics of Structures: Theory and Analysis

Outline of parts 1 and 2

HOMEWORK 4: MATH 265: SOLUTIONS. y p = cos(ω 0t) 9 ω 2 0

Wave Phenomena Physics 15c. Lecture 2 Damped Oscillators Driven Oscillators

Linear Spring Oscillator with Two Different Types of Damping

Topic 5 Notes Jeremy Orloff. 5 Homogeneous, linear, constant coefficient differential equations

5.6 Unforced Mechanical Vibrations

Vibrations Qualifying Exam Study Material

Math 215/255 Final Exam (Dec 2005)

M A : Ordinary Differential Equations

4. Complex Oscillations

MA 266 Review Topics - Exam # 2 (updated)

0.1 Infinite Limits 2 = (2 1) =. Finally, we conclude that

APPM 2360: Midterm 3 July 12, 2013.

Math 2250 Maple Project 2, March Tacoma Narrows. CLASSTIME VERSION A-E, F-K, L-R, S-Z

Chaotic motion. Phys 420/580 Lecture 10

LAB 11: FREE, DAMPED, AND FORCED OSCILLATIONS

Section 4.9; Section 5.6. June 30, Free Mechanical Vibrations/Couple Mass-Spring System

Section 3.4. Second Order Nonhomogeneous. The corresponding homogeneous equation. is called the reduced equation of (N).

Dynamics of structures

Chapter 15 - Oscillations

Physics 2310 Lab #3 Driven Harmonic Oscillator

Section 3.4. Second Order Nonhomogeneous. The corresponding homogeneous equation

Physics 141, Lecture 7. Outline. Course Information. Course information: Homework set # 3 Exam # 1. Quiz. Continuation of the discussion of Chapter 4.

Solutions to the Homework Replaces Section 3.7, 3.8

18.03SC Practice Problems 14

Mathematical Physics

ODE. Philippe Rukimbira. Department of Mathematics Florida International University PR (FIU) MAP / 92

Finish up undetermined coefficients (on the midterm, on WW3) Physics applications - mass springs (not on midterm, on WW4)

BSc/MSci EXAMINATION. Vibrations and Waves. Date: 4 th May, Time: 14:30-17:00

Chapter 1. Harmonic Oscillator. 1.1 Energy Analysis

Solutions to the Homework Replaces Section 3.7, 3.8

L = 1 2 a(q) q2 V (q).

Chapter 2 Second-Order Linear ODEs

KEELE UNIVERSITY PHYSICS/ASTROPHYSICS MODULE PHY OSCILLATIONS AND WAVES PRACTICE EXAM

Oscillatory Motion. Simple pendulum: linear Hooke s Law restoring force for small angular deviations. Oscillatory solution

SECOND ORDER ODE S. 1. A second order differential equation is an equation of the form. F (x, y, y, y ) = 0.

Oscillatory Motion. Simple pendulum: linear Hooke s Law restoring force for small angular deviations. small angle approximation. Oscillatory solution

Math 240: Spring-mass Systems

Math 216 Second Midterm 19 March, 2018

Forced Response - Particular Solution x p (t)

Second Order Linear ODEs, Part II

Chapter 5 Design. D. J. Inman 1/51 Mechanical Engineering at Virginia Tech

Transcription:

Definition of Pure Resonance Typical Pure Resonance Graphic Pure and Practical Resonance in Forced Vibrations Pure Resonance Explained by Undetermined Coefficients The Wine Glass Experiment. Glass at 5mm Rim Excursion, before Breakage. Soldiers Breaking Cadence. Millenium Bridge over The Thames. Definition of Practical Resonance Boundedness of the Homogeneous Solution. Transient Solution. Undetermined Coefficients Steady-State Solution and x p Practical Resonance. Visualization of Practical Resonance.

Definition of Pure Resonance The notion of pure resonance in the differential equation (1) x (t) + ω 2 0 x(t) = F 0 cos(ωt) is the existence of a solution that is unbounded as t. We already know that for ω ω 0, the general solution of (1) is the sum of two harmonic oscillations, hence it is bounded. Equation (1) for ω = ω 0 has by the method of undetermined coefficients the unbounded oscillatory solution x(t) = F 0 2ω 0 t sin(ω 0 t). Pure resonance occurs exactly when the natural internal frequency ω 0 matches the natural external frequency ω, in which case all solutions of the differential equation are unbounded.

Typical Pure Resonance Graphic In Figure 1, pure resonance is illustrated for x (t) + 16x(t) = 8 cos 4t, which in (1) corresponds to ω = ω 0 = 4 and F 0 = 8. x t Figure 1. Pure resonance. Graphed are the envelope curves x = ±t and the solution x(t) = t sin 4t of the equation x (t) + 16x(t) = 8 cos ωt, where ω = 4.

Pure Resonance Explained by Undetermined Coefficients An initial trial solution of is x (t) + 16x(t) = 8 cos ωt x = d 1 cos ωt + d 2 sin ωt. The homogeneous solution x h = c 1 cos 4t + c 2 sin 4t considered in the correction rule has duplicate terms exactly when the natural frequencies match: ω = 4. Then the final trial solution is classified as follows. ω 4 = x(t) = d 1 cos ωt + d 2 sin ωt is bounded, no resonance, ω = 4 = x(t) = t(d 1 cos ωt + d 2 sin ωt) is unbounded, pure resonance. Even before the undetermined coefficients d 1, d 2 are evaluated, we can decide that unbounded solutions occur exactly when frequency matching ω = 4 occurs, because of the amplitude factor t. If ω 4, then x p (t) is a pure harmonic oscillation, which implies it is bounded. If ω = 4, then x p (t) equals a time varying amplitude Ct times a pure harmonic oscillation, hence it is unbounded.

The Wine Glass Experiment Equation x (t) + ω 2 0 x(t) = F 0 cos(ωt) is advertised as the basis for a physics experiment in the Public Television Annenberg CPB Project, called the wine glass experiment. Cal Tech physicist Goodstein, in front of an audience of physics students, equips a lab table with a frequency generator, an amplifier and an audio speaker. The valuable wine glass is replaced by a glass lab beaker. The frequency generator is tuned to the natural frequency of the glass beaker (ω ω 0 ), then the volume knob on the amplifier is suddenly turned up (F 0 adjusted larger), whereupon the sound waves emitted from the speaker break the glass beaker. The 29-minute CPB video can be viewed in Windows Media Player, no cost, in Waves #17 at the link below. The wine glass experiment is the last 4 minutes of the film. http://www.free-ed.net/free-ed/science/physics/ The glass itself will vibrate at a certain frequency, as can be determined experimentally by pinging the glass rim. This vibration operates within elastic limits of the glass and the glass will not break under these circumstances. A physical explanation for the breakage is that an incoming sound wave from the speaker is timed to add to the glass rim excursion. After enough amplitude additions [from 1mm initially to eventually 5mm], the glass rim moves beyond the elastic limit and the glass breaks. The explanation implies that the external frequency from the speaker has to match the natural frequency of the glass. But there is more to it: the glass has some natural damping that nullifies feeble attempts to increase the glass rim amplitude. The physicist uses to great advantage this natural damping to tune the external frequency to the glass. The reason for turning up the volume on the amplifier is to nullify the damping effects of the glass. The amplitude additions then build rapidly as sound wave energy gets stored in the glass, seen as rim vibrations, and the glass breaks.

Glass at 5mm Rim Excursion, before Breakage Figure 2. Still image from a wine glass experiment video. The source for the still image is the wine glass experiment Quicktime video found at Blaze Labs, http://www.blazelabs.com/f-p-glass.asp

Soldiers Breaking Cadence The collapse of the Broughton bridge near Manchester, England in 1831 is blamed for the now standard practise of breaking cadence when soldiers cross a bridge. Bridges like the Broughton bridge have many natural low frequencies of vibration, so it is possible for a column of soldiers to vibrate the bridge at one of the bridge s natural frequencies. The bridge locks onto the frequency while the soldiers continue to add to the excursions with every step, causing larger and larger bridge oscillations.

Millenium Bridge over The Thames In 2000, the Millenium Bridge was opened in London, and closed shortly thereafter, due to unusual vibrations caused by the foot traffic. It was re-opened in 2002 after spending an additional 5 million pounds on bridge re-design. Read about the bridge here: http://en.wikipedia.org/wiki/millennium_bridge_(london). Figure 3. London s Millenium Suspension Footbridge over the Thames.

Definition of Practical Resonance The notion of pure resonance is easy to understand both mathematically and physically, because frequency matching characterizes the event. This ideal situation never happens in the physical world, because damping is always present. In the presence of damping c > 0, it will be established below that only bounded solutions exist for the forced springmass system (2) mx (t) + cx (t) + kx(t) = F 0 cos ωt. Our intuition about resonance seems to vaporize in the presence of damping effects. But not completely. Most would agree that the undamped intuition is correct when the damping effects are nearly zero. Practical resonance is said to occur when the external frequency ω has been tuned to produce the largest possible solution (a more precise definition appears below). It will be shown that this happens for the condition (3) ω = k/m c 2 /(2m 2 ), k/m c 2 /(2m 2 ) > 0. Pure resonance ω = ω 0 k/m is the limiting case obtained by setting the damping constant c to zero in condition (3). This strange but predictable interaction exists between the damping constant c and the size of solutions, relative to the external frequency ω, even though all solutions remain bounded.

Boundedness of the Homogeneous Solution The decomposition of x(t) into homogeneous solution x h (t) and particular solution x p (t) gives some intuition into the complex relationship between the input frequency ω and the size of the solution x(t). The homogeneous solution. For positive damping, c > 0, equation (2) has homogeneous solution x h (t) = c 1 x 1 (t) + c 2 x 2 (t) where according to the recipe the basis elements x 1 and x 2 are given in terms of the roots of the characteristic equation mr 2 +cr+k = 0, as classified by the discriminant D = c 2 4mk, as follows: Case 1, D > 0 Case 2, D = 0 Case 3, D < 0 Symbols α, β satisfy β > 0 and α = c/(2m) < 0. x 1 = e r1t, x 2 = e r2t with r 1 and r 2 negative. x 1 = e r1t, x 2 = te r1t with r 1 negative. x 1 = e αt cos βt, x 2 = e αt sin βt. It follows that x h (t) contains a negative exponential factor, regardless of the positive values of m, c, k. Then y h (t) is bounded.

Transient Solution A solution x(t) is called a transient solution provided it satisfies the relation lim t x(t) = 0. The conclusion: The homogeneous solution x h (t) of the equation mx (t) + cx (t) + kx(t) = 0 is a transient solution for all positive values of m, c, k. A transient solution graph x(t) for large t lies atop the axis x = 0, as in Figure 4, because lim t x(t) = 0. 2 x 0 Figure 4. Transient oscillatory solution x = 2e t (cos t + sin t) of the differential equation x + 2x + 2x = 0. t 6π

Undetermined Coefficients The method of undetermined coefficients applied to mx + cx + kx = F 0 cos ωt gives a trial solution of the form x(t) = A cos ωt + B sin ωt with coefficients A, B satisfying the equations (4) (k mω 2 )A + (cω)b = F 0, ( cω)a + (k mω 2 )B = 0. Solving (4) with Cramer s rule or elimination produces the solution (5) A = (k mω 2 )F 0 (k mω 2 ) 2 + (cω) 2, B = cωf 0 (k mω 2 ) 2 + (cω) 2.

Steady-State Solution and x p The steady state solution, periodic of period 2π/ω, is given by (6) x p (t) = = F 0 ( (k mω 2 ) cos ωt + (cω) sin ωt ) (k mω 2 ) 2 + (cω) 2 F 0 (k mω2 ) 2 + (cω) 2 cos(ωt α), where α is defined by the phase amplitude relations (7) C cos α = k mω 2, C sin α = cω, C = F 0 / (k mω 2 ) 2 + (cω) 2. The terminology steady state refers to that part xss(t) of the solution x(t) that remains when the transient portion is removed, that is, when all terms containing negative exponentials are removed. As a result, for large T, the graphs of x(t) and xss(t) on t T are the same. This feature of xss(t) allows us to find its graph directly from the graph of x(t). We say that xss(t) is observable, because it is the solution visible in the graph after the transients (negative exponential terms) die out.

Practical Resonance Practical resonance is said to occur when the external frequency ω has been tuned to produce the largest possible steady state amplitude. Mathematically, this happens exactly when the amplitude function C = C(ω) defined in (7) has a maximum. If a maximum exists on 0 < ω <, then C (ω) = 0 at the maximum. The power rule implies (8) C (ω) = F 0 2 2(k mω 2 )( 2mω) + 2c 2 ω ((k mω 2 ) 2 + (cω) 2 ) 3/2 = ω ( 2mk c 2 2m 2 ω 2) C(ω) 3 F 2 0 If 2km c 2 0, then C (ω) does not vanish for 0 < ω < and hence there is no maximum. If 2km c 2 > 0, then 2km c 2 2m 2 ω 2 = 0 has exactly one root ω = k/m c 2 /(2m 2 ) in 0 < ω < and by C( ) = 0 it follows that C(ω) is a maximum. Practical resonance for mx (t) + cx (t) + kx(t) = F 0 cos ωt occurs precisely when the external frequency ω is tuned to ω = k/m c2 /(2m 2 ) and k/m c 2 /(2m 2 ) > 0.

Visualization of Practical Resonance In Figure 5, the amplitude of the steady state periodic solution is graphed against the external natural frequency ω, for the differential equation x + cx + 26x = 10 cos ωt and damping constants c = 1, 2, 3. The practical resonance condition is ω = 26 c2 /2. As c increases from 1 to 3, the maximum point (ω, C(ω)) satisfies a monotonicity condition: both ω and C(ω) decrease as c increases. The maxima for the three curves in the figure occur at ω = 25.5, 24, 21.5. Pure resonance occurs when c = 0 and ω = 26. C c = 1 2 c = 3 Figure 5. Practical resonance for x + cx + 26x = 10 cos ωt: amplitude C = 10/ (26 ω 2 ) 2 + (cω) 2 versus external frequency ω for c = 1, 2, 3. ω