Aaron Boyd's Assignment 8: Difference between revisions

From Class Wiki
Jump to navigation Jump to search
No edit summary
No edit summary
 
(6 intermediate revisions by the same user not shown)
Line 1: Line 1:
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 \theta<sub>0</sub>. Find a function to determine the angle at any time t.
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 <math>\theta<sub>0</sub></math>. Find a function to determine the angle at any time t.
The summation of forces yields
The summation of forces yields
<math>\begin{align}
<math>\begin{align}
Line 5: Line 5:
F_y &= T\cos(\theta)-mg = 0
F_y &= T\cos(\theta)-mg = 0
\end{align}</math>
\end{align}</math>





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


now:




Line 20: Line 18:


<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}
\ddot\theta - (g/L)\sin(\theta) = 0\\
\\
\text{Now we take the laplace transform of this.}\\
\text{Unfortunately the laplace transform of that is horrible. Long, and impossible to solve.}\\
\text{So we use the approximation } \\
\\
\sin(\theta) = \theta
\\
\text{ where } \theta \text{ is small. }\\
\\
\text{(I tried to leave } sin(\theta) \text{ in the equation.}\\
\text{After 4 hours and many wolframalpha.com timeouts I gave up) }\\
\\
\text{with this new equation 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\\
\\
\\
\text{solving for } \theta(t) \text{ we get} \\
\\
\theta(t) = s*\frac{\theta_0}{(\frac{-g}{L}+S^2)}\\
\\
\\
\text{now we take the inverse laplace transform of that which yields }\\
\\
\\
\theta(t) = \theta_0cos(t\sqrt(\frac{g}{L}))\\
\end{align}
</math>


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


You can solve for the same thing from the cartesian coordinates. Taking:


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)


<math>\begin{align}
with this new equation we get:
F_x &= T\sin(\theta) = 0\\

\text { and recognizing } T = mg\\

\end{align}</math>
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*\theta<sub>0</sub>/((-g/L)+S^2)


now we take the inverse laplace transform of that which yields




you can arrive at the same answer
\theta(t) = cosh(t*(g/L)^(1/2))

Latest revision as of 11:16, 19 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 . 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.



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


You can solve for the same thing from the cartesian coordinates. Taking:



you can arrive at the same answer