2
$\begingroup$

Hello I am having trouble trying to find the correct model for this coupled spring system. The scenario is the following we have: Ceiling - Spring - Mass(1) - Spring(2) - Mass(2) - Spring (3) - Mass(3) End.

I came up with the following system of differential equations in the 2nd order to model this problem.

$x_1^{''}=[-k_1x_1-k_2(x_2-x_1)-k_3(x_3-x_2)]/m_1$

$x_2^{''}=[-k_2(x_2-x_1)-k_3(x_3-x_2)]/m_2$

$x_3^{''}=-k_3(x_3-x_2)/m_3$

Is this the correct model? Afterwards I am trying to linearize these equations into 6 differential equations that I can input in matlab and plot the position of each spring.

So I linearized them and obtained the following:

$y_1^{'}=y_2$

$y_2^{'}=(-k_1y_1-k_2(y_3-y_1)-k_3(y_5-y_3)/m_1$

$y_3^{'}=y_4$

$y_4^{'}=(-k_2(y_3-y_1)-k_3(y_5-y_3)/m_2$

$y_5^{'}=y_6$

$y_6^{'}=(-k_3(y_5-y_3))/m_3$

I am not sure if this is correct or not. When I plot them in matlab I dont get a sinusoidal wave. A big plus if you guys can tell me how I could animate this system in matlab so that I can see the change in position in all three of the springs.

$\endgroup$
  • $\begingroup$ You aren't expecting simple sinusoids except in a few special cases. $\endgroup$ – dmckee Aug 7 '14 at 2:53
  • $\begingroup$ What do you mean? So the model is correct? $\endgroup$ – adam Aug 7 '14 at 3:02
  • 1
    $\begingroup$ I have no idea if the model is correct or not, but not finding simple sinusoids does not, in and of itself, point to a bug. There are a few regular modes, but finding them is an eigenvalue problem and I am unsure if you know that term and what it implies. $\endgroup$ – dmckee Aug 7 '14 at 3:44
  • $\begingroup$ I know what eigenvalues are but ususally spring mass systems eigenvalues are complex. $\endgroup$ – adam Aug 7 '14 at 3:47
  • $\begingroup$ IIRC, this example is solved in Taylor's Classical Mechanics. $\endgroup$ – jinawee Aug 7 '14 at 8:50
3
$\begingroup$

From the free body diagram you must have

$$ \begin{align} m_1 \ddot{x}_1 &= F_1 - F_2 \\ m_2 \ddot{x}_2 &= F_2 - F_3 \\ m_2 \ddot{x}_3 &= F_3 \end {align} $$

with the spring forces defined as

$$ \begin{align} F_1 & = -k_1 x_1 \\ F_2 & = -k_2 (x_2-x_1) \\ F_3 & = -k_3 (x_3-x_1) \end{align} $$

The above is combined as

$$ \begin{bmatrix} m_1 & 0 & 0 \\ 0& m_2 & 0 \\ 0 & 0 & m_3 \end{bmatrix} \begin{pmatrix} \ddot{x}_1 \\ \ddot{x}_2 \\ \ddot{x}_3 \end{pmatrix} = - \begin{bmatrix} k_1 + k_2 & -k_2 & 0 \\ -k_2 & k_2 + k_3 & -k_3 \\ 0 & -k_3 & k_3 \end{bmatrix} \begin{pmatrix} {x}_1 \\ {x}_2 \\ {x}_3 \end{pmatrix}$$

Which I think matches your equations (you have to check).

To make an ODE out of this you need a state vector

$$ y = \begin{pmatrix} x_1 \\ x_2 \\ x_3 \\ v_1 \\ v_2 \\ v_3 \end{pmatrix} $$

and its derivative

$$ \dot{y} = A\,y $$ $$ \begin{pmatrix} \dot{x}_1 \\ \dot{x}_2 \\ \dot{x}_3 \\ \dot{v}_1 \\ \dot{v}_2 \\ \dot{v}_3 \end{pmatrix} = \begin{bmatrix} 0 & 0 & 0 & 1 & 0 & 0 \\ 0 & 0 & 0 & 0 & 1 & 0 \\ 0 & 0 & 0 & 0 & 0 & 1 \\ - \frac{k_1+k_2}{m_1} & \frac{k_2}{m_1} & 0 & 0 & 0 & 0 \\ \frac{k_2}{m_2} & - \frac{k_2+k_3}{m_2} & \frac{k_3}{m_2} & 0 & 0 & 0 \\ 0 & \frac{k_3}{m_3} & - \frac{k_3}{m_3} & 0 & 0 & 0 \end{bmatrix} \begin{pmatrix} x_1 \\ x_2 \\ x_3 \\ v_1 \\ v_2 \\ v_3 \end{pmatrix}$$

As long as $\ddot{x}_i \propto - x_i $ there would be a harmonic response. If you not seeing this, then there is something wrong on how you are using ode45().

$\endgroup$

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Not the answer you're looking for? Browse other questions tagged or ask your own question.