Aaron Boyd's Assignment 8: Difference between revisions
Aaron.boyd (talk | contribs) (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>. …') |
Aaron.boyd (talk | contribs) No edit summary |
||
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 |
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. |
||
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 |
|||
\end{align}</math> |
|||
F<sub>x</sub> = T*sin(A) |
|||
F<sub>y</sub> = T*cos(A) - mg = 0 |
|||
Line 13: | Line 13: | ||
<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> |
|||
⚫ | |||
⚫ | |||
<math>\begin{align} |
|||
⚫ | |||
⚫ | |||
\end{align}</math> |
|||
sin( |
sin(\theta)*mg = mL\theta" |
||
Line 30: | Line 31: | ||
\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( |
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( |
(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: |
with this new equation we get: |
||
g* |
g*\theta(t)/L +s^2*\theta(t) - s*\theta(0) - \theta'(0) = 0 |
||
we know that |
we know that \theta(0) = \theta0 and \theta'(0) = 0 |
||
solving for |
solving for \theta(t) we get |
||
\theta(t) = s*\theta<sub>0</sub>/((-g/L)+S^2) |
|||
Line 53: | Line 54: | ||
\theta(t) = cosh(t*(g/L)^(1/2)) |
Revision as of 10:11, 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:
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))