振り子 (4回目)

2024/1/3(水)
振り子 (4回目)
 
二重振り子(Double pendulum)
 
連立微分方程式の導出
(ニュートンの運動方程式)
 
 
■ 導出
▼ 定義
ℓi :紐の長さ(m)
mi :質点の質量(kg)
θi:鉛直下方向からの角度(rad)
Ti :張力
g  :重力加速度
 
 
▼ 位置r
xi :支点からの右方向の変位(m)
yi :支点からの下方向の変位(m)
r1 = (x1, y1) = (ℓ1sinθ1, ℓ1cosθ1)
r2 = (x2, y2)
= (x1 + ℓ2sinθ2, y1 + ℓ2cosθ2)
 
 
▼ 速度v
v1 = r(・)1 = (x(・)1, y(・)1) = (ℓ1θ(・)1cosθ1, -ℓ1θ(・)1sinθ1)
v2 = r(・)2 = (x(・)2, y(・)2)
= (x(・)1 + ℓ2θ(・)2cosθ2, y(・)1 - ℓ2θ(・)2sinθ2)
 
 
▼ 加速度a
a1 = v(・)1 = (x(･･)1, y(･･)1) = ℓ1(d/dt)(θ(・)1cosθ1, -θ(・)1sinθ1)
= ℓ1(θ(･･)1cosθ1 - θ(・)12sinθ1, -θ(･･)1sinθ1 - θ(・)12cosθ1)
a2 = v(・)2 = (x(･･)2, y(･･)2)
= (d/dt)(x(・)1 + ℓ2θ(・)2cosθ2, y(・)1 - ℓ2θ(・)2sinθ2)
= (x(･･)1 + ℓ2θ(･･)2cosθ2 - ℓ2θ(・)22sinθ2,
   y(･･)1 - ℓ2θ(･･)2sinθ2 - ℓ2θ(・)22cosθ2)
 
 
▼ ニュートンの運動方程式
(x(･･)1, y(･･)1)
= ℓ1(θ(･･)1cosθ1 - θ(・)12sinθ1, -θ(･･)1sinθ1 - θ(・)12cosθ1)
(x(･･)2, y(･･)2)
= (x(･･)1 + ℓ2θ(･･)2cosθ2 - ℓ2θ(・)22sinθ2,
   y(･･)1 - ℓ2θ(･･)2sinθ2 - ℓ2θ(・)22cosθ2)
 
ma = Fより
 
m1x(･･)1 = -T1sinθ1 + T2sinθ2 
m1y(･･)1 = -T1cosθ1 + T2cosθ2 + m1g
 
m2x(･･)2 = -T2sinθ2 
m2y(･･)2 = -T2cosθ2 + m2g
 
▼ 運動方程式(整理)
T2sinθ2 = -m2x(･･)2 
T2cosθ2 = -m2y(･･)2 + m2g
を
m1x(･･)1 = -T1sinθ1 + T2sinθ2 
m1y(･･)1 = -T1cosθ1 + T2cosθ2 + m1g
に代入
m1x(･･)1 + m2x(･･)2 = -T1sinθ1 
m1y(･･)1 + m2y(･･)2 = -T1cosθ1 + m2g + m1g
さらに
(m1x(･･)1 + m2x(･･)2)cosθ1 = -T1sinθ1cosθ1  
(m1y(･･)1 + m2y(･･)2)sinθ1 = -T1sinθ1cosθ1 + (m2 + m1)gsinθ1
差を取ると
(m1y(･･)1 + m2y(･･)2)sinθ1 - (m1x(･･)1 + m2x(･･)2)cosθ1 
= (m2 + m1)gsinθ1 
(m1y(･･)1 + m2y(･･)2 - m1g - m2g)sinθ1 = (m1x(･･)1 + m2x(･･)2)cosθ1 
 
また
 
m2x(･･)2 = -T2sinθ2 
m2y(･･)2 = -T2cosθ2 + m2g
より
m2x(･･)2cosθ2 = -T2sinθ2cosθ2 
m2y(･･)2sinθ2 = -T2sinθ2cosθ2 + m2gsinθ2
差を取ると
m2y(･･)2sinθ2 - m2x(･･)2cosθ2 = m2gsinθ2
(m2y(･･)2 - m2g)sinθ2 = m2x(･･)2cosθ2 
よって
 
(m1y(･･)1 + m2y(･･)2 - m1g - m2g)sinθ1 = (m1x(･･)1 + m2x(･･)2)cosθ1 
(m2y(･･)2 - m2g)sinθ2 = m2x(･･)2cosθ2 
となる
 
 
▼ 運動方程式(変形1)
(m1y(･･)1 + m2y(･･)2 - m1g - m2g)sinθ1 = (m1x(･･)1 + m2x(･･)2)cosθ1 
 
に
(x(･･)1, y(･･)1)
= ℓ1(θ(･･)1cosθ1 - θ(・)12sinθ1, -θ(･･)1sinθ1 - θ(・)12cosθ1)
(x(･･)2, y(･･)2)
= (x(･･)1 + ℓ2θ(･･)2cosθ2 - ℓ2θ(・)22sinθ2,
   y(･･)1 - ℓ2θ(･･)2sinθ2 - ℓ2θ(・)22cosθ2)
を代入
 
{-m1ℓ1(θ(･･)1sinθ1 + θ(・)12cosθ1)
- m2ℓ1(θ(･･)1sinθ1 + θ(・)12cosθ1)
- m2ℓ2(θ(･･)2sinθ2 + θ(・)22cosθ2) - m1g - m2g}sinθ1 
= {m1ℓ1(θ(･･)1cosθ1 - θ(・)12sinθ1 
+ m2(ℓ1θ(･･)1cosθ1 - ℓ1θ(・)12sinθ1 
+ ℓ2θ(･･)2cosθ2 - ℓ2θ(・)22sinθ2)}cosθ1 
 
- m1ℓ1θ(･･)1sin2θ1 - m1ℓ1θ(・)12cosθ1sinθ1 
- m2ℓ1θ(･･)1sin2θ1 - m2ℓ1θ(・)12cosθ1sinθ1 
- m2ℓ2θ(･･)2sinθ2sinθ1 - m2ℓ2θ(・)22cosθ2sinθ1 
- m1gsinθ1 - m2gsinθ1 
- m1ℓ1θ(･･)1cos2θ1 + m1ℓ1θ(・)12sinθ1cosθ1 
- m2ℓ1θ(･･)1cos2θ1 + m2ℓ1θ(・)12sinθ1cosθ1
- m2ℓ2θ(･･)2cosθ2cosθ1 + m2ℓ2θ(・)22sinθ2cosθ1 
= 0
 
m1ℓ1θ(･･)1sin2θ1 + m1ℓ1θ(･･)1cos2θ1 
+ m2ℓ1θ(･･)1cos2θ1 + m2ℓ1θ(･･)1sin2θ1 
+ m2ℓ2θ(･･)2cosθ1cosθ2 + m2ℓ2θ(･･)2sinθ1sinθ2 
+ m1ℓ1θ(・)12cosθ1sinθ1 - m1ℓ1θ(・)12sinθ1cosθ1 
+ m2ℓ1θ(・)12cosθ1sinθ1 - m2ℓ1θ(・)12sinθ1cosθ1
+ m2ℓ2θ(・)22cosθ2sinθ1 - m2ℓ2θ(・)22sinθ2cosθ1 
+ m1gsinθ1 + m2gsinθ1 
= 0
 
(m1 + m2)ℓ1θ(･･)1 
+ m2ℓ2θ(･･)2(cosθ1cosθ2 + sinθ1sinθ2)
+ m1ℓ1θ(・)12(sinθ1cosθ1 - cosθ1sinθ1)
+ m2ℓ1θ(・)12(sinθ1cosθ1 - cosθ1sinθ1)
+ m2ℓ2θ(・)22(sinθ1cosθ2 - cosθ1sinθ2)
+ (m1g + m2g)sinθ1 
= 0
 
(m1 + m2)ℓ1θ(･･)1 + m2ℓ2θ(･･)2cos(θ1 - θ2)
+ m2ℓ2θ(・)22sin(θ1 - θ2) + (m1g + m2g)sinθ1 
= 0
 
(m1 + m2)ℓ 1θ(･･)1 
+ m2ℓ2{θ(･･)2cos(θ1 - θ2) + θ(・)22sin(θ1 - θ2)}
= -(m1 + m2)gsinθ1  … ①
 
▼ 運動方程式(変形2)
(m2y(･･)2 - m2g)sinθ2 = m2x(･･)2cosθ2 
に
(x(･･)1, y(･･)1)
= ℓ1(θ(･･)1cosθ1 - θ(・)12sinθ1, -θ(･･)1sinθ1 - θ(・)12cosθ1)
(x(･･)2, y(･･)2)
= (x(･･)1 + ℓ2θ(･･)2cosθ2 - ℓ2θ(・)22sinθ2,
   y(･･)1 - ℓ2θ(･･)2sinθ2 - ℓ2θ(・)22cosθ2)
を代入
 
{m2(-ℓ1θ(･･)1sinθ1 - ℓ1θ(・)12cosθ1 
- ℓ2θ(･･)2sinθ2 - ℓ2θ(・)22cosθ2) - m2g}sinθ2 
= m2(ℓ1θ(･･)1cosθ1 - ℓ1θ(・)12sinθ1 
+ ℓ2θ(･･)2cosθ2 - ℓ2θ(・)22sinθ2)cosθ2 
 
- m2ℓ1θ(･･)1sinθ1sinθ2 - m2ℓ1θ(・)12cosθ1sinθ2 
- m2ℓ2θ(･･)2sin2θ2 - m2ℓ2θ(・)22cosθ2sinθ2 - m2gsinθ2  
- m2ℓ1θ(･･)1cosθ1cosθ2 + m2ℓ1θ(・)12sinθ1cosθ2  
- m2ℓ2θ(･･)2cos2θ2 + m2ℓ2θ(・)22sinθ2cosθ2 = 0
 
m2ℓ1θ(･･)1sinθ1sinθ2 + m2ℓ1θ(･･)1cosθ1cosθ2 
+ m2ℓ2θ(･･)2sin2θ2 + m2ℓ2θ(･･)2cos2θ2 
+ m2ℓ1θ(・)12cosθ1sinθ2 - m2ℓ1θ(・)12sinθ1cosθ2  
+ m2ℓ2θ(・)22cosθ2sinθ2 - m2ℓ2θ(・)22sinθ2cosθ2 
+ m2gsinθ2 = 0
 
m2ℓ1θ(･･)1(cosθ1cosθ2 + sinθ1sinθ2) + m2ℓ2θ(･･)2 
+ m2ℓ1θ(・)12(cosθ1sinθ2 - sinθ1cosθ2)
+ m2ℓ2θ(・)22(cosθ2sinθ2 - sinθ2cosθ2)
+ m2gsinθ2 = 0
 
m2ℓ1θ(･･)1cos(θ1 - θ2) + m2ℓ2θ(･･)2 
- m2ℓ1θ(・)12sin(θ1 - θ2) + m2gsinθ2 = 0
 
ℓ2θ(･･)2 + ℓ1{θ(･･)1cos(θ1 - θ2) - θ(・)12sin(θ1 - θ2)}
= - gsinθ2  … ②
 
 
 
■ 結果
▼ 定義
ℓi :紐の長さ(m)
mi :質点の質量(kg)
θi:鉛直下方向からの角度(rad)
g  :重力加速度
xi :支点からの右方向の変位(m)
yi :支点からの下方向の変位(m)
 
 
▼ 連立微分方程式
(m1 + m2)ℓ 1θ(･･)1 
+ m2ℓ2{θ(･･)2cos(θ1 - θ2) + θ(・)22sin(θ1 - θ2)}
= -(m1 + m2)gsinθ1  … ①
 
ℓ2θ(･･)2 + ℓ1{θ(･･)1cos(θ1 - θ2) - θ(・)12sin(θ1 - θ2)}
= - gsinθ2  … ②
 
 
 

このブログの人気の投稿

NEWS

2024/12/27(金 ) N88-BASICでハノイの塔 (2回目)(再帰なし) を公開しました    2024/12/25(水) VL,NL,XL-BASIC(ver～28w3) の?LBasic01.txt修正版公開しました    2024/12/24(火) VL,NL,XL-BASIC(ver～28w3)  を公開しました    2024/12/21(土 ) ハノイの塔 (3回目)  を公開しました
続きを読む

N88-BASICでゲーム (1回目)

イメージ
2021/8/16(月) N88-BASICでゲーム  (1回目) by ULproject   球同士の衝突と跳返り   大文字ベクトル、小文字スカラー 文字の後の、 'は衝突後、数字は球番号   球同士の反発係数、衝突面に水平 ,垂直u,e 球の中心位置 P,半径r,質量m,速度V 衝突面の球 1方向の法線N = (P 1 -P 2 )/|P 1 -P 2 |   球同士の衝突は、 |P 1 -P 2 | ≦ r 1 +r 2   なので √{(P 1 -P 2 )･(P 1 -P 2 )} ≦ r 1 +r 2   (P 1 -P 2 )･(P 1 -P 2 ) ≦ (r 1 +r 2 )(r 1 +r 2 ) で判断できます。   ここで、ベクトルの計算ですが、 P(i) (i=0,1,2) をPのx,y,z成分と すると、 P = P 1  - P 2 は、 P(i) = P1(i) - P2(i) (i=0,1,2) P･Pは ΣP(i)P(i) (i=0,1,2) となります   BASICでは、 P(0) = P1(0) - P2(0) P(1) = P1(1) - P2(1) P(2) = P1(2) - P2(2) D = P(0)*P(0)+P(1)*P(1)+P(2)*P(2) R = R1 + R2 IF D <= R*R THEN 重なり又は接触 となります。   反発させるときは、球が接触している 間は再び反発しないようにしないと 球が離れなくなります。   次に反発についてですが、衝突面に 水平な成分は、球の回転に関係します が、ここでは球の回転は考えずに 垂直成分と同じように反発係数で 減速させる事にします。   衝突面に垂直な成分 (v 1 は球 1の速度、v 1 'は衝突後) 運動量 p = mv 運動量保存則 p 1 +p 2 =p 1 '+ p 2 'より p 1 +p 2  = m 1 v 1 '+ m 2 v 2 ' 反発係数 e = -(v 1 '- v 2 ')/(v 1  - v 2 )より ev 1  - ev 2  = -v 1 '+ v ...
続きを読む