M x t = K x F t x t = A x M 1 F t. M x t = K x cos t G 0. x t = A x cos t F 0

Similar documents
Math Week 15: Section 7.4, mass-spring systems. These are notes for Monday. There will also be course review notes for Tuesday, posted later.

Thus the force is proportional but opposite to the displacement away from equilibrium.

Lecture 23 Damped Motion

t is a basis for the solution space to this system, then the matrix having these solutions as columns, t x 1 t, x 2 t,... x n t x 2 t...

Viscous Damping Summary Sheet No Damping Case: Damped behaviour depends on the relative size of ω o and b/2m 3 Cases: 1.

Mon Apr 9 EP 7.6 Convolutions and Laplace transforms. Announcements: Warm-up Exercise:

Hint: There's a table of particular solutions at the end of today's notes.

2.1 Harmonic excitation of undamped systems

Some Basic Information about M-S-D Systems

Chapter 9 Sinusoidal Steady State Analysis

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.

Week 1 Lecture 2 Problems 2, 5. What if something oscillates with no obvious spring? What is ω? (problem set problem)

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

Homework 2 Solutions

Math Week 14 April 16-20: sections first order systems of linear differential equations; 7.4 mass-spring systems.

Announcements: Warm-up Exercise:

Concourse Math Spring 2012 Worked Examples: Matrix Methods for Solving Systems of 1st Order Linear Differential Equations

Module 2 F c i k c s la l w a s o s f dif di fusi s o i n

Fourier Series & The Fourier Transform. Joseph Fourier, our hero. Lord Kelvin on Fourier s theorem. What do we want from the Fourier Transform?

Physics 240: Worksheet 16 Name

Two Coupled Oscillators / Normal Modes

Chapter 2. First Order Scalar Equations

Chapter 6. Systems of First Order Linear Differential Equations

WEEK-3 Recitation PHYS 131. of the projectile s velocity remains constant throughout the motion, since the acceleration a x

Chapter Three Systems of Linear Differential Equations

SIGNALS AND SYSTEMS LABORATORY 8: State Variable Feedback Control Systems

Math 2142 Exam 1 Review Problems. x 2 + f (0) 3! for the 3rd Taylor polynomial at x = 0. To calculate the various quantities:

KINEMATICS IN ONE DIMENSION

Chapter 7: Solving Trig Equations

Math Wednesday March 3, , 4.3: First order systems of Differential Equations Why you should expect existence and uniqueness for the IVP

Wall. x(t) f(t) x(t = 0) = x 0, t=0. which describes the motion of the mass in absence of any external forcing.

dt = C exp (3 ln t 4 ). t 4 W = C exp ( ln(4 t) 3) = C(4 t) 3.

System of Linear Differential Equations

Introduction to Mechanical Vibrations and Structural Dynamics

On Measuring Pro-Poor Growth. 1. On Various Ways of Measuring Pro-Poor Growth: A Short Review of the Literature

Math 333 Problem Set #2 Solution 14 February 2003

e a s a f t dt f t dt = p = p. t = a

Second-Order Differential Equations

Challenge Problems. DIS 203 and 210. March 6, (e 2) k. k(k + 2). k=1. f(x) = k(k + 2) = 1 x k

ECE 2100 Circuit Analysis

Echocardiography Project and Finite Fourier Series

10. State Space Methods

Problem set 6: Solutions Math 207A, Fall x 0 2 x

Math Spring 2015 PRACTICE FINAL EXAM (modified from Math 2280 final exam, April 29, 2011)

INDEX. Transient analysis 1 Initial Conditions 1

A note on diagonalization of integral quadratic forms modulo p m

(1) (2) Differentiation of (1) and then substitution of (3) leads to. Therefore, we will simply consider the second-order linear system given by (4)

Traveling Waves. Chapter Introduction

Lecture 20: Riccati Equations and Least Squares Feedback Control

t, x 2 t above, using input-output modeling. Assume solute concentration is uniform in each tank. If x 1 0 = b 1

FITTING EQUATIONS TO DATA

Y 0.4Y 0.45Y Y to a proper ARMA specification.

Linear Motion, Speed & Velocity

LAPLACE TRANSFORM AND TRANSFER FUNCTION

Electrical Circuits. 1. Circuit Laws. Tools Used in Lab 13 Series Circuits Damped Vibrations: Energy Van der Pol Circuit

Chapter 7 Response of First-order RL and RC Circuits

Explaining Total Factor Productivity. Ulrich Kohli University of Geneva December 2015

Section 3.5 Nonhomogeneous Equations; Method of Undetermined Coefficients

5.1 - Logarithms and Their Properties

KEY. Math 334 Midterm I Fall 2008 sections 001 and 003 Instructor: Scott Glasgow

HOMEWORK # 2: MATH 211, SPRING Note: This is the last solution set where I will describe the MATLAB I used to make my pictures.

ME 391 Mechanical Engineering Analysis

non -negative cone Population dynamics motivates the study of linear models whose coefficient matrices are non-negative or positive.

Lecture 2-1 Kinematics in One Dimension Displacement, Velocity and Acceleration Everything in the world is moving. Nothing stays still.

EECE 301 Signals & Systems Prof. Mark Fowler

Section 7.4 Modeling Changing Amplitude and Midline

Notes 04 largely plagiarized by %khc

Review - Quiz # 1. 1 g(y) dy = f(x) dx. y x. = u, so that y = xu and dy. dx (Sometimes you may want to use the substitution x y

Math 334 Fall 2011 Homework 11 Solutions

Conservation of Momentum. The purpose of this experiment is to verify the conservation of momentum in two dimensions.

Sections 2.2 & 2.3 Limit of a Function and Limit Laws

4.1 - Logarithms and Their Properties

EXERCISES FOR SECTION 1.5

Section 4.4 Logarithmic Properties

d 1 = c 1 b 2 - b 1 c 2 d 2 = c 1 b 3 - b 1 c 3

AP Chemistry--Chapter 12: Chemical Kinetics

Physics 101 Fall 2006: Exam #1- PROBLEM #1

ECE 2100 Circuit Analysis

Finish reading Chapter 2 of Spivak, rereading earlier sections as necessary. handout and fill in some missing details!

Problem Set 7-7. dv V ln V = kt + C. 20. Assume that df/dt still equals = F RF. df dr = =

) were both constant and we brought them from under the integral.

Solutions of Sample Problems for Third In-Class Exam Math 246, Spring 2011, Professor David Levermore

5.2. The Natural Logarithm. Solution

( ) = b n ( t) n " (2.111) or a system with many states to be considered, solving these equations isn t. = k U I ( t,t 0 )! ( t 0 ) (2.

This is an example to show you how SMath can calculate the movement of kinematic mechanisms.

( ) a system of differential equations with continuous parametrization ( T = R + These look like, respectively:

Inventory Analysis and Management. Multi-Period Stochastic Models: Optimality of (s, S) Policy for K-Convex Objective Functions

MA 214 Calculus IV (Spring 2016) Section 2. Homework Assignment 1 Solutions

ME 452 Fourier Series and Fourier Transform

2.3 SCHRÖDINGER AND HEISENBERG REPRESENTATIONS

x y θ = 31.8 = 48.0 N. a 3.00 m/s

Math 334 Test 1 KEY Spring 2010 Section: 001. Instructor: Scott Glasgow Dates: May 10 and 11.

Introduction to Physical Oceanography Homework 5 - Solutions

= ( ) ) or a system of differential equations with continuous parametrization (T = R

!!"#"$%&#'()!"#&'(*%)+,&',-)./0)1-*23)

4.5 Constant Acceleration

1. VELOCITY AND ACCELERATION

1 Widrow-Hoff Algorithm

Problem Set 5. Graduate Macro II, Spring 2017 The University of Notre Dame Professor Sims

ln 2 1 ln y x c y C x

Transcription:

Forced oscillaions (sill undaped): If he forcing is sinusoidal, M = K F = A M F M = K cos G wih F = M G = A cos F Fro he fundaenal heore for linear ransforaions we now ha he general soluion o his inhoogeneous linear proble is of he for = P H, and we've been discussing how o find he hoogeneous soluions H As long as he driving frequency is NOT one of he naural frequencies, we don' epec resonance; he ehod of undeerined coefficiens predics here should be a paricular soluion of he for P = cos c where he vecor c is wha we need o find Eercise ) Subsiue he guess P = cos c ino he DE syse = A cos F o find a ari algebra forula for c = c Noice ha his forula aes sense precisely when is NOT one of he naural frequencies of he syse Soluion: c = A I F Noe, ari inverse eiss precisely if is no an eigenvalue

Eercise ) Coninuing wih he configuraion fro Monday's noes, bu now for an inhoogeneous forced proble, le =, and force he second ass sinusoidally: = cos We now fro previous wor ha he naural frequencies are =, = and ha H = C cos C cos Find he forula for P, as on he preceding page Noice ha his seady periodic soluion blows up as or (If we don' have ie o wor his by hand, we ay sip direcly o he echnology chec on he ne page Bu since we have quic forulas for inverses of by arices, his is definiely a copuaion we could do by hand) Soluion: As long as,, he general soluion = P H is given by = cos C cos C cos Inerpreaion as far as inferred pracical resonance for slighly daped probles: If here was even a sall aoun of daping, he hoogeneous soluion would acually be ransien (i would be eponenially decaying and oscillaing - underdaped) There would sill be a sinusoidal paricular soluion, which would have a forula close o our paricular soluion, he firs er above, as long as, (There would also be a relaively saller sin d er as well) So we can infer he pracical resonance behavior for differen values wih sligh daping, by looing a he size of he c er for he undaped problesee ne page for visualizaions 6

resar : wih LinearAlgebra : A Mari,,,,, : F Vecor, : Iden IdeniyMari : c c ; A wih plos : wih LinearAlgebra : Iden F : # he forula we wored ou by hand plo Nor c,, =, agniude =, color = blac, ile = `iicing pracical resonance` ; # Nor(c( ),) is he agniude of he c( ) vecor agniude iicing pracical resonance 8 6 () plo Nor c Pi T,, T = 5, agniude = 5, color = blac, ile = `pracical resonance as funcion of forcing period` ; pracical resonance as funcion of forcing period 5 agniude 5 5 5 T

There are srong connecions beween our discussion here and he odeling of how earhquaes can shae buildings: As i urns ou, for our physics lab springs, he odes and frequencies are alos idenical:

An ineresing shae-able deonsraion: hp://wwwyouubeco/wach?v=m_jokahzm Below is a discussion of how o odel he unforced "hree-sory" building shown shaing in he video above, fro which we can see which odes will be ecied There is also a "wo-sory" building odel in he video, and is ari and eigendaa follow Here's a scheaic of he hree-sory building: For he unforced (hoogeneous) proble, he acceleraions of he hree assive floors (he op one is he roof) above ground and of ass, are given by = Noe he value in he las diagonal enry of he ari This is because is easuring displaceens for he op floor (roof), which has nohing above i The "" is jus he linearizaion proporionaliy facor, and depends on he ension in he walls, and he heigh beween floors, ec, as discussed on he previous page Here is eigendaa for he unscaled ari = For he scaled ari you'd have he sae eigenvecors, bu he eigenvalues would all be uliplied by he scaling facor and he naural frequencies would all be scaled by Syeric arices lie ours (ie ari equals is ranspose) are always diagonalizable wih real eigenvalues and eigenvecorsand you can choose he eigenvecors o be uually perpendicular This is called he "Specral Theore for syeric arices" and is an iporan fac in any ah and science applicaionsyou can read abou i here: hp://enwiipedia org/wii/syeric_ari) If we ell Maple ha our ari is syeric i will no confuse us wih unsiplified nubers and vecors ha ay loo cople raher han real wih LinearAlgebra : A Mari,,,,,,,,,, ; # I used a leas one decial value so Maple would evaluae in floaing poin A := (5) Digis 5 : # 5 digis should be fine, for our decial approiaions

eigendaa Eigenvecors Mari A, shape = syeric : # o ae advanage of he # specral heore labdas eigendaa : #eigenvalues evecors eigendaa : #corresponding eigenvecors - for fundaenal odes oegas ap sqr, labdas ; # naural angular frequencies f evalf Pi : periods ap f, oegas ; #naural periods eigenvecors ap evalf, evecors ; # ge digis down o 5 89 oegas := 7 5 periods := 87 586 8 eigenvecors := 59 7698 799 7698 799 59 799 59 7698 B Mari,,,,, : eigendaa Eigenvecors Mari B, shape = syeric : # o ae advanage of he # specral heore labdas eigendaa : #eigenvalues evecors eigendaa : #corresponding eigenvecors - for fundaenal odes oegas ap sqr, labdas ; # naural angular frequencies f evalf Pi : periods ap f, oegas ; #naural periods eigenvecors ap evalf, evecors ; # ge digis down o 5 oegas := 68 68 (6) periods := 88 66 eigenvecors := 8565 557 557 8565 (7) Eercise ) Inerpre he daa above, in ers of he naural odes for he shaing building In he youube video he firs ode o appear is he slow and dangerous "sloshing ode", where all hree floors oscillae in phase, wih apliude raios : 59 : 7 fro he firs o he hird floor Wha's he second ode ha ges ecied? The hird ode? (They don' show he hird ode in he video) Rear) All of he ideas we've discussed in secion 5 also apply o olecular vibraions The eigendaa in hese cases is relaed o he "specru" of ligh frequencies ha correspond o he naural fundaenal odes for olecular vibraions

Mah 8- Fri Mar 56 Mari eponenials and linear syses: The analogy beween firs order syses of linear differenial equaions (Chaper 5) and scalar linear differenial equaions (Chaper ) is uch sronger han you ay have epeced This will becoe especially clear on Monday, when we sudy secion 57 Definiion Consider he linear syse of differenial equaions for : = A where A n n is a consan ari as usual If,, n is a basis for he soluion space o his syse, hen he ari having hese soluions as coluns, n is called a Fundaenal Mari (FM) o his syse of differenial equaions Noice ha his equivalen o saying ha X = solves X = A X X nonsingular ie inverible (jus loo colun by colun) Noice ha a FM is jus he Wronsian ari for a soluion space basis Eaple page 5 = : = 5: general soluion Possible FM: general soluion: y = y = = 5 6 6 = c e c = = e e v = e e v = c e 5 e 5 e 5 e 5 e 5 c c