How do you find the position function x(t) if you suppose that the mass in a mass-spring-dashpot system with m=25, c=10 and k=226 is set in motion with x(0)=20 and x'(0)=41?

1 Answer
Aug 9, 2016

x = (21.1062 cos(3t)+ 15.0737sin(3t))e^{-t/5}-1.10619

Explanation:

In those considerations that follow, the positive axis is oriented from down-up for forces, velocities and displacements.

Supposing the mass m vertically placed over the spring-dashpot, the mass is actuated by

f-m g = m ddot x

where f is the force exerted by the spring-dashpot.

The spring-dashpot reacts according to

-f = kx + c dot x

Joining both equations we have

m ddot x + c dot x + k x = -mg

This is a linear non-homogeneous differential equation whose solution is given by

x = x_h + x_p

where x_h is the homogeneous solution

m ddot x_h + c dot x_h + k x_h = 0

and

x_p is a particular solution for

m ddot x_p + c dot x_p + k x_p = -mg

In this case x_p = -(mg)/k

The homogeneous equation has the general solution

x_h = e^{lambda t}

Substituting this generic solution we get at

(m lambda^2+c lambda + k)e^{lambda t} = 0

The feasible lambda are given by

lambda=(-c pm sqrt(c^2-4mk))/(2m). For our case we have

lambda = (-10 pm sqrt(100-4 xx 25 xx 226))/(50) = -1/5 pm 3i

so

x_h = e^{-t/5} e^{pm i 3 t} = e^{-t/5}(cos(3t)pm i sin(3t))

when lambda is a complex conjugate pair, for each pair the solution proposal is as follows.

x_h = (C_1 cos(3t)+ C_2sin(3t))e^{-t/5}

(This can be demostrated but requires a lot of algebra.)

The general solution is given by

x = (C_1 cos(3t)+ C_2sin(3t))e^{-t/5}-(m g)/k The constants C_1,C_2 are determined according to the movement initial conditions. So

{(x(0) = C_1-(m g)/k = 20), (dot x(0) = -C_1/5+3C_2=41) :}

From this system, we obtain the constants

C_1 = (m g)/k + 20, C_2 = (m g)/(15k) + 15 assuming g = 10 we have

C_1 = 21.1062, C_2 = 15.0737

So finally

x = (21.1062 cos(3t)+ 15.0737sin(3t))e^{-t/5}-1.10619