Aaron Boyd's Assignment 8: Difference between revisions

From Class Wiki
Jump to navigation Jump to search
No edit summary
No edit summary
Line 20: Line 20:


<math>\begin{align}
<math>\begin{align}
Now since F_r &= 0 we can ignore it and look only at F_\theta.\\
\text{Now since } F_r &= 0 \text{ we can ignore it and look only at } F_\theta.\\
Since we know F_\theta &= maL and ma &= mLa". We can conclude
\text{ Since we know } F_\theta &= maL \text{ and } ma = mLa. \text{ We can conclude}
\end{align}</math>
\end{align}</math>


<math>\begin{align}


sin(\theta)*mg = mL\theta"
\sin(\theta)*mg = mL\ddot\theta
\end{align}</math>



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


<math>\begin{align}

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

\text{Now we take the laplace transform of this.}\\

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.
\text{Unfortunately the laplace transform of that is horrible. Long, and impossible to solve.}\\
\text{So we use the approximation } \\
(I tried to leave sin(\theta) in the equation. After 4 hours and many wolframalpha.com timeouts I gave up)
\\

\sin(\theta) = \theta
with this new equation we get:
\\

\text{ where } \theta \text{ is small. }\\

\\
g*\theta(t)/L +s^2*\theta(t) - s*\theta(0) - \theta'(0) = 0
\text{(I tried to leave } sin(\theta) \text{ in the equation.}\\

\text{After 4 hours and many wolframalpha.com timeouts I gave up) }\\
we know that \theta(0) = \theta0 and \theta'(0) = 0
\\

\text{with this new equation we get:}\\

\\
solving for \theta(t) we get
g*\frac{\theta(t)}{L +s^2}*\theta(t) - s\theta(0) - \dot\theta(0) = 0\\

\\

\text{we know that } \theta(0) = \theta_0 \text{ and } \dot\theta(0) = 0\\
\theta(t) = s*\theta<sub>0</sub>/((-g/L)+S^2)
\\

\\

\text{solving for } \theta(t) \text{ we get} \\
now we take the inverse laplace transform of that which yields
\\

\theta(t) = s*\frac{\theta_0}{(\frac{-g}{L}+S^2)}\\

\\
\theta(t) = cosh(t*(g/L)^(1/2))
\\
\text{now we take the inverse laplace transform of that which yields }\\
\\
\\
\theta(t) = cosh(t\sqrt(\frac{g}{L})\\
\end{align}</math>

Revision as of 10:40, 1 November 2010

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:



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