この記事は、高校物理(物理・物理基礎)で単振動や振り子を学んでいる人、大学初年次の物理学実験で「単振り子による重力加速度 gg の測定」や「ばね定数と単振動」のレポート課題に取り組んでいる人、実験データの誤差要因(有限振幅・空気抵抗・ばねの直列並列接続)を正しく考察・計算したい人向けです。

「単振り子の周期公式 T=2πL/gT = 2\pi\sqrt{L/g} で、なぜおもりの重さは関係ないのか?」「実験レポートで振れ角を大きくしすぎたとき、ボルダの補正公式をどう使えばいいのか?」「ばねを直列につないだときと並列につないだときで、周期はどう変わる?」「鉛直ばね振り子で重力があるのに、なぜ水平ばねと同じ周期公式が使えるのか?」「速度が変位より90°進むとはどういう意味か?」。

単振動と振り子の力学でつまずく最大の理由は、「微小角近似が成り立つ前提条件」と「復元力の比例関係 F=−kxF = -kx」を数式と図解で整理しないまま、T=2π…T = 2\pi\sqrt{\dots} の公式だけを丸暗記しようとすることにあります。

この記事では、単振り子とばね振り子の周期公式の導出から、大振幅時のボルダ補正早見表、実験データから重力加速度 gg を逆算する実践手計算、ばねの直列・並列合成、単振動の変位・速度・加速度の位相差(90°/180°)、力学的エネルギー保存則まで、ステップ順にわかりやすく解説します。


まず確認:実験と計算ツールで同じ条件を使う

確認項目 本文・ツールでの扱い レポートに残す内容
1周期 出発位置へ同じ進行方向で戻る1往復 往復回数N、合計時間t、T=t/N
長さL 支点からおもりの重心まで 糸のみか、有効長かを区別
角度 画面は度、展開式にはradで代入 初期角と採用した補正次数
順算モード 1+θ₀²/16+11θ₀⁴/3072 打ち切り近似で厳密解ではない
重力逆算モード 1+θ₀²/16でT₀を復元 順算とは補正次数が異なる
必要長モード 微小角周期T₀から長さを計算 大振幅の周期指定ではない

単振動・振り子計算ツールの重力逆算で、L=1 m、測定周期2.008 s、角度10°を入力すると、T₀≈2.0042 s、g≈9.8284 m/s²です。これは本文の架空の測定データによる検算例で、特定の場所の重力測定結果ではありません。

20往復40.16秒を20で割って周期2.008秒、角度10度で補正して微小角周期約2.0042秒
測定周期と近似補正後の周期を分けて記録します。図を選択して拡大できます。

このモデルは点質量・伸びない軽い糸・摩擦なしの理想化です。大角度では打ち切り近似の差も増えます。画面のアニメーションはモデル理解用で、実機の損失や安全性の評価ではありません。

結論:単振動・振り子の重要公式早見表

単振動(Simple Harmonic Motion: SHM)の本質は、物体の変位 xx に比例する中心向きの復元力 F=−KxF = -K x が働く運動であるという点にあります。このとき運動方程式 md2xdt2=−Kxm \frac{d^2x}{dt^2} = -K x より、角振動数は ω=K/m\omega = \sqrt{K/m}、周期は T=2πω=2πmKT = \frac{2\pi}{\omega} = 2\pi\sqrt{\frac{m}{K}} で決定されます。

単振り子の公式(微小角近似):
周期: T₀ = 2π √(L/g) (振れ角 θ₀ ≪ 1 rad のとき、おもりの質量 m に無関係)
有限振幅(ボルダ補正): T ≈ T₀ [1 + (1/16)θ₀²] (θ₀ はラジアン単位)
重力加速度の逆算: g = 4π² L / T₀²
必要長さの逆算: L = g T₀² / (4π²) (秒振子 T = 2.0 s のとき L ≈ 99.4 cm)

ばね振り子の公式:
周期: T = 2π √(m/k) (水平・鉛直共通)
並列合成(硬くなる): k_eff = k₁ + k₂ ⇒ T_parallel = 2π √[m / (k₁ + k₂)]
直列合成(柔らかくなる): 1/k_eff = 1/k₁ + 1/k₂ ⇒ k_eff = (k₁k₂) / (k₁ + k₂)
鉛直つりあいの伸び: x₀ = mg / k (この位置が振動中心)

単振動の一般解・力学的エネルギー:
変位: x(t) = A cos(ωt + φ)
速度: v(t) = −Aω sin(ωt + φ) (変位より 90° 進む、最大値 v_max = Aω)
加速度: a(t) = −Aω² cos(ωt + φ) = −ω² x (変位と 180° 逆位相、最大値 a_max = Aω²)
力学的エネルギー: E = K + U = ½mv² + ½kx² = ½kA² = ½mv_max² (常に一定)

振動系 周期公式 T 復元力 F 振動数 f / 角振動数 ω 特徴と適用条件
単振り子(小振幅) T₀ = 2π √(L/g) F ≈ −(mg/L)x ω = √(g/L), f = 1/T₀ 質量 m に非依存。振れ角 θ₀ ≤ 5° で高精度(誤差0.05%以下)。
単振り子(大振幅) T ≈ T₀(1 + θ₀²/16) F = −mg sin θ 実効角振動数は低下 振幅が大きくなると周期は延びる。30° で約 +1.7% の遅れ。
単一ばね振り子 T = 2π √(m/k) F = −kx ω = √(k/m), f = 1/T 重力の有無(水平・鉛直・斜面)で周期は不変。
ばね並列接続 T = 2π √[m/(k₁+k₂)] F = −(k₁+k₂)x ω = √[(k₁+k₂)/m] 同じばね2本なら剛性2倍で周期は 1/√2 ≈ 0.707 倍に短縮。
ばね直列接続 T = 2π √[m(k₁+k₂)/(k₁k₂)] F = −k_eff x ω = √[k₁k₂/{m(k₁+k₂)}] 同じばね2本なら周期は √2 ≈ 1.414 倍に延長。

1. 単振り子の運動方程式・周期公式と微小角近似

単振り子とは、伸び縮みせず質量の無視できる長さ LL の糸の先端に、質量 mm のおもりを取り付けて鉛直面内で往復運動させた理想的な系です。

     支点 O
      | \
      |  \  糸の長さ L
      |   \ 
      |θ   \
鉛直線 |     ● おもり (質量 m)
      |    ↙ 接線方向復元力: F = -mg sin θ
      |   ↓ 重力: mg

運動方程式と復元力の導出

鉛直線からの傾き角を θ\theta(ラジアン)、最下点を基準とする円弧上の変位を s=Lθs = L\theta とします。おもりに働く力は「糸の張力 SS」と「鉛直下向きの重力 mgmg」の2つです。

糸に垂直な円弧の接線方向について運動方程式を立てると:

md2sdt2=−mgsin⁡θm \frac{d^2s}{dt^2} = -mg \sin\theta

ここで s=Lθs = L\theta より d2sdt2=Ld2θdt2\frac{d^2s}{dt^2} = L \frac{d^2\theta}{dt^2} であるため、両辺を mLmL で割ると次の厳密な非線形微分方程式が得られます。

d2θdt2+gLsin⁡θ=0\frac{d^2\theta}{dt^2} + \frac{g}{L} \sin\theta = 0

微小角近似(sin⁡θ≈θ\sin\theta \approx \theta)による線形化

この方程式には sin⁡θ\sin\theta という非線形項が含まれており、そのままでは厳密な単振動になりません。しかし、角度 θ\theta が十分小さいとき、テイラー展開(マクローリン展開):

sin⁡θ=θ−θ33!+θ55!−…\sin\theta = \theta - \frac{\theta^3}{3!} + \frac{\theta^5}{5!} - \dots

の第1項のみを採用し、sin⁡θ≈θ\sin\theta \approx \theta と近似(微小角近似)することができます。これにより運動方程式は標準的な単振動の方程式に変形されます。

d2θdt2+gLθ=0\frac{d^2\theta}{dt^2} + \frac{g}{L} \theta = 0

ここで ω2=gL\omega^2 = \frac{g}{L} と置くと、角振動数は ω=gL\omega = \sqrt{\frac{g}{L}} となり、周期 T0T_0 は次のように導出されます。

T0=2πω=2πLgT_0 = \frac{2\pi}{\omega} = 2\pi \sqrt{\frac{L}{g}}

おもりの質量 mm に依存しない理由

周期の式 T0=2πL/gT_0 = 2\pi\sqrt{L/g} には、おもりの質量 mm が一切含まれていません。
その理由は、おもりの動きにくさを決める「慣性質量 minertialm_{\text{inertial}}」と、重力による引く力を決める「重力質量 mgravitationalm_{\text{gravitational}}」が物理的に等価(等価原理)であるため、運動方程式の両辺で質量が完全に相殺されるからです。

支点からおもりの重心までの有効長が同じで、点質量・質量のない糸・摩擦なしというモデルの条件が成り立つなら、質量を変えても周期は同じです。実際の球の大きさや糸の伸びは別に評価します。


2. 大振幅時の誤差とボルダの補正公式

高校物理の教科書では「振幅が小さければ周期は一定(振り子の等時性)」と教わりますが、大学の実験や精密測定では**「どこまでが微小角なのか」「振幅が大きくなると何秒ズレるのか」**が決定的に重要になります。

ボルダ(Borda)の展開補正公式

振幅角 θ0\theta_0 が大きくなると、復元力 −mgsin⁡θ-mg\sin\theta は線形値 −mgθ-mg\theta よりも小さくなる(sin⁡θ<θ\sin\theta < \theta)ため、おもりの戻る勢いが弱まり、周期は微小角近似値 T0T_0 よりも長くなります。

厳密な周期は第1種完全楕円積分 K(k)K(k) を用いて表されますが、実用的にはフランスの天文学者ジャン=シャルル・ド・ボルダ(Jean-Charles de Borda)が導いた次の展開近似式が使われます。

T=T0(1+116θ02+113072θ04+… )T = T_0 \left( 1 + \frac{1}{16}\theta_0^2 + \frac{11}{3072}\theta_0^4 + \dots \right)

ここで θ0\theta_0 は**度(°)ではなくラジアン(rad)**で代入します(θrad=θdeg×π180\theta_{\text{rad}} = \theta_{\text{deg}} \times \frac{\pi}{180})。

振れ角別・誤差早見表

振幅角 θ0\theta_0 ごとに、微小角近似値 T0T_0 に対する周期の延び(誤差率)を計算した表です。

振幅角 θ₀ (度) 角度 (rad) ボルダ補正係数 (1 + θ₀²/16) 周期の増加率 (%) 必要精度との比較
2.0° 0.0349 rad 1.000076 +0.008% 必要な測定精度より十分小さいかを確認。
5.0° 0.0873 rad 1.000476 +0.048% 0.05%程度の差を許容できるかで判断。
10.0° 0.1745 rad 1.001904 +0.19% 有効数字3〜4桁の測定では補正を考慮すべき領域。
15.0° 0.2618 rad 1.004284 +0.43% 0.43%より小さい不確かさを目指すなら補正を検討。
20.0° 0.3491 rad 1.007615 +0.76% 重力加速度 g の算出値に約 1.5% の過小誤差を生む。
30.0° 0.5236 rad 1.017135 +1.71% 約1.7%の周期差。必要精度に合わせて高次補正を検討。

この表は1+θ₀²/16までの打ち切り近似です。30°でθ₀⁴項も加えると係数は約1.017404になります。合格・不合格の共通基準ではなく、指定された必要精度と比較します。

レポート作成の要点:
周期 TT が 1% 長くなると、逆算する重力加速度 g∝1/T2g \propto 1/T^2 は約 2% 小さく見積もられてしまいます。
実験で振幅を 15°〜20° 程度振ってしまった場合は、考察欄で「微小角近似の破れによる系統誤差」としてボルダ補正を適用した修正値を併記すると、採用した近似とその影響を読み手が確認できます。


3. 物理実験での重力加速度 g 測定手順・誤差評価とレポート考察例

理系大学1年次の基礎物理実験で最も頻出の「単振り子による重力加速度 gg の測定」の実践手順をまとめます。

【測定系】
支点 ──── 糸の長さ L (ノギス・巻き尺で測定) ──── おもり (半径 r を足して重心まで L_eff = L + r)

なぜ「20〜50回」の往復時間を測るのか?

手動ストップウォッチでは、20回や50回など複数往復の合計時間を測り、T=tN/NT=t_N/N とすると開始・停止の時間ずれの相対的な影響を小さくできます。回数は装置・減衰・授業の指定に合わせます。光電ゲートや動画なら1周期を測る方法もあり、「必ず20回以上」という規則ではありません。

例えば合計時間の不確かさを0.15 sと仮定すると、2 sの1周期では7.5%、40 sの20周期では0.375%です。これは同じ絶対的な時間ずれを仮定した計算例で、人の反応時間の測定結果ではありません。回数を増やしても、長さのずれ・振幅の変化・数え間違いが自動で小さくなるわけではありません。同じ位置・同じ進行方向を1周期の基準にし、繰り返し測定も行います。

実験データからの手計算と逆算例

  • 測定条件:
    • 糸の長さ+おもり重心までの有効長 L=1.000 mL = 1.000\,\text{m}
    • 振幅角 θ0=10.0∘\theta_0 = 10.0^\circ(=0.1745 rad= 0.1745\,\text{rad})
    • 20回往復の測定時間 t20=40.16 st_{20} = 40.16\,\text{s}
    • 測定周期 Tmeas=40.1620=2.0080 sT_{\text{meas}} = \frac{40.16}{20} = 2.0080\,\text{s}

ステップ1:ボルダ補正による微小角周期 T0T_0 の復元

T0=Tmeas1+116θ02=2.00801+116(0.1745)2=2.00801.00190≈2.0042 sT_0 = \frac{T_{\text{meas}}}{1 + \frac{1}{16}\theta_0^2} = \frac{2.0080}{1 + \frac{1}{16}(0.1745)^2} = \frac{2.0080}{1.00190} \approx 2.0042\,\text{s}

ステップ2:重力加速度 gg の計算

g=4π2LT02=4×(3.14159)2×1.000(2.0042)2=39.47844.0168≈9.828 m/s2g = \frac{4\pi^2 L}{T_0^2} = \frac{4 \times (3.14159)^2 \times 1.000}{(2.0042)^2} = \frac{39.4784}{4.0168} \approx 9.828\,\text{m/s}^2

(もしボルダ補正を無視して Tmeas=2.0080 sT_{\text{meas}} = 2.0080\,\text{s} でそのまま計算すると g=39.4784/(2.0080)2=9.791 m/s2g = 39.4784 / (2.0080)^2 = 9.791\,\text{m/s}^2 となり、補正を入れた値とは約 0.037 m/s² 異なります。標準重力加速度は約束された基準値で、その場所の実測値ではないため、近い方が正解とは判断できません)


4. ばね振り子の周期公式・直列並列合成と鉛直つりあいの位置シフト

ばねの先端に質量 mm のおもりを取り付けたばね振り子では、フックの法則 F=−kxF = -kx に基づく単振動が生じます。

基本周期公式

ばね定数を kk(N/m)、おもりの質量を mm(kg)とすると、運動方程式は md2xdt2=−kxm \frac{d^2x}{dt^2} = -kx です。
角振動数は ω=km\omega = \sqrt{\frac{k}{m}}、周期は次のようになります。

T=2πmkT = 2\pi \sqrt{\frac{m}{k}}

ばね振り子の周期は、理想的な線形ばねで、ばね質量・摩擦を無視できる範囲では、振幅 AA にも一定の重力加速度 gg にも依存しません。

なぜ鉛直ばね振り子でも周期は水平と同じなのか?

「水平ばね振り子はおもりを横に引っ張るだけだが、鉛直に吊るすと常に下向きの重力 mgmg が加わる。それなのに、なぜ周期公式 T=2πm/kT = 2\pi\sqrt{m/k} は変わらないのか?」という疑問は、物理学習者が最も頻繁に抱く論点です。

【自然長からの変位 y】
天井 ─── (自然長) ─── [伸び x₀ = mg/k (つりあい点)] ─── (任意の変位 y) ─── ●
復元力: F_total = mg - k(x₀ + y) = mg - kx₀ - ky = -ky
  • 証明:
    1. 自然長から質量 mm のおもりを吊るすと、重力とおもりの弾性力がつりあう位置で静止します。このつりあいの伸びを x0x_0 とおくと、kx0=mg  ⟺  x0=mgkk x_0 = mg \iff x_0 = \frac{mg}{k} です。
    2. この**つりあいの位置を新たな原点(y=0y = 0)**にとり、鉛直下向きに yy だけ変位した瞬間のおもりに働く合力 FF を求めます。
    3. ばねの全体の伸びは x0+yx_0 + y ですから: F=mg−k(x0+y)=mg−kx0−kyF = mg - k(x_0 + y) = mg - kx_0 - ky
    4. つりあいの条件 mg−kx0=0mg - kx_0 = 0 を代入すると、定数項が完全に消去されます。 F=−kyF = -ky
    5. したがって、運動方程式は md2ydt2=−kym \frac{d^2y}{dt^2} = -ky となり、水平ばねの式と全く同じになります。

結論: 重力の効果は「振動の中心をつりあいの位置 x0=mg/kx_0 = mg/k だけ下方にシフトさせる」だけであり、復元力の傾き(ばね定数 kk)には一切影響を与えないため、周期は水平時と完全に一致します。

ばねの並列接続と直列接続の合成ばね定数

複数のばねを組み合わせた場合の合成ばね定数 keffk_{\text{eff}} の導出です。

【並列接続:剛性2倍】           【直列接続:剛性半分】
   天井                           天井
  ┌──┴──┐                          │ ばね k₁
 ⌇ k₁   ⌇ k₂                      │ (伸び x₁)
  └──┬──┘                          ├──
     ● おもり                      │ ばね k₂
                                   │ (伸び x₂)
                                   ● おもり

1. 並列接続(Parallel)

おもりを xx だけ変位させると、2本のばねはどちらも同じ量 xx だけ伸びます。おもりに働く合力はそれぞれの復元力の和になります。 F=−(k1x+k2x)=−(k1+k2)xF = -(k_1 x + k_2 x) = -(k_1 + k_2) x したがって、合成ばね定数は単純な足し算になります。 kparallel=k1+k2k_{\text{parallel}} = k_1 + k_2 ばねが硬くなるため、周期は短縮されます。同じばね2本(k1=k2=kk_1 = k_2 = k)なら keff=2kk_{\text{eff}} = 2k となり、周期は Tparallel=T2≈0.707TT_{\text{parallel}} = \frac{T}{\sqrt{2}} \approx 0.707 T になります。

2. 直列接続(Series)

2本のばねを直列につないでおもりを引くと、ばねの各部分にかかる張力 FF は等しくなります。全体の伸び xx は各ばねの伸びの合計です。 x=x1+x2=Fk1+Fk2=F(1k1+1k2)x = x_1 + x_2 = \frac{F}{k_1} + \frac{F}{k_2} = F \left( \frac{1}{k_1} + \frac{1}{k_2} \right) したがって、逆数の和として合成ばね定数が求まります。 1kseries=1k1+1k2  ⟺  kseries=k1k2k1+k2\frac{1}{k_{\text{series}}} = \frac{1}{k_1} + \frac{1}{k_2} \iff k_{\text{series}} = \frac{k_1 k_2}{k_1 + k_2} ばね全体として柔らかくなるため、周期は延びます。同じばね2本なら keff=k2k_{\text{eff}} = \frac{k}{2} となり、周期は Tseries=2T≈1.414TT_{\text{series}} = \sqrt{2} T \approx 1.414 T になります。

ばねを切るとどうなるか?(切断の物理):
長さ LL、ばね定数 kk のばねをちょうど半分(L/2L/2)に切ると、ばね定数は 2k2k(2倍硬く) になります。
なぜなら、もとのばねは「半分サイズのばね2本が直列につながったもの」とみなせるからです(1/k=1/(2k)+1/(2k)1/k = 1/(2k) + 1/(2k))。半分の長さにすると、同じ力で引っ張ったときの伸びが半分になるため、剛性は2倍になります。


5. 単振動の変位・速度・加速度の位相差と力学的エネルギー保存

単振動の一般解から、変位・速度・加速度の数学的関係と力学的エネルギー保存のメカニズムを整理します。

変位・速度・加速度の微分関係

初期位相を ϕ\phi、振幅を AA とすると、単振動の基本解は以下のようになります。

  1. 変位 x(t)x(t): x(t)=Acos⁡(ωt+ϕ)x(t) = A \cos(\omega t + \phi)
  2. 速度 v(t)v(t)(変位の1階微分): v(t)=dxdt=−Aωsin⁡(ωt+ϕ)=Aωcos⁡(ωt+ϕ+π2)v(t) = \frac{dx}{dt} = -A\omega \sin(\omega t + \phi) = A\omega \cos\left(\omega t + \phi + \frac{\pi}{2}\right) 速度の最大値(振幅中心通過時)は vmax⁡=Aωv_{\max} = A\omega です。位相角に +π2+\frac{\pi}{2}(+90∘+90^\circ)が加わっているため、**「速度は変位より 90° 位相が進む」**といいます。
  3. 加速度 a(t)a(t)(速度の1階微分): a(t)=dvdt=−Aω2cos⁡(ωt+ϕ)=−ω2x(t)=Aω2cos⁡(ωt+ϕ+π)a(t) = \frac{dv}{dt} = -A\omega^2 \cos(\omega t + \phi) = -\omega^2 x(t) = A\omega^2 \cos(\omega t + \phi + \pi) 加速度の最大値(折り返し点・最大変位時)は amax⁡=Aω2a_{\max} = A\omega^2 です。変位に対してちょうど符号が反転(−x-x)しており、位相角が +π+\pi(+180∘+180^\circ)ずれているため、**「加速度は変位と 180° 逆位相」**になります。

位相差 90° と 180° の直観的イメージ

  • 最大変位点(x=+Ax = +A、折り返し点): 物体は一瞬静止するため速度はゼロ(v=0v = 0)ですが、ばねは最も強く引き戻そうとするため、負の向きに最大加速度(a=−amax⁡a = -a_{\max})が発生します。
  • 振動中心通過時(x=0x = 0): ばねの伸び縮みがないため力も加速度もゼロ(a=0a = 0)ですが、これまでに加速されてきた結果、速度は最大(∣v∣=vmax⁡|v| = v_{\max})に達します。

力学的エネルギー保存則の数式証明

おもりの運動エネルギー KK と、ばねの弾性位置エネルギー UU をそれぞれ計算します。

K=12mv(t)2=12m(−Aωsin⁡(ωt+ϕ))2=12mω2A2sin⁡2(ωt+ϕ)K = \frac{1}{2}mv(t)^2 = \frac{1}{2}m (-A\omega \sin(\omega t + \phi))^2 = \frac{1}{2}m \omega^2 A^2 \sin^2(\omega t + \phi)

ここで ω2=km  ⟺  mω2=k\omega^2 = \frac{k}{m} \iff m\omega^2 = k ですから:

K=12kA2sin⁡2(ωt+ϕ)K = \frac{1}{2}kA^2 \sin^2(\omega t + \phi)

一方、弾性位置エネルギー UU は:

U=12kx(t)2=12k(Acos⁡(ωt+ϕ))2=12kA2cos⁡2(ωt+ϕ)U = \frac{1}{2}kx(t)^2 = \frac{1}{2}k (A \cos(\omega t + \phi))^2 = \frac{1}{2}kA^2 \cos^2(\omega t + \phi)

両者の合計(全力学的エネルギー EE)をとると:

E=K+U=12kA2(sin⁡2(ωt+ϕ)+cos⁡2(ωt+ϕ))E = K + U = \frac{1}{2}kA^2 \left( \sin^2(\omega t + \phi) + \cos^2(\omega t + \phi) \right)

三角関数の基本関係式 sin⁡2θ+cos⁡2θ=1\sin^2\theta + \cos^2\theta = 1 により、括弧の中身は常に 1 になります。

E=12kA2=12mvmax⁡2=一定E = \frac{1}{2}kA^2 = \frac{1}{2}m v_{\max}^2 = \text{一定}

物理的結論: 運動エネルギー KK と位置エネルギー UU は、周期 T/2T/2(運動自体の半分の時間)ごとに互いに姿を変えて行き来しますが、その総和 EE は時間に一切依存せず完全に一定に保たれます。


6. 実践計算例題 5選

理論式を定着させるための、具体的かつ典型的な5つの計算問題です。

例題1:秒振子(時計標準の振り子)の必要長

問題: 周期がちょうど T=2.00 sT = 2.00\,\text{s}(片道 1.00 s でカチッと進む秒振子)となる単振り子を、標準重力加速度 g=9.80665 m/s2g = 9.80665\,\text{m/s}^2 の場所で作る場合、必要な糸の長さ LL は何 cm か?

  • 計算: 周期公式 T=2πL/gT = 2\pi\sqrt{L/g} を LL について解きます。 L=gT24π2=9.80665×(2.00)24×π2=9.80665×4.0039.4784=39.226639.4784≈0.9936 m=99.36 cmL = \frac{g T^2}{4\pi^2} = \frac{9.80665 \times (2.00)^2}{4 \times \pi^2} = \frac{9.80665 \times 4.00}{39.4784} = \frac{39.2266}{39.4784} \approx 0.9936\,\text{m} = 99.36\,\text{cm}
  • 解答: 約 99.4 cm(微小角・理想振り子モデルでの計算値です)。

例題2:長い理想振り子の周期

問題: 仮想的な長い振り子の有効長が L=67.0 mL = 67.0\,\text{m}、おもりの質量が m=28.0 kgm = 28.0\,\text{kg} です。この問題での重力加速度を g=9.809 m/s2g = 9.809\,\text{m/s}^2 とするとき、この振り子の周期 TT は何秒か?また、おもりを 56 kg に倍増させた場合、周期はどう変化するか?

  • 計算: T=2πLg=2π67.09.809=2π6.8305≈2×3.14159×2.6135≈16.42 sT = 2\pi \sqrt{\frac{L}{g}} = 2\pi \sqrt{\frac{67.0}{9.809}} = 2\pi \sqrt{6.8305} \approx 2 \times 3.14159 \times 2.6135 \approx 16.42\,\text{s}
  • 解答: 周期は約 16.4 秒。おもりの質量を 56 kg に変えてもこの理想モデルでは周期は変化しません(質量非依存性)。

例題3:実験室のばね振り子の最大速度と最大加速度

問題: ばね定数 k=25.0 N/mk = 25.0\,\text{N/m} のばねに質量 m=250 gm = 250\,\text{g}(=0.250 kg= 0.250\,\text{kg})のおもりを取り付け、振幅 A=5.0 cmA = 5.0\,\text{cm}(=0.050 m= 0.050\,\text{m})で水平方向に振動させた。周期 TT、最大速度 vmax⁡v_{\max}、最大加速度 amax⁡a_{\max} を求めよ。

  • 計算:
    1. 角振動数:ω=km=25.00.250=100=10.0 rad/s\omega = \sqrt{\frac{k}{m}} = \sqrt{\frac{25.0}{0.250}} = \sqrt{100} = 10.0\,\text{rad/s}
    2. 周期:T=2πω=2×3.1415910.0≈0.628 sT = \frac{2\pi}{\omega} = \frac{2 \times 3.14159}{10.0} \approx 0.628\,\text{s}
    3. 最大速度:vmax⁡=Aω=0.050×10.0=0.500 m/sv_{\max} = A\omega = 0.050 \times 10.0 = 0.500\,\text{m/s}
    4. 最大加速度:amax⁡=Aω2=0.050×(10.0)2=5.00 m/s2a_{\max} = A\omega^2 = 0.050 \times (10.0)^2 = 5.00\,\text{m/s}^2
  • 解答: 周期 0.628 s、最大速度 0.500 m/s、最大加速度 5.00 m/s²。

例題4:ばねの並列接続 vs 直列接続の比較

問題: ばね定数 k1=30.0 N/mk_1 = 30.0\,\text{N/m} と k2=60.0 N/mk_2 = 60.0\,\text{N/m} の2本のばねがある。質量 m=0.500 kgm = 0.500\,\text{kg} のおもりを (a) 並列に接続した場合、(b) 直列に接続した場合の合成ばね定数と周期をそれぞれ求めよ。

  • 計算:
    • (a) 並列接続: kparallel=k1+k2=30.0+60.0=90.0 N/mk_{\text{parallel}} = k_1 + k_2 = 30.0 + 60.0 = 90.0\,\text{N/m} Tparallel=2π0.50090.0=2π0.005556≈2π×0.07454≈0.468 sT_{\text{parallel}} = 2\pi \sqrt{\frac{0.500}{90.0}} = 2\pi \sqrt{0.005556} \approx 2\pi \times 0.07454 \approx 0.468\,\text{s}
    • (b) 直列接続: kseries=k1k2k1+k2=30.0×60.030.0+60.0=180090.0=20.0 N/mk_{\text{series}} = \frac{k_1 k_2}{k_1 + k_2} = \frac{30.0 \times 60.0}{30.0 + 60.0} = \frac{1800}{90.0} = 20.0\,\text{N/m} Tseries=2π0.50020.0=2π0.0250≈2π×0.15811≈0.993 sT_{\text{series}} = 2\pi \sqrt{\frac{0.500}{20.0}} = 2\pi \sqrt{0.0250} \approx 2\pi \times 0.15811 \approx 0.993\,\text{s}
  • 解答: 並列時は k=90.0 N/mk = 90.0\,\text{N/m}、周期 0.468 s。
    直列時は k=20.0 N/mk = 20.0\,\text{N/m}、周期 0.993 s(直列にすると約2.1倍周期が長くなります)。

例題5:鉛直ばね振り子のつりあいの伸びと力学的エネルギー

問題: ばね定数 k=50.0 N/mk = 50.0\,\text{N/m} の鉛直ばねにおもり m=0.400 kgm = 0.400\,\text{kg} を吊るした。重力加速度 g=9.80 m/s2g = 9.80\,\text{m/s}^2 とする。

  1. つりあいの位置でのばねの伸び x0x_0 は何 cm か?
  2. つりあいの位置からさらに下へ A=6.0 cmA = 6.0\,\text{cm}(=0.060 m= 0.060\,\text{m})引っ張って静かに放したとき、振動の力学的エネルギー EE は何 J か?
  • 計算:
    1. つりあいの伸び: x0=mgk=0.400×9.8050.0=3.9250.0=0.0784 m=7.84 cmx_0 = \frac{mg}{k} = \frac{0.400 \times 9.80}{50.0} = \frac{3.92}{50.0} = 0.0784\,\text{m} = 7.84\,\text{cm}
    2. 単振動の力学的エネルギー(つりあいの位置を基準とするエネルギー): E=12kA2=12×50.0×(0.060)2=25.0×0.00360=0.0900 J=90.0 mJE = \frac{1}{2}kA^2 = \frac{1}{2} \times 50.0 \times (0.060)^2 = 25.0 \times 0.00360 = 0.0900\,\text{J} = 90.0\,\text{mJ}
  • 解答: つりあいの伸び 7.84 cm、力学的エネルギー 0.090 J (90.0 mJ)。

7. 無料計算ツール(Tools Hub)での検算・レポート作成手順

手計算による導出やレポートの数値を、ブラウザ上でリアルタイムに検証・作図できる 単振動・振り子計算ツール の活用方法です。

[単振動・振り子計算ツールの画面構成]
┌────────────────────────────────────────────────────────┐
│ [単振り子タブ]  [ばね振り子タブ]  [波形・エネルギー保存タブ] │
├──────────────────────────┬─────────────────────────────┤
│ ・目標選択(周期/重力/長さ) │ ・リアルタイム往復シミュレーション  │
│ ・長さ L (m / cm)         │ ・速度ベクトル表示           │
│ ・振れ角 θ (度)           │ ・周期 T、ボルダ補正 T_corr   │
│ ・重力プリセット (地球/月/火星)│ ・ステップ別導出数式          │
├──────────────────────────┴─────────────────────────────┤
│ [ レポート用テキストをコピー ]  [ 初期値にリセット ]            │
└────────────────────────────────────────────────────────┘

1. 単振り子の重力逆算を検算する場合

  1. ツール上部のタブから**「単振り子」**を選択します。
  2. 計算目標を**「実験周期から重力加速度 g を逆算」**に切り替えます。
  3. 有効長 LL(例: 1.00 m)と、測定した1周期 TT(例: 2.006 s)、振れ角 θ0\theta_0(例: 5.0°)を入力します。
  4. 出力エリアに、ボルダ補正から近似的に復元した微小角周期 T0T_0 と、逆算された重力加速度 g≈9.8200 m/s2g \approx 9.8200\,\text{m/s}^2 が途中式付きで即座に表示されます。

2. ばねの直列・並列合成を比較する場合

  1. **「ばね振り子」**タブをクリックします。
  2. 接続方式で**「並列(Parallel)」または「直列(Series)」**を選択し、それぞれのばね定数 k1,k2k_1, k_2 とおもり質量 mm を入力します。
  3. 合成ばね定数 keffk_{\text{eff}}、周期 TT、鉛直つりあいの伸び x0=mg/kx_0 = mg/k、最大速度 vmax⁡v_{\max} が一瞬で計算されます。
  4. 画面右側のシミュレーターで、ばねコイルとおもりブロックが実際に伸縮運動する様子を確認できます。

3. レポート用テキストのワンクリックコピー

計算結果の下部にある**「レポート用テキストをコピー」**ボタンを押すと、入力パラメータ、合成値、周期、最大速度、全エネルギー、ステップ別導出式が整形された Markdown テキストとしてクリップボードにコピーされ、実験レポートや課題提出ファイルにそのまま貼り付けることができます。


8. よくある疑問・つまずきポイント FAQ

Q1. 単振り子の糸に質量がある場合、周期はどう変わりますか?

糸を一様な剛体の棒として近似でき、小角で揺れるモデルでは、糸の質量 mstringm_{\text{string}} が無視できない場合、系は単純な質点ではなく**「剛体振り子(物理振り子)」**として扱わなければなりません。
支点のまわりの慣性モーメントが I=mbobL2+13mstringL2I = m_{\text{bob}} L^2 + \frac{1}{3} m_{\text{string}} L^2 となり、周期公式は T=2πI/(mtotalgdCM)T = 2\pi\sqrt{I / (m_{\text{total}} g d_{\text{CM}})} で計算されます。糸の質量が加わると、重心位置 dCMd_{\text{CM}} がやや支点寄りに上がる一方、慣性モーメントも増加するため、通常の単振り子よりも周期がわずかに変化します。

Q2. 振り子をエレベーターの中に持ち込むと周期はどうなりますか?

エレベーターが上向きに加速度 α\alpha で加速しているとき、慣性力 −mα-m\alpha が下向きに加わるため、見かけの重力加速度は g′=g+αg' = g + \alpha に増加します。したがって、周期は T=2πL/(g+α)T = 2\pi\sqrt{L / (g + \alpha)} となり、**通常より短く(速く)**なります。
逆に、エレベーターが下向きに加速度 α\alpha で加速すると見かけの重力は g−αg - \alpha に減少し、周期は長くなります。ワイヤーが切れて自由落下(α=g\alpha = g)すると見かけの重力はゼロになり、振り子は往復運動しなくなります(周期無限大)。

Q3. ばね振り子を宇宙ステーション(無重力)に持っていくと振動しますか?

全く問題なく同じ周期で振動します。
ばね振り子の復元力はばね自身の弾性力 F=−kxF = -kx のみによって生じており、重力に依存しません。宇宙ステーション内(無重力環境)でも、おもりの質量を mm、ばね定数を kk とし、線形性・取付条件が同じなら、理想モデルでは地上と同じ T=2πm/kT = 2\pi\sqrt{m/k} の周期で振動します。実際に宇宙飛行士の体重測定器(質量測定器)には、このばね振り子の周期から質量を逆算する原理が応用されています。

Q4. 単振動で「速度が変位より90°進む」とはグラフ上でどういうことですか?

時間 tt を横軸にとった波形グラフを見ると、変位 x(t)x(t) が山のピーク(最大変位 +A+A)に達するより 1/4 周期(時間で T/4T/4、位相角で 90∘90^\circ または π/2\pi/2)早いタイミングで、速度 v(t)v(t) のピークが現れることを意味します。変位がピークのときは折り返し点なので速度はゼロですが、その 1/4 周期前に中心を通過した瞬間に速度が正の最大値をとっていたため、時間的に速度のピークが先行して現れます。

Q5. 減衰振動(空気抵抗がある場合)では周期はどうなりますか?

空気抵抗や摩擦などの速度比例抵抗力 Fdrag=−γvF_{\text{drag}} = -\gamma v が働く実際の実験では、振動は次第に小さくなる減衰振動(Damped Oscillation)になります。
不足減衰(γ/(2m)<ω0\gamma/(2m) < \omega_0)で振動が続く場合、このときの減衰角振動数は ω′=ω02−(γ/2m)2\omega' = \sqrt{\omega_0^2 - (\gamma / 2m)^2} となり、抵抗のない理想的な角振動数 ω0\omega_0 よりも小さくなります。したがって、空気抵抗があると周期 T′=2π/ω′T' = 2\pi / \omega' は理想値より長くなります(弱い減衰なら差は小さくなります)。臨界減衰・過減衰では往復振動が続かないため、この周期式は使えません。


9. 参考文献・関連ツール

信頼できる学術・一次資料

Tools Hub の関連計算ツール

測定の不確かさと再利用のチェックリスト

長さLと周期T₀が独立な標準不確かさを持ち、角度補正の不確かさを無視できる場合、一次近似では u(g)/g≈[u(L)/L]2+[2u(T0)/T0]2u(g)/g \approx \sqrt{[u(L)/L]^2+[2u(T_0)/T_0]^2} です。長さのずれ、測定時刻のばらつき、振幅の減衰、補正の打ち切りを別々に考察します。相関や角度の不確かさが無視できない場合はこの簡略式だけでは足りません。

  1. ツールの計算結果をコピーし、実測値・仮定・単位・補正次数をレポートに追記します。コピー文だけを実験結果の証拠にしないでください。
  2. ツール名付近の「お気に入り」で登録すると、同じブラウザのマイツールから再訪できます。入力データは保存されません。
  3. 次の測定では同じ設定かを確認し、結果を前回の記録と比較します。

2026年10月24日に本文の検算例と公開実装のモード・補正次数を照合しました。参考資料は理論の根拠であり、本サイトの校正や専門家による実験装置の認証を意味しません。