既定値は L=1.0 m、θ₀=45°、g=9.81 m/s²、m=1.0 kg。振り子アニメーションは厳密周期 T に基づき実時間で 1 周期を 1 振動として表示し、最大振幅位置（±θ₀）を点線で示します。質量 m は周期に影響しませんが位置エネルギーには寄与します。
固定点から長さ L で吊るされた質量 m の振り子。青球が現在角、左右の灰色点線が最大振幅 ±θ₀ の停止位置。再生中は厳密周期 T で実時間振動。
横軸 振幅 θ₀ (°) [0–175]、縦軸 周期誤差 (%)。青曲線は楕円積分による厳密誤差、橙点線は級数展開 1 項目 θ₀²/16。黄マーカーが現在の動作点 (θ₀, 誤差)。θ₀→180° で発散。
振り子の周期は小振幅近似と厳密解で大きく異なります。それぞれを並べて掲載します。
小角度近似周期（等時性）：
$$T_{0} = 2\pi\sqrt{\dfrac{L}{g}}$$
厳密周期（第1種完全楕円積分 K(k)）：
$$T = 4\sqrt{\dfrac{L}{g}}\,K\!\left(\sin\dfrac{\theta_{0}}{2}\right),\quad K(k)=\int_{0}^{\pi/2}\dfrac{d\varphi}{\sqrt{1-k^{2}\sin^{2}\varphi}}$$
級数展開と位置エネルギー：
$$T \approx T_{0}\left(1+\dfrac{\theta_{0}^{2}}{16}+\dfrac{11\theta_{0}^{4}}{3072}+\dfrac{173\theta_{0}^{6}}{737280}+\cdots\right),\qquad E = m g L (1-\cos\theta_{0})$$
$L$ は振り子の長さ [m]、$\theta_{0}$ は振幅 [rad]、$g$ は重力加速度 [m/s²]、$m$ は質量 [kg]。$T_{0}$ は質量・振幅に依存しない（等時性）が、$T$ は $\theta_{0}$ が大きいほど長くなり、$\theta_{0}\to\pi$ で発散します。位置エネルギー $E$ は最高点（振幅 $\theta_{0}$）での値です。
大振幅振り子 シミュレーターとは
🙋
既定値 θ₀=45° で T₀=2.006 s、T=2.086 s、誤差 4.00% って出てます。たった 45° で 4% も違うんですか？「振り子は等時性」って習ったので意外です。
🎓
いい疑問だね。「等時性」が成り立つのは sinθ≈θ と近似できる小振幅領域（だいたい θ₀<10° くらい）。θ₀=45° では sinθ が θ より小さく、復元力が弱くなる分だけ周期が伸びる。級数展開 T≈T₀(1+θ₀²/16+...) の第 1 項だけで θ₀=45°(0.785 rad) なら +(0.785)²/16=3.85% と既に近い値が出る。本ツールは AGM 法で完全楕円積分 K(k) を直接計算して厳密値 4.00% を返している。ガリレオが 1602 年に「等時性」を発見できたのは、振幅 5° 程度の小さな揺れを観察したから成立した、いわば近似的な経験則だったんだ。
🙋
スライダーで θ₀=175° まで上げると、誤差が >1000% になって振り子が画面のほぼ真上で止まりかけてます。これは何が起きてるんですか？
🎓
θ₀→180° で K(1)=∞ つまり周期そのものが発散する。物理的には「真上の不安定平衡点」に到達するまで無限時間かかることを意味する——エネルギー保存則で θ=180° では運動エネルギー=0 となる瞬間、わずかなずれが指数関数的に成長してそこから先に行けない（または行ったら元に戻れない）。位相空間（θ-ω 平面）ではこの軌道が「分離線（separatrix）」と呼ばれ、閉じた振動軌道と回転軌道（逆さ回り）を分ける境界線になる。鞍点（saddle point）の周りで時間が「ゆっくり進む」古典的な例だね。
🙋
青曲線（厳密）と橙点線（級数 1 項）のズレが θ₀=60°くらいから急に開きますね。級数って何項取れば実用的になりますか？
🎓
良い観察。級数 T≈T₀(1+θ₀²/16+11θ₀⁴/3072+173θ₀⁶/737280+22931θ₀⁸/1321205760+...) の収束半径は θ₀=π。θ₀=45° なら 2 項で 0.005% 精度、4 項で機械精度。だが θ₀=120°(2.09 rad) では 1 項目 27.4% に対し厳密 31.7% でズレが大きく、5〜6 項必要。実務的には θ₀<60° なら 2〜3 項展開で十分、それ以上は AGM 法による楕円積分が計算量・精度の両面で優位。本ツールは安定性を優先して全域で AGM を使い、30 回反復で機械精度に達している。
🙋
質量 m を変えると位置エネルギー E が変わりますが、周期 T₀・T は変わらないですね。なぜ振り子の周期は質量に依存しないんですか？
🎓
運動方程式 m·L·θ'' = -m·g·sinθ から両辺 m を消去すると θ'' = -(g/L)sinθ となり、質量 m が完全に消える——これは「慣性質量と重力質量が等しい」（等価原理）から来る性質で、ガリレオが「異なる質量の球が同時に落下する」と示したのと同じ理屈。アインシュタインの一般相対論の出発点にもなった原理だ。一方で位置エネルギー E=mgL(1−cosθ₀) は m に比例するから、エネルギー貯蔵デバイス（重力エネルギー貯蔵）の設計では質量を増やすことに意味がある。スイス Energy Vault 社は 35 トンのコンクリートブロックを高所に上げ下げして MWh 級の電力を貯蔵するシステムを商業化しているよ。
物理モデルと主要な数式
長さ $L$ の糸で質量 $m$ を吊るした剛体振り子（理想振り子）の運動方程式は、ニュートンの第二法則から $m L \ddot{\theta} = -m g \sin\theta$、すなわち $\ddot{\theta} = -(g/L)\sin\theta$ となります。小振幅近似 $\sin\theta\approx\theta$ では $\ddot{\theta}+(g/L)\theta=0$ の調和振動子方程式となり、固有角振動数 $\omega_{0}=\sqrt{g/L}$、周期 $T_{0}=2\pi/\omega_{0}=2\pi\sqrt{L/g}$ が得られます。質量・振幅に依存しない等時性が成立します。
大振幅では $\sin\theta$ を保持したまま、エネルギー保存則 $\tfrac{1}{2}L^{2}\dot{\theta}^{2}+gL(1-\cos\theta)=gL(1-\cos\theta_{0})$ から $\dot{\theta}=\pm\sqrt{(2g/L)(\cos\theta-\cos\theta_{0})}$。これを 0 から $\theta_{0}$ まで積分し 4 倍すると周期：
$$T = 4\int_{0}^{\theta_{0}}\dfrac{d\theta}{\sqrt{(2g/L)(\cos\theta-\cos\theta_{0})}} = 4\sqrt{\dfrac{L}{g}}\,K\!\left(\sin\dfrac{\theta_{0}}{2}\right)$$
ここで $K(k)=\int_{0}^{\pi/2}d\varphi/\sqrt{1-k^{2}\sin^{2}\varphi}$ は第1種完全楕円積分です。$k=\sin(\theta_{0}/2)\to 1$ で $K(k)\to\infty$ となり、周期が発散します。級数展開で表すと $T=T_{0}(1+\theta_{0}^{2}/16+11\theta_{0}^{4}/3072+173\theta_{0}^{6}/737280+\cdots)$ で、第1項 $\theta_{0}^{2}/16$ が代表的な「大振幅補正」です。最高点での位置エネルギーは $E=mgL(1-\cos\theta_{0})$ で、質量に比例しますが周期には現れません。本ツールは AGM 法（arithmetic-geometric mean）で $K(k)$ を機械精度で計算します。
実世界での応用
振り子時計の精度設計： ホイヘンス（1656年）以来の振り子時計は、脱進機（escapement）で「振幅をほぼ一定に保つ」ことで等時性を実現します。Riefler や Shortt の精密天文時計は振幅 1.5°〜2° に厳密制御され、周期誤差は θ²/16 から (0.026)²/16=4.3×10⁻⁵=0.0043% に抑え込まれ、1 日 1/100 秒の精度を達成しました。本ツールに θ₀=2° を入れると誤差 0.0076% でこれを再現できます。GPS や原子時計が登場するまでの 300 年間、最も正確な時間標準だった理由がここにあります。
地震応答解析と TMD（同調質量ダンパー）： 建物の 1 次固有振動モードは振り子に似た形で、超高層ビルの頂部に振り子型 TMD を設置して逆位相で揺らし主構造の揺れを吸収します。台北 101（660 トン振り子）、上海中心大廈（1000 トン）など。設計では「想定地震動の振幅で TMD 振り子が周期 4% ずれて主構造との同調を失う」リスクを評価し、本ツールの大振幅補正がそのまま指針になります。CAE では非線形動的解析（Newmark-β や中央差分法）で sinθ をそのまま積分します。
サーカディアン・体内時計の生物物理モデル： 哺乳類の生体時計は転写・翻訳負帰還ループの非線形振動子で、振幅依存的に周期が変動する「等時性の破れ」を持ちます。光刺激で振幅を変えると周期が 0.5〜1 時間シフトする現象が、本ツールの大振幅振り子と同じ「非線形振り子モデル」で記述されます。Hodgkin-Huxley 型ニューロン発火モデルでも同じ数学構造が現れ、てんかん発作の理論解析に応用されています。
重力測定と地球の形： 振り子重力計（Kater 振り子、1818 年）は $g=4\pi^{2}L/T^{2}$ で重力加速度を 0.001% 精度で測定し、20 世紀前半まで地球の重力分布図作成の標準器でした。本ツールで $L=1$ m、$T=2.006$ s から逆算すると $g=4\pi^{2}\cdot 1/2.006^{2}=9.81$ m/s² が復元できます。現在は超伝導重力計や絶対重力計（自由落下式）が主流ですが、振り子重力計のアイデアは地下空洞・鉱床探査や地殻変動モニタリングに今も継承されています。
よくある誤解と注意点
最も多い誤解は、「振り子の周期は常に T=2π√(L/g) で振幅に依存しない」 です。これは小振幅 (θ₀≪1 rad) の極限でのみ正しい近似で、$\theta_{0}=10°$ で 0.19%、$45°$ で 4.00%、$90°$ で 18.1%、$120°$ で 32%、$170°$ で 290% と急速に発散します。本ツールの誤差曲線で「等時性が実用上成り立つのは θ₀<10° 程度」がひと目で分かります。物理教科書の「振り子の等時性」は前提として「微小振動の極限」を必ず明記しているので、応用時には振幅の桁を確認してください。
次に多いのが、「質量を重くすれば周期が短くなる（または長くなる）」 という誤解です。理想振り子では運動方程式から質量 $m$ が完全に消えるため、周期は質量に依存しません（等価原理）。本ツールで $m$ スライダーを動かすと位置エネルギー $E$ だけが変化し、$T_{0}$ と $T$ は固定であることが確認できます。ただし「実物の振り子」では空気抵抗が質量効率に比例して大きくなったり（小球は減衰しやすい）、糸の質量を無視できなかったり（実体振り子モデルが必要）するため、極限的にこの法則が破れる場合があります。
最後に、「楕円積分は閉じた形で表せないので近似計算しかできない」 という思い込みです。確かに初等関数では書けませんが、AGM 法（算術幾何平均）を使うと反復回数 $n$ につき桁数が 2 倍ずつ増える「二次収束」が得られ、30 回反復で機械精度（10⁻¹⁵）に達します。Python の scipy.special.ellipk や C 言語の boost::math::ellint_1 も内部は AGM 法です。本ツールも 50 回反復で θ₀=174.999° まで安定計算できます。級数展開は θ₀<60° までは実用的ですが、大振幅では多項必要で AGM 法が圧倒的に高速です。
よくある質問
大振幅振り子の周期はなぜ振幅に依存するのですか？
運動方程式 d²θ/dt² = -(g/L)sinθ には非線形項 sinθ が含まれます。小振幅では sinθ≈θ と近似でき周期 T₀=2π√(L/g) が振幅に依存しない「等時性」が成り立ちますが、大振幅では sinθ < θ となり復元力が相対的に弱まるため、周期が長くなります。厳密には T=4√(L/g)·K(sin(θ₀/2)) と第1種完全楕円積分 K(k) で表され、本ツールに既定値 L=1.0 m、θ₀=45°、g=9.81 m/s² を入れると T₀=2.006 s、T=2.086 s、誤差 4.00% が得られます。
完全楕円積分 K(k) はどのように計算しますか？
第1種完全楕円積分 K(k)=∫₀^(π/2) dφ/√(1-k²sin²φ) は閉じた初等関数で表せませんが、算術幾何平均 (AGM) 法を使うと数回の反復で機械精度に達します。a₀=1、b₀=√(1-k²) から始め、a_{n+1}=(a_n+b_n)/2、b_{n+1}=√(a_n·b_n) を繰り返し、収束値 M(1,√(1-k²)) を得て K=π/(2M)。本ツールは AGM 法で 50 回反復し、θ₀=175° まで正確に計算します。級数展開 T≈T₀(1+θ₀²/16+11θ₀⁴/3072+...) も使えますが、θ₀ が大きいと収束が遅くなります。
振幅 θ₀=90° と 175° では周期はどう変わりますか？
L=1.0 m、g=9.81 m/s² で θ₀=90° のとき K(sin45°)=1.8541、T=4·0.3193·1.8541=2.369 s、誤差 18.1%。θ₀=175° では k=sin(87.5°)≒0.9990 で K≒5.435、T≒6.943 s、誤差 246%（周期が約 3.5 倍）になります。θ₀→180° では K(1)=∞ で周期が発散し、これは「不安定平衡点」（真上）に到達するのに無限時間かかることを意味します。位相空間ではこの極限が「分離線（separatrix）」として現れます。
級数展開と楕円積分のどちらを使うべきですか？
θ₀ が小さい（< 30° 程度）なら級数展開 T≈T₀(1+θ₀²/16+11θ₀⁴/3072+173θ₀⁶/737280+...) で 1〜2 項取れば 0.01% 精度が得られます。θ₀=45° で 2 項展開すると T/T₀≒1.0386 で誤差 0.005%、厳密値とほぼ一致。一方 θ₀>90° では級数収束が遅くなり項を多数必要とするため、AGM 法による楕円積分の直接計算が高速・高精度で推奨されます。本ツールは内部で AGM 法を採用し、θ₀=1〜175° の全域で精密に厳密周期を計算します。