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<sub>x</sub> = T*sin(A)
        F_y &= T\cos(\theta)-mg = 0
\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:
F<sub>r</sub> = T- cos(A)*mg = 0
F<sub>A</sub> = sin(A)*mg = maL




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
<math>\begin{align}
        F_r &= T - mg\cos(\theta) = 0\\
        F_\theta &= \sin(\theta)mg = maL
\end{align}</math>




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


sin(A)*mg = mLA"
<math>\begin{align}


\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.}\\
(I tried to leave sin(A) in the equation. After 4 hours and many wolframalpha.com timeouts I gave up)
\text{So we use the approximation } \\
 
\\
with this new equation we get:
\sin(\theta) = \theta
 
\\
 
\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.}\\
we know that A(0) = A0 and A'(0) = 0
\text{After 4 hours and many wolframalpha.com timeouts I gave up) }\\
 
\\
 
\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>




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




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




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

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 θ<sub>0</sub>. Find a function to determine the angle at any time t. The summation of forces yields Fx=Tsin(θ)Fy=Tcos(θ)mg=0


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


Fr=Tmgcos(θ)=0Fθ=sin(θ)mg=maL


Now since Fr=0 we can ignore it and look only at Fθ. Since we know Fθ=maL and ma=mLa. We can conclude

sin(θ)*mg=mLθ¨

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

θ¨(g/L)sin(θ)=0Now we take the laplace transform of this.Unfortunately the laplace transform of that is horrible. Long, and impossible to solve.So we use the approximation sin(θ)=θ where θ is small. (I tried to leave sin(θ) in the equation.After 4 hours and many wolframalpha.com timeouts I gave up) with this new equation we get:g*θ(t)L+s2*θ(t)sθ(0)θ˙(0)=0we know that θ(0)=θ0 and θ˙(0)=0solving for θ(t) we getθ(t)=s*θ0(gL+S2)now we take the inverse laplace transform of that which yields θ(t)=θ0cos(t(gL))


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


Fx=Tsin(θ)=0 and recognizing T=mg


you can arrive at the same answer