Coupled Horizontal Spring Mass Oscillator: Difference between revisions

From Class Wiki
Jump to navigation Jump to search
Line 120: Line 120:


<math>\tilde{z}=\tilde{T}\tilde{x}\,</math>
<math>\tilde{z}=\tilde{T}\tilde{x}\,</math>
<math>\tilde{x}=\tilde{T}^{-1}\tilde{z}\,</math>
<math>\text {Where }\,</math>
<math>\tilde{T}^{-1}=\,</math>

Revision as of 19:44, 9 December 2009

Coupled Oscillator Spring Mass Oscillator: State Space

Problem Statement

Two 4 Kg Weights are suspended between two walls. They are connected by a spring between them with a spring constant k2. They are connected to the walls by two springs k1 and k3 with k1=k3. m1 is a distance x1 form m2 and m2 is x2 from the wall.


Solution

Things we know

m1=5kg

m2=5kg

k1=50Nm

k2=100Nm

k3=50Nm

So now that we have are problem we need to start setting up the equations we need to solve it.

x1˙=x1˙

x1¨+k1+k2m1x1k2m1x2=0

x2˙=x2˙

x2¨+k3+k2m2x2k2m2x1=0

Now we take these equations and put them in a state space model.

[x1˙x1¨x2˙x2¨] = [0100(k1+k2)m10k1m100001k1m20(k1+k2)m20][x1x˙1x2x˙2]+[0]

Now we make the appropriate numerical substitutions.


[x1˙x1¨x2˙x2¨] = [01001505050500001505015050][x1x˙1x2x˙2]+[0]



[x1˙x1¨x2˙x2¨]=[01003001000001100300][x1x˙1x2x˙2]+[0]

So using maple i was able to obtain the eigenvalues and eigenvectors.

Eigenvalues.


λ1=210 λ2=210 λ3=25 λ4=25

Eigenvectors.

K1=[12(10)12(10)],K2=[12(10)12(10)],K3=[12(5)12(5)],K4=[12(5)12(5)]

So then the answer is...

x=c1[12(10)12(10)]e210+c2[12(10)12(10)]e2*210+c3[12(5)12(5)]e3*25+c4[12(5)12(5)]e4*25


Solve with the Matrix exponential

So first we need to know what the matrix exponential equation looks like.

it is...


x~=eA~tx(0)~

Where A is a matrix

Also 

z~=T~x~

x~=T~1z~

Where 

T~1=