Aaron Boyd's Assignment 8

From Class Wiki
Revision as of 10:11, 1 November 2010 by Aaron.boyd (talk | contribs)
Jump to navigation Jump to search

I decided to use laplace transforms to solve a pendulum equation. A pendulum with a weight of mass m and a massless rod length L is released from an initial angle \theta0. Find a function to determine the angle at any time t. The summation of forces yields


Polar coordinates may be easier to use, lets try that.

now:



Failed to parse (unknown function "\begin{align}"): {\displaystyle \begin{align} Now since F_r &= 0 we can ignore it and look only at F_\theta.\\ Since we know F_\theta &= maL and ma &= mLa". We can conclude \end{align}}


sin(\theta)*mg = mL\theta"


canceling the common mass term and rearranging a bit we get.


\theta" - sin(\theta)(g/L) = 0


now we take the laplace transform of this. Unfortunately the laplace transform of that is horrible. Long, complicated, and nearly impossible to solve. So we use the approximation sin(\theta) = \theta where \theta is small. (I tried to leave sin(\theta) in the equation. After 4 hours and many wolframalpha.com timeouts I gave up)

with this new equation we get:


g*\theta(t)/L +s^2*\theta(t) - s*\theta(0) - \theta'(0) = 0

we know that \theta(0) = \theta0 and \theta'(0) = 0


solving for \theta(t) we get


\theta(t) = s*\theta0/((-g/L)+S^2)


now we take the inverse laplace transform of that which yields


\theta(t) = cosh(t*(g/L)^(1/2))