質量のない2mの竿の両端にm1とm2の質量を取り付け、中央を支点とします。
この装置の2階の微分方程式の数値解をグラフに求めました。
m1とm2の差が1kgの場合のは、周期Tは約5秒前後になります。
より正確には、
m1:2kgとm2:1kgの場合はT=約4.5秒
3kgと2kgの場合はT=約5秒
4kgと3kgの場合はT=約5.5秒
尚、m1単体の場合は、つまりm2=0の場合は、公式通りT=2秒となります。
T=2π(L/g)^0.5=2π(1/9.8)^0.5=2
Q1)上記の数値解のT=約5秒等は正しいでしょうか?どうでしょうか?
例えば、2秒となる筈だとか
以上、宜しくお願いします。
No.3ベストアンサー
- 回答日時:
>Tを計算しました所下記となり、T=8秒をなりました。
各パラメータにどのような値を代入したのか分かりませんが、
>T=2π(L/g)^0.5=2π(1/9.8)^0.5=2
これの((m1+m2)/(m1-m2))^0.5倍なので、例えばm1=2 kg, m2=1 kgの場合を考えているであれば√3倍の3.5秒くらいになるはずです。
>Q2)振れ幅が微小とは、どれくらいでしょうか?
単振り子の周期が
>T=2π(L/g)^0.5
になるという導出の過程では、sinθ≒θとし単振動で近似していますので、この近似が有効である時です。
どの程度のずれを許容するかは、考えたい精度によって変わり得るので一概には言えません。
まぁ、θ^2/6が0.1を超える (θが45°くらいを超える)ようだと微小とはほとんどの場合に微小とは言えないでしょうし、
θ^2/6が0.01を超えない(θが15°くらいを超えない)のであれば現実的には微小といっていいことが多いのでは。
>なお、私は90度でやりました
「単振り子 厳密解」などで検索すると、いろいろ情報が出てくるかと思いますが、
最大の振れ幅が90°であるのなら、周期は単振動だと思った時の1.2倍くらいになるようです。
冒頭の質問に
>尚、m1単体の場合は、つまりm2=0の場合は、公式通りT=2秒となります。
とありますが、正しく単振り子の周期が計算できているのであれば、2.4秒前後になるべきで、「(単振動の)公式通り」になったのはむしろ正しくありません。(0.4秒の差を気にしていないという事なのかもしれませんが)
「m2=0の場合と比べると、周期が((m1+m2)/(m1-m2))^0.5倍になる」という部分に振れ幅は関係ないので、最大の振れ幅を90°で考えたのなら、この2.4秒の((m1+m2)/(m1-m2))^0.5倍になるかな。
Octaveの計算環境(シミレーション)
でTの計算していますので、m2=0でやり、8が出たと思います。
確認のためm1 = 2、m2 = 1では下記となりました。
>> m1 m1 = 2
>> m2 m2 = 1
>> T=2*pi*(L*(m1+m2)/g*(m1-m2))^0.5 T = 3.4764
グラフの角度曲線を適当に読んでいますので下記は誤差が大きと
思います。
m1 = 2、m2 = 1では、T=約4.5
m1 = 3、m2 = 2では、T=約5
m1 = 4、m2 = 3では、T=約5.5
今回の質問の目的は、2階の微分方程式が正しいか否かでした。
%==================================
% m1*L^2*dxdt(2,Θ) + m2*L^2*dxdt(2, Θ+angle)=(-c*dxdt(1)+( -m1Lg*sin(Θ)-m2Lg*sin(Θ+angle))
% (m1+m2)*L^2*dxdt(2,Θ) =(-c*dxdt(1)+(-sm1Lg*sin(Θ)-m2Lg*sin(Θ+angle))
感じで、微分方程式を作成、最終的にODE45には、数値解を求めています。
===================
function dx = f_pendulumA(t,x,L, g, m1, m2, c, N,D)
%L:竿の長さの半分
%g:重力加速度
%m1:x軸の右側の質量
%m2:x軸の左側の質量
%c=0.05:摩擦係数
%N,D=8 竿は中央で曲なし、πの180度
%angle=pi*N/D; %竿の曲がり
ml2=(m1+m2)*L^2;
dx(1) = x(2); %角速度、x(1):
angle=pi*N/D; %竿の曲
dx(2) = ( -c*x(2)+(-m1*L*g*sin(x(1))-m2*L*g*sin(x(1)+angle)) )/ml2;
end
===================
最終目標は、地球ゴマが床に平行に歳差運動をして、重心に
gが懸かりますが、水平を保ちます。
初速ゼロから始まり、所定の速度に達しますが、その動作を
微分方程式で記述したいわけです。
以上、また別途この件の質問があれば、お教え頂けないでしょうか?
No.2
- 回答日時:
重りm、長さL、重力加速度gの単振り子から、
m→m1+m2
L→L
g→(m1-m2)/(m1+m2) g
と置き換えたものがお考えの系に対応する事になっていると思うので、振れ幅が微小のときの周期は
T=2π(L(m1+m2)/g(m1-m2))^0.5
になるかと思いますので、お示しの計算結果はずれていそうですね。
微小でないのなら振れ幅で考えていたら上の式からずれるのでそのせいでかもしれませんが。
お世話になります。
Tを計算しました所下記となり、T=8秒をなりました。
Q2)振れ幅が微小とは、どれくらいでしょうか?
お教え頂ければやって見たいと思います。
なお、私は90度でやりました。
>> T=2*pi*(L*(m1+m2)/g*(m1-m2))^0.5
T = 8.02835969261797e+00
お手数ですが宜しくお願いします。
No.1
- 回答日時:
両端に取り付ける質量によって「回転運動の慣性モーメント」が変わるので、角加速度が変わります。
質量差による「トルク」が同じでも、「慣性モーメント」が異なれば「回転運動のしかた」要するに「角加速度」は異なります。
「慣性モーメント」はご存じですよね?
https://ja.wikipedia.org/wiki/%E6%85%A3%E6%80%A7 …
http://www.buturigaku.net/main01/RigidBody/Rigid …
https://www.sit.ac.jp/user/konishi/JPN/L_Support …
お探しのQ&Aが見つからない時は、教えて!gooで質問しましょう!
おすすめ情報
- ・漫画をレンタルでお得に読める!
- ・街中で見かけて「グッときた人」の思い出
- ・「一気に最後まで読んだ」本、教えて下さい!
- ・幼稚園時代「何組」でしたか?
- ・激凹みから立ち直る方法
- ・1つだけ過去を変えられるとしたら?
- ・【あるあるbot連動企画】あるあるbotに投稿したけど採用されなかったあるある募集
- ・【あるあるbot連動企画】フォロワー20万人のアカウントであなたのあるあるを披露してみませんか?
- ・映画のエンドロール観る派?観ない派?
- ・海外旅行から帰ってきたら、まず何を食べる?
- ・誕生日にもらった意外なもの
- ・天使と悪魔選手権
- ・ちょっと先の未来クイズ第2問
- ・【大喜利】【投稿~9/7】 ロボットの住む世界で流行ってる罰ゲームとは?
- ・推しミネラルウォーターはありますか?
- ・都道府県穴埋めゲーム
- ・この人頭いいなと思ったエピソード
- ・準・究極の選択
- ・ゆるやかでぃべーと タイムマシンを破壊すべきか。
- ・歩いた自慢大会
- ・許せない心理テスト
- ・字面がカッコいい英単語
- ・これ何て呼びますか Part2
- ・人生で一番思い出に残ってる靴
- ・ゆるやかでぃべーと すべての高校生はアルバイトをするべきだ。
- ・初めて自分の家と他人の家が違う、と意識した時
- ・単二電池
- ・チョコミントアイス
デイリーランキングこのカテゴリの人気デイリーQ&Aランキング
マンスリーランキングこのカテゴリの人気マンスリーQ&Aランキング
-
電池式振り子時計について
-
振り子時計の振り子はなぜ止ま...
-
物理のエッセンスの発展問題で...
-
力学の問題です。水平なレール...
-
力学的エネルギーの保存・弾性...
-
単振り子の振動を止める
-
振り子の長さと周期の関係のグ...
-
難問題の系統とその解き方、例...
-
振り子の長さをy かかる時間を...
-
単振り子
-
スマホ振り子ってdヘルスケア、...
-
物理の単振動のx=A sinωt とか...
-
ガリレオの振り子が止まる理由
-
高校物理で質問です。一定の加...
-
三ヶ月半の乳児連れで特急しお...
-
2重振り子がカオス運動を引き起...
-
イッツ・ア・スモールワールド...
-
単振動と慣性力
-
ふりこ 糸の長さ 一往復の時間
-
剛体振り子の周期
おすすめ情報