Published on

06-Jan-2016View

26Download

0

Embed Size (px)

DESCRIPTION

Sect. 6.6: Damped, Driven Pendulum. Consider a plane pendulum subject to an an applied torque N & subject to damping by the viscosity of the medium (say, air) in which it swings . Mass m . See figure: The angle of swing is allowed in the entire range: - - PowerPoint PPT Presentation

Transcript

Sect. 6.6: Damped, Driven PendulumConsider a plane pendulum subject to an an applied torque N & subject to damping by the viscosity of the medium (say, air) in which it swings. Mass m. See figure:

The angle of swing is allowed in the entire range: - The usual assumption of a massless rod. Length = R

First, a simple pendulum, no applied torque. Newtons 2nd Law or Lagranges equations give the equation of motion: Ngrav = I = I(d2/dt2), I = mR2 Ngrav = - mgR sin mR2(d2/dt2) + mgRsin = 0 (1)NOTE: The motion is clearly NOT simple harmonic!For small angles sin & (1) becomes: mR2(d2/dt2) + mgR = 0 Or: (d2/dt2) + (g/R) = 0 (2) (2) is clearly simple harmonic with frequency 0 = (g/R)

Now, add an applied torque = NNo damping for now. We still have gravitational torqueNgrav = - mgRsin Newtons 2nd Law or Lagranges equations give the equation of motion: N + Ngrav = mR2(d2/dt2) (3) There is a critical angle c at which the applied torque & the gravitational torque exactly balance (for which the pendulum is at static equilibrium!). This is given byN = - Ngrav = mgR sinc (4) At = c have static equilib. Both angular velocity & angular acceleration vanish: (d/dt) = 0 & (d2/dt2) =0

At critical angle c for staticequilib the applied torque is N = mgRsinc (4)(4) Obviously, the greater the applied torque N, the larger the critical angle c.(4) There is a critical torque Nc for which the critical anglec = () & the pendulum is horizontal; c = () is obviouslythe maximum critical angle for static equilibrium. From (4), this is: Nc = mgR (5)If the applied torque N > Nc, then N > mgRsin for all angles Thependulum can rotate over the top & goes from 0 to 2. The motionthen is rotation, not oscillation. The rotational motion still satisfies theequation of motion: N+ Ngrav= mR2(d2/dt2) (3)Solution to (3) will give time dependent angular velocity = (d/dt).

Its useful to remember all this when now considering the motion of the driven, damped pendulum. Assume, as usual, that the torque due to the frictionalor damping force is angular velocity = (d/dt), with damping constant equal to the viscosity : Nfr - Newtons 2nd Law gives the equation of motion: N + Nfr + Ngrav = mR2(d2/dt2) (3) Putting in the forms for Nfr & Ngrav & putting them on the other side of the equation gives: N = mR2(d2/dt2) + (d/dt) + mgRsin (6)

N = mR2(d2/dt2) + (d/dt) + mgRsin (6) Define a critical angular speed the angular speed at which the damping torque Nfr = exactly equals the critical torque for static equilibrium Nc = mgR:c Nc-1 = mgR-1 (7)Write (6) in dimensionless form: Divide by Nc & use (7): (N/Nc) = (0)-2(d2/dt2) + (c)-1(d/dt) + sin (8)Recall: 0 = (g/R) = harmonic oscillation frequency for the small angle problem. (8) is the equation we will now work with in several special cases. There is rich, varied, complex behavior which can be obtained from it! This includes the possibility of chaotic behavior! For discussion of conditions under which (8) can lead to chaos, see the text by Marion, Ch. 4.

(N/Nc) = (0)-2(d2/dt2) + (c)-1(d/dt) + sin (8) where: Nc = mgR, 0 = (g/R) , c = mgR-1 (8): A nonlinear, inhomogeneous, 2nd order differential equation! No general theory of nonlinear differential equations exists! We can get the solution in special cases or else we need numerical solutions. Solution to the inhomogeneous equation = solution to homogeneous equation + a particular solution to the inhomogeneous equation. Due to damping, the 2nd term on right side (8), the solution to the homogeneous equation will exponentially decay with time (transient solution!): exp(-t/c) = exp[- t(mgR)-1] We look at solutions to (8) at long times (t) >> (mgR) so the transient solution has decayed away & only the particular solution to The inhomogeneous equation remains Dynamic Steady State

(N/Nc) = (0)-2(d2/dt2) + (c)-1(d/dt) + sin (8) where: Nc = mgR, 0 = (g/R) , c = mgR-1 Now, a qualitative discussion of several special cases of solutions to (8) under Dynamic Steady State conditions:1. Low applied torques: N Nc (the applied torque is less than the critical torque for static equilibrium at c = ()): In this case, after the transient dies away, the pendulum will eventually go to a Static Steady State, at which applied torque & gravitational torque exactly balance (so the pendulum is at static equilibrium!). We had this earlier. It is given by N = - Ngrav = mgRsinc = Ncsinc For constant N, after long enough times, the pendulum will stop at angle c !

(N/Nc) = (0)-2(d2/dt2) + (c)-1(d/dt) + sin (8) Special cases of solutions to (8) in the Dynamic Steady State:2. Undamped Motion ( = 0) & constant applied torque N: Its best in this case to go back to Newtons 2nd Law: Nnet = N + Ngrav = N - Ncsin = mR2(d2/dt2)a. The total torque, Nnet is dependent! It takes on special values at 4 angles: = 0Nnet = N = ()Nnet = N - Nc = Nnet = N = 3(), Nnet = N + Nc .

(N/Nc) = (0)-2(d2/dt2) + (c)-1(d/dt) + sin (8) Special cases of solutions to (8) in Dynamic Steady State:2. Undamped Motion ( = 0) & constant applied torque N: Nnet = N + Ngrav = N - Ncsin = mR2(d2/dt2)b. If N > Nc, the motion is continuously accelerated rotation! The pendulum gains energy as time goes on. The angular speed increases with time, but with periodic fluctuations every cycle. The time averaged angular speed: = increases continually with time. See figure (which is with damping).

(N/Nc) = (0)-2(d2/dt2) + (c)-1(d/dt) + sin (8) Special cases of solutions to (8) in Dynamic Steady State:3. Damping is present, small c NcThe angular speed increases until, at some time t, the damping term approaches the applied torque: (d/dt) N. At this time, average angular speed approaches the limiting value: L as in the figure. After this time, the angular acceleration fluctuates about a zero average value: = 0. In this state, the pendulum undergoes quasi-staticmotion: It rotates with angular speed which has periodicfluctuations about its average L, but remains close to it.

(N/Nc) = (0)-2(d2/dt2) + (c)-1(d/dt) + sin (8) Special cases of solutions to (8) in Dynamic Steady State:3. Damping is present, small c NcMore on quasi-static motion: Rotates with near L,fluctuating periodically about it. In this state, the angularacceleration fluctuates about a 0 average value: = 0. Approximately: (N/Nc) (c)-1(d/dt) + sin (8) (8): Solve analytically! Student exercise. Then, time average & find: = 0, N < Nc = c[(N/Nc)2 - 1] , N > Nc c(N/Nc) , N >> Nc. See figure:

- (N/Nc) = (0)-2(d2/dt2) + (c)-1(d/dt) + sin (8) Special cases of solutions to (8) in Dynamic Steady State:3. Damping is present, small c
(N/Nc) = (0)-2(d2/dt2) + (c)-1(d/dt) + sin (8) Special cases of solutions to (8) in Dynamic Steady State:3. Damping is present, small c > 2Nc c(N/Nc) >> c, the solution to (8) approximately gives (t) + sin(t) ( & constants) as for point B in the figures.

(N/Nc) = (0)-2(d2/dt2) + (c)-1(d/dt) + sin (8) Special cases of solutions to (8) in Dynamic Steady State:4. Negligible (or small) damping ( 0, c >> 0):a. If we wait long enough, we can get Static Steady State, at which the applied torque & the gravitational torque exactly balance (static equilibrium!) Given byN = - Ngrav = mgRsin c = Nc sinc Nc = 0, N Nc (a)b. Also, the solution just discussed, in which the applied torque balances the time average damping torque also applies for all N. c(N/Nc), all N (b) See figure:

(N/Nc) = (0)-2(d2/dt2) + (c)-1(d/dt) + sin (8) Special cases of solutions to (8) in Dynamic Steady State:4. Negligible (or small) damping ( 0, c >> 0): = 0, N Nc (a) c(N/Nc), all N (b)Discussion: See figure. The system hashysteresis (differing behavior for increasing & decreasing N): Starting from N = 0, as N increases for N Nc, the pendulum is always stopped atangles c given by N = Ncsinc (a) gives = 0. When N reaches Nc, JUMPSdiscontinuously: (b) gives = c(N=Nc). Also, increases linearly as c(N/Nc) as N increases > Nc. Fordecreasing N, (b) applies, all N & decreases linearly to 0.

(N/Nc) = (0)-2(d2/dt2) + (c)-1(d/dt) + sin (8) Special cases of solutions to (8) in Dynamic Steady State:5. General case: Response for c > 0 What is the response for intermediate conditions where c 0? To answer this, we must solve the general equation (8). Must do this numerically. Results for N vs. for c = 20 are shown here Increasing N: Get = 0 until N = Nc. At N = Nc, JUMPS discontinuously to = 2c. Then, increases non-linearly with N until c >> 0 is reached & it becomeslinear. Decreasing N: Hysteresis with =0 for N = Nc < Nc