透過 MCP 連接 →

輸入計算

數學公式

廣告

結果

阻尼類型
under-damped
x(0) = 由 x0 處靜止釋放
阻尼角頻率 w_d 4.898979 rad/s
自然週期 T0 1.256637 s
繪製時間範圍(4 個週期) 5.026548 s
時間步長 dt 0.100531 s
第一步的位移 x 0.884163
時間 t(秒) 位移 x(t)
0 1
0.100531 0.884163
0.201062 0.591278
0.301593 0.219229
0.402124 -0.134186
0.502655 -0.393158
0.603186 -0.516853
0.703717 -0.502188
0.804248 -0.377651
0.904779 -0.191163
1.00531 0.004199
1.105841 0.163038
1.206372 0.256408
1.306903 0.275309
1.407434 0.229069
1.507964 0.139915
1.608495 0.035644
1.709026 -0.057623
1.809557 -0.120989
1.910088 -0.1457
2.010619 -0.133257
2.11115 -0.093217
2.211681 -0.039579
2.312212 0.013159
2.412743 0.053343
2.513274 0.074253
2.613805 0.074714
2.714336 0.058347
2.814867 0.031886
2.915398 0.003098
3.015929 -0.021141
3.11646 -0.036227
3.216991 -0.040447
3.317522 -0.034833
3.418053 -0.022436
3.518584 -0.007269
3.619115 0.006765
3.719646 0.016726
3.820177 0.021125
3.920708 0.019988
4.021239 0.014579
4.12177 0.006888
4.222301 -0.000946
4.322831 -0.007141
4.423362 -0.010608
4.523893 -0.011065
4.624424 -0.00896
4.724955 -0.00523
4.825486 -0.001007
4.926017 0.002672
5.026548 0.005082

這個計算器能做什麼

本工具用於計算一維阻尼諧振子的位移 \(x(t)\),假設系統從初始位移 \(\text{x}_0\) 處靜止釋放。它求解標準的(以質量歸一化的)運動方程式,並將位移在四個自然週期內逐點列表,讓你清楚看到系統如何逐漸回到平衡位置。同時,工具還會自動判別系統屬於欠阻尼、臨界阻尼或過阻尼。

展示欠阻尼、臨界阻尼與過阻尼衰減的三條阻尼振盪曲線
三種阻尼狀態:欠阻尼會振盪並衰減,臨界阻尼無超調地最快回歸,過阻尼則緩慢回歸。

主導方程式

運動遵循線性常微分方程式 $$\frac{d^2x}{dt^2} + 2k\frac{dx}{dt} + \omega_0^2 x = 0$$ 其中 \(\omega_0\) 為無阻尼角頻率,\(k\) 為阻力(阻尼)係數(兩者單位皆為 \(1/\text{s}\))。在初始條件 \(x(0) = \text{x}_0\) 與 \(\frac{dx}{dt}(0) = 0\) 之下,封閉解的形式取決於 \(k\) 與 \(\omega_0\) 的大小關係。

當 \(k < \omega_0\) 時,系統為欠阻尼,會以較低的阻尼角頻率 \(\omega_d = \sqrt{\omega_0^2 - k^2}\) 振盪,且振幅依 \(e^{-kt}\) 衰減。當 \(k = \omega_0\) 時,系統為臨界阻尼,能在不振盪的前提下以最快速度回到靜止:$$x(t) = \text{x}_0\left(1 + \omega_0\, t\right) e^{-\omega_0\, t}$$ 當 \(k > \omega_0\) 時,系統為過阻尼,會緩慢爬回平衡而不發生振盪。

質量-彈簧-阻尼器系統示意圖
物理模型:帶阻尼器的彈簧上的質量塊,遵循運動方程式。

使用方式

輸入無阻尼角頻率 \(\omega_0\)(必須大於零)、阻尼係數 \(k\)(大於或等於零,\(k = 0\) 表示純粹無阻尼運動)、初始位移 \(\text{x}_0\),以及數據表的時間分割數。自然週期為 \(T_0 = \frac{2\pi}{\omega_0}\);表格範圍涵蓋 \(4 T_0\),並以相等的步長 \(dt = \frac{\text{timeSpan}}{\text{divisions}}\) 切分,因此會產生 \(\text{divisions}+1\) 列資料。

實際範例

當 \(\omega_0 = 5\)、\(k = 1\)、\(\text{x}_0 = 1\)、分割數為 \(50\) 時,系統屬於欠阻尼,阻尼角頻率 \(\omega_d = \sqrt{25 - 1} = 4.89898\ \text{rad/s}\)。自然週期為 \(1.256637\) 秒,總時間範圍為 \(5.026548\) 秒,\(dt = 0.100531\) 秒。在 \(t = 0\) 時位移為 \(1\);在第一步 \(t = 0.100531\) 秒時,位移約為 \(0.884153\)。

常見問題

阻尼係數 \(k\) 代表什麼?它是以質量歸一化的「半阻尼」項;每單位質量所受的阻力等於 \(2k\) 乘以速度。

如果 \(k\) 剛好等於 \(\omega_0\) 會怎樣?在這一點上,欠阻尼與過阻尼的公式會出現可去奇異點,因此只要 \(k\) 落在極小的容差範圍內接近 \(\omega_0\),工具就會改用臨界阻尼公式。

為什麼剛好取四個週期?四個自然週期已足以呈現完整的衰減包絡,同時又能讓表格保持精簡、易於閱讀。

最後更新: