Skip to main content

State space model of vertical two-mass system

Short answer

Define states x1=y1, x2=dy1/dt, x3=y2, x4=dy2/dt. The state-space model is x' = Ax + Bu with A = [0 1 0 0; -K/M1 -B/M1 0 B/M1; 0 0 0 1; 0 B/M2 0 -B/M2], B = [0; 1/M1; 0; 0], y = [0 0 1 0]x.

Consider a mechanical system shown in figure below. The external force u(t) is input to the system and displacement y(t) of the mass is the output. Obtain the state space representation of the system. [Figure: vertical system; mass M1 hangs from a fixed support through spring K, with external force u(t) applied to M1; mass M2 hangs below M1 through damper B; y(t) is the downward displacement of M2; a fixed surface is drawn below M2]

Let y1y_1 be the downward displacement of M1M_1 and y=y2y = y_2 that of M2M_2, both measured from the static equilibrium position, so gravity and the static spring force cancel and do not appear.

   //////////  support
       |
       K  (spring)
       |
   +-------+
   |  M1   | <-- u(t), y1
   +-------+
       |
       B  (damper)
       |
   +-------+
   |  M2   |  y(t) = y2
   +-------+

Equations of motion

For M1M_1: spring force Ky1Ky_1 and damper force B(y˙1−y˙2)B(\dot{y}_1 - \dot{y}_2) oppose its motion:

M1y¨1+B(y˙1−y˙2)+Ky1=u(t)M_1\ddot{y}_1 + B(\dot{y}_1 - \dot{y}_2) + Ky_1 = u(t)

For M2M_2: only the damper acts on it:

M2y¨2+B(y˙2−y˙1)=0M_2\ddot{y}_2 + B(\dot{y}_2 - \dot{y}_1) = 0

State variables

x1=y1,x2=y˙1,x3=y2,x4=y˙2x_1 = y_1,\quad x_2 = \dot{y}_1,\quad x_3 = y_2,\quad x_4 = \dot{y}_2 x˙1=x2x˙2=−KM1x1−BM1x2+BM1x4+1M1ux˙3=x4x˙4=BM2x2−BM2x4\begin{aligned} \dot{x}_1 &= x_2 \\ \dot{x}_2 &= -\frac{K}{M_1}x_1 - \frac{B}{M_1}x_2 + \frac{B}{M_1}x_4 + \frac{1}{M_1}u \\ \dot{x}_3 &= x_4 \\ \dot{x}_4 &= \frac{B}{M_2}x_2 - \frac{B}{M_2}x_4 \end{aligned}

State-space representation

[x˙1x˙2x˙3x˙4]=[0100−KM1−BM10BM100010BM20−BM2][x1x2x3x4]+[01M100]u(t)\begin{bmatrix} \dot{x}_1 \\ \dot{x}_2 \\ \dot{x}_3 \\ \dot{x}_4 \end{bmatrix} = \begin{bmatrix} 0 & 1 & 0 & 0 \\ -\dfrac{K}{M_1} & -\dfrac{B}{M_1} & 0 & \dfrac{B}{M_1} \\ 0 & 0 & 0 & 1 \\ 0 & \dfrac{B}{M_2} & 0 & -\dfrac{B}{M_2} \end{bmatrix}\begin{bmatrix} x_1 \\ x_2 \\ x_3 \\ x_4 \end{bmatrix} + \begin{bmatrix} 0 \\ \dfrac{1}{M_1} \\ 0 \\ 0 \end{bmatrix}u(t) y(t)=[0010][x1x2x3x4]y(t) = \begin{bmatrix} 0 & 0 & 1 & 0 \end{bmatrix}\begin{bmatrix} x_1 \\ x_2 \\ x_3 \\ x_4 \end{bmatrix}

y2y_2 itself is needed as a state only because it is the output; the dynamics depend on y1y_1, y˙1\dot{y}_1 and y˙2\dot{y}_2.

✓ Answer checked by Claude Opus 5.5 (high effort)

All Control Systems chapter 8 questions →