活塞運動方程

本篇以數學公式說明，從連桿聯結到曲柄的非偏移活塞運動（在內燃機中）；並且就這些運動方程式如何導出，附一示例圖。

曲軸幾何

活塞销、曲柄销和曲柄中心几何布局示意图

定義

$l$ 活塞桿長度（活塞銷與曲柄銷之間的距離）

$r$ 曲柄半徑（曲柄銷與曲柄中心之間的距離，即半行程）

$A$ 曲軸轉角（從缸膛中心線開始）

$x$ 活塞銷位置（從曲軸中心沿缸膛中心線向上）

$v$ 活塞銷速度（從曲軸中心沿缸膛中心線向上）

$a$ 活塞銷加速度（從曲軸中心沿缸膛中心線向上）

$\omega$ 曲柄角速度

角速度

所述曲軸的角速度，即是發動機每分鐘轉數（RPM）：

$\omega =\frac{2\pi \cdot \mathrm{R}\mathrm{P}\mathrm{M}}{60}$

三角關係

如圖所示，曲柄銷，曲柄中心和活塞銷形成角NOP。根據餘弦定理可以看出：

$l^2=r^2+x^2-2\cdot r\cdot x\cdot \cos A$

關於角度位置的方程（角度域）

以下公式描述活塞相對於曲軸轉角的往復運動。下面顯示了這些方程的示例圖。

位置

相對於曲軸轉角的位置（通過重新排列三角關係）：

$l^2-r^2=x^2-2\cdot r\cdot x\cdot \cos A$

$l^2-r^2=x^2-2\cdot r\cdot x\cdot \cos A+r^2[({\cos }^{2}A+{\sin }^{2}A)-1]$

$l^2-r^2+r^2-r^2{\sin }^{2}A=x^2-2\cdot r\cdot x\cdot \cos A+r^2{\cos }^{2}A$

$l^2-r^2{\sin }^{2}A=(x-r\cdot \cos A{)}^2$

$x-r\cdot \cos A=\sqrt{l^2-r^2{\sin }^{2}A}$

$x=r\cdot \cos A+\sqrt{l^2-r^2{\sin }^{2}A}$

速度

相對於曲軸轉角的速度（採用一階導數，使用链式法则）：

$\begin{array}{ccc}{x}^{\prime }& =& \frac{dx}{dA}\\ & =& -r\sin A+\frac{(\frac{1}{2}).(-2).r^2\sin A\cos A}{\sqrt{l^2-r^2{\sin }^{2}A}}\\ & =& -r\sin A-\frac{r^2\sin A\cos A}{\sqrt{l^2-r^2{\sin }^{2}A}}\end{array}$

加速

關於曲軸轉角的加速度（採用二階導數，使用链式法则和商法則）：

$\begin{array}{ccc}{x}^{″}& =& \frac{d^{2}x}{d{A}^2}\\ & =& -r\cos A-\frac{r^2{\cos }^{2}A}{\sqrt{l^2-r^2{\sin }^{2}A}}-\frac{-r^2{\sin }^{2}A}{\sqrt{l^2-r^2{\sin }^{2}A}}-\frac{r^2\sin A\cos A.(-\frac{1}{2})\cdot (-2).r^2\sin A\cos A}{{\left(\sqrt{l^2-r^2{\sin }^{2}A}\right)}^3}\\ & =& -r\cos A-\frac{r^2({\cos }^{2}A-{\sin }^{2}A)}{\sqrt{l^2-r^2{\sin }^{2}A}}-\frac{r^4{\sin }^{2}A{\cos }^{2}A}{{\left(\sqrt{l^2-r^2{\sin }^{2}A}\right)}^3}\end{array}$

時間方程（時域）

角速度導數

如果角速度是恆定的，那麼

$A=\omega t$

並適用以下關係：

$\frac{dA}{dt}=\omega$

$\frac{d^{2}A}{d{t}^2}=0$

從角度域轉換到時域

下面的等式描述了活塞相對於時間的往復運動。如果時間域是必需的，而不是角度域，先用 ω 取代在等式噸，然後擴展為角速度如下：

位置

單純相對於時間的位置即是：

$x$

速度

相對於時間的速度（使用链式法则）：

$\begin{array}{ccc}v& =& \frac{dx}{dt}\\ & =& \frac{dx}{dA}\cdot \frac{dA}{dt}\\ & =& \frac{dx}{dA}\cdot \omega \\ & =& x^{\prime }\cdot \omega \end{array}$

加速

相對於時間的加速度（使用链式法则和乘積法則以及角速度導數）：

$\begin{array}{ccc}a& =& \frac{d^{2}x}{d{t}^2}\\ & =& \frac{d}{dt}\frac{dx}{dt}\\ & =& \frac{d}{dt}(\frac{dx}{dA}\cdot \frac{dA}{dt})\\ & =& \frac{d}{dt}(\frac{dx}{dA})\cdot \frac{dA}{dt}+\frac{dx}{dA}\cdot \frac{d}{dt}(\frac{dA}{dt})\\ & =& \frac{d}{dA}(\frac{dx}{dA})\cdot (\frac{dA}{dt}{)}^2+\frac{dx}{dA}\cdot \frac{d^{2}A}{d{t}^2}\\ & =& \frac{d^{2}x}{d{A}^2}\cdot (\frac{dA}{dt}{)}^2+\frac{dx}{dA}\cdot \frac{d^{2}A}{d{t}^2}\\ & =& \frac{d^{2}x}{d{A}^2}\cdot {\omega }^2+\frac{dx}{dA}\cdot 0\\ & =& x^{″}\cdot {\omega }^2\end{array}$

速度最大值/最小值

加速過零點

速度的最大值和最小值都沒有在曲柄角度發生（A）加或減90°。速度最大值和最小值出現在取決於桿長（l）和半行程（r）的曲柄角上，並且對應於加速度為零（橫過水平軸）的曲柄角。

曲柄桿角度不正確

當曲柄與桿成直角時，速度最大值和最小值不一定會發生。有反例證明當曲柄角度正確時速度最大值/最小值出現的觀點。

示例

對於 6 英吋的連桿長和 2 英吋的曲柄半徑，通過數值求解加速度過零點，可確定速度最大值/最小值在曲柄角±73.17615°處。然後使用三角正弦定律，發現曲柄角為88.21738°，桿垂直角為18.60647°。在這個例子中，顯然地曲柄和連桿之間的角度並不是直角。

總結三角形的角度 88.21738°+ 18.60647°+ 73.17615°給出 180.000°。這個反例就反駁了“速度最大值/最小值是發生在當曲柄與桿成直角時”的說法。

活塞運動示例圖

曲線圖顯示了相對於曲軸轉角的 x，x'，x" 的各種半衝程，其中 L =長度（l）和 R =半行程（r）：

The vertical axis units are inches for position, [inches/rad] for velocity, [inches/rad²] for acceleration.
The horizontal axis units are crank angle degrees.

Template:Clr 杆长和曲柄半径与上图相同的活塞动画：

活塞各半冲程动画

參考文獻

• -{R|http://www.epi-eng.com/piston_engine_technology/piston_motion_basics.htm}-

延伸閱讀

• 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.

外部链接

• epi-eng Piston Motion
• codecogs Velocity and Acceleration of a Piston
• animated engines Four Stroke Engine
• desmos interactive crank animation
• networcs D & T Mechanisms - Interactive Tools for Teachers
• mecamedia piston motion animation
• youtube Rotating chevy 350 short block.
• youtube 3D animation of a V8 ENGINE
• youtube Inside a V8 Engine at Idle Speed