Dynamics

Energy and Work

Work, Forces and Moments

Work done by a force $\vec{F}$.

\[\begin{aligned} W &= \int_{\vec{r}_1}^{\vec{r}_2} \vec{F} \cdot d\vec{r} \\ & = \int_{t_1}^{t_2} \vec{F} \cdot \dot{\vec{r}}dt \end{aligned}\]

Work done by a constant force $\vec{F}$.

\[ W = \vec{F} \cdot \Delta \vec{r} \]

The change in position is $\Delta \vec{r} = \vec{r}_2 - \vec{r}_1$.

Compute the work definition for a constant force by pulling $\vec{F}$ out of the integral, so the resulting trivial integral is just the difference between the endpoints.

Work done by a constraint force.

\[ W = 0 \]

The work $W$ is by a force of constraint force $\vec{F}$.

By definition a constraint force $\vec{F}$ does not permit movement in the direction that it applies. This means that the constraint force and the displacement are always orthogonal and the dot product between orthogonal vectors is zero.

Work done by a moment $M$.

\[\begin{aligned} W &= \int_{\theta_1}^{\theta_2} M \, d\theta \\ & = \int_{t_1}^{t_2} M \dot\theta \, dt \\ &= M \, \Delta \theta \text{ (for constant $M$)} \end{aligned}\]

The rotation angle $\theta$ is measured around the same axis $\vec{a}$ about which the moment $M$ is applied.

Work and Energy

Conservative vs non-conservative forces

A conservative force is one such that the work done by the force only depends on the initial and final positions of the object (i.e independent of the path taken by the object). Examples of conservative forces are gravity and spring force, where only the initial and final displacements are needed to calculate the work done by those forces. Examples of a non-conservative force are any type of drag forces (e.g friction, air resistance).

The diagram below illustrates different paths and the work done by a conservative forces vs. a non-conservative force.

If the force moving an object from point $A$ to point $B$ was conservative, the work done by that force in all paths is equal. If the force was non-conservative, the work done by that force is not equal, and it depends on the path length.

Potential Energy

Force and potential energy.

\[\begin{aligned} \vec{F} &= -\nabla V \\ & = -\frac{\partial V}{\partial x}\hat{\imath} -\frac{\partial V}{\partial y}\hat{\jmath} - \frac{\partial V}{\partial z}\hat{k} \end{aligned}\]

Beginning with #ren-ec and generalizing it to one-dimensional displacement in the y-direction, the equation becomes: \[\begin{aligned} W &= \vec{F} \cdot \Delta{\vec{y}} \end{aligned}\]

Assuming the only force applied is the force $\vec{F}$, then the system gains potential energy, so the amount of work done is equal to $\Delta{V}$. \[\begin{aligned} \Delta{V} &= \vec{F} \cdot \Delta{\vec{y}} \end{aligned}\]

Rearranging: \[\begin{aligned} F &= \frac{\Delta{V}}{\Delta{y}} \end{aligned}\]

Since we are interested in the force the field exerts on the system, due to Newton's Third Law, the force is equal and opposite. Therefore, the above result becomes: \[\begin{aligned} F &= -\frac{\Delta{V}}{\Delta{y}} \end{aligned}\]

If the change in potential energy and displacement is infinitesimally small, this becomes a derivative: \[\begin{aligned} F &= -\frac{dV}{dy} \end{aligned}\]

Generalizing this result to 3 dimensions, we use the gradient to obtain the following result: \[\begin{aligned} \vec{F} &= -\nabla V \\ & = -\frac{\partial V}{\partial x}\hat{\imath} -\frac{\partial V}{\partial y}\hat{\jmath} - \frac{\partial V}{\partial z}\hat{k} \end{aligned}\]

Example Problem: Gravitational potential energy.

Consider an object of mass $m$ in the air under the influence of gravity. Derive the potential energy $V_g$ of the object.

Using #ren-efp, we can derive the gravitational potential energy. The force due to gravity is: \[\begin{aligned} \vec{F} &= -mg\hat{\jmath} \end{aligned}\]

Since gravity only acts in one direction, we can simplify #ren-efp to: \[\begin{aligned} F &= -\frac{dV}{dy} \\ -mg &= -\frac{dV}{dy} \\ mg \, dy &= dV \\ \end{aligned}\]

Integrating both sides: \[\begin{aligned} \int_{V_1}^{V_2} dU &= \int_{y_1}^{y_2} mg \, dy \\ \Delta V & = mg\Delta y \end{aligned}\]

We can simplify the result by assuming $U_1$ is zero, and $y_1$ is at the reference zero-potential height. Therefore, the above becomes: \[\begin{aligned} V_g & = mgy = mgh \end{aligned}\]

The above result also proves that gravity is a conservative force, as it only depends on the initial and final position of the object.

Example Problem: Potential energy in a spring

Consider an object of mass $m$ attached to a spring of stiffness $k$. Derive the potential energy $V_s$ stored in the spring.

Using #ren-efp, we can derive the potential energy in a spring. From the above diagram, the spring force is: \[\begin{aligned} \vec{F} &= -kx\hat{\imath} \end{aligned}\]

Since the restoring force in a spring only acts in one direction, we can simplify #ren-efp to: \[\begin{aligned} F &= -\frac{dV}{dx} \\ -kx &= -\frac{dV}{dx} \\ kx \, dx &= dV \\ \end{aligned}\]

Integrating both sides: \[\begin{aligned} \int_{V_1}^{V_2} dU &= \int_{x_1}^{x_2} kx \, dx \\ \Delta V & = \frac{1}{2}k\Delta x^2 = \frac{1}{2}k(x^2_2 - x^2_1) \end{aligned}\]

We can simplify the result by assuming $U_1$ is zero, and $x_1$ is at the reference zero-potential displacement. Therefore, the above becomes: \[\begin{aligned} V & = \frac{1}{2}kx^2 \end{aligned}\]

The above result also proves that the restoring force in a spring is a conservative force, as it only depends on the initial and final position of the object.

Warning: Work and energy are not vectors.

Work and energy are scalar quantities, as seen by the dot product that results in a scalar value. Positive and negative work/energy values represent either a gain or loss in energy in the system. For instance, in thermodynamics, work done by the system is negative, and work done on the system is positive. In mechanics, positive work represents an object speeding up (putting energy into the system), and negative work represents an object slowing down (taking energy away from the system).

Gravitational potential energy with uniform gravitational acceleration $g$.

\[V_g = mgh_C\]

The distance $h_C$ is the height of the center of mass of the body above the reference height.

Potential energy stored in a spring with stiffness $k$.

\[V_s = \frac{1}{2}kx_C^2\]

The distance $x_C$ is the displacement of the center of mass of the body relative to the reference displacement.

Work-potential energy theorem.

\[W = -\Delta V\]

Beginning with the work definition: \[\begin{aligned} W = \vec{F} \cdot \Delta \vec{r} \end{aligned}\]

Assuming one-dimensional motion, and using the force-potential energy relationship: \[\begin{aligned} W & = -\frac{\Delta V}{\Delta x}\Delta x\\ & = -\Delta V \end{aligned}\]

Kinetic Energy

Kinetic energy of a point mass.

\[T = \frac{1}{2}mv^2\]

We begin with Newton's 2nd Law assuming $\vec{F}$ is the total external force applied: \[\begin{aligned} \vec{F} &= m\vec{a} \\ & = m\frac{d\vec{v}}{dt} \end{aligned}\]

Using #ren-efp to relate force and energy: \[\begin{aligned} \frac{dU(\vec{r})}{d\vec{r}} & = m\frac{d\vec{v}}{dt} \end{aligned}\]

Note that the negative is not found here, as we're interested in the force exerted on the field, not by the field (see derivation in #ren-efp).

Rearranging then yields: \[\begin{aligned} dU & = m\frac{d\vec{v}}{dt} d\vec{r} \end{aligned}\]

We can then do the following: \[\begin{aligned} dU & = m\frac{d\vec{r}}{dt} d\vec{v} \\ & = m\vec{v} \, d\vec{v} \\ \end{aligned}\]

Integrating both sides: \[\begin{aligned} \int_{U_1}^{U_2} dU & = \int_{\vec{v_1}}^{\vec{v_2}}m\vec{v} \, d\vec{v} \\ \Delta U & = \frac{1}{2}m\left[\vec{v} \cdot \vec{v}\right]_{\vec{v_1}}^{\vec{v_2}} \\ \end{aligned}\]

Assuming an initial velocity of zero, and zero potential energy (i.e the object is at the reference height): \[\begin{aligned} U & = \frac{1}{2}m\left(\vec{v} \cdot \vec{v}\right) \end{aligned}\]

Using conservation of energy, and the fact that a vector dotted with itself equals its magnitude squared (see #rvi-eg), then: \[\begin{aligned} T & = \frac{1}{2}mv^2\\ \end{aligned}\]

Another perhaps less intuitive way is to begin with the work definition #ren-ef: \[\begin{aligned} W &= \int_{\vec{r}_1}^{\vec{r}_2} \vec{F} \cdot d\vec{r} \\ & = \int_{t_1}^{t_2} \vec{F} \cdot \dot{\vec{r}}dt \\ & = \int_{t_1}^{t_2} \vec{F} \cdot \vec{v} \, dt \end{aligned}\]

Using Newton's 2nd Law assuming $\vec{F}$ is the total external force applied: \[\begin{aligned} W & = \int_{t_1}^{t_2} m\vec{a} \cdot \vec{v} \, dt \end{aligned}\]

Note that this integrand looks similar to the following: \[\begin{aligned} \frac{d}{dt}(\vec{u} \cdot \vec{u}) = \dot{\vec{u}} \cdot \vec{u} + \vec{u} \cdot \dot{\vec{u}} = 2\vec{u} \cdot \dot{\vec{u}} \end{aligned}\]

By the fundamental theorem of calculus, this must mean that: \[\begin{aligned} W = \frac{1}{2}m\left(\vec{v} \cdot \vec{v}\right) \end{aligned}\]

And we've arrived at the last step of the previous derivation.

Kinetic energy of an arbitrary body.

\[T = \iiint_{\mathcal{B}} \frac{1}{2} \rho v^2 \, dV\]

Kinetic energy of a rigid body about the center of mass $C$.

\[T = \frac{1}{2} m v_C^2 + \frac{1}{2} I_{C,\hat\omega} \omega^2\]

Using #ren-ea below, in the case that the point $Q$ is the center of mass $C$ we see that $\vec{r}_{QC} = \vec{r}_{CC} = 0$ and so the middle term in #ren-ea is eliminated, leaving the result.

Kinetic energy of a rigid body about the instantaneous center $M$.

\[T = \frac{1}{2} I_{M,\hat\omega} \omega^2\]

Using #ren-ea below, in the case that the point $Q$ is the instantaneous center $M$ we have that $\vec{v}_Q = \vec{v}_M = 0$ and so the first two terms are eliminated.

Kinetic energy of a rigid body about an arbitrary body point $Q$.

\[ T = \frac{1}{2} m v_Q^2 + m \vec{v}_Q \cdot \left( \vec\omega \times \vec{r}_{QC} \right) + \frac{1}{2} I_{Q,\hat\omega} \omega^2 \]

We start with the general expression #rem-eb: \[\begin{aligned} T &= \iiint_{\mathcal{B}} \frac{1}{2} \rho v_P^2 \, dV, \end{aligned}\] where we integrate over the body with a location $P$. We choose a point $Q$ fixed to the body and use #rkg-er to express the velocity of $P$ in terms of $\vec{v}_Q$ and $\vec\omega$, giving \[\begin{aligned} T &= \iiint_{\mathcal{B}} \frac{1}{2} \rho \|\vec{v}_Q + \vec\omega \times \vec{r}_{QP}\|^2 \, dV \\ &= \iiint_{\mathcal{B}} \frac{1}{2} \rho \left(\vec{v}_Q + \vec\omega \times \vec{r}_{QP}\right) \cdot \left(\vec{v}_Q + \vec\omega \times \vec{r}_{QP}\right) \, dV \\ &= \iiint_{\mathcal{B}} \frac{1}{2} \rho \left( \vec{v}_Q \cdot \vec{v}_Q + 2 \vec{v}_Q \cdot (\vec\omega \times \vec{r}_{QP}) + (\vec\omega \times \vec{r}_{QP}) \cdot (\vec\omega \times \vec{r}_{QP}) \right) \, dV. \end{aligned}\] We next use the fact that $\vec{v}_Q$ and $\vec\omega$ do not depend on the integration point within the body to pull them outside of the integrals, which results in \[\begin{aligned} T &= \frac{1}{2} \|\vec{v}_Q\|^2 \iiint_{\mathcal{B}} \rho \, dV + \vec{v}_Q \cdot \bigg( \vec\omega \times \iiint_{\mathcal{B}} \rho \vec{r}_{QP} \, dV \bigg) + \iiint_{\mathcal{B}} \frac{1}{2} \rho (\vec\omega \times \vec{r}_{QP}) \cdot (\vec\omega \times \vec{r}_{QP}) \, dV \\ &= \frac{1}{2} m v_Q^2 + \vec{v}_Q \cdot \left( \vec\omega \times m \vec{r}_{QC} \right) + \iiint_{\mathcal{B}} \frac{1}{2} \rho \| \vec\omega \times \vec{r}_{QP} \|^2 \, dV \\ &= \frac{1}{2} m v_Q^2 + m \vec{v}_Q \cdot \left( \vec\omega \times \vec{r}_{QC} \right) + \iiint_{\mathcal{B}} \frac{1}{2} \rho \omega^2 r_{QP}^2 \sin^2\theta \, dV \\ &= \frac{1}{2} m v_Q^2 + m \vec{v}_Q \cdot \left( \vec\omega \times \vec{r}_{QC} \right) + \frac{1}{2} \omega^2 \iiint_{\mathcal{B}} \rho (r_{QP} \sin\theta)^2 \, dV, \end{aligned}\] where $\theta$ is the angle between $\vec\omega$ and $\vec{r}_{QP}$. But $r_{QP}\sin\theta$ is simply the orthogonal distance to point $P$ from the line through $Q$ in direction $\vec\omega$, so from #rem-ei we see that the final integral above is the moment of inertia $I_{Q,\hat\omega}$ about the $\hat\omega$ axis through point $Q$.

Work-kinetic energy theorem.

\[W = \Delta T\]

Skipping the derivation shown in #ren-ep, taking the last step with the integral: \[\begin{aligned} \int dW & = \int_{\vec{v_1}}^{\vec{v_2}}m\vec{v} \, d\vec{v} \\ W & = \frac{1}{2}m\left[\vec{v} \cdot \vec{v}\right]_{\vec{v_1}}^{\vec{v_2}} \\ & = \Delta T \end{aligned}\]

Work-energy theorem.

\[W = \Delta T + \Delta V\]

where $W$ is the work done by non-conservative forces. If non-conservative forces are not present, then it becomes conservation of energy.

Work and Power

Power transferred by a force $\vec{F}$.

\[P = \vec{F} \cdot \vec{v}\]

The velocity $\vec{v}$ is the velocity of the point where $\vec{F}$ is applied.

Work in terms of power.

\[W = \int_{t_1}^{t_2} P \, dt = \int_{t_1}^{t_2} \vec{F} \cdot \vec{v} \, dt\]

The velocity $\vec{v}$ is the velocity of the point where $\vec{F}$ is applied.

We start from the expression #ren-ef for work as the integral of incremental work $\vec{F} \cdot d\vec{r}$. We change the variable of integration from $\vec{r}$ to $t$, and we substitute $d\vec{r} = \vec{v} dt$ which is simply a rearrangement of the definition #rkv-ev of velocity as $\vec{v} = \frac{d\vec{r}}{dt}$.