Skip to main content

Section 7.5 Center of Mass Motion

Subsection 7.5.1 Center of Mass Velocity, Acceleration, and Momentum

Center of mass of an object provides an important way to follow the overall motion of the object. We have already seen the formula for the position of the center of mass. For a body that is a collection of discrete masses \(m_1\) at \(\vec r_1\text{,}\) \(m_2\) at \(\vec r_2\text{,}\) \(\cdots\text{,}\) \(m_N\) at \(\vec r_N\text{,}\) the center of mass will be at \(\vec R_\text{CM}\text{,}\) given by

\begin{equation} \vec R_\text{CM} = \frac{1}{M}\left( m_1\vec r_1 + m_2\vec r_2 + \cdots + m_N\vec r_N \right),\label{eq-cm-position-N-particles}\tag{7.5.1} \end{equation}

where total mass \(M=m_1+m_2+\cdots+m_N\text{.}\) Now, studying the rate at which position \(\vec R_\text{CM}\) changes will give us the velocity of center of mass, \(\vec V_{\text{CM}}\text{.}\)

\begin{equation} \vec V_{\text{CM}} = \dfrac{m_1\vec v_1 + m_2\vec v_2 + \cdots + m_N\vec v_N }{M}.\label{eq-cm-velocity-non-calculus}\tag{7.5.2} \end{equation}

Similarly, when we look at the rate which \(\vec V_{\text{CM}}\) changes, we get the acceleration of the center of mass, \(\vec A_{\text{CM}}\text{.}\)

\begin{equation} \vec A_{\text{CM}} = \dfrac{m_1\vec a_1 + m_2\vec a_2 + \cdots + m_N\vec a_N }{M}. \label{eq-cm-acceleration-non-calculus}\tag{7.5.3} \end{equation}

The total momentum of all particles will be

\begin{equation*} \vec p_\text{total} = m_1\vec v_1 + m_2\vec v_2 + \cdots + m_N\vec v_N. \end{equation*}

Now, this is just \(M \vec V_{\text{CM}}\) as you can see from Eq. (7.5.2) when you multiply it by \(M\text{.}\)

\begin{equation*} M \vec V_{\text{CM}} = m_1\vec v_1 + m_2\vec v_2 + \cdots + m_N\vec v_N \end{equation*}

We call the product \(M \vec V_{\text{CM}}\) the center of mass momentum, \(\vec P_{\text{CM}}\text{,}\) which is actually total momentum of all particles in the body.

\begin{equation} \vec P_{\text{CM}} = M \vec V_{\text{CM}} = \vec p_1 + \vec p_2 + \cdots + \vec p_N = \vec p_\text{total}. \tag{7.5.4} \end{equation}

Subsection 7.5.2 Center of Mass Motion - Constant Mass Case

Let us consider a system with \(N \) parts of masses \(m_1 \text{,}\) \(m_2 \text{,}\) \(\cdots \text{,}\) \(m_N \text{.}\) Let \(\vec F_1 \text{,}\) \(\vec F_2 \text{,}\) \(\cdots \text{,}\) \(\vec F_N \) be forces acting on those parts. For the sake of simplicity, let us look at the constant mass system, i.e., where masss of each part is fixed in time.

Now, we write Newton's equation of motion for each part for constant mass.

\begin{align*} \amp \vec F_1 = m_1 \vec a_1 \\ \amp \vec F_2 = m_2 \vec a_2 \\ \amp \cdots \\ \amp \vec F_N = m_N \vec a_N \end{align*}

Summing them, we find that the left side can be written as sum over all the forces from objects outside the \(N \) parts making up the system, which we call the external forces, and a sum over all the forces that one part applies on another part.

\begin{equation*} \text{Left Side } = \vec F_1 + \vec F_2 + \cdots + \vec F_N = \vec F_{\text{net}}^{\text{ext}} + \vec F_{\text{net}}^{\text{int}} \end{equation*}

By Newton's third law, internal forces, say by part 1 on 2 and by part 2 on 1 come in pairs whose sum is zero. Therefore,

\begin{equation*} \vec F_{\text{net}}^{\text{int}} = 0. \end{equation*}

That gives us the following simplified equation after summing the equations of all \(N \) parts.

\begin{equation*} F_{\text{net}}^{\text{ext}} = m_ 1\vec a_1 + m_2\vec a_2 + \cdots + m_N\vec a_N, \end{equation*}

whose right side is just the total mass \(M \) times the center-of-mass acceleration \(\vec A_{\text{CM}} \text{.}\)

\begin{equation} F_{\text{net}}^{\text{ext}} = M \vec A_{\text{CM}}.\label{eq-cm-eom-fext-ma}\tag{7.5.5} \end{equation}

It is interesting that the overall motion of a malti-part system can be captured essentially by the motion of a single (fictitious) particle of total mass placed at the special point of center of mass. This is the equation we solve whenever we study the overall translation motion of finite-size objects since they are inherently mutiparticle systems.

General Case (Calculus):

We have seen before that the general case of mass varying as well as constant mass systems is given in terms of rate of change of momentum.

\begin{equation*} \vec F = \frac{d\vec p}{dt}. \end{equation*}

Now, if you apply this to the \(N\)-particle system, you will get

\begin{equation*} F_{\text{net}}^{\text{ext}} = \frac{d\vec P_\text{CM}}{dt} = \frac{d\vec p_\text{total}}{dt}. \end{equation*}

Similarly, the integrated form in terms of impulse and change in momentum will be

\begin{equation} \vec J_{\text{net}}^{\text{ext}}= \int F_{\text{net}}^{\text{ext}} dt = \vec p_{\text{total},f} - \vec p_{\text{total},i}\tag{7.5.6} \end{equation}

Subsection 7.5.3 (Calculus) Center of Mass Velocity, Acceleration, and Momentum

(You can safely skip this section)

To study the motion of the center of mass (CM) of a multi-part system, we need to define velocity and acceleration of the CM. We use Calculus to do that. Recall that velocity is defined as the time derivative of position and acceleration the time derivative of velocity. We do the same with CM-position in Eq. (7.5.1).

\begin{align} \vec V_{\text{CM}} \amp = \dfrac{1}{M}\dfrac{d}{dt}\left( m_1\vec r_1 + m_2\vec r_2 + \cdots + m_N\vec r_N \right), \label{eq-cm-velocity}\tag{7.5.7}\\ \amp =\dfrac{m_1\vec v_1 + m_2\vec v_2 + \cdots + m_N\vec v_N }{M}, \tag{7.5.8} \end{align}

where \(\vec v_1 \text{,}\) \(\vec v_2 \text{,}\) \(\cdots \text{,}\) \(\vec v_N \text{,}\) are the velocities of the individual parts, or more precisely, the velocities of the centers of masses of the parts. Therefore, we find that the total momentum, i.e., the sum of the momenta of parts, is equal to the total mass times the momentum of the center of mass.

\begin{equation} m_1\vec v_1 + m_2\vec v_2 + \cdots + m_N\vec v_N = M\vec V_{\text{CM}} \equiv \vec P_{\text{CM}}. \label{eq-cm-momentum-equals-total-momentum}\tag{7.5.9} \end{equation}

We call the \(M\vec V_{\text{CM}} \) the center-of-mass momentum, \(\vec P_{\text{CM}}\text{.}\) Writing momentum of the parts as \(\vec p_1\text{,}\) \(\vec p_2\text{,}\) etc, we get

\begin{equation} \vec p_{\text{tot}} = \vec p_1 + \vec p_2 + \cdots = \vec P_{\text{CM}}. \label{eq-cm-momentum-equals-total-momentum-2}\tag{7.5.10} \end{equation}

Taking a derivative of \(\vec V_{\text{CM}}\) with respect to \(t \) will give the acceleration of the CM of a system whose mass is constant in time,

\begin{align} \vec A_{\text{CM}} \amp = \dfrac{1}{M}\dfrac{d}{dt}\left( m_1\vec v_1 + m_2\vec v_2 + \cdots + m_N\vec v_N \right), \label{eq-cm-acceleration}\tag{7.5.11}\\ \amp = \dfrac{m_1\vec a_1 + m_2\vec a_2 + \cdots + m_N\vec a_N }{M}, \tag{7.5.12} \end{align}

where \(\vec a_1 \text{,}\) \(\vec a_2 \text{,}\) \(\cdots \text{,}\) \(\vec a_N \text{,}\) are the accelerations of the individual parts.

Subsection 7.5.4 (Calculus) Center of Mass Motion - General Case

(You can safely skip this section)

For the general case of \(N \) part system, we will have to start with equations of motion of individual parts in the momentum form.

\begin{align*} \amp \vec F_1 = \dfrac{d\vec p_1}{dt} \\ \amp \vec F_2 = \dfrac{d\vec p_2}{dt} \\ \amp \cdots \\ \amp \vec F_N = \dfrac{d\vec p_N}{dt} \end{align*}

Now, summing them, we get the left side to equal the \(F_{\text{net}}^{\text{ext}} \) as above, and the right side equals the rate of change of total momentum.

\begin{equation*} \vec F_{\text{net}}^{\text{ext}} = \dfrac{d}{dt} \left( \vec p_1 + \vec p_2 + \cdots + \vec p_N \right). \end{equation*}

But by Eq. (7.5.10), we can replace the sum of momenta of the parts by the center-of-mass momentum, or the total momentum.

\begin{equation} \vec F_{\text{net}}^{\text{ext}} = \dfrac{d \vec P_{\text{CM}}}{dt} = \dfrac{d \vec p_{\text{tot}}}{dt},\label{eq-cm-eom-fext-dpdt}\tag{7.5.13} \end{equation}

where \(\vec P_{\text{CM}} = M \vec V_{\text{CM}}\text{.}\) This equation takes care of even those situations in which the total mass of the system is changing, regardless of which part that may be gaining or losing mass. Of course, if a system is not gaining or losing mass, we would use much simpler analysis based on Eq. (7.5.5).

Eq. (7.5.13) can be integrated to yield the following important relation between the change in total momentum and impulse by external forces.

\begin{equation} \vec p_{\text{tot},f} - \vec p_{\text{tot},i} = \vec J_{\text{net}}^{\text{ext}},\tag{7.5.14} \end{equation}

where

\begin{equation} \vec J_{\text{net}}^{\text{ext}} = \int_{t_i}^{t_f}\, \vec F_{\text{net}}^{\text{ext}}\, dt.\tag{7.5.15} \end{equation}

As an example consider the motion of pieces of an exploded shell in free-fall.

Although the pieces fly away in various directions, the CM falls as if no explosion had taken place since during the explosion all forces are internal to the system. The motion of CM is that of a single particle with total mass of the system. The fictitious particle at the CM falls freely with acceleration equal to g, the acceleration due to gravity, until one of the pieces hits the ground.

At the instant first piece hits the ground, there will be a new upward impulse on the collective body. That is, the CM's free fall motion will end.

Suppose you are trying to practice a vertical jump. You crouch down and push against the ground.

The magnitude of the normal force from the ground on you varies with time as shown in Figure 7.5.3 of \(F\) versus \(t\text{.}\) The initial normal force is just your weight. As you push down, the normal force increases. You leave the ground at the point normal force is zero.

Figure 7.5.3. Figure for Checkpoint 7.5.2.

As a net result of normal force from the ground, and weight, your center of mass has acceleration upward. What will be your speed at the instant you are no longer in contact with the ground?

Hint

Add the area under the curve to get impulse from the normal force. Then, subtract the impulse from weight since it is in opposite direction. Use the net impulse to deduce the change in momentum of the center of mass.

Answer

\(1.47\text{ m/s}.\)

Solution

Let positive \(y\) axis be pointed up. We will implement the following law for the center of mass since CM momentum is equal to \(p_\text{tot}\text{.}\)

\begin{equation*} J_y^{\text{net}} = \Delta p_y^{CM}. \end{equation*}

To get the net impulse, we need to vectorially add up impulses from all forces on the body. We have two forces on the body here - the weight and the normal force. From \(t=0\) point in Figure 7.5.3, we note that magnitude of weight is \(500\text{ N}\text{.}\) Since weight has negative \(y\) component, its \(J_y\) will be negative.

\begin{equation*} J_y^{\text{weight}} = - 500\text{ N}\times 0.200\text{ s} = -100\text{ N.s}. \end{equation*}

The impulse from the normal force is the area under the given plot.

\begin{align*} J_y^{\text{normal}} \amp = \left( 500\times 0.2 - \dfrac{500\times 0.02}{2}\right) \\ \amp \ \ \ \ + \dfrac{1000\times 0.08}{2} + 1000\times 0.02 + \dfrac{1000 \times 0.04}{2}, \\ \amp = 175\text{ N.s}, \end{align*}

which is obtained by slicing the area into triangles and rectangle. Therefore,

\begin{equation*} J_y^{\text{net}} = 175\text{ N.s} - 100\text{ N.s} = 75\text{ N.s}. \end{equation*}

The mass of the person is obtained from the weight.

\begin{equation*} m = \dfrac{500\text{ N}}{9.81\text{ m/s}^2} = 51\text{ kg}. \end{equation*}

Initially, the momentum was zero. Therefore, the velocity at the end is

\begin{equation*} v_{\text{cm}} = \dfrac{75\text{ N.s}}{51\text{ kg}} = 1.47\text{ m/s}. \end{equation*}

A \(78.0\text{-kg}\) athlete jumps straight up so that his CM moves by a height of \(1.2\text{ m}\) after he leaves the ground. Find the magnitude of impulse from the ground.

Hint

Impulse will equal \(mv_0\) where \(v_0\) is the velocity of CM just before the athlete leaves ground.

Answer

\(378.3\text{ N.s}\text{.}\)

Solution

The jump is due to push by the athlete on the ground which causes a normal force that exceeds his weight and gives upward momentum when he leaves the ground. This force from the ground during the jumping process is responsible for imparting initial upward momentum.

Let \(v_0\) be the velocity pointed upwards at the instant athlete leaves the ground and \(m\) mass of the athlete. Then, the change in momentum from zero to \(mv_0\) will equal the the impulse by ground. Let us work with \(y\) axis pointed up. Then, \(y\)-component of impulse will be

\begin{equation*} J_y = \Delta p_y = mv_0., \end{equation*}

The CM starts with velocity \(v_0\) and has constant acceleration under the external force of gravity so that its velocity is zero at maximum height \(h\text{.}\) Using constant acceleration of the motion of CM along \(y\) axis we get

\begin{equation*} 0 - v_0^2 = - 2 g h\ \ \longrightarrow\ \ v_0 = \sqrt{2gh}. \end{equation*}

Hence, impulse by ground is

\begin{equation*} J = m\sqrt{2gh} = 78\sqrt{2\times 1.2\times 9.8} = 378.3\text{ N.s}. \end{equation*}

A diver jumps off a board with initial speed \(12\text{ m/s}\) at an angle \(60^\circ\) with the horizontal direction. She falls a total height of \(8\text{ m}\) before entering water. Find the speed of her center of mass when she enters water?

Hint

The CM moves as if it were a particle in free fall.

Answer

\(17.4\text{ m/s}\text{.}\)

Solution

Since the only force on her body after she takes off is her weight. The CM will obey the free fall equations. Therefore, we start with components of her initial velocity.

\begin{align*} v_{ix} \amp = 12\,\cos\,60^{\circ} = 6\text{ m/s}, \\ v_{iy} \amp = 12\,\sin\,60^{\circ} = 10.4\text{ m/s}. \end{align*}

The final \(v_x\) is same as the initial, and the final \(v_y\) is obtained from

\begin{equation*} v_{fy} = - \sqrt{ v_{iy}^2 - 2 g \Delta y}, \end{equation*}

where \(\Delta y = - 8\text{ m}\text{,}\) which is negative since the final location is below the initial point, and \(v_{fy}\) is negative since it is pointed down while positive \(y\) axis is pointed up.

\begin{equation*} v_{fy} = -\sqrt{ 10.4^2 + 2\times9.81\times 8} = -16.3\text{ m/s}. \end{equation*}

From \(v_{fx}\) and \(v_{fy}\) we get the speed just before she enters water.

\begin{equation*} v_f = \sqrt{ v_{fx}^2 + v_{fy}^2 } = \sqrt{6^2 + 16.3^2} = 17.4\text{ m/s}. \end{equation*}

Two blocks of mass \(m\) each is attached by a spring of spring constant \(k\) and unstretched length \(l_0\text{.}\) The blocks are placed on a horizontal frrictionless track.

Figure 7.5.7.

While holding one of the blocks, say B in the figure, the other block is pushed in so that spring compresses by \(D\text{,}\) and then everything is released from rest at \(t=0\text{.}\) Use \(x\)-axis so that at \(t=0\text{,}\) \(x_A=D\) and \(x_B=l_0\text{.}\)

Find the position of the center of mass, \(x_\text{cm}\text{,}\) at an arbitrary instant \(t\text{.}\)

Hint

Since there is no external force, \(x_\text{cm}\) will remain fixed in time. You can show this by a detailed calculation.

Answer

\(x_\text{cm} (t)= \frac{D+l_0}{2}\text{.}\)

Solution

Since there is no external force, \(x_\text{cm}\) will remain fixed in time. Since the center of mass is at \(x=(D+l_0)/2\) at \(t=0\text{,}\) it will remain there. This does not mean A and B do not move - they do, but their CM remains at the same point in space. We can also show this by a detailed calculation as follows.

Figure 7.5.8 shows situation at an arbitrary instant when the spring is stretched. The spring forces on the two blocks are shown.

Figure 7.5.8.

The \(x\)-equations of motion for the two blocks will be

\begin{align} \amp m_A\frac{d^2x_A}{dt^2} = k (x_B - x_A - l_0) \label{eq-two-blocks-atttached-by-a-spring-eom-A}\tag{7.5.16}\\ \amp m_B\frac{d^2x_B}{dt^2} = -k (x_B - x_A - l_0) \label{eq-two-blocks-atttached-by-a-spring-eom-B}\tag{7.5.17} \end{align}

Here \(m_A=m_B=m\) and \(x_B \gt x_A\text{.}\) Since center of mass is give by

\begin{equation*} x_\text{cm} = \frac{1}{2}\left( x_A + x_B \right), \end{equation*}

we can take two derivatives with respect to time to get

\begin{equation*} \frac{d^2x_\text{cm}}{dt^2}= \frac{1}{2}\left( \frac{d^2x_A}{dt^2} + \frac{d^2x_B}{dt^2} \right). \end{equation*}

Using the equations of motions for \(x_A\) and \(x_B\text{,}\) we immediately get

\begin{equation*} \frac{d^2x_\text{cm}}{dt^2} = 0. \end{equation*}

This has a simple solution.

\begin{equation} x_\text{cm}(t) = a t + b,\label{eq-wo-blocks-atttached-by-a-spring-cm-gen-sol}\tag{7.5.18} \end{equation}

where \(a\) and \(b\) will be fixed from the initial conditions.

\begin{equation*} x_A=D,\ \ x_B = l_0,\ \ v_A = 0 ,\ \ v_B=0. \end{equation*}

Writing Eq. (7.5.18) in terms of \(x_A\) and \(x_B\) and also writing their derivative we have

\begin{align*} \amp \frac{1}{2}\left( x_A + x_B \right) = a t + b\\ \amp \frac{1}{2}\left(v_A + v_B \right) = a \end{align*}

Evaluating these at \(t=0\) we get

\begin{equation*} a = 0,\ \ b = \frac{D+l_0}{2}. \end{equation*}

Therefore, we have the following answer for the center of mass at arbitrary instant \(t\text{.}\)

\begin{equation} x_\text{cm} = \frac{D+l_0}{2}.\label{eq-two-blocks-atttached-by-a-spring-eom-cm-ans}\tag{7.5.19} \end{equation}

That is, CM remains fixed at \(x= (D+l_0)/2\text{.}\)

In Checkpoint 7.5.6, the two blocks moved freely on the frictionless horizontal track.

Now, suppose you were pushing block A, you kept block B against a wall as shown in Figure 7.5.10.

Figure 7.5.10.

After compressing the spring by amount \(D\text{,}\) you release both blocks at rest. What will be the motion of the center of mass if \(l_0\) is the natural length of the spring?

Hint

First solve the situation when block B is in contact with the wall. Then solve the situation after B loses contact with the wall. When in contact, set up equation of motion for the CM and another equation of motion for B. Also, use the relation between \(x_A\) and \(x_\text{cm}\text{.}\)

Answer

Solution

Overall strategy:

The motion of the center of mass will be controlled by external forces on the two-block (and spring) system. When block B is in contact with the wall, the net external force on the two-block system comes from the normal force of the wall, which will vary with the varying length of the spring. When block B is no longer in contact with the wall, the normal force from the wall will be zero, which will make net external force zero on the two block system. Thus, after block B is no longer in contact with the wall, the center of mass will just move with a constant velocity.

Block B in Contact with the Wall:

Figure 7.5.11 shows situation at an arbitrary instant when block B is still in contact with the wall. The equation of motion of the center of mass will be

\begin{equation*} M \frac{d^2x_\text{cm}}{dt^2} = - N(t). \end{equation*}

Figure 7.5.11.

We need an expression for \(N(t)\) to solve this equation. We get that by looking at force on B as shown in Figure 7.5.12. Balancing spring and wall forces on B gives

\begin{equation*} N(t) = k (l_0 - x_A). \end{equation*}

Figure 7.5.12.

We also note that when B is at the wall \(x_\text{cm}\) and \(x_A\) are related by

\begin{equation*} x_\text{cm} = \frac{m x_A + m l_0}{2m} = \frac{x_A + l_0}{2}. \end{equation*}

Using the three equations above, we get the following equation of motion for \(x_\text{cm}\text{.}\)

\begin{equation*} \frac{d^2x_\text{cm}}{dt^2} = \frac{k}{m}\left(x_\text{cm} - l_0\right). \end{equation*}

This has the following general solution.

\begin{equation} x_\text{cm}(t) = l_0 + A\, e^{\sqrt{k/m}\, t} + B\, e^{-\sqrt{k/m}\, t}.\label{eq-block-against-wall-cm-pos}\tag{7.5.20} \end{equation}

The velocity of CM will be

\begin{equation} v_\text{cm}(t) = \sqrt{k/m}\,A\, e^{\sqrt{k/m}\, t}- \sqrt{k/m}\,B\, e^{-\sqrt{k/m}\, t}.\label{eq-block-against-wall-cm-vel}\tag{7.5.21} \end{equation}

Using the initial conditions on A and B, we have the following intial conditions on CM.

\begin{align*} \amp x_\text{cm}(t) = \frac{x_A (0)+ l_0}{2} = \frac{D+ l_0}{2}.\\ \amp v_\text{cm}(t) = \frac{v_A (0)}{2} = 0. \end{align*}

Using these in Eqs. (7.5.20) and (7.5.21), we find

\begin{equation*} A = B = \frac{D-l_0}{4}. \end{equation*}

Therefore, the position of CM and its velocity are

\begin{align} \amp x_\text{cm}(t) = l_0 - \left(\frac{l_0-D}{2} \right)\,\cosh(\sqrt{k/m}\,t). \label{eq-block-against-wall-cm-pos-fin}\tag{7.5.22}\\ \amp v_\text{cm}(t) = -\frac{(l_0-D)\sqrt{k/m}}{2}\,\sinh(\sqrt{k/m}\,t).\label{eq-block-against-wall-cm-vel-fin}\tag{7.5.23} \end{align}

Block B not in Contact with the Wall:

At the instant, block B loses contact with the wall, the CM will be at \(l_0/2\text{.}\) At that instant onward the velocity of the CM will remain unchanged since there are net external force on the two-block system is zero now. We can find the time when this happens from Eq. (7.5.22). Let this instant be at \(t=T\text{.}\)

\begin{equation*} \frac{l_0}{2} = l_0 - \left(\frac{l_0-D}{2} \right)\,\cosh(\sqrt{k/m}\,T) \end{equation*}

Therefore,

\begin{equation*} T = \sqrt{\frac{m}{k}}\, \cosh^{-1}\left( \frac{l_0}{l_0 - D} \right). \end{equation*}

From Eq. (7.5.23), the velocity of CM after this instant will be constat at

\begin{align*} v_\text{cm}(t>T)\amp = - \frac{(l_0-D)\sqrt{k/m}}{2} \,\sinh(\sqrt{k/m}\,T)\\ \amp = -\frac{(l_0-D)\sqrt{k/m}}{2}\,\sinh\left[ \cosh^{-1}\left( \frac{l_0}{l_0 - D} \right) \right]. \end{align*}

This expression can be simplified by noting that

\begin{equation*} \sinh(\cosh^{-1}(x)) = \sqrt{x^2 - 1}. \end{equation*}

To prove this try \(y=\cosh^{-1}(x)\) and solve for \(e^y = x + \sqrt{x^2-1}\) and \(e^{-y} = x - \sqrt{x^2-1}\text{.}\) Using this, we simplify CM velocity for \(t\gt T\) to be

\begin{equation*} v_\text{cm}(t>T) = -\frac{1}{2}\sqrt{ \frac{k}{m} (2l_0 D - D^2) }. \end{equation*}

The velocity is \(x\)-component and pointed towards negative \(x\)-axis. That is why there is a negative sign in the front.