この記事は個人で行った実験結果と考察をもとに書かれています.エンタメとしてご覧ください.
はじめに
この実験は,「第19回全国物理コンテスト物理チャレンジ2023」の「第1チャレンジ2023実験レポート課題」の実験課題に基づき行われています.
振り子の周期を、振れ角を変えて調べてみよう
単振り子の周期を測定する実験では振れ角を小さくして行います が、振れ角を大きくしていくと周期はどうなるでしょうか。振り子の 振れ角と周期の関係について、実験を行って調べてみましょう。 周期を正確に測る工夫、減衰の小さい振り子の作製、振れ角と周期 の関係についての考察などを期待しています。
物理チャレンジについてはこちらから:
全国物理コンテスト 物理チャレンジ!
概要
参考文献[1]によるとによると,単振り子の周期は「振れが小さいとき,周期は振幅によらず,糸の長さと重力加速度の大きさだけで決まる.これを振り子の当時性という.」とある.つまり,
単振り子の振幅が微小のとき,周期は,
のように表されるが,振幅が一般の角のとき,周期は楕円積分を用いると,
と表される.ただし,は二重階乗.これを利用して以下の実験を行った.
実験1
上記の周期の厳密解にどれほどの正確性があるのかを,振幅と弦の長さを変えて,カメラで単振り子が振れる様子を撮影し,周期を測定することによって調べた.
実験2
振幅を変えて,周期の変化について,カメラで単振り子が振れる様子を撮影し,周期を測定する事によって調べた.
実験2から,周期と振幅の関係について,周期が初期の振幅の六次関数として表されることを発見し,その理由を実験1と理論的考察から推測した.
また,空気抵抗と慣性モーメントについての考察も行い,単振り子の測定した周期と理論的な周期のずれについて徹底的に検証を行った.
実験の目的
今回の実験の目的を以下のように設定した.
- 単振り子の周期の厳密解にどれほどの正確性があるのか調べる.
- 周期と振幅の関係を調べる.
- 単振り子の周期において,空気抵抗と慣性モーメントが与える影響を調べる.
実験手法
実験に使用したもの
実験道具(実験自体に使用したもの)
単振り子(島津理化機械株式会社製 TN‐3 形,ケニス株式会社製 No.110-0210 単振り
子 SP ).金属球:スチールボール1 5/8(直径:2a = 413mm,質量:m = 286g).弦:
ピアノ線(太さ:0.35mm).
単振り子(金属球を使用した)
単振り子(弦を使用した)
鉄製スタンド
分度器
C型クランプ(鉄製スタンドを机に固定するのに使用)
3.5mメジャー
段ボールとコピー用紙A4で作成した角度確認用分度器
角度確認用分度器
ペン
卓上ボール盤(段ボールに穴をあけるのに使用)
測定道具(実験の測定及び解析に使用したもの)
- iPhoneXR(撮影用)
- 三脚
- パソコン(ストップウォッチ用と解析用)
実験原理
以下,単振り子の弦の長さを(この節では,単振り子のおもりを質点として扱うが,実験においては,おもりの半径を加えた長さを弦の長さとする.),振幅を(は振り子が振れ始めてからの経過時間の関数),周期を,小球の質量を,重力加速度を,最下点の位置を,点から円弧に沿った変位をと表す.
振幅が非常に小さいとき,単振り子の周期は以下のように表される.
なお,周期がこのように表される理由は以下のようである.
円弧の接線方向の運動方程式を解く.
であるので,
が十分に小さいとき,として考えられるので,
この二階微分方程式を解くと,
ただし,は任意定数.そしてこの式は単振動を表すので,角速度は,となり,周期は,
となる.
では,を任意の大きさの角として,厳密解を求めるとどうなるのか.以下,簡略化のためにの時間一階微分を,時間二階微分をと記す.式の両辺にをかけて,右辺を左辺に移行する.
よって,は定数となるので,定数と置く.(はEnergyの頭文字である.実際,位置エネルギーの基準点を最下点に置くと,振り子の持つ力学的エネルギー(位置エネルギーと運動エネルギーの和)は,となり,上の式と対応している.)
のとき,だとすると,
よって,
いま,からまでの運動を考えると,
さらに,からまでの運動にかかる時間はであるので,
ここで,を満たすに変数変換すると,
となり,第一種完全楕円積分(楕円積分の節を参照)の形となる.よって楕円積分の節で求めた結果を用いて計算すると,
となる。ただし,は二重階乗.
なお、簡略化のため、関数を以下のように定める.
Pythonを利用した級数の計算方法
前節からの厳密解が分かったが,無限級数の形で表されているため,実際に計算しつくすことは不可能である.そのため,できるだけ理論値の正確性を向上させるために,プログラミング言語Pythonを利用して,の数値計算を行う.以下のコードを使用する.ただし定数aは角度であり,常々設定する.
使用したソースコード
これによって得られた,の値は以下の表の通りである
得られた値
一個の角度あたり、300~600秒くらいの計算時間がかかった.
実験に使用した値
実験手順
大まかな実験手順は以下の通りである.実験ごとの細かい解析手順についてはその下で示す.
- 弦の長さ(金属球の中心から弦についたリングと支持棒の接点までの長さ)を調べる対象に合わせて測り,固定.
- 実験装置を図のようにセットする.
- 鉄製スタンド横に設置したパソコンのストップウォッチをスタートさせておく.
- 実験者Aが振り子の位置を調節し,実験者BがiPhoneのカメラから振り子の位置を確認し,適切な位置を指示する.
- 実験者BがiPhoneのスローカメラ[2]をスタートさせ,実験者Aが振り子をそっと放し,振り子が5回以上振れたときに実験者Bがカメラを止める.
- 撮影した動画を確認し,振り子が静止した状態で振り子の中心を通るように鉛直方向に引いた線(線とする.)を利用して,初めの時間(振り子が線を初めて通過したときの時間.)と終わりの時間(振り子が5回振れたときに線を通過した時間.)を計測することによって,その差から5回振れる間にかかった時間()を測定する.
実験装置
実験1
実験1では,目的1のため,以下の表のような場合について実験を行う.それぞれの場合について,5回ずつ行う.その後の手順は以下の通りである.
実験1の実験対象
- 式について,として,変形すると,
となる. - また,式について,変形すると,
となる. - 測定した周期から,上の式により,重力加速度を算出する.ここでは,式から算出した重力加速度を,式から算出した重力加速度をと記す.(添え字のはexactの頭文字,はapproximateの頭文字.)
- 実際の重力加速度と算出した重力加速度,の相対誤差から,厳密解の正確性を考察する.
実験2
実験2では,目的2のため,以下の表のような場合について実験を行う.それぞれの場合について,5回ずつ行う.(なお,実験1で計測済みの値についてはその値を使用した.)その後の手順は以下の通りである.
実験2の実験対象
- 計測した周期から,縦軸が周期,横軸が初期の振幅のグラフを作成する.
- 近似曲線をExcelの機能を用いて作成.
- 振幅と周期の関係について考察する.
誤差を少なくするための実験方法の工夫
の測定誤差について
測定したの値が実際の値とずれていたら,もちろんそれだけの誤差が生じる.よって,測定したの値と実際の値との測定誤差を小さくするために,スローカメラを利用し,できるだけ誤差が生じないようにした.
の測定誤差について
測定したが計算に用いた値とずれていたら,もちろんそれだけの誤差が生じる.よって,による測定誤差を小さくするために,金属球の直径とリングの大きさを差し引いた値でピアノ線の長さを測定した.
円錐振り子による誤差について
単振り子が円錐振り子になってしまえば,もちろんそれだけの誤差が生じる.よって,円錐振り子による誤差を小さくするために,実験者Aが振り子の振れる面からずれないように,真横から手を放した.
減衰による誤差について
単振り子が振れるうちに,減衰し,周期の測定誤差が生じる可能性がある.しかし,実験を進める過程で,5回振れる間ならほとんど無視していいほど減衰も少ないとわかった.(5回振れてすぐ後のの最大角を見ると,初期の振幅とほぼ変わらないとわかった.)よって,減衰による誤差を小さくするために,周期を測定する際,単振り子が5回振れるときの時間を測定した.
金属球の回転による誤差
ピアノ線がねじれたままの状態で金属球を放すと,金属球が回転しながら振れて,周期の測定誤差が生じる可能性がある.よって,手を放した状態で金属球の回転が止まったところから,初期の振幅を合わせた.
カメラの固定位置による誤差について
カメラの画面から実験者Bは初期の振幅を合わせるわけだが,カメラの位置により,実際は違う振幅に合わせられる可能性がある.よって,カメラの固定位置による誤差を小さくするために,支持棒の奥行きが見えなくなるように,つまりカメラの画面に支持棒の側面の円の部分しか見えないようにカメラを固定した.
スタンドの棒の揺れによる誤差
鉄製スタンドの棒が単振り子の弦に引っ張られることにより,横に揺れ,単振り子の周期に測定に影響することが想定されたので,鉄製スタンドをC型クランプにより机に固定し,スタンドの棒が揺れないようにした.
実験結果
実験1の結果
以下,実験1の結果を示す.縦軸は実際の重力加速度と算出した重力加速度,の相対誤差,横軸は初期の振幅である.青色でプロットしてある点が、計測した値である.
実験1 相対誤差(=0.25m)
実験1 相対誤差(=0.25m)
実験1 相対誤差(=0.50m)
実験1 相対誤差(=0.50m)
実験1 相対誤差(=1.0m)
実験1 相対誤差(=1.0m)
実験2の結果
以下,実験2の結果を示す.縦軸が周期,横軸が初期の振幅である.青色でプロットしある点が計測した周期,緑色の線が線形近似した直線,オレンジ色の線が6次関数で多項式近似した曲線である.
実験2 振幅と周期の関係()
考察
周期の厳密解の妥当性について
まず,実験1について考察する.実験1から作成したグラフより,厳密解(式)で算出した重力加速度は実際の重力加速度との相対誤差がどの弦の長さにおいても以内に抑えられている.一方,近似解(式)で算出した重力加速度は実際の重力加速度との相対誤差がどの長さにおいても初期の振幅が大きくなるにつれ大きくなり,においてはとなっている.しかし,近似解も初期の振幅が小さいときには厳密解と同程度の正確性を持っていると言える.
振幅と周期の関係について
次は,実験2について考察する.実験2から作成した図15より,周期は初期の振幅の6次関数のように表現できることを発見した.ここで,厳密解(式)の式についてもう一度見てみると,
となり,の6次関数となっている.では,との関係はどうなっているだろうか.
との関係
図16は,青色の点でとの関係をプロットしてある.また,オレンジ色の直線はプロットした青色の点の線形近似である.すると,の一次関数的にが表現できていることがわかる.
実験1から,厳密解(式)がある程度の妥当性を持っていることがわかっているから,つまり,周期が振幅の6次関数のように表現できることが正しいと言える.ということは,無限級数の6次の項までの影響がそれほど大きいということだろうか.実際,以下のPythonのコードを利用して計算してみると,
使用したソースコード
の6次の項までの値
となる.ただし,表の左から3つ目の「におけるの6次の項までの値の割合」とは,表の左から2つ目の「 の6次の項までの値」を「」で割り,100を掛けた値である.
重力加速度の真の値とのずれについて
考察の周期の厳密解の妥当性についての節で,厳密解(式)で算出した重力加速度は,実際の重力加速度との相対誤差が1に抑えられていると述べた.しかし,やはり,誤差が生じてしまっている.なぜだろうか.その可能性は主に5つ考えられる.
の測定誤差により誤差が生じた可能性
の測定誤差はかなり大きな影響を与える.実際には,程度の誤差が生じていた可能性が疑われる(スローカメラで測定した際,一コマ程度だとわかった.).改善方法としては,撮影に使用したカメラを高性能なものにして,パソコンの編集機能を利用する方法である.ただ,今回予算的に行うことができなかった.今後の改善点である.
ピアノ線のねじれ及びリングと支持棒の接触部分の干渉により誤差が生じた可能性
今回実験に使用した,ピアノ線はかなり古いものであったので,少しねじれており,実験を行う際は延ばしてから使用したが,残っていた小さなねじれから誤差が生じた可能性がある.また,弦に取り付けたリングは,までなら,支持棒に干渉せずにおもりが振れることができる.しかし,以上においては,リングが支持棒に干渉してしまう場合があった.これにより,誤差が生じた可能性がある.
実際の重力加速度がずれている可能性
実験場所の重力加速度を参考文献[11],[14]を使用して設定したが,実験場所が地上5階であったので高さをとしたが,これが実際とは違い,誤差が生じた可能性がある.この議論については,後で詳しく考察するが,結果的には,ほとんど無視していいほどの誤差であると結論付ける.
空気抵抗により誤差が生じた可能性
厳密解を導く過程で,速度に比例する抵抗力を一切無視して,計算を行った.しかし,実際は空気抵抗により,測定した周期に誤差が生じていた可能性がある.これについては,後で詳しく考察する.
単振り子の慣性モーメントにより誤差が生じた可能性
厳密解を導く過程で,おもりを単なる質点として,計算を行った.しかし,実際は単振り子の慣性モーメントにより,測定した周期に誤差が生じていた可能性がある.これについては,後で詳しく考察する.
国土地理院の重力値測定サービスと単振り子による重力加速度の測定について
参考文献[12]によると,「国土地理院は、全国で重力測量を実施し、国際的な基準に整合した重力値を提供しています。重力値は、正確な標高を決めるために必要であるほか、地下の活断層や資源の探査、計量機器の較正などに広く利用されています。日本の標高の基準は、測量法で平均海面と定められており、この平均海面を陸域にまで仮想的に延長した面をジオイドといいます。国土地理院では、重力測量や水準測量の結果などから、ジオイド・モデルを作成し、提供しています。ジオイド・モデルは、衛星測位で標高を求めるために使用されています。」とある.地上での測定方法について詳しく調べてみると,「重力の測定は、絶対的な重力値を正確に測定する「絶対重力測定[3]」と、相対的な重力の差を求める「相対重力測定[4]」に分けられます。国土地理院では、全国に重力の基準となる3種類の重力点(基準重力点、一等重力点、二等重力点)を設置しており、基準重力点の重力を絶対重力測定で、一等重力点及び二等重力点の重力を基準重力点からの相対重力測定でそれぞれ測定しています。」とある.また,それによって得られる絶対的な重力加速度は,「有効数字9桁」とある.つまり,相対的な重力加速度としては,「有効数字7桁」となる.
一方,今回,実験1で周期の厳密解により測定した重力加速度のとの相対誤差は,弦の長さ,初期の振幅においては,程度に抑えることができた.つまり,算出した重力加速度は,有効数字2桁程度までなら信頼できると言える.実際には,重力値は,大気圧や潮汐力などによっても変化するため,国土地理院の測定方法ほどの精度とは言えないが,ある程度の精度を持って,単振り子の周期の測定により重力加速度を算出することができる.
また,国土地理院による重力値測定は,「測定する高さを高くすると、重力値の7桁目が3程度減少する。」とあるので,真の値とのずれについての可能性3つ目「実際の重力加速度がずれている可能性」は,ほぼ否定できる.
空気抵抗が与える影響について
空気抵抗が与える影響はどれほどのものなのだろうか.単振り子のおもりに物体の速度に比例する抵抗力がはたらくとする.そのとき,空気抵抗は,Stokesの抵抗法則(詳しくは,参考文献[7]を参照.)より,
である.ただし,おもりの半径を,空気の粘性率をとする.この空気抵抗による項を単振り子における運動方程式に加えると,
ここで,が十分小さいとすると,より,上の式は,
となる.これは,二階の線形常微分方程式となっているので,一般解は,
と書ける.ただし,は定数であり,は特性方程式の解である.のとき,を求めると,
となり,よって,一般解は,
となる.しかも,実験から,単振り子は加減衰ではなく,減衰振動であるので,である.この条件の下,Eulerの公式より,解は下のように書き換えられる.
ここで,は任意定数.よって,周期は,
となる.
この式により導き出された空気抵抗を考慮した理論的な周期をと記す.(添え字のはresistanceの頭文字.)また,近似解(式)の周期をと記す.この式を実験を行った状況に適用する.参考文献[7]によると,室温(約20℃),圧力1気圧において,空気の粘性率はである.また,,である.単位に気を付けて適用する.実際の周期との差について見ると,以下の表のようになる.表から,空気抵抗を考慮せず初期の振幅を微小とした近似解と空気抵抗を考慮し初期の振幅を微小とした解の差は,周期の測定誤差であるよりも小さい.
また,初期の振幅を一般の角だとして,空気抵抗を考慮した解を導出すれば,より厳密解(式)に近づくと推測される.
空気抵抗を考慮した周期について
つまり,空気抵抗を考慮したところで,それによって生まれる差はほとんど無視していいほどのものであるということである.すなわち,真の値とのずれについての可能性4つ目「空気抵抗により誤差が生じた可能性」は,ほぼ否定できる.
単振り子の慣性モーメントが与える影響について
単振り子の慣性モーメントが周期に与える影響はどれほどのものなのだろうか.おもりの質量を,おもりの半径をとすると,おもり(球)の中心周りの慣性モーメントは,
となる.(補足を参照)よって,平行軸の定理(補足を参照)により,単振り子の慣性モーメントは,
である.ここで,は弦の長さ(支点からおもりの中心までの長さ.)である.振り子が受ける力のモーメントは,振り子がだけ振れた位置では,
となるので,が十分小さいとすると,より,運動方程式は,
と書ける.これを整理すると,
となり,これは単振動の方程式であるから,角振動数をとすれば,
であるから,周期は,
となる.
この式により導き出された単振り子の慣性モーメントを考慮した理論的な周期をと記す.[5]この式を実験を行った状況に適用する.,である.単位に気を付けて適用する.実際の周期との差について見ると,以下の表のようになる.表から,単振り子の慣性モーメントを考慮せず初期の振幅を微小とした近似解と単振り子の慣性モーメントを考慮し初期の振幅を微小とした解の差は,周期の測定誤差であるよりも小さい.
また,初期の振幅を一般の角だとして,単振り子の慣性モーメントを考慮した解を導出すれば,より厳密解(式)に近づくと推測される.
単振り子の慣性モーメントを考慮した周期について
つまり,単振り子の慣性モーメントを考慮したところで,それによって生まれる差はほとんど無視していいほどのものであるということである.(しかし,今回の実験とは違って,おもりの半径が弦の長さに比べて,ほぼ同じくらいの大きさとなるとさすがに無視できないほどの差が生まれるはずである.)すなわち,真の値とのずれについての可能性5つ目「単振り子の慣性モーメントにより誤差が生じた可能性」は,ほぼ否定できる.
結論
今回,2つの実験と考察から,以下,4つの事がわかった.
結論
- 単振り子の周期について,
がある程度の妥当性(実際の重力加速度との相対誤差1以内)を持って成り立つ. - 単振り子の周期の測定による重力加速度の算出はある程度の信頼性(有効数字2桁)を持つ.
- 単振り子の周期は初期の振幅の6次関数として表現することができる.
- 単振り子の周期において,空気抵抗と慣性モーメントは,ほぼ影響しない.
補足
補足では,楕円積分の概要とその応用及び慣性モーメントに関する捕捉事項について述べる.
楕円積分
一般に,という形の積分のことを楕円積分(Elliptic integral)という.その中でも,第一種完全楕円積分と第二種完全楕円積分と呼ばれるものを導入する.
楕円積分(ルジャンドルの標準形)
- 第一種完全楕円積分
- 第二種完全楕円積分
ただし,はパラメーター.
ここで,について考えると,この積分は解析的に解くことはできない.そのため解析的な積分ではなく,級数展開によってアプローチする.ここで,を利用すると,
となる.ただし,は二重階乗.この級数はの範囲で一様収束するので,和と積分の順序を入れ替えることが許される.よって,
では,以下のWallisの公式(参考文献[10]を参照)を利用した.
楕円の周を求める
まず,で表される楕円の周の長さを求める.曲線の微小要素の長さは,三平方の定理からと表せるので,に留意すると,
となり,第一種楕円積分の時と同様に級数展開して,
Wallisの公式を適用すると,
となる.ただし,.実際,とすると,となるので,円周の長さに帰着する.
単振り子の周期の近似解と厳密解の差について
周期の近似解(式)と厳密解(式)の相対差(近似解と厳密解の比から1を引き,100を掛けた.)は,図4のソースコードを利用して求めると以下の表のようになる.
近似解と厳密解の相対差
明らかに,角度が大きくになるにつれ厳密解と近似解の差が大きくなるのがわかる.
おもり(球)の中心周りの慣性モーメントについて
まず,半径,高さの一様な円柱について考える.この円柱について,円柱の軸に関する慣性モーメントを計算する.円柱の軸を軸とする円柱座標系をとる.円柱の密度をとすれば,軸に関する慣性モーメントは,
に等しい.円柱の全質量をとすれば,であるから,
となる.
ここで,半径の球体の中心周りの慣性モーメントについて考える.
1つの直径を選んでそれを軸とする.球を軸に垂直な薄い円盤に分割したとすると,求める慣性モーメントはこれらの円盤の慣性モーメントの総和である.
密度をとすれば,の位置にある厚さの円盤の質量はであるから,式により,求める慣性モーメントは,
となる.ここでは,球体の全質量を使った.
平行軸の定理
単振り子のおもりの慣性モーメントを計算するために以下の定理を導入する.
平行軸の定理
任意の直線に関する物体の慣性モーメントを,その物体の重心を通ってに平行な直線に関する慣性モーメントをとする.また,物体の質量を,との距離をとすると,である.
を原点,を軸とする座標系をとる.慣性モーメントの定義から,
である.ただし,は物体の密度,はその微小体積要素を表わす.
式の右辺第1項は定義によりに等しい.一方,この物体の重心の座標をとすると,第2項の積分は,
であるが,にとってあるからこれはである.第3項はに等しい.したがって,
となる.
付録
pdfはこちらから
楕円積分と Python を用いた単振り子の周期の検証とExcel を用いた振幅との関係についての考察
[3]: 参考文献[13]によると,国土地理院は絶対重力計の絶対重力計FG5と量子型絶対重力計AGQを保有しており,絶対重力計FG5は,「 落下槽と呼ばれる筒状の部分の内部が真空に保たれており、この落下槽の中で特殊な鏡を落下させ、落下距離をレーザーで、落下に要した時間を原子時計で計測し、重力値を求めます。」とあるのに対し,量子型絶対重力計AGQは,「原子干渉を用いた次世代型の測定装置で、真空槽内でルビジウム原子を自由落下させて重力値を求めます。実用化された量子型絶対重力計としては、我が国で初めて導入されました。FG5と異なり落体を自由落下させるための機械的な構造がないため、系統誤差の低減が期待出来るほか、取扱いが容易で可搬性に優れています。また、屋外での運用が可能であり、今後幅広い活躍が期待されます。」とある.↩ [4]: 参考文献[13]によると,国土地理院では,ラコスト重力計(LaCoste&Romberg gravimeter)を使用しており,ラコスト重力計は,「バネに「おもり」を吊り下げると、バネの伸び(弾性力)と重力が釣り合ったところで静止します。 重力の変化によりバネの伸び方も変わるので、その差を計ることによって重力の差を求めます。絶対重力計に比べ小型で運搬が容易ですが、測定精度は10μGalで絶対重力計より1桁落ちます。相対重力計でも、繰り返し観測を行うことにより重力の時間変化を求めることができます。」とある.↩ [5]: 添え字のは物理学者ジャン=シャルル・ド・ボルダ(Jean-Charles, chevalier de Borda,1733-1799)のBordaの頭文字.一般に,剛体球のおもりによる振り子をBordaの振り子という.↩