Piston motion equations
|This article does not cite any references or sources. (December 2009)|
The motion of a non-offset piston connected to a crank through a connecting rod (as would be found in internal combustion engines), can be expressed through several mathematical equations. This article shows how these motion equations are derived, and shows an example graph.
- 1 Crankshaft geometry
- 2 Equations with respect to angular position (Angle Domain)
- 3 Equations with respect to time (Time Domain)
- 4 Velocity maxima/minima
- 5 Example graph of piston motion
- 6 See also
- 7 References
- 8 Further reading
- 9 External links
l = rod length (distance between piston pin and crank pin)
r = crank radius (distance between crank pin and crank center, i.e. half stroke)
A = crank angle (from cylinder bore centerline at TDC)
x = piston pin position (upward from crank center along cylinder bore centerline)
v = piston pin velocity (upward from crank center along cylinder bore centerline)
a = piston pin acceleration (upward from crank center along cylinder bore centerline)
ω = crank angular velocity in rad/s
Equations with respect to angular position (Angle Domain)
The equations that follow describe the reciprocating motion of the piston with respect to crank angle.
Example graphs of these equations are shown below.
Position with respect to crank angle (by rearranging the triangle relation):
Equations with respect to time (Time Domain)
Angular velocity derivatives
If angular velocity is constant, then
and the following relations apply:
Converting from Angle Domain to Time Domain
The equations that follow describe the reciprocating motion of the piston with respect to time.
Position with respect to time is simply:
Scaling for angular velocity
You can see that x is unscaled, x' is scaled by ω, and x" is scaled by ω².
To convert x' from velocity vs angle [inch/rad] to velocity vs time [inch/s] multiply x' by ω [rad/s].
To convert x" from acceleration vs angle [inch/rad²] to acceleration vs time [inch/s²] multiply x" by ω² [rad²/s²].
Note that dimensional analysis shows that the units are consistent.
Equation for diameter of a piston
R2= (F2)(R1)^2 / F1
Acceleration zero crossings
The velocity maxima and minima do not occur at crank angles (A) of plus or minus 90°.
The velocity maxima and minima occur at crank angles that depend on rod length (l) and half stroke (r),
and correspond to the crank angles where the acceleration is zero (crossing the horizontal axis).
Crank-Rod angle not right angled
The velocity maxima and minima do not necessarily occur when the crank makes a right angle with the rod.
Counter-examples exist to disprove the idea that velocity maxima/minima occur when crank-rod angle is right angled.
For rod length 6" and crank radius 2", numerically solving the acceleration zero-crossings finds the velocity maxima/minima to be at crank angles of ±73.17615°.
Then, using the triangle sine law, it is found that the crank-rod angle is 88.21738° and the rod-vertical angle is 18.60647°.
Clearly, in this example, the angle between the crank and the rod is not a right angle.
(Sanity check, summing the angles of the triangle 88.21738° + 18.60647° + 73.17615° gives 180.00000°)
A single counter-example is sufficient to disprove the statement "velocity maxima/minima occur when crank makes a right angle with rod".
Example graph of piston motion
The graph shows x, x', x" with respect to crank angle for various half strokes, where L = rod length (l) and R = half stroke (r):
Pistons motion animation with same rod length and crank radius values in graph above :
- John Benjamin Heywood, Internal Combustion Engine Fundamentals, McGraw Hill, 1989.
- Charles Fayette Taylor, The Internal Combustion Engine in Theory and Practice, Vol. 1 & 2, 2nd Edition, MIT Press 1985.