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)


with this new equation we get:
<math>\begin{align}
 
        F_x &= T\sin(\theta) = 0\\
 
\text { and recognizing } T = mg\\
g*\theta(t)/L +s^2*\theta(t) - s*\theta(0) - \theta'(0) = 0
\end{align}</math>
 
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




\theta(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