Aaron Boyd's Assignment 8: Difference between revisions

From Class Wiki
Jump to navigation Jump to search
No edit summary
No edit summary
Line 20: Line 20:


<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}
\theta" - sin(\theta)(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(\theta) = \theta where \theta is small.  
\text{Unfortunately the laplace transform of that is horrible. Long, and impossible to solve.}\\
(I tried to leave sin(\theta) 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*\theta(t)/L +s^2*\theta(t) - s*\theta(0) - \theta'(0) = 0
\\
 
\text{(I tried to leave } sin(\theta) \text{ in the equation.}\\
we know that \theta(0) = \theta0 and \theta'(0) = 0
\text{After 4 hours and many wolframalpha.com timeouts I gave up) }\\
 
\\
 
\text{with this new equation we get:}\\
solving for \theta(t) we get  
\\
 
g*\frac{\theta(t)}{L +s^2}*\theta(t) - s\theta(0) - \dot\theta(0) = 0\\
 
\\
\theta(t) = s*\theta<sub>0</sub>/((-g/L)+S^2)
\text{we know that } \theta(0) = \theta_0 \text{ and } \dot\theta(0) = 0\\
 
\\
 
\\
now we take the inverse laplace transform of that which yields
\text{solving for } \theta(t) \text{ we get} \\
 
\\
 
\theta(t) = s*\frac{\theta_0}{(\frac{-g}{L}+S^2)}\\
\theta(t) = cosh(t*(g/L)^(1/2))
\\
\\
\text{now we take the inverse laplace transform of that which yields }\\
\\
\\
\theta(t) = cosh(t\sqrt(\frac{g}{L})\\
\end{align}</math>

Revision as of 11:40, 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 Fx=Tsin(θ)Fy=Tcos(θ)mg=0


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

now:


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)=cosh(t(gL)