CodeCogs

An analysis of the velocity and acceleration of a piston

You're viewing an older version of this page (#3610). View the current version.

In order to define the velocity and acceleration of a piston, consider the mechanism in Figure 1, where the crank $OC$ is driven with the uniform angular velocity $\omega$ . Also, let $r$ be the crank radius, $l$ the rod length, $x$ the position of the piston pin from the crank center, $\phi$ the $\angle {CPO}$ angle, and $\theta$ the $\angle {COP}$ angle (the crank angle). Such a mechanism in motion can be seen in Figure 2, where the crankshaft is depicted with red, and the pistons with gray (to see the animation click on the thumbnail).

From Figure 1 it can be seen that:

$x = r \cos \theta + l \cos \phi$

(3)

It can also be noted from Figure 1 that:

$r \sin \theta = l \sin \phi$

(4)

Squaring equation (#2) gives:

$r^2 \sin^2 \theta = l^2 \sin^2 \phi$

(5)

which can also be written as:

$r^2 \sin^2 \theta = l^2 (1- \cos^2 \phi)$

(6)

This eventually leads to:

$\cos \phi = \sqrt{1-\frac{\sin^2 \theta}{n^2}}$

(7)

where $\displaystyle n = \frac{l}{r}$ .

By using this expression of $\cos \phi$ in equation (#1), we obtain:

$x = r \cos \theta + l \sqrt{1-\frac{\sin^2 \theta}{n^2}}$

(8)

which can also be written as:

$x = r \left( \cos \theta + n \sqrt{1-\frac{\sin^2 \theta}{n^2}} \right)$

(9)

and, furthermore, as:

$x = r \left( \cos \theta + \sqrt{n^2 - \sin^2 \theta \right)$

(10)

In order to calculate the piston velocity, we differentiate (#8) with respect to time, when we get:

$v = \frac{dx}{dt} ( = \dot x ) = -r \omega \left( \sin \theta + \frac{\sin 2\theta}{2\sqrt{n^2 - \sin^2 \theta}} \right)$

(11)

The piston acceleration can then be calculated by differentiating again with respect to time, when we obtain:

$a = \frac{dv}{dt} ( = \ddot x ) = -r \omega^2 \left[ \cos \theta + \frac{n^2 \cos 2\theta + \sin^4 \theta}{(n^2 - \sin^2 \theta)^\frac{3}{2}} \right]$

(12)

Note that in these equations the positive direction of velocity and acceleration is away from the crankshaft.

In normal situations $\sin^2 \theta$ and $\sin^4 \theta$ can be neglected in comparison with $n^2$ , and, thus, the above equations can be reduced to:

$v = -r \omega \left( \sin \theta + \frac{\sin 2\theta}{2n} \right)$

(13)

$a = -r \omega^2 \left( \cos \theta + \frac{\cos 2\theta}{n} \right)$

(14)

Example 1 [imperial]

Problem

In the crank and slotted lever mechanism diagramed in Figure E1, the crank $OP$ is driven at a uniform speed of $\omega$ radians per second. If $OL$ is the perpendicular from $O$ on $XQ$ , $XQ$ being the center-line of the slotted lever, prove that the angular acceleration of the slotted lever is given by:

$\ddot \phi = \omega^2 \frac{OL}{XP} \left( \frac{2PL}{XP} - 1 \right)$

Also, find the magnitude and direction of the acceleration of the point $Q$ when the crank angle is $\theta = 60^\circ$ , given that $OP = 3\; inch$ , $OX = 10\; inch$ , and $XQ = 18\; inch$ .

Workings

From the diagram we can write that:

$OL = OX \sin \phi$

(15)

$OL = OP \sin (\theta - \phi)$

(16)

from which:

$OX \sin \phi = OP \sin (\theta - \phi)$

(17)

By differentiating (#3) with respect to time, we get:

$OX \dot \phi \cos \phi = OP (\omega - \dot \phi) \cos (\theta - \phi)$

(18)

which leads to:

$\dot \phi [ OX \cos \phi + OP \cos (\theta - \phi) ] = OP \omega \cos (\theta - \phi)$

(19)

From Figure E1 we can write that:

$XL = OX \cos \phi$

(20)

$PL = OP \cos (\theta - \phi)$

(21)

Thus, equation (#5) becomes:

$XP \dot \phi = PL \omega$

(22)

or:

$\dot \phi = \frac{PL}{XP} \omega$

(23)

By differentiating (#5) with respect to time, we obtain:

$\ddot \phi [ OX \cos \phi + OP \cos(\theta - \phi)] - \dot \phi^2 OX \sin \phi - \dot \phi (\omega - \phi) OP \sin (\theta - \phi) = -\omega (\omega - \dot \phi) OP \sin (\theta - \phi)$

(24)

By substituting $\dot \phi$ in equation (#10) with the expression from (#9), we can write:

$XP \ddot \phi - OL \frac{PL^2}{XP^2} \omega^2 - OL \frac{PL}{XP} \left( 1 - \frac{PL}{XP} \right) \omega^2 = -OL \left( 1 - \frac{PL}{XP} \right) \omega^2$

(25)

or:

$\ddot \phi = \omega^2 \frac{OL}{XP} \left( \frac{2PL}{XP} - 1 \right) \; \blacksquare$

(26)

By replacing the numerical values in equation (#3), we can write that:

$10 \sin \phi = 3 \sin (60 - \phi)$

(27)

As:

$3 \sin (60 - \phi) = 3(\sin 60 \cos \phi - \sin \phi \cos 60)$

(28)

equation (#13) becomes:

$10 \sin \phi = 3\frac{\sqrt{3}}{2} \cos \phi - \frac{3}{2} \sin \phi$

(29)

which gives:

$\tan \phi = \frac{2.59}{11.5}$

(30)

Thus, we get that:

$\phi = 12^\circ 32'$

(31)

and, by using this value, we also obtain that:

$OL = 10 \sin \phi = 2.20 \; inch$

(32)

$PL = 3 \cos (60-\phi) = 2.035 \; inch$

(33)

$XP = 10 \cos \phi + PL = 11.79 \; inch$

(34)

Furthermore, by replacing numerical values in (#9) and (#12), we respectively get that:

$\dot \phi = \frac{2.035}{11.79} \cdot 10 = 1.75 \; \frac{rad}{s}$

(35)

$\ddot \phi = 100 \cdot \frac{2.2}{11.79} \cdot \left( \frac{2 \cdot 2.035}{11.79} - 1 \right) = -12.2 \; \frac{rad}{s^2}$

(36)

In order to calculate the acceleration of $Q$ , $a$ , we are going to split it in a tangential component, $a_T$ , perpendicular to $XQ$ , and a radial component, $a_r$ , in the direction of $XQ$ (see Figure E1). Also, let $\alpha$ be the angle which defines the direction of the acceleration. We can write that: