Aaron Boyd's Assignment 8: Difference between revisions

From Class Wiki
Jump to navigation Jump to search
(Created page with '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 A<sub>0</sub>. …')
 
No edit summary
 
(7 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 A<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}

F_x &= T\sin(\theta)\\

F_y &= T\cos(\theta)-mg = 0
F<sub>x</sub> = T*sin(A)
\end{align}</math>
F<sub>y</sub> = T*cos(A) - mg = 0




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


now:




<math>\begin{align}
F<sub>r</sub> = T- cos(A)*mg = 0
F_r &= T - mg\cos(\theta) = 0\\
F_\theta &= \sin(\theta)mg = maL
\end{align}</math>


F<sub>A</sub> = sin(A)*mg = maL


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


<math>\begin{align}
now since F<sub>r</sub> = 0 we can ignore it and look only at FA.

Since we know F<sub>A</sub> = maL and ma = mLa". We can conclude



sin(A)*mg = mLA"


\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}

A" - sin(A)(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(A) = A where A 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(A) 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*A(t)/L +s^2*A(t) - s*A(0) - A'(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 A(0) = A0 and A'(0) = 0
\\

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

\\
solving for A(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\\
\\
\\
\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>




You can solve for the same thing from the cartesian coordinates. Taking:
A(t) = s*A0/((-g/L)+S^2)




<math>\begin{align}
now we take the inverse laplace transform of that which yields
F_x &= T\sin(\theta) = 0\\
\text { and recognizing } T = mg\\
\end{align}</math>




you can arrive at the same answer
A(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