5
エンタメ議論
文献あり

【自由研究】単振り子の周期の厳密解

1463
0

この記事は個人で行った実験結果と考察をもとに書かれています.エンタメとしてご覧ください.

はじめに

この実験は,「第19回全国物理コンテスト物理チャレンジ2023」の「第1チャレンジ2023実験レポート課題」の実験課題に基づき行われています.

振り子の周期を、振れ角を変えて調べてみよう
単振り子の周期を測定する実験では振れ角を小さくして行います が、振れ角を大きくしていくと周期はどうなるでしょうか。振り子の 振れ角と周期の関係について、実験を行って調べてみましょう。 周期を正確に測る工夫、減衰の小さい振り子の作製、振れ角と周期 の関係についての考察などを期待しています。

物理チャレンジについてはこちらから: 全国物理コンテスト 物理チャレンジ!

概要

参考文献[1]によるとによると,単振り子の周期は「振れが小さいとき,周期は振幅によらず,糸の長さと重力加速度の大きさだけで決まる.これを振り子の当時性という.」とある.つまり,
単振り子の振幅が微小のとき,周期は,
T=2πlg
のように表されるが,振幅が一般の角のとき,周期は楕円積分を用いると,
T=2πlg{n=0[(2n1)!!(2n)!!]2n2nθ02}
と表される.ただし,!!は二重階乗.これを利用して以下の実験を行った.

実験1

上記の周期の厳密解にどれほどの正確性があるのかを,振幅θ0[°]と弦の長さl[m]を変えて,カメラで単振り子が振れる様子を撮影し,周期を測定することによって調べた.

実験2

振幅θ0[°]を変えて,周期の変化について,カメラで単振り子が振れる様子を撮影し,周期を測定する事によって調べた.
実験2から,周期と振幅の関係について,周期が初期の振幅の六次関数として表されることを発見し,その理由を実験1と理論的考察から推測した.
また,空気抵抗と慣性モーメントについての考察も行い,単振り子の測定した周期と理論的な周期のずれについて徹底的に検証を行った.

実験の目的

今回の実験の目的を以下のように設定した.

  • 単振り子の周期の厳密解にどれほどの正確性があるのか調べる.
  • 周期と振幅の関係を調べる.
  • 単振り子の周期において,空気抵抗と慣性モーメントが与える影響を調べる.

実験手法

実験に使用したもの

実験道具(実験自体に使用したもの)
  • 単振り子(島津理化機械株式会社製 TN‐3 形,ケニス株式会社製 No.110-0210 単振り
    子 SP ).金属球:スチールボール1 5/8(直径:2a = 413mm,質量:m = 286g).弦:
    ピアノ線(太さ:0.35mm).
    単振り子(金属球を使用した) 単振り子(金属球を使用した)
    単振り子(弦を使用した) 単振り子(弦を使用した)

  • 鉄製スタンド

  • 分度器

  • C型クランプ(鉄製スタンドを机に固定するのに使用)

  • 3.5mメジャー

  • 段ボールとコピー用紙A4で作成した角度確認用分度器
    角度確認用分度器 角度確認用分度器

  • ペン

  • 卓上ボール盤(段ボールに穴をあけるのに使用)

測定道具(実験の測定及び解析に使用したもの)
  • iPhoneXR(撮影用)
  • 三脚
  • パソコン(ストップウォッチ用と解析用)

実験原理

以下,単振り子の弦の長さをl[m](この節では,単振り子のおもりを質点として扱うが,実験においては,おもりの半径を加えた長さを弦の長さlとする.),振幅をθ[rad]θは振り子が振れ始めてからの経過時間t[s]の関数),周期をT[s],小球の質量をm[kg],重力加速度をg[m/s2],最下点の位置をOO点から円弧に沿った変位をx[m]と表す.
振幅が非常に小さいとき,単振り子の周期は以下のように表される.
(1)T=2πlg

なお,周期がこのように表される理由は以下のようである.

円弧の接線方向の運動方程式を解く.
md2xdt2=mgsinθd2xdt2=gsinθ
x=lθであるので,
(2)ld2θdt2=gsinθd2θdt2=glsinθ
θが十分に小さいとき,sinθθとして考えられるので,
d2θdt2=glθ
この二階微分方程式を解くと,
θ(t)=Csin(glt+ϕ)
ただし,C,ϕは任意定数.そしてこの式は単振動を表すので,角速度ωは,ω=glとなり,周期Tは,
T=2πω=2πgl=2πlg
となる.

では,θを任意の大きさの角として,厳密解を求めるとどうなるのか.以下,簡略化のためにθの時間一階微分をθ˙,時間二階微分をθ¨と記す.式(2)の両辺にθ˙をかけて,右辺を左辺に移行する.
θ¨θ˙+glθ˙sinθ=0ddt[12θ˙2glcosθ]=0
よって,12θ˙2glcosθは定数となるので,定数Eと置く.(EはEnergyの頭文字である.実際,位置エネルギーの基準点を最下点に置くと,振り子の持つ力学的エネルギー(位置エネルギーと運動エネルギーの和)は,12ml2θ˙2+mgl(1cosθ)となり,上の式と対応している.)
θ(0)=θ0のとき,θ˙=0だとすると,
E=glcosθ0
よって,
12θ˙2gl(cosθcosθ0)=0θ˙2=2gl(cosθcosθ0)
いま,θ=θ0からθ=0までの運動を考えると,
θ˙=2gl(cosθcosθ0)
さらに,θ=θ0からθ=0までの運動にかかる時間はT4であるので,
T=40T4dt=4θ0012gl(cosθcosθ0)dθ(θ˙=dθdtθ˙dt=dθ)=22lg0θ01cosθcosθ0dθ
ここで,sinθ2=sinθ02sinϕを満たすϕに変数変換すると,
T=4lg0π2dϕ1sin2θ02sin2ϕ
となり,第一種完全楕円積分(楕円積分の節を参照)の形となる.よって楕円積分の節で求めた結果を用いて計算すると,
T=4lg0π211sin2θ02sin2ϕdϕ=4lgK(sinθ02)(3)=2πlg{n=0[(2n1)!!(2n)!!]2sin2nθ02}
となる。ただし,!!は二重階乗.

なお、簡略化のため、関数σ(θ0)を以下のように定める.
σ(θ0)=n=0[(2n1)!!(2n)!!]2sin2nθ02

Pythonを利用した級数の計算方法

前節からTの厳密解(3)が分かったが,無限級数の形で表されているため,実際に計算しつくすことは不可能である.そのため,できるだけ理論値の正確性を向上させるために,プログラミング言語Pythonを利用して,Tの数値計算を行う.以下のコードを使用する.ただし定数aは角度であり,常々設定する.
使用したソースコード 使用したソースコード
これによって得られた,σ(θ0)の値は以下の表の通りである
得られた値 得られた値
一個の角度あたり、300~600秒くらいの計算時間がかかった.

実験に使用した値

  • 測定場所における重力加速度 gs=9.796895m/s2[1]
  • 円周率 π=3.1415926535

実験手順

大まかな実験手順は以下の通りである.実験ごとの細かい解析手順についてはその下で示す.

  1. 弦の長さ(金属球の中心から弦についたリングと支持棒の接点までの長さ)を調べる対象に合わせて測り,固定.
  2. 実験装置を図のようにセットする.
  3. 鉄製スタンド横に設置したパソコンのストップウォッチをスタートさせておく.
  4. 実験者Aが振り子の位置を調節し,実験者BがiPhoneのカメラから振り子の位置を確認し,適切な位置を指示する.
  5. 実験者BがiPhoneのスローカメラ[2]をスタートさせ,実験者Aが振り子をそっと放し,振り子が5回以上振れたときに実験者Bがカメラを止める.
  6. 撮影した動画を確認し,振り子が静止した状態で振り子の中心を通るように鉛直方向に引いた線(線Iとする.)を利用して,初めの時間(振り子が線Iを初めて通過したときの時間.)と終わりの時間(振り子が5回振れたときに線Iを通過した時間.)を計測することによって,その差から5回振れる間にかかった時間(5T)を測定する.
    実験装置 実験装置
実験1

実験1では,目的1のため,以下の表のような場合について実験を行う.それぞれの場合について,5回ずつ行う.その後の手順は以下の通りである.
実験1の実験対象 実験1の実験対象

  1. (3)について,σ(θ0)=σとして,変形すると,
    (4)T=2πlgσ(θ0)g=4π2lσ2T2
    となる.
  2. また,式(1)について,変形すると,
    (5)T=2πlgg=4π2lT2
    となる.
  3. 測定した周期Tから,上の式により,重力加速度を算出する.ここでは,式(4)から算出した重力加速度をge,式(5)から算出した重力加速度をgaと記す.(添え字のeはexactの頭文字,aはapproximateの頭文字.)
  4. 実際の重力加速度gsと算出した重力加速度gegaの相対誤差から,厳密解の正確性を考察する.
実験2

実験2では,目的2のため,以下の表のような場合について実験を行う.それぞれの場合について,5回ずつ行う.(なお,実験1で計測済みの値についてはその値を使用した.)その後の手順は以下の通りである.
実験2の実験対象 実験2の実験対象

  1. 計測した周期から,縦軸が周期,横軸が初期の振幅のグラフを作成する.
  2. 近似曲線をExcelの機能を用いて作成.
  3. 振幅と周期の関係について考察する.

誤差を少なくするための実験方法の工夫

Tの測定誤差について

測定したTの値が実際の値とずれていたら,もちろんそれだけの誤差が生じる.よって,測定したTの値と実際の値との測定誤差を小さくするために,スローカメラを利用し,できるだけ誤差が生じないようにした.

lの測定誤差について

測定したlが計算に用いた値とずれていたら,もちろんそれだけの誤差が生じる.よって,lによる測定誤差を小さくするために,金属球の直径とリングの大きさを差し引いた値でピアノ線の長さを測定した.

円錐振り子による誤差について

単振り子が円錐振り子になってしまえば,もちろんそれだけの誤差が生じる.よって,円錐振り子による誤差を小さくするために,実験者Aが振り子の振れる面からずれないように,真横から手を放した.

減衰による誤差について

単振り子が振れるうちに,減衰し,周期の測定誤差が生じる可能性がある.しかし,実験を進める過程で,5回振れる間ならほとんど無視していいほど減衰も少ないとわかった.(5回振れてすぐ後のθの最大角を見ると,初期の振幅とほぼ変わらないとわかった.)よって,減衰による誤差を小さくするために,周期を測定する際,単振り子が5回振れるときの時間を測定した.

金属球の回転による誤差

ピアノ線がねじれたままの状態で金属球を放すと,金属球が回転しながら振れて,周期の測定誤差が生じる可能性がある.よって,手を放した状態で金属球の回転が止まったところから,初期の振幅を合わせた.

カメラの固定位置による誤差について

カメラの画面から実験者Bは初期の振幅を合わせるわけだが,カメラの位置により,実際は違う振幅に合わせられる可能性がある.よって,カメラの固定位置による誤差を小さくするために,支持棒の奥行きが見えなくなるように,つまりカメラの画面に支持棒の側面の円の部分しか見えないようにカメラを固定した.

スタンドの棒の揺れによる誤差

鉄製スタンドの棒が単振り子の弦に引っ張られることにより,横に揺れ,単振り子の周期に測定に影響することが想定されたので,鉄製スタンドをC型クランプにより机に固定し,スタンドの棒が揺れないようにした.

実験結果

実験1の結果

以下,実験1の結果を示す.縦軸は実際の重力加速度gsと算出した重力加速度gegaの相対誤差,横軸は初期の振幅である.青色でプロットしてある点が、計測した値である.

実験1 !FORMULA[88][36409691][0] 相対誤差(!FORMULA[89][37980][0]=0.25m) 実験1 ge 相対誤差(l=0.25m)

実験1 !FORMULA[90][36409567][0] 相対誤差(!FORMULA[91][37980][0]=0.25m) 実験1 ga 相対誤差(l=0.25m)

実験1 !FORMULA[92][36409691][0] 相対誤差(!FORMULA[93][37980][0]=0.50m) 実験1 ge 相対誤差(l=0.50m)

実験1 !FORMULA[94][36409567][0] 相対誤差(!FORMULA[95][37980][0]=0.50m) 実験1 ga 相対誤差(l=0.50m)

実験1 !FORMULA[96][36409691][0] 相対誤差(!FORMULA[97][37980][0]=1.0m) 実験1 ge 相対誤差(l=1.0m)

実験1 !FORMULA[98][36409567][0] 相対誤差(!FORMULA[99][37980][0]=1.0m) 実験1 ga 相対誤差(l=1.0m)

実験2の結果

以下,実験2の結果を示す.縦軸が周期,横軸が初期の振幅である.青色でプロットしある点が計測した周期,緑色の線が線形近似した直線,オレンジ色の線が6次関数で多項式近似した曲線である.
実験2 振幅と周期の関係(!FORMULA[100][1470830325][0]) 実験2 振幅と周期の関係(l=1.0m

考察

周期の厳密解の妥当性について

まず,実験1について考察する.実験1から作成したグラフより,厳密解(式(3))で算出した重力加速度geは実際の重力加速度gsとの相対誤差がどの弦の長さにおいても1%以内に抑えられている.一方,近似解(式(1))で算出した重力加速度gaは実際の重力加速度gsとの相対誤差がどの長さにおいても初期の振幅が大きくなるにつれ大きくなり,θ0=45°においては8%となっている.しかし,近似解も初期の振幅が小さいときには厳密解と同程度の正確性を持っていると言える.

振幅と周期の関係について

次は,実験2について考察する.実験2から作成した図15より,周期Tは初期の振幅θ0の6次関数のように表現できることを発見した.ここで,厳密解(式(3))の式についてもう一度見てみると,
T=2πlg{n=0[(2n1)!!(2n)!!]2sin2nθ02}=2πlg{1+14sin2θ02+964sin4θ02+2252304sin6θ02+}
となり,sinθ02の6次関数となっている.では,θ0sinθ02の関係はどうなっているだろうか.
!FORMULA[114][1782773758][0]と!FORMULA[115][-1053334291][0]の関係 θsinθ02の関係
図16は,青色の点でθ0sinθ02の関係をプロットしてある.また,オレンジ色の直線はプロットした青色の点の線形近似である.すると,θ0の一次関数的にsinθ02が表現できていることがわかる.
実験1から,厳密解(式3)がある程度の妥当性を持っていることがわかっているから,つまり,周期が振幅の6次関数のように表現できることが正しいと言える.ということは,無限級数の6次の項までの影響がそれほど大きいということだろうか.実際,以下のPythonのコードを利用して計算してみると,
使用したソースコード 使用したソースコード
!FORMULA[120][1146214121][0]の6次の項までの値 σ(θ0)の6次の項までの値
となる.ただし,表の左から3つ目の「σ(θ0)におけるσ(θ0)の6次の項までの値の割合[%]」とは,表の左から2つ目の「 σ(θ0)の6次の項までの値」を「σ(θ0)」で割り,100を掛けた値である.

重力加速度の真の値とのずれについて

考察の周期の厳密解の妥当性についての節で,厳密解(式(3))で算出した重力加速度は,実際の重力加速度との相対誤差が1%に抑えられていると述べた.しかし,やはり,誤差が生じてしまっている.なぜだろうか.その可能性は主に5つ考えられる.

Tの測定誤差により誤差が生じた可能性

Tの測定誤差はかなり大きな影響を与える.実際には,ΔT=±0.015s程度の誤差が生じていた可能性が疑われる(スローカメラで測定した際,一コマ±0.015s程度だとわかった.).改善方法としては,撮影に使用したカメラを高性能なものにして,パソコンの編集機能を利用する方法である.ただ,今回予算的に行うことができなかった.今後の改善点である.

ピアノ線のねじれ及びリングと支持棒の接触部分の干渉により誤差が生じた可能性

今回実験に使用した,ピアノ線はかなり古いものであったので,少しねじれており,実験を行う際は延ばしてから使用したが,残っていた小さなねじれから誤差が生じた可能性がある.また,弦に取り付けたリングは,θ0=40°までなら,支持棒に干渉せずにおもりが振れることができる.しかし,θ0=40°以上においては,リングが支持棒に干渉してしまう場合があった.これにより,誤差が生じた可能性がある.

実際の重力加速度がずれている可能性

実験場所の重力加速度gsを参考文献[11][14]を使用して設定したが,実験場所が地上5階であったので高さを15mとしたが,これが実際とは違い,誤差が生じた可能性がある.この議論については,後で詳しく考察するが,結果的には,ほとんど無視していいほどの誤差であると結論付ける.

空気抵抗により誤差が生じた可能性

厳密解を導く過程で,速度に比例する抵抗力を一切無視して,計算を行った.しかし,実際は空気抵抗により,測定した周期に誤差が生じていた可能性がある.これについては,後で詳しく考察する.

単振り子の慣性モーメントにより誤差が生じた可能性

厳密解を導く過程で,おもりを単なる質点として,計算を行った.しかし,実際は単振り子の慣性モーメントにより,測定した周期に誤差が生じていた可能性がある.これについては,後で詳しく考察する.

国土地理院の重力値測定サービスと単振り子による重力加速度の測定について

参考文献[12]によると,「国土地理院は、全国で重力測量を実施し、国際的な基準に整合した重力値を提供しています。重力値は、正確な標高を決めるために必要であるほか、地下の活断層や資源の探査、計量機器の較正などに広く利用されています。日本の標高の基準は、測量法で平均海面と定められており、この平均海面を陸域にまで仮想的に延長した面をジオイドといいます。国土地理院では、重力測量や水準測量の結果などから、ジオイド・モデルを作成し、提供しています。ジオイド・モデルは、衛星測位で標高を求めるために使用されています。」とある.地上での測定方法について詳しく調べてみると,「重力の測定は、絶対的な重力値を正確に測定する「絶対重力測定[3]」と、相対的な重力の差を求める「相対重力測定[4]」に分けられます。国土地理院では、全国に重力の基準となる3種類の重力点(基準重力点、一等重力点、二等重力点)を設置しており、基準重力点の重力を絶対重力測定で、一等重力点及び二等重力点の重力を基準重力点からの相対重力測定でそれぞれ測定しています。」とある.また,それによって得られる絶対的な重力加速度は,「有効数字9桁」とある.つまり,相対的な重力加速度としては,「有効数字7桁」となる.

一方,今回,実験1で周期の厳密解により測定した重力加速度gegsとの相対誤差は,弦の長さl=0.25cm,初期の振幅θ0=5.0°においては,0.20%程度に抑えることができた.つまり,算出した重力加速度geは,有効数字2桁程度までなら信頼できると言える.実際には,重力値は,大気圧や潮汐力などによっても変化するため,国土地理院の測定方法ほどの精度とは言えないが,ある程度の精度を持って,単振り子の周期の測定により重力加速度を算出することができる.

また,国土地理院による重力値測定は,「測定する高さを1m高くすると、重力値の7桁目が3程度減少する。」とあるので,真の値とのずれについての可能性3つ目「実際の重力加速度がずれている可能性」は,ほぼ否定できる.

空気抵抗が与える影響について

空気抵抗が与える影響はどれほどのものなのだろうか.単振り子のおもりに物体の速度に比例する抵抗力がはたらくとする.そのとき,空気抵抗fvは,Stokesの抵抗法則(詳しくは,参考文献[7]を参照.)より,
fv=6πηav
である.ただし,おもりの半径をa,空気の粘性率をηとする.この空気抵抗による項を単振り子における運動方程式に加えると,
md2xdt2=6πηadxdtmgsinθmd2xdt2+6πηadxdt+mgsinθ=0
ここで,θが十分小さいとすると,sinθxlより,上の式は,
md2xdt2+6πηadxdt+mgxl=0d2xdt2+6πηamdxdt+glx=0
となる.これは,二階の線形常微分方程式となっているので,一般解は,
x=C1eλ1t+C2eλ2t
と書ける.ただし,C1,C2は定数であり,λ1,λ2は特性方程式λ2+6πηamλ+gl=0の解である.λ1>λ2のとき,λ1,λ2を求めると,
λ1=3πηam+(3πηam)2glλ2=3πηam(3πηam)2gl
となり,よって,一般解は,
x=C1exp(3πηamt)exp(t(3πηam)2gl)+C2exp(3πηamt)exp(t(3πηam)2gl)
となる.しかも,実験から,単振り子は加減衰ではなく,減衰振動であるので,(3πηam)2gl<0である.この条件の下,Eulerの公式eiθ=cosθ+isinθより,解は下のように書き換えられる.
x=Cexp(3πηamt)sin(tgl(3πηam)2)
ここで,Cは任意定数.よって,周期は,
(6)T=2πlg(1lg(3πηam)2)1
となる.
この式により導き出された空気抵抗を考慮した理論的な周期をTrと記す.(添え字のrはresistanceの頭文字.)また,近似解(式(1))の周期をTaと記す.この式を実験を行った状況に適用する.参考文献[7]によると,室温(約20℃),圧力1気圧において,空気の粘性率はη=1.809g/(cms)である.また,2a=413mmm=286gである.単位に気を付けて適用する.実際の周期との差について見ると,以下の表のようになる.表から,空気抵抗を考慮せず初期の振幅を微小とした近似解Taと空気抵抗を考慮し初期の振幅を微小とした解Trの差は,周期の測定誤差であるΔT=±0.015sよりも小さい.
また,初期の振幅を一般の角だとして,空気抵抗を考慮した解を導出すれば,より厳密解(式(3))に近づくと推測される.
空気抵抗を考慮した周期について 空気抵抗を考慮した周期について
つまり,空気抵抗を考慮したところで,それによって生まれる差はほとんど無視していいほどのものであるということである.すなわち,真の値とのずれについての可能性4つ目「空気抵抗により誤差が生じた可能性」は,ほぼ否定できる.

単振り子の慣性モーメントが与える影響について

単振り子の慣性モーメントが周期に与える影響はどれほどのものなのだろうか.おもりの質量をm,おもりの半径をaとすると,おもり(球)の中心周りの慣性モーメントIwは,
Iw=25ma2
となる.(補足を参照)よって,平行軸の定理(補足を参照)により,単振り子の慣性モーメントIは,
I=ml2+25ma2
である.ここで,lは弦の長さ(支点からおもりの中心までの長さ.)である.振り子が受ける力のモーメントNは,振り子がθだけ振れた位置では,
N=mglsinθ
となるので,θが十分小さいとすると,sinθθより,運動方程式は,
m(l2+25a2)θ¨=mglθ
と書ける.これを整理すると,
θ¨=gl11+25a2l2θ
となり,これは単振動の方程式であるから,角振動数をωとすれば,
ω2=gl11+25a2l2
であるから,周期Tは,
T=2πω=2πlg(1+25a2l2)
となる.

この式により導き出された単振り子の慣性モーメントを考慮した理論的な周期をTBと記す.[5]この式を実験を行った状況に適用する.2a=413mmm=286gである.単位に気を付けて適用する.実際の周期との差について見ると,以下の表のようになる.表から,単振り子の慣性モーメントを考慮せず初期の振幅を微小とした近似解Taと単振り子の慣性モーメントを考慮し初期の振幅を微小とした解TBの差は,周期の測定誤差であるΔT=±0.015sよりも小さい.
また,初期の振幅を一般の角だとして,単振り子の慣性モーメントを考慮した解を導出すれば,より厳密解(式3)に近づくと推測される.
単振り子の慣性モーメントを考慮した周期について 単振り子の慣性モーメントを考慮した周期について
つまり,単振り子の慣性モーメントを考慮したところで,それによって生まれる差はほとんど無視していいほどのものであるということである.(しかし,今回の実験とは違って,おもりの半径が弦の長さに比べて,ほぼ同じくらいの大きさとなるとさすがに無視できないほどの差が生まれるはずである.)すなわち,真の値とのずれについての可能性5つ目「単振り子の慣性モーメントにより誤差が生じた可能性」は,ほぼ否定できる.

結論

今回,2つの実験と考察から,以下,4つの事がわかった.

結論
  1. 単振り子の周期について,
    T=2πlg{n=0[(2n1)!!(2n)!!]2sin2nθ02}
    がある程度の妥当性(実際の重力加速度との相対誤差1%以内)を持って成り立つ.
  2. 単振り子の周期の測定による重力加速度の算出はある程度の信頼性(有効数字2桁)を持つ.
  3. 単振り子の周期は初期の振幅の6次関数として表現することができる.
  4. 単振り子の周期において,空気抵抗と慣性モーメントは,ほぼ影響しない.

補足

補足では,楕円積分の概要とその応用及び慣性モーメントに関する捕捉事項について述べる.

楕円積分

一般に,cxf(x,p(x))dx(c=const.,pは3次か4次の多項式,fは有理式)という形の積分のことを楕円積分(Elliptic integral)という.その中でも,第一種完全楕円積分第二種完全楕円積分と呼ばれるものを導入する.

楕円積分(ルジャンドルの標準形)
  • 第一種完全楕円積分
    K(m)=0π2dθ1m2sin2θ
  • 第二種完全楕円積分
    E(m)=0π21m2sin2θdθ
    ただし,m(0m1)はパラメーター.

ここで,K(m)について考えると,この積分は解析的に解くことはできない.そのため解析的な積分ではなく,級数展開によってアプローチする.ここで,11x=n=0(2n1)!!(2n)!!xnを利用すると,
(1m2sin2θ)12=1+12m2sin2θ+38m4sin4θ+1548m6sin6θ+=n=0(2n1)!!(2n)!!m2nsin2nθ
となる.ただし,!!は二重階乗.この級数はm<1の範囲で一様収束するので,和と積分の順序を入れ替えることが許される.よって,
0π2(1m2sin2θ)12dθ=0π2n=0(2n1)!!(2n)!!m2nsin2nθdθ=n=0(2n1)!!(2n)!!m2n0π2sin2nθdθ(7)=π2{n=0[(2n1)!!(2n)!!]2m2n}
(7)では,以下のWallisの公式(参考文献[10]を参照)を利用した.
0π2sin2nθdθ=(2n1)!!(2n)!!π2

楕円の周を求める

まず,x=asinθ,y=bcosθ(a>b)で表される楕円の周の長さを求める.曲線の微小要素の長さdsは,三平方の定理からds=dy2+dx2と表せるので,dx=acosθdθ,dy=bsinθdθに留意すると,
L=ds=dy2+dx2=40π2b2sin2θ+a2cos2θdθ=4a0π21sin2θ+b2a2sin2θdθ=4a0π21m2sin2θdθ   (0m=a2b2a1)
となり,第一種楕円積分の時と同様に級数展開して,
(1m2sin2θ)12=112m2sin2θ18m4sin4θ348m6sin6θ...=1n=1(2n1)!!(2n)!!m2n2n1sin2nθ
Wallisの公式を適用すると,
L=4a0π21m2sin2θdθ=4a0π2{1n=1(2n1)!!(2n)!!m2n2n1sin2nθdθ}=4a0π2dθn=1(2n1)!!(2n)!!m2n2n14a0π2sin2nθdθ=2aπ{1n=1[(2n1)!!(2n)!!]2m2n2n1}
となる.ただし,m=a2b2a.実際,a=b=rとすると,m=0となるので,円周の長さL=2πrに帰着する.

単振り子の周期の近似解と厳密解の差について

周期の近似解(式(1))と厳密解(式(3))の相対差(近似解と厳密解の比から1を引き,100を掛けた.)は,図4のソースコードを利用して求めると以下の表のようになる.
近似解と厳密解の相対差 近似解と厳密解の相対差
明らかに,角度が大きくになるにつれ厳密解と近似解の差が大きくなるのがわかる.

おもり(球)の中心周りの慣性モーメントについて

まず,半径a,高さhの一様な円柱について考える.この円柱について,円柱の軸に関する慣性モーメントを計算する.円柱の軸をz軸とする円柱座標系(r,θ,z)をとる.円柱の密度をρとすれば,z軸に関する慣性モーメントIzは,
Iz=ρr2dV=ρ0ar3dr02πdθ0hdz=π2ρa4h
に等しい.円柱の全質量をMとすれば,M=πρa2hであるから,
(8)Iz=12Ma2
となる.
ここで,半径aの球体の中心周りの慣性モーメントについて考える.
1つの直径を選んでそれをz軸とする.球をZ軸に垂直な薄い円盤に分割したとすると,求める慣性モーメントはこれらの円盤の慣性モーメントの総和である.
密度をρとすれば,zの位置にある厚さdzの円盤の質量はπρ(a2z2)dzであるから,式(8)により,求める慣性モーメントIwは,
Iw=aa12πρ(a2z2)dz(a2z2)=12πρaa(a2z2)2dz=815πρa5=25ma2
となる.ここでは,球体の全質量m=ρ43πa3を使った.

平行軸の定理

単振り子のおもりの慣性モーメントを計算するために以下の定理を導入する.

平行軸の定理

任意の直線lに関する物体の慣性モーメントをI,その物体の重心Gを通ってlに平行な直線lGに関する慣性モーメントをIGとする.また,物体の質量をMllGの距離をhとすると,I=IG+Mh2である.

Gを原点,lGz軸とする座標系をとる.慣性モーメントの定義から,
I=ρ((xh)2+y2)dV(9)=ρ(x2+y2)dV2hρxdV+h2ρdV
である.ただし,ρは物体の密度,dVはその微小体積要素を表わす.

(9)の右辺第1項は定義によりIGに等しい.一方,この物体の重心のx座標をxGとすると,第2項の積分は,
ρxdV=MxG
であるが,xG=0にとってあるからこれは0である.第3項はMh2に等しい.したがって,
I=IG+Mh2
となる.

付録

pdfはこちらから

楕円積分と Python を用いた単振り子の周期の検証とExcel を用いた振幅との関係についての考察




[1]: 国土地理院の重力地測定サービス(参考文献[11])を利用した.

[2]: 1080p HD/240fps

[3]: 参考文献[13]によると,国土地理院は絶対重力計の絶対重力計FG5と量子型絶対重力計AGQを保有しており,絶対重力計FG5は,「 落下槽と呼ばれる筒状の部分の内部が真空に保たれており、この落下槽の中で特殊な鏡を落下させ、落下距離をレーザーで、落下に要した時間を原子時計で計測し、重力値を求めます。」とあるのに対し,量子型絶対重力計AGQは,「原子干渉を用いた次世代型の測定装置で、真空槽内でルビジウム原子を自由落下させて重力値を求めます。実用化された量子型絶対重力計としては、我が国で初めて導入されました。FG5と異なり落体を自由落下させるための機械的な構造がないため、系統誤差の低減が期待出来るほか、取扱いが容易で可搬性に優れています。また、屋外での運用が可能であり、今後幅広い活躍が期待されます。」とある.

[4]: 参考文献[13]によると,国土地理院では,ラコスト重力計(LaCoste&Romberg gravimeter)を使用しており,ラコスト重力計は,「バネに「おもり」を吊り下げると、バネの伸び(弾性力)と重力が釣り合ったところで静止します。 重力の変化によりバネの伸び方も変わるので、その差を計ることによって重力の差を求めます。絶対重力計に比べ小型で運搬が容易ですが、測定精度は10μGalで絶対重力計より1桁落ちます。相対重力計でも、繰り返し観測を行うことにより重力の時間変化を求めることができます。」とある.

[5]: 添え字のBは物理学者ジャン=シャルル・ド・ボルダ(Jean-Charles, chevalier de Borda,1733-1799)のBordaの頭文字.一般に,剛体球のおもりによる振り子をBordaの振り子という.

参考文献

[1]
國友正和,今和泉卓也,井上邦雄,黒田楯彦,河本敏郎,小林雅之,萩野浩一,田原輝夫,深 津晋,橋本道雄,牧島一夫,増渕哲夫,漆原次郎, 文部科学省検定済教科書 物理, 数研出版, 2022
[2]
兵頭俊夫, 考える力学 第 2 版, 学術図書出版会, 2021
[3]
今井功,高見穎郎,高木隆司,吉澤徴,下村裕, セミナーライブラリ物理学 2 演習 力学 (新訂版), サイエンス社, 2021
[4]
塩野正明,松村敬治, 小学校理科における「振り子の運動」の実験指導と誤差の扱いについて, 西南学院大学人間科学論集, 2011, 107-121
[7]
巽友正, 流体力学, 培風館, 2022
[8]
杉山弘,遠藤剛,新井隆景, 流体力学 第2版, 森北出版株式会社, 2021
投稿日:20231220
更新日:2024110
OptHub AI Competition

この記事を高評価した人

高評価したユーザはいません

この記事に送られたバッジ

バッジはありません。
バッチを贈って投稿者を応援しよう

バッチを贈ると投稿者に現金やAmazonのギフトカードが還元されます。

投稿者

PRESENT
PRESENT
10
3528
‘06 高3 | JPhO'22,’23 / 科甲J’19,科甲’22,’23 / JKA小中’16,’17 記事の一覧は下のリンクから

コメント

他の人のコメント

コメントはありません。
読み込み中...
読み込み中
  1. はじめに
  2. 概要
  3. 実験の目的
  4. 実験手法
  5. 実験に使用したもの
  6. 実験原理
  7. Pythonを利用した級数の計算方法
  8. 実験に使用した値
  9. 実験手順
  10. 誤差を少なくするための実験方法の工夫
  11. 実験結果
  12. 実験1の結果
  13. 実験2の結果
  14. 考察
  15. 周期の厳密解の妥当性について
  16. 振幅と周期の関係について
  17. 重力加速度の真の値とのずれについて
  18. 国土地理院の重力値測定サービスと単振り子による重力加速度の測定について
  19. 空気抵抗が与える影響について
  20. 単振り子の慣性モーメントが与える影響について
  21. 結論
  22. 補足
  23. 楕円積分
  24. 楕円の周を求める
  25. 単振り子の周期の近似解と厳密解の差について
  26. おもり(球)の中心周りの慣性モーメントについて
  27. 平行軸の定理
  28. 付録
  29. pdfはこちらから
  30. 参考文献