- 追加された行はこの色です。
- 削除された行はこの色です。
- 力学基礎 へ行く。
#author("2021-06-08T18:13:36+09:00","","")
#author("2021-06-22T18:44:30+09:00","","")
* 波についての章が教科書にないが [#c78d359a]
調和振動子は２階の線型方程式で、減衰振動はその自然な拡張であることを学んだ。調和振動子が複数個ありそれらが結合している場合の運動の問題が更に自然な力学の問題の拡張となっている。複数個でなく無限個の拡張となるとそれは物理の問題の重要な分野である波動現象となる。この題材の歴史的な意義も含めて以下の動画をダウンロードし見てみよう。この分野の古典的名著は「バークレー物理学コース 波動」だと個人的には思います。
https://natsci.kyokyo-u.ac.jp/~takasima/MU/Waves.avi
手っ取り早く波の概念を掴むにはシミュレーションも便利な手段です。今どきは以下のようなページを利用するのも良いかもしれないです。
https://qiita.com/kuroitu/items/5d6054ddf13e68720c6e
* 減衰振動の時の微妙な位相差 [#g7568f71]
通常単振動なら初速を与えない限り、引っ張って静かに離すだけではコサイン関数内の位相差は生まれない。ところが減衰振動の場合はどうもそうでないらしい。減衰力から自動的に位相差を生み出すことになるようだ。素粒子でも粒子が転換していくときに減衰振動となるが、いつも位相差のパラメータを入れているのが不思議だったが、たしかに位相差は減衰力を通じたパラメータになっているらしい。
* 力学、電磁気学でアニメーションで可視化 [#o2c8181a]
qiita, javascript, canvasで検索すると以下のページが出てきた。最後のリンクで地球の運動のアニメーションがある。南半球から見てるようになってる。
https://qiita.com/kyrieleison/items/a3ebf7c55295c3e7d8f0
またまたパソコンをひらって来てapache2を起動して、上の最後のリンクのページの演習をやってみた。
#ref(developer.tar.gz)
sudo su -
cd /var/www/html
tar zxvf developer.tar.gz
exit
このあと自分のところでブラウザのlocalhostからそれぞれのところを見る。結局、java, flashはもう使われなくなってる。animeやマウスを使うものはjavascript。
windowsでもメモリがでかくなってるので、apache2を起動しっぱなしでOKでないか。以下のページでやってみた。広告をクリックしないように。
https://webkaru.net/php/windows-apache-install/
* 火星軌道がどんなものかを見てみよう [#i252ba4c]
#ref(ellipsMars.py)
中のコードをみると、角度のデータを３６０個作成して、半径のデータも同じだけ作り、プロットさせているみたい。地球の軌道も同様に計算して、プロットする。
kepler's 3 law のビデオでは、火星の１年を使って、地球軌道を推定している。太陽の方位とその方位に対する火星の方位のなす角で地球の位置を示す交点を決めているように見える。そのあと火星軌道は合を使って火星方位を決めた後、火星の１年後の地球の位置で太陽の方位をきめ、火星の方位からその最初の火星の方位との交点を決めるというやり方となっている。次にまた衝の位置に移る。同じことを繰り返して、火星軌道を求めているように見える。
火星軌道に火星をアニメーションで示す時は角度を増やしていけば良いが、長半径の平方根をかけて半径の２乗で割ってやると良いようだ。金星も同時に示して、「5金星は8太陽」のドレスデン・コデックスが示せればよいかも。金星の衝の周期が8/5太陽年となるという意味。これで太陽が８回まわる間に金星は１３回まわることになる。このあたりはコスミックフロントのマヤ天文学で示されていた。
#ref(planet.java)
#ref(planet.html)
javac -encoding UTF-8 planet.java
appletviewer planet.html
youtubeで火星軌道を描くのがあったので、その表を探してみた。Mars opposition protractorで検索したところ宿題を集めたところがあり、そこにチコブラーエの観測らしき表があった。初めの行で地球の場所を決め火星の方向に線を引く、次の行でも同じことをする。角度はすべて春分点からの角度になっている。分度器の中心を太陽に合わせるか、地球に合わせる。二つの行は火星の1年間隔なので火星は同じ場所にいる。表は1度刻みなので荒っぽい。チコの観測はその60分の1の精度だったらしい。別のところにもこんなのがあった。
youtubeで火星軌道を描くのがあったので、その表を探してみた。Mars opposition protractorで検索したところ宿題を集めたところがあり、そこにチコブラーエの観測らしき表があった。
https://web2.utc.edu/~Johnny-Cantrell/cp/labs/cplb1445.pdf
初めの行で地球の場所を決め火星の方向に線を引く、次の行でも同じことをする。角度はすべて春分点からの角度になっている。分度器の中心を太陽に合わせるか、地球に合わせる。二つの行は火星の1年間隔なので火星は同じ場所にいる。表は1度刻みなので荒っぽい。チコの観測はその60分の1の精度だったらしい。別のところにもこんなのがあった。
https://faculty.humanities.uci.edu/bjbecker/ExploringtheCosmos/planetdistances.html
* 位置エネルギーとエネルギーの保存 [#iefc2de2]
自分にわかりやすいようにスライドを作ってみました。y軸は地表から鉛直上向き。宇宙の彼方に物を持ち上げるときも同じ。
#ref(potential.pdf)
* 第4回の課題 [#tfc3e986]
教科書を読んでいくと、落下運動では速度に比例する減衰力、速度の2乗に比例する減衰力の話がある。上向きを正としている関係で式が成立している。速度の2乗の式は下向きに運動しているので上向きの抵抗力となっているあたりに注意する。上向きに高速運動する場合は加速してしまい間違いとなる。
今回は、タコマ橋の崩落について適当にyoutubeで検索してもらい、カルマン渦や強制振動の問題を考察する。演習としてはanaconda3を使って強制振動のシミュレーションをやってみてください。anaconda3のプロンプトから
conda install matplotlib
conda install scipy
idle
pipでなく、condaコマンドでパッケージを入れたほうが良いらしい。そのあと以下のページからプログラムをコピペして実行してみてください。
https://watlab-blog.com/2019/06/14/mck-force-vibration/
anaconda3の使い方については京都大学の喜多さん、pythonについての演習のpdfを見てください。コードで分かりにくい部分は変数の型が明示されない点。solはnparrayの方程式の解のsolに対して、位置と速度のリストの位置部分をmatplotのグラフ(t,ｙ)のy値に設定するところのようである。sol[:,0]は行に対してスライスで全部、列に対して第0列。
昔、試験問題を作った時、ばねと重りの問題で下向きを正としたときもあった。その時はつり合いから始まるので重力は式に出てこない。その時も速度に比例する抵抗力なら負の符号でよく、復元力も負の符号でよい。式としては横向きと同じになる。強制振動は天井をゆする感じになる。
ここでodeintについて調べてみたところ、以下のページがあり、自由落下の例について実行してみると、odeintの出力がnumpyの行列表現であることがわかる。
https://qiita.com/Sunset_Yuhi/items/d938718ad277eeab746c
* 第3回の課題はファイルにした [#w9855c83]
一応関連するビデオをアップしてみた。最後のバイロンの詩のところは飛ばしてよい。古語など使ってわかりません。
https://natsci.kyokyo-u.ac.jp/~takasima/MU/AppleMoon.avi
* 授業開始までの課題2 [#y398bd41]
慣性の法則のビデオを見るのを2番目の課題とします。prevail(overcome), ponder(think deeply),、projectile(投射体）、ploy（策略）、recant（取り消す）など辞書を引きながら見てみるとよいでしょう。
https://natsci.kyokyo-u.ac.jp/~takasima/MU/Inertia.avi
* 授業開始までの課題１ [#d87c2577]
30年ほど前に作られたCaltecのGoodstein教授の力学基礎の講義
https://natsci.kyokyo-u.ac.jp/~takasima/MU/Introduction.avi
よくわからないところがでてくるはず。海上保安隊とビールを飲みながら漂う船のエピソードは連絡を受けた地点から海流を考えて、救助に向かう話。ベクトルの足し算の話題の具体例だったりする。30回分の力学の授業のビデオのイントロであるため。
* フーコーの振り子が右にずれていく [#h0d1e729]
最初は左にずれるがそのうち右に行く。コンパイルしてappletviewerで実行。
#ref(turnTable.java)
#ref(turnTable.html)
* 剛体の勉強法 [#e970ddb7]
剛体の慣性モーメントを勉強するのに良いのは坂道ころりん問題。角運動量とも関係するが、モーメントを勉強するのに良いのが壁たてパッタッン問題。摩擦が関係するので、坂道を滑り始めるタンジェントが摩擦係数となるのに関係付けるのに上の角度で考えるのが良い。まずは上側が摩擦なしとして、タンジェントが摩擦係数の２倍で滑り始める。上が同じ摩擦ならそれより少し角度が付く。(1-mu^2)で割ったぶん。
* 減衰振動を方程式で表現する [#sb6b6c87]
振動現象は複素数で表現する。
plot [-6:6] [-1.5:1.5] x-x**3/3!+x**5/5!-x**7/7!+x**9/9!
*二重振り子再び [#cd9a8a66]
クリストフさんのページの説明がmathmaticaなんで、同じような計算は下のページにある。微分の関係を分離計算する必要があるので多少面倒になる。
https://www.aihara.co.jp/~taiji/pendula-equations/present-node2.html
* javaの演習 [#tf4cca37]
http://www.ita.uni-heidelberg.de/~chfeder/applets.shtml?lang=en
* 中性子星の合体 [#e0b2d1b8]
最近の最大の話題は中性子星の合体
https://www.ligo.org/detections/GW170817.php
10月にはマルチメッセンジャー論文というのがアストロ・フィジカル・ジャーナルから刊行された。
https://www.ligo.caltech.edu/page/detection-companion-papers
Dan GreenのParticle Physicsのp437に10kmの干渉計で50Aのずれが観測される例が示されている。この場合は10e-13のスケールでの変化となる。10e-21でも観測できるというのは本当なのか？->色々見てもこの差は埋まらない。無理やり振幅の変動を見ているらしい。位相が8桁分の1ずれるのは分かるものらしい。
* 宇宙の臨界密度の計算 [#sd90b4c2]
仮想的な質量中心からの脱出速度をハッブル定数に取り、普通の力学での脱出速度から臨界質量が計算できる。普通の物質が５％ほどだが残りの９５％がダークマターとダークエネルギーとするとちょうど観測量が臨界質量にほぼ一致するらしい。ウィキペディアでは平坦性問題にそう書いてある。
* よくわかる初等力学のjava [#v3e46b8a]
最近のブラウザはjavaに厳しい。ソースをコピーしてjdkを使ってコンパイルする。macだとappletviewerがあるのでnkf -w xxx.java > xxxU.javaとしてxxxU.javaのクラス名をxxxUにする。そのあとコンパイルして、xxxU.htmlを作って実行。UbuntuやSL6でもjdk1.8で同様に動作することを確認した。windows7でもjdk8_91というのができてbinにはappletviewerが入っていた。windows10もできる。コントロールパネル->ユーザー->環境変数の設定でPATHにProgram Files\Java\jdkxxx\binを含める。Java3dというのもあるが入れるのが面倒そう。Macだとjreが/Library以下に有ったりする。現状OpenGLがらみで動作しない。でもJOGLというのがあるそうだ。
nkf -g xxx.java //漢字コードを確認
nkf -w xxx.java > xxxU.java //UTFに変換(Mac Linux, windows はS-JIS)
vi xxxU.java // クラス名を一致させる 初めの一行のみ変更
javac xxxU.java
appletviewer xxxU.html //ウェブのhtmlのソースをコピーして作成
* 質量比1:3のスーパーボール団子 [#e34a50a3]
重い方のスーパーボールに針金をさして、軽い方はボール盤で穴を開けておき、針金を摘んで落とすと、重い方は原理的に静止して、すべてのエネルギーが軽い方に行くというのは意外と知られていないかも。
* キャベンディッシュの実験 [#o2e4b9ac]
http://fnorio.com/0006Chavendish/Chavendish.htm
原理はわかったので、昔からある装置を動かそうとしたら、つり線が切れた。古い機械なのでつり線は交換不能。島津のは32万円。
*重力場内のエネルギー保存則からわかること [#u24a0b74]
ローレンス クラウス教授の３回目の講義で、ハッブル定数と宇宙の質量分布を使って、宇宙の運命を予想する話があった。膨張宇宙の説明などは目から鱗の説明で、面白かった。秒速200km/sの銀河の公転速度から、銀河の質量が導ける話もあった。
'A Universe from nothing'で検索すると動画が見れる。
*単振動の数値計算 [#k12e04cd]
ニュートンが見いだした運動の法則を単振動に適用して、計算機で数値解を求めてみましょう。gfortranかg77を使ってデータファイルを作ります。その後、gnuplotで結果を見てみます。筑波大学の原子核理論のHPにあったプログラムです。
#ref(springEuler.f)
gfortran -o springEuler springEuler.f
./springEuler
gnuplot
plot 'Newton_Euler.dat', 'Newton_Euler.dat' using 1:3
Runge Kutta法を用いると精度が上がります。同じようにプロットを作成しましょう。
#ref(springRunge.f)
* 才差運動 [#se4b1853]
力学的世界観とは自分を地球の外に持っていく視点を持つこと。例えば「地球はでっかいジャイロスコープ」なんてことがピンとくる人になること。おもちゃの「地球ごま」などもそのたぐい。おもちゃを地球に例えるわけです。そしてその原理を角運動量の保存という物理法則に求めることができればあなたも立派な物理愛好家。
* 2重振り子の運動 [#g73dd1b9]
この運動方程式は一般化座標を使って解くことができる。変分法を使うなどして、棒の長さがそれぞれ&math(\ell_1);と&math(\ell_2);で質量が&math(m_1);と&math(m_2);、鉛直となす角が&math(\theta_1);と&math(\theta_2);の振り子の運動エネルギー&math(T);を&math(h_1=\ell_1/2);、&math(h_2=\ell_2/2);として求めると
#math(\frac{1}{2}m_1 h_1^2\dot{\theta}_1^2+\frac{1}{2}m_2\{\ell_1^2\dot{\theta}_1^2+h_2^2\dot{\theta}_2^2+2\ell_1 h_2\dot{\theta}_1 \dot{\theta}_2\cos(\theta_1-\theta_2)\})
となる。位置エネルギー&math(P);は幾何学的に求まる。不思議なことに最小作用の原理を適用するのは&math(T+P);でなく&math(T-P);である。これから方程式を作って、ルンゲクッタ法で時間発展を計算する。
* 2重振り子の作り方 [#a0c0263c]
http://www2.kek.jp/engineer/oho/giken/procedng/paper/met005.pdf
ベアリングはミニ四駆の四隅につけるガードベアリングをヨドバシで買った。アルミの棒に穴をあけて取り付ければOKなはず。M2.6のタップセットをアマゾンで注文した。
二重振り子のラグランジアンは変分の考え方で、運動エネルギー引くポテンシャルエネルギーを計算して作る。作った後の時間発展は、ルンゲクッタ法の良い利用例となる。やってみよう。
* 浮力の実験が力学でベクトルを学ぶ時、有効 [#h0ef2e5d]
物理天秤というものがある。水の中の重さの向きを変える実験が行える。
![[PukiWiki]](image/pukiwiki.png "[PukiWiki]")
