Mutual Gravitational Interaction

Three material points of masses m 1 , m 2 , m 3 m_{1},m_{2},m_{3} are at the vertices of an equilateral triangle of side d d . The system is rotating in free space in such a way that under the mutual gravitational interaction of the three particles, the system is neither expanding nor contracting. Find angular velocity ω \omega of this rotation.

Universal gravitational constant is G G .

Details and Assumptions
1) m 1 = 1 , m 2 = 2 , m 3 = 3 , d = 4 , G = 5 m_{1}=1,m_{2}=2, m_{3}=3, d=4, G=5


The answer is 0.684.

This section requires Javascript.
You are seeing this because something didn't load right. We suggest you, (a) try refreshing the page, (b) enabling javascript if it is disabled on your browser and, finally, (c) loading the non-javascript version of this page . We're sorry about the hassle.

1 solution

Mark Hennings
Sep 16, 2020

Suppose that the position vectors of the masses m 1 , m 2 , m 3 m_1,m_2,m_3 are a , b , c \mathbf{a},\mathbf{b},\mathbf{c} , respectively. The gravitational attraction on the the mass m 1 m_1 at a \mathbf{a} is thus G m 1 m 2 d 3 ( b a ) + G m 1 m 3 d 3 ( c a ) = G m 1 d 3 ( ( m 2 + m 3 ) a + m 2 b + m 3 c ) = G m 1 d 3 ( m 1 + m 2 + m 3 ) ( p a ) \begin{aligned} \frac{Gm_1m_2}{d^3}(\mathbf{b}-\mathbf{a}) + \frac{Gm_1m_3}{d^3}(\mathbf{c}-\mathbf{a}) & = \; \frac{Gm_1}{d^3}\big(-(m_2+m_3)\mathbf{a} + m_2\mathbf{b} + m_3\mathbf{c}\big) \\ & = \; \frac{Gm_1}{d^3}(m_1+m_2+m_3)(\mathbf{p} - \mathbf{a}) \end{aligned} where p = 1 m 1 + m 2 + m 3 ( m 1 a + m 2 b + m 3 c ) \mathbf{p} \; = \; \frac{1}{m_1+m_2+m_3}(m_1\mathbf{a}+m_2\mathbf{b}+m_3\mathbf{c}) In other words, the gravitational attraction of the other two masses results in a force on m 1 m_1 of magnitude G m 1 d 3 ( m 1 + m 2 + m 3 ) p a \frac{Gm_1}{d^3}(m_1+m_2+m_3)|\mathbf{p}-\mathbf{a}| directed towards the point P P with position vector p \mathbf{p} . Of course, P P is the centre of mass of the system. In the absence of external forces, P P remains constant throughout the motion. Similar calculations can be performed for the other two masses. with similar results.

The particles orbit about the centre of mass P P with angular velocity ω \omega , where G m 1 d 3 ( m 1 + m 2 + m 3 ) p a = m 1 p a ω 2 \frac{Gm_1}{d^3}(m_1+m_2+m_3)|\mathbf{p}-\mathbf{a}| \; = \; m_1|\mathbf{p}-\mathbf{a}|\omega^2 and hence ω 2 = G ( m 1 + m 2 + m 3 ) d 3 \omega^2 = \frac{G(m_1+m_2+m_3)}{d^3} In this case we obtain ω = 1 8 30 = 0.6845631969 \omega = \tfrac18\sqrt{30} = \boxed{0.6845631969} .

@Mark Hennings Yeah, very nice solution.Thanks.

Talulah Riley - 8 months, 4 weeks ago

@Mark Hennings if you don't mind , what is the concept behind writing answer upto 10 decimal places. ?

Talulah Riley - 8 months, 4 weeks ago

Log in to reply

None, except you asked for a decimal solution. Another time you could say that ω = 1 a b \omega = \tfrac{1}{a}\sqrt{b} where a , b a,b are positive integers and b b is square-free, and ask (for example) for a + b = 38 a+b=38 . I think this is a nicer way to handle exact answers.

Mark Hennings - 8 months, 4 weeks ago

@Mark Hennings can you help me in this problem
Thanks in advance

Talulah Riley - 8 months, 4 weeks ago

Log in to reply

Suppose that the normal reaction between the cylinder and the wedge is N N , and that the normal reaction between the wedge and the slope is R R . Let the friction force between the cylinder and the wedge be F F . Suppose that the acceleration of the wedge is a a , and that the angular acceleration of the cylinder is α \alpha . The wedge has acceleration a ( cos θ sin θ ) a\binom{-\cos\theta}{-\sin\theta} Since the cylinder rolls on the wedge, the acceleration of the point of contact of the cylinder with the wedge must be the same as that of the wedge. Thus the centre of mass of the cylinder has acceleration a ( cos θ sin θ ) + r α ( 1 0 ) a\binom{-\cos\theta}{-\sin\theta} + r\alpha\binom{1}{0} where r r is the radius of the cylinder. The equations of motion of the centre of mass of the cylinder and the wedge are. respectively ( F N m g ) = m a ( cos θ sin θ ) + m r α ( 1 0 ) ( F R sin θ R cos θ N M g ) = M a ( cos θ sin θ ) \begin{aligned} \binom{-F}{N - mg} & = \; ma\binom{-\cos\theta}{-\sin\theta} + mr\alpha\binom{1}{0} \\ \binom{F-R\sin\theta}{R\cos\theta - N - Mg} & = \; Ma\binom{-\cos\theta}{-\sin\theta} \end{aligned} and the angular equation of motion of the cylinder is F r = m r 2 α Fr \; = \; mr^2\alpha Thus we deduce that F = m r α F = mr\alpha and hence that r α = 1 2 a cos θ r\alpha = \tfrac12a\cos\theta , and, finally a = 2 ( M + m ) g sin θ 2 ( M + m ) m cos 2 θ a \; = \; \frac{2(M+m)g\sin\theta}{2(M+m) - m\cos^2\theta}

Mark Hennings - 8 months, 4 weeks ago

Log in to reply

@Mark Hennings Thank you so much for such a fast response.

Talulah Riley - 8 months, 4 weeks ago

@Foolish Learner You can see the solution above. The problem is very clear. And please delete the report.

Talulah Riley - 8 months, 4 weeks ago

@Mark Hennings What a nice solution.

Krishna Karthik - 8 months, 4 weeks ago

So, in all honestly, I could not solve this analytically. I did solve this numerically by trial and error by plotting the distance between two masses, varying ω \omega at each iteration, which should be constant at the appropriate value of ω \omega . The trajectory and successive distance between masses is as shown below:

Ideally, the distance should be exactly 4 4 but that is not so. Also, if I simulate for a large enough time, I see that the particles diverge.

Please do share your thoughts on this, if any. @Mark Hennings

I could not have concluded that the motion is circular based on this trajectory plot.

Karan Chatrath - 8 months, 4 weeks ago

Log in to reply

The value of ω \omega used here is after seeing the closed-form expression. In my iterations, I was further off the mark, so the trajectories diverged earlier. I think this also has to do with the numeric complexities. The system solution is very sensitive to the parameter ω \omega

Karan Chatrath - 8 months, 4 weeks ago

Log in to reply

In the study of the three-body problem, and in astronomy, there are things called Greek and Trojan orbits. They are precisely the positions that are possible for an object to orbit stably, while forming an equilateral triangle with (say) Jupiter and the Sun (as the two largest bodies in the solar system), although the Sun is often viewed as begin stationary, in which case the Trojans are points on the same circular orbit as Jupiter, but 6 0 60^\circ ahead or behind it. Other positions are not stable. I guess your numerical techniques are discovering the instability of not being in the Trojan.

Figure 3 here gives a nice animation of this solution.

Mark Hennings - 8 months, 4 weeks ago

Log in to reply

@Mark Hennings Thank you for your insights. I will take a deeper dive into this article when I can. It looks really interesting.

So if I understand correctly, my numeric approximation of ω \omega is what causes the solution to not be in the Trojan? And the system is indeed very sensitive to this parameter?

Karan Chatrath - 8 months, 4 weeks ago

Log in to reply

@Karan Chatrath That is why you find moonlets in the Trojan of a planet - it is only place they can stay for long...

Mark Hennings - 8 months, 4 weeks ago

Log in to reply

@Mark Hennings @Mark Hennings Exactly.

Talulah Riley - 8 months, 4 weeks ago

@Mark Hennings Sir can you please help me in this problem, whenever you will be free. Thanks in advance.

Talulah Riley - 8 months, 4 weeks ago

Log in to reply

Let x x be the extension of the string, so that x ¨ \ddot{x} is the horizontal acceleration of the centre of mass of the wheel, as well as the downward accelerations of the two masses. Let the tension in the top section of the string be T 1 T_1 , and the tension in the bottom section be T 2 T_2 . Let the angular velocity of the wheel be ω \omega . Let the friction force acting on the wheel from the table be F F . The equations of motion of the centre of the wheel, A A and B B are F + T 1 k x = m x ¨ T 2 + m g T 1 = m x ¨ 3 m g T 2 = 3 m x ¨ \begin{aligned} F + T_1 - kx & = \; m\ddot{x} \\ T_2 + mg - T_1 & = \; m\ddot{x} \\ 3mg - T_2 & = \; 3m\ddot{x} \end{aligned} respectively, while the angular equation of motion of the wheel is F r = m r 2 ω ˙ -Fr \; = \; mr^2\dot{\omega} where r r is the radius of the wheel. Since the wheel rolls on the table, we deduce that x ˙ r ω = 0 \dot{x} - r\omega = 0 , so that r ω = x ˙ r\omega = \dot{x} , and hence F = m x ¨ F = -m\ddot{x} . Adding the first three equations gives F + 4 m g k x = 5 m x ¨ 6 m x ¨ + k x = 4 m g \begin{aligned} F + 4mg - kx & = \; 5m\ddot{x} \\ 6m\ddot{x} + kx & = \; 4mg \end{aligned} and hence, since the system starts from rest with x = 0 x=0 , we deduce that x = 4 m g k [ 1 cos ( t k 6 m ) ] x \; = \; \frac{4mg}{k}\Big[1 - \cos\left(t\sqrt{\tfrac{k}{6m}}\right)\Big] As might be suspected, the system oscillates. Then T 2 = 3 m ( g x ¨ ) = m g [ 3 2 cos ( t k 6 m ) ] T_2 \; =\; 3m(g - \ddot{x}) \; = \; mg\Big[3 - 2\cos\left(t\sqrt{\tfrac{k}{6m}}\right)\Big] making the maximum value of T 2 T_2 equal to 5 m g 5mg .

Mark Hennings - 8 months, 4 weeks ago

Log in to reply

@Mark Hennings What does 1 st line mean, I didn't understand? Can you elaborate more.

Talulah Riley - 8 months, 4 weeks ago

Log in to reply

@Talulah Riley @Mark Hennings Thank you so much sir for the solution.
I understand the whole now.
Your mechanics skill are next level, awesome and very nice.
I have posted a note on Differential equations, and I am not able to mention you. Can you have a look on that.
Thanks in advance.

Talulah Riley - 8 months, 4 weeks ago

@Talulah Riley Sorry. My mistake. Your comment applied to a report on the 3-body question, not the supplementary question you asked me.

Mark Hennings - 8 months, 4 weeks ago

Log in to reply

@Mark Hennings @Mark Hennings Exactly,you may thought that when I said “see the above solution, I was showing your solution of this problem not the doubt which you have cleared me.
Cheers

Talulah Riley - 8 months, 4 weeks ago

@Mark Hennings @Mark Hennings Sir from now,whenever I ask you a doubt, would you like to see my attempt ???

Talulah Riley - 8 months, 4 weeks ago

Log in to reply

@Talulah Riley I am happy to receive questions - I cannot always promise to have the time to answer them!

Mark Hennings - 8 months, 4 weeks ago

Log in to reply

@Mark Hennings @Mark Hennings Ok No problem.
Sir in which place should I ask you doubt, by making a note of a problem or in any question in which you have posted a solution.

Talulah Riley - 8 months, 4 weeks ago

Log in to reply

@Talulah Riley Mention me in the problem itself - it makes it easier to find again!

Mark Hennings - 8 months, 4 weeks ago

Log in to reply

@Mark Hennings @Mark Hennings Thank you so much.
If you don't mind, generally in which time in a day you are most probably active on Brilliant?

Talulah Riley - 8 months, 4 weeks ago

0 pending reports

×

Problem Loading...

Note Loading...

Set Loading...