Analytical Dynamics of Variable-Mass System
Analytical Dynamics of Variable-Mass System
Analytical Dynamics of Variable-Mass System
A Presentation on
Prosanto Chokroborty
To learn, to share
Community of Physics
To learn, to share
Community of Physics
1 1
𝜎= σ𝑗 𝒗𝑗 ⋅ 𝑚ሶ𝑗 𝒗𝑗 = 𝑚ሶ𝑗 𝑁𝑖𝑗𝑘 𝑞𝑖ሶ 𝑞𝑘ሶ
2 2
Here we are considering a body with N number of ducts where mass can be absorbed or can be expelled.
𝑚ሶ 𝑗 is the rate of absorbing or expelling mass, 𝒗𝑗 is the velocity of absorbed or expelled mass at the jth duct.
To learn, to share
Community of Physics
𝑑 𝑚𝒗
=𝒇+𝝌 (1)
𝑑𝑡
𝑑 𝑰⋅𝝎
=𝒍+𝝃 (2)
𝑑𝑡
𝑑 𝑚𝑗 𝒗𝑗
= 𝒇𝑗 + 𝝌𝑗 (3)
𝑑𝑡
To learn, to share
Community of Physics
𝑑 𝑑 𝑑
(𝑚𝒗) ⋅ 𝛿𝒓 + 𝑑𝑡 (𝑰 ⋅ 𝝎) ⋅ 𝛿𝜽 + σ𝑗 𝑑𝑡 (𝑚𝑗 𝒗𝑗 ) ⋅ 𝛿𝒓𝑗 = 𝒇 ⋅ 𝛿𝒓 + 𝒍 ⋅ 𝛿𝜽 + 𝒇𝑗 ⋅ 𝛿𝒓𝑗 + 𝝌 ⋅ 𝛿𝒓 + 𝝃 ⋅ 𝛿𝜽 + 𝒇𝑗 ⋅ 𝒓𝑗 (4)
𝑑𝑡
We are considering a system where total virtual work by constraint forces is zero and
𝜕𝒓 𝜕𝜽 𝜕𝒓𝑗
𝒇 ⋅ 𝛿𝒓 + 𝒍 ⋅ 𝛿𝜽 + 𝒇𝑗 ⋅ 𝛿𝒓𝑗 = 𝒇⋅ +𝒍⋅ + 𝒇𝑗 ⋅ 𝛿𝑞𝑘
𝜕𝑞k 𝜕𝑞𝑘 𝜕𝑞𝑘
= Q k 𝛿𝑞𝑘 (5)
𝑑 𝜕𝒓 𝑑 𝜕𝜽 𝑑 𝜕𝒓𝑗
(𝑑𝑡 (𝑚𝒗) ⋅ 𝜕𝑞 + 𝑑𝑡 (𝑰 ⋅ 𝝎) ⋅ 𝜕𝑞𝑘
+ σ𝑗 𝑑𝑡
(𝑚𝑗 𝒗𝑗 ) ⋅ 𝜕𝑞 − 𝑄𝑘 ) 𝛿𝑞_𝑘 = 0 (6)
𝑘 𝑘
To learn, to share
Community of Physics
Now,
𝑑 𝑚𝒗 𝜕𝒓 𝑑 𝜕𝒓 𝑑 𝜕𝒓
⋅ = (𝑚𝒗 ⋅ ) − ( ) ⋅ 𝑚𝒗
𝑑𝑡 𝜕𝑞𝑘 𝑑𝑡 𝜕𝑞𝑘 𝑑𝑡 𝜕𝑞𝑘
𝑑 𝜕𝒗 𝜕𝒗
= 𝑚𝒗 ⋅ ሶ − ⋅ 𝑚𝒗
𝑑𝑡 𝜕𝑞𝑘 𝜕𝑞𝑘
ሶ 𝑡𝑟𝑎𝑛𝑠𝑙𝑎𝑡𝑖𝑜𝑛𝑎𝑙
𝑑 𝜕𝑆 ሶ 𝑡𝑟𝑎𝑛𝑠𝑙𝑎𝑡𝑖𝑜𝑛𝑎𝑙
𝜕𝑆
= − (7)
𝑑𝑡 𝜕ሶ 𝑞𝑘ሶ ሶ 𝑘
𝜕𝑞
1
Where, 𝑆𝑡𝑟𝑎𝑛𝑠𝑙𝑎𝑡𝑖𝑜𝑛𝑎𝑙 = 𝒗 ⋅ 𝑚𝒗
2
Similarly,
𝑑 𝜕𝜽 ሶ 𝑟𝑜𝑡𝑎𝑡𝑖𝑜𝑛𝑎𝑙
𝑑 𝜕𝑆 ሶ 𝑟𝑜𝑡𝑎𝑡𝑖𝑜𝑛𝑎𝑙
𝜕𝑆
𝑰⋅𝝎 ⋅ = − (8)
𝑑𝑡 𝜕𝑞𝑘 𝑑𝑡 𝜕ሶ 𝑞𝑘ሶ ሶ 𝑘
𝜕𝑞
1
Where, 𝑆𝑟𝑜𝑡𝑎𝑡𝑖𝑜𝑛𝑎𝑙 = 𝝎 ⋅ 𝑰 ⋅ 𝝎
2
To learn, to share
Community of Physics
Now final term can be dealt with in this way,
𝑑 𝜕𝒓𝑗 𝜕𝒓𝑗
σ𝑗 𝑚𝑗 𝒗𝑗 ⋅ 𝜕𝑞 = σ𝑗 𝑚ሶ 𝑗 𝒗𝑗 + 𝑚𝑗 𝒂𝑗 ⋅ 𝜕𝑞
𝑑𝑡 𝑘 𝑘
𝜕𝒓𝑗
= σ𝑗 𝑚ሶ 𝑗 𝒗𝑗 ⋅
𝜕𝑞𝑘
𝜕𝒗𝑗
= σ𝑗 𝑚ሶ 𝑗 𝒗𝑗 ⋅ 𝜕𝑞 ሶ
𝑘
𝜕𝜎ሶ
= 𝜕ሶ 𝑞𝑘ሶ
Here, the second order effect of 𝑚𝑗 𝒂𝑗 was assumed to be negligible.
Now accumulating all the terms we can rewrite (4),
ሶ 𝑡𝑟𝑎𝑛𝑠𝑙𝑎𝑡𝑖𝑜𝑛𝑎𝑙
𝑑 𝜕𝑆 ሶ 𝑡𝑟𝑎𝑛𝑠𝑙𝑎𝑡𝑖𝑜𝑛𝑎𝑙 𝑑 𝜕𝑆
𝜕𝑆 ሶ 𝑟𝑜𝑡𝑎𝑡𝑖𝑜𝑛𝑎𝑙 ሶ 𝑟𝑜𝑡𝑎𝑡𝑖𝑜𝑛𝑎𝑙
𝜕𝑆 𝜕𝜎 ሶ
− + − + = 𝑄𝑘
𝑑𝑡 𝜕ሶ 𝑞𝑘ሶ ሶ 𝑘
𝜕𝑞 𝑑𝑡 𝜕ሶ 𝑞𝑘ሶ ሶ 𝑘
𝜕𝑞 𝜕ሶ 𝑞𝑘ሶ
𝑑 𝜕𝑆 ሶ 𝜕𝑆ሶ 𝜕𝜎ሶ
− + = 𝑄𝑘 (9)
𝑑𝑡 𝜕ሶ 𝑞𝑘ሶ ሶ 𝑘
𝜕𝑞 𝜕ሶ 𝑞𝑘ሶ
To learn, to share
Community of Physics
Rocket Equation:
1 1
𝑆 = 2 𝒗 ⋅ 𝑚𝒗 = 2 𝑚𝑥ሶ 2 (i)
A rocket has only one duct to expel mass. If the relative speed of expelled mass
is 𝑢 ,
𝒗1 = 𝑥ሶ − 𝑢 𝒏 ෝ (ii)
So, the Meshcherskii function,
1
𝜎 = 2 𝒗1 ⋅ 𝑚ሶ 1 𝒗1
1 2
= 2 𝑚ሶ 1 𝑥ሶ − 𝑢 (iii)
To learn, to share
Community of Physics
To learn, to share
Community of Physics
Thank You!
To learn, to share