MIT - 16.20 Fall, 2002 Unit 20 Solutions for Single Spring-Mass Systems Paul A. Lagace, Ph.D. Professor of Aeronautics & Astronautics and Engineering Systems Paul A. Lagace ? 2001 qt MIT - 16.20 Fall, 2002 Return to the simplest system: the single spring-mass… This is a one degree-of-freedom system with the governing equation: mq ˙˙ + k q = F First consider … Free Vibration ? Set F = 0 resulting in: mq ˙˙ + k q = 0 The solution to this is the homogeneous solution to the general equation. For an Ordinary Differential Equation of this form, know that the solution is of the form: () = e pt ? mp 2 e pt + k e pt = 0 2 ? mp + k = 0 (in order to hold for all t) Paul A. Lagace ? 2001 Unit 20 - 2 qt it it MIT - 16.20 Fall, 2002 2 k ? p = ? m ? p =± i k m where: i = ? 1 k m = ω = natural frequency of single degree-of-freedom system [rad/sec] Important concept that natural frequency = stiffness mass So have the equation: () = C 1 e + ω + C 2 e ? ω Paul A. Lagace ? 2001 Unit 20 - 3 qt MIT - 16.20 Fall, 2002 from mathematics, know this is: qt () = C 1 ′ sin ω t + C 2 ′ cos ω t general solution Now use the Initial Conditions: @ t = 0 q = q 0 ? C 2 ′ = q 0 @ t = 0 q ˙ = q ˙ 0 ? C 1 ′ = q ˙ 0 ω This results in: () = q ˙ 0 sin ω t + q 0 cos ω t ω with: ω = k m This is the basic, unforced response of the system So if one gives the system an initial displacement A and then lets go: q 0 = A q ˙ 0 = 0 Paul A. Lagace ? 2001 Unit 20 - 4 qt MIT - 16.20 Fall, 2002 The response is: () = A cos ω t Figure 20.1 Basic unforced dynamic response of single spring-mass system But, generally systems have a force, so need to consider: Paul A. Lagace ? 2001 Unit 20 - 5 MIT - 16.20 Fall, 2002 Forced Vibration The homogeneous solution is still valid, but must add a particular solution The simplest case here is a constant load with time … Figure 20.2 Representation of constant applied load with time (think of the load applied suddenly ? step function at t = 0) The governing equation is: mq ˙˙ + k q = F 0 The particular solution has no time dependence since the force has no time dependence: F 0 q particular = k Paul A. Lagace ? 2001 Unit 20 - 6 qt MIT - 16.20 Fall, 2002 Now use the homogeneous solution with this to get the total solution: F () = C 1 sin ω t + C 2 cos ω t + 0 k The Initial Conditions are: @ t = 0 q = 0 q ˙ = 0 0 q () = 0 ? C = ? F 0 2 k q 0 ˙ () = 0 ? C 1 = 0 So the final solution is: F q = 0 ( 1 ? cos ω t ) k with ω = k m Plotting this: Paul A. Lagace ? 2001 Unit 20 - 7 MIT - 16.20 Fall, 2002 Figure 20.3 Ideal dynamic response of single spring-mass system to constant force Dynamic response Static response Note that: Dynamic response = 2 x static response “dynamic magnification factor ” - will be larger when considering stresses over their static values Know this doesn ’ t really happen (i.e. response does not continue forever) What has been left out? DAMPING Paul A. Lagace ? 2001 Unit 20 - 8 MIT - 16.20 Fall, 2002 Actual behavior would be … Figure 20.4 Actual dynamic response (with damping) of single spring- mass system to constant force rate of damping dependent on magnitude of damper (c) Have considered a simple case. But, in general forces are not simple steps. Consider the next “level”… The Unit Impulse Paul A. Lagace ? 2001 Unit 20 - 9 MIT - 16.20 Fall, 2002 An impulse occurs at time t = τ . Such a force is represented by the Dirac delta function : Ft () = δ ( t ? τ ) where: δ ( t ? τ ) = 0 @ t ≠ τ δ ( t ? τ ) →∞ @ t = τ and: ∞ ?∞ δ ( t ? τ ) d t = 1 ∫ Recall: force x time = impulse Hence, ∞ ?∞ gt ∫ () δ ( t ? τ ) d t = 0 everywhere but at t = τ ∞ τ = g () ∫ ?∞ δ ( t ? τ ) d t now a constant = 1 with regard to time Paul A. Lagace ? 2001 Unit 20 - 1 0 τττ MIT - 16.20 Fall, 2002 ∞ () δ ( t ? τ ) d t = g () ? ∫ ?∞ gt τ Use this force function in the governing equation: mq ˙˙ + k q = δ ( t ? τ ) Figure 20.5 Representation of unit impulse force at time τ Use Initial Condition of system at rest: @ t < τ q = 0 q ˙ = 0 Paul A. Lagace ? 2001 Unit 20 - 1 1 ? ? MIT - 16.20 Fall, 2002 To get the result, integrate the equation over the regime where the force is nonzero (t = τ to τ + ε ) 2 + τε τε + + ? ∫ τ τε m dq dt + k ∫ τ q dt = ∫ τ δ ( t ? τ ) d t dt 2 use trapezoidal rule Trapezoidal Rule : B 1 ∫ A qd t = 2 [ q ( B ) + q ( A ) ] × ( t B ? t A ) This gives: ? m ? ? ? ? dq ? ? ? ?? dq ?? ?? + k 1 [ q ( τ + ε ) + q () ] τ ε = 1 ? ? dt ? τε ? dt ? τ ? 2 + = 0 goes to zero = 0 from Initial from since ε→ 0 Condition Initial Condition All this results in: Paul A. Lagace ? 2001 Unit 20 - 1 2 qt MIT - 16.20 Fall, 2002 dq dt m ? ? ? ? ? ? = + τε 1 q () = + τε 0 A unit impulse at time τ is thus a free vibration problem with an initial velocity (equal in magnitude to the inverse of the mass) Note : Units are consistent since the integral of δ (t - τ ) is a force x time. So: mass × length 2 × time ? dq ? Force × time ( time ) = = ? ? dt ? ? mass mass = length = [ velocity ] time So use the homogeneous solution: () = C 1 sin ω t + C 2 cos ω t but: Paul A. Lagace ? 2001 Unit 20 - 1 3 qt qt MIT - 16.20 Fall, 2002 ? shift the time { t ’ = 0 @ t = τ ? (t - τ )} ? use: 1 qt = τ ) = ˙ ( m qt = τ ) = 0 ( Thus: () = C 1 sin ω ( t ? τ ) + C 2 cos ω ( t ? τ ) Initial Conditions give: qt = τ ) = 0 ? C 2 = 0 ( 1 1 1 qt = τ ) = ? = C 1 ω ? C 1 = m ω ˙ ( m m Thus: 1 () = sin ω ( t ? τ ) m ω with: ω = k m Paul A. Lagace ? 2001 Unit 20 - 1 4 τττ qt MIT - 16.20 Fall, 2002 Figure 20.6 Dynamic response of single spring-mass system to unit force impulse at time τ Summarizing … the response to a unit impulse at t = τ is: 1 () = sin ω ( t ? τ ) for t ≥ τ m ω 0 for t ≤ τ For an impulse of arbitrary magnitude I o , multiply the solution by I o Paul A. Lagace ? 2001 Unit 20 - 1 5 MIT - 16.20 Fall, 2002 For convenience, write this particular response as: h (t - τ ) = response to unit impulse Now, how does one progress to an … Arbitrary Force Consider some arbitrary force with time: Figure 20.7 Representation of arbitrary force with time Break the force into/represent it as a series of impulses: Paul A. Lagace ? 2001 Unit 20 - 1 6 MIT - 16.20 Fall, 2002 Figure 20.8 Representation of arbitrary force with time as series of impulses over any interval ?τ : Impulse = F ( τ ) δ (t - τ ) The response to any particular impulse is: qt ( ) = F ( τ ) ? τ h ( t ? τ ) where h(t - τ ) is characteristic response of system to an impulse. So the total response is the summation of the responses to all the impulses. In the limit, there are infinite impulses leading to the integral: t qt τ () = ∫ 0 F () h ( t ? τ ) d τ Paul A. Lagace ? 2001 Unit 20 - 1 7 MIT - 16.20 Fall, 2002 This is known as: Duhamel ’s integral or The convolution integral or The linear superposition integral This is a general case. For the particular single spring-mass system: 1 t () = m ω ∫ 0 F () sin ω ( t ? τ ) d τ qt τ response to arbitrary F(t) Additional initial conditions (velocity and displacement) can be added to this through the homogeneous solution to give: 0 () = m 1 ω ∫ 0 t F () sin ω ( t ? τ ) d τ + q ˙ () sin ω t qt τ ω 0 + q () cos ω t Of all the arbitrary forces, there is one form of particular interest: Paul A. Lagace ? 2001 Unit 20 - 1 8 MIT - 16.20 Fall, 2002 Sinuisoidal Force This is an important basic input ? motors ? helicopters ? other cyclical cases Use the integral previously developed for the case of the basic force: Ft () = F sin ? t o (use capital ? to differentiate from system natural frequency) ? - forcing frequency ω - system (response) natural frequency Paul A. Lagace ? 2001 Unit 20 - 1 9 MIT - 16.20 Fall, 2002 Figure 20.9 Representation of sinusoidal forcing function The solution is thus: t qt m ω ∫ 0 sin ? τ sin ω ( t ? τ ) d τ ( ) = F 0 One can perform this integral to get the solution or Go back to the original Ordinary Differential Equation: mq ˙˙ + k q = F o sin ? t Paul A. Lagace ? 2001 Unit 20 - 2 0 MIT - 16.20 Fall, 2002 The overall solution is: qt () = C 1 sin ω t + C 2 cos ω t + q particular homogeneous solution Can see: q particular = C 3 sin ? t Plugging this in the governing Ordinary Differential Equation: ? m ? 2 C 3 sin ? t + k C 3 sin ? t = F o sin ? t F F o o = ? C 3 = ? m ? 2 + k k ? 1 ? ? 2 m ? ? k ? use: F o ω = k m ? C = 3 ? ? 2 ? k ? 1 ? ω 2 ? ? ? This gives the overall solution of: Paul A. Lagace ? 2001 Unit 20 - 2 1 MIT - 16.20 Fall, 2002 F qt () = C sin ω t + C cos ω t + o ? 2 ? sin ? t ? 1 2 Starting transient with natural frequency (this will eventually die out due to damping) k ? 1 ? ω 2 ? ? ? Steady state response following frequency of forcing function Now, using the Initial Conditions: q (0) = 0 ? C 2 = 0 F ? q ˙ (0) = 0 ? ω C 1 + ? o ? 2 ? = 0 k ? 1 ? ω 2 ?? ? ? F o ? C 1 = ? ? ω ? 2 ? k ? 1 ? ω 2 ? ? ? Thus, the final solution is: Paul A. Lagace ? 2001 Unit 20 - 2 2 MIT - 16.20 Fall, 2002 qt F k t o () sin sin = ? ? ? ? ? ? ? ? + ? ? ? ? 1 2 2 ? ? ? ω ω ω t starting steady transient state Notes: F o ? = static response k 1 ? ? 2 ? = Dynamic Magnification Factor (DMF) ? ? 1 ? ω 2 ? ? ? ? For low ? , response approximately static ? For high ? , response goes to zero (waves!) ? For medium ? … .DMF varies. Medium depends on ω which is system structural response ? At ω = ? , DMF → ∞ This is known as “Resonanc e ” Paul A. Lagace ? 2001 Unit 20 - 2 3 MIT - 16.20 Fall, 2002 Figure 20.10 Dynamic Magnification Factor versus ratio of frequency of forcing function to system natural frequency waves Thus, if one excites a system at/near its natural frequency, very large responses result (damping and nonlinearities keep the response from going to infinity) Paul A. Lagace ? 2001 Unit 20 - 2 4 MIT - 16.20 Fall, 2002 Examples: Flutter (window blinds in wind!) Rule of Thumb : stay away from natural frequency( ies ) of system. Further Note : an arbitrary force can be broken down into a sinusoidal series. If any of the force components of importance have a frequency near a system natural frequency ? Trouble ! The next task is to take the basic points learned for this single degree-of-function system and extend them to a multiple degree-of-function system. But first need to see how one can accomplish this … Paul A. Lagace ? 2001 Unit 20 - 2 5