Aaron Boyd's Assignment 8: Difference between revisions

From Class Wiki
Jump to navigation Jump to search
No edit summary
No edit summary
 
(5 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 31: Line 29:
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}
<math> \begin{align}
\ddot\theta - (g/L)\sin(\theta) = 0\\
\ddot\theta - (g/L)\sin(\theta) = 0\\
\\
\\
Line 60: Line 58:
\\
\\
\\
\\
\theta(t) = cosh(t\sqrt(\frac{g}{L})\\
\theta(t) = \theta_0cos(t\sqrt(\frac{g}{L}))\\
\end{align}
</math>
 
 
You can solve for the same thing from the cartesian coordinates. Taking:
 
 
<math>\begin{align}
        F_x &= T\sin(\theta) = 0\\
\text { and recognizing } T = mg\\
\end{align}</math>
\end{align}</math>
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