1 / 18

Chapter 3

Chapter 3. Simple Harmonic Motion. 3.1 Simple Harmonic Motion. Force on the pendulum. Equation of motion (EoM). Analytical solution:. constants determined by initial conditions. The period of the SHO is given by. Euler method to solve the EoM.

mimis
Download Presentation

Chapter 3

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript


  1. Chapter 3 Simple Harmonic Motion

  2. 3.1 Simple Harmonic Motion Force on the pendulum Equation of motion (EoM) Analytical solution: constants determined by initial conditions. The period of the SHO is given by

  3. Euler method to solve the EoM Run t for a few period (T) to see the oscillatory behavior Set Dt, q0, w0, l, g, iterate ti+1=ti + Dt sample code 3.1.1

  4. Total energy of SHO You may use the Mathematica command Series[Cos[x],{x,0,2}]to obtain the Taylor series expansion of cosq at q = 0. This is amount to making the assumption that q  0.

  5. Euler method fails in this case The oscillation amplitude increases with time  Euler method is not suitable for oscillatory motion (but may be good for projectile motion which is non-periodic) The peculiar effect of the numerical solution can also be investigated from the energy of the oscillator when Euler method is applied: E≈ (1/2)m l2 [w2 + (g/l)q 2] (with q  0) Descretise the equation by applying Euler method, Ei+1=Ei+ (1/2)mgl2[wi2 + (g/l)qi2](Dt)2  E increases with i  this is an undesirable result.

  6. Euler vs. Euler-Cromer method In the Euler method: wi+1 is calculated based on previous values wi and qi So is qi+1 is calculated based on previous values wi and qi • In the Euler-Cromer method: • wi+1 is calculated based on previous values wi and qi • but qi+1 is calculated based on present values wi+1 and previous value qi sample code 3.2.2

  7. Why Euler-Cromer works better? Because it conserver energy over each complete period of motion.

  8. 3.2 Forced SHM Drag force on a moving object, fd = - kv For a pendulum, instantaneous velocity v = wl = l (dq/dt) Hence, fd = - kl (dq/dt). The net force on the forced pendulum along the tangential direction l fd - kl (dq/dt). r w E.o.M

  9. Analytical solution Underdamped regime (small damping). Still oscillate, but amplitude decay slowly over many period before dying totally. Overdamped regime (very large damping), decay slowly over several period before dying totally. q is dominated by exponential term. Critically damped regime, intermediate between under- and overdamping case.

  10. Cromer-Euler descretisation sample code 3.2.3

  11. Numerical solution of the damped oscillator for q = 10,5,1, l=1.0 m This is generated with NDSolve sample code 3.2.4

  12. Adding driving force WD frequency of the applied force - kl (dq/dt) + FDsin(WDt)

  13. Analytical solution Resonance happens when

  14. Euler-Kromer descretisation sample code 3.2.5

  15. Second order Runge-Kutta method Consider a generic second order differential equation. It can be numerically solved using second order Runge-Kutta method. First, split the second order DE into two first order parts:

  16. Algorithm Set initial conditions: Calculate Calculate Calculate Calculate

  17. RK2 for SHO

  18. Use RK2 to calculate the total energy of SHO We have seen previously that energy as calculated using Euler method is unstable. Now try to calculate E again, using the values of wiandqias calculated using RK2. You will have to do this in the HA3. Ei+1

More Related