syst. eq
state var
syst. sol
lapl. transf
stability
plots
glossary

 

Laplace Transforms

Case 1:
Solution using Laplace Transform

Mz'' + Bz' + Kz = KU(t)
ML{z''} + BL{z'} + KL{z} = KL{1}
z'(0) = z(0) = 0

Therefore:

Ms''Z(s) + BsZ(s) + KZ(s) = K/s
Z(s)[ Msē + Bs + K ] = K/s
Z(s) = K/s/( Msē + Bs + K)

The roots of (Msē + Bs + K) are calculated above.

r1 = - 0.14 + 0.87i
r2 = - 0.14 - 0.87i

Z(s) = K/s(s-r1)(s-r2)
or
Z(s) = K/s(s + 0.14 - 0.87i )*( s + 0.14 + 0.87i )
(K = 2000)

Using partial fraction expansion:

Z(s) = C1/s + C2/(s + 0.14 - 0.87i ) + C3/(s + 0.14 + 0.87i )


C1=sZ(s)|s=0= 0
C2 =(s + 0.14 + 0.87i )Z(s)| s=-0.14-0.87i
C3 =(s + 0.14 - 0.87i )Z(s)|s=-0.14+0.87i

C1 = K/(r1*r2) = 1
C2 = K/s(s-r2) = k/r1(r1-r2) = -1287.8 + 207.24i
C3 = K/s(s-r1)=K/r2(r2-r1) = -1287.8 - 207.24i

Write C1 and C2 in exponential form.:
C1 = 1304exp(0.16i)
C2 = 1304exp(-0.16i)

Z(s) = 2576/s + [1304exp(-0.16i)] / (s + 0.14 - 0.87i ) + [1304exp(0.16i)] / (s + 0.14 + 0.87i )

z(t) = Inverse Laplace Transform of Z(s)

z(t) = 1+ 0.98exp(-0.14t)* cos(0.87t - 0.16)

Previous
Next


Next topic:

Previous topics:

All graphs are included in the above pages, to see graphs ONLY (no solutions) go here.

 


Home | History | About us | System Analysis | Sitemap | Email
This is an educational website designed and maintenaned by Julia Maslov