Calculating the Position of Center of Gravity for a Slender Rod Satellite

AI Thread Summary
The discussion focuses on calculating the center of gravity (CoG) for a slender rod satellite in a circular orbit, establishing the relationship between the center of gravity (ρ_G) and the center of mass (ρ_C) using the formula ρ_G = ρ_C (1 + l²/4ρ_C²)^(1/4). Participants clarify the definitions of the position vectors and the orientation of the rod, confirming it is perpendicular to the position vector. The integral setup for calculating gravitational force is debated, with suggestions to use angle-based relationships for accuracy due to the rod's orientation. The conversation highlights the importance of considering directionality in the force calculations, ultimately leading to a successful evaluation of the integral. The discussion concludes with acknowledgment of the complexity of the problem and the effectiveness of the proposed methods.
sun18
Messages
16
Reaction score
1

Homework Statement


Consider a slender rod satellite in a circular orbit. Show that the position of the center of gravity, ##\rho_G##, can be written in terms of the position of the center of mass, ##\rho_C##, as: ##\rho_G = \rho_C (1+\frac{l^2}{4\rho_C^2})^\frac{1}{4} ##, where ## l## is the length of the rod.

I'm not sure how to best describe this, but the rod is oriented like this, perpendicular to the position vector:
|
|
##\rho_C##--------------------##\rho##-----------------------M
|
|

Homework Equations


I know the total force on the rod acts on the CoG:
## F = \frac{\mu mM}{\rho_G^2} ##

m is the total mass of the bar, M is the mass of the attracting body, ##\rho## is a position vector.
I need to set this force equal to an integral, which I believe goes something like this:

## F = \int_{0}^{m} \frac{\mu M dm}{\rho^2} ##

The Attempt at a Solution


I changed to a length integral (half of the bar multiplied by 2), and have this:
##F = 2 \int_{\rho_C}^{\sqrt{\rho_C^2+l^2/4}} \frac{\mu Mmd \rho }{l \rho^2} ##

using the fact that ##dm = \frac{m}{l} d\rho##

I'm not confident that this is correct, as I can't seem to get the right answer. Something seems off about this integral, but I'm not sure what.
 
Physics news on Phys.org
Can you define the terms?
ρ = position vector = distance from the Earth's center of mass M to any arbitrary point?
ρC = distance from the Earth's center of mass M to the cylinder's center of mass m?
ρG = distance from the Earth's center of mass M to the cylinder's center of gravity?

Are you sure the cylinder is oriented perpendicular to the position vector?
 
Last edited:
You are correct on all three accounts:
ρ is the position vector from the Earth's center (mass M), to the infinitesimal mass element dm.
ρc is the distance from Earth's center to the rod's (mass m) center of mass
ρG is the distance from Earth's center to the rod's center of gravity.

And yes, the rod is perpendicular to the position vector. The first part of the problem had the rod aligned "vertically" with respect to the surface of the earth. In this case, it is perpendicular.
 
OK, then it looks like what you are doing is taking ρ from ρC to one end of the rod in your integral and then using symmetry to double it.

I am not sure that you can do this since this would give you the same center of gravity as a parallel rod with length = √(ρC2+(L/2)2) -ρC. Maybe it works out to be the same.

Could you instead take the angle θ and create an integral from θ=0 to L/(2ρC) ?
 
Hm I hadn't thought about using the angle. I will give that a try tomorrow and report back. Thanks for the responses by the way, I really appreciate it.
 
sun18 said:
##dm = \frac{m}{l} d\rho##
Doesn't look right to me. Working with the angle should make it clearer.
paisiello2 said:
what you are doing is taking ρ from ρC to one end of the rod in your integral and then using symmetry to double it.
I think that's ok.
 
dm=m/l dρ
haruspex said:
Doesn't look right to me. Working with the angle should make it clearer.
Seems OK to me. m/l is the density of the rod.
 
paisiello2 said:
dm=m/l dρ

m/l is the density of the rod.
Sure, but ρ is the distance from Earth's centre, not a distance along the rod.
 
But dρ is along the length of the rod.
 
  • #10
paisiello2 said:
But dρ is along the length of the rod.
Eh?! dρ is a change in distance from Earth's centre as we move along an element of the rod. That will not be the same as the length of the element.
 
  • #11
Yes, of course, you're right. I was thinking for a moment the rod was parallel.

That's probably why he was getting the wrong answer.
 
  • #12
Hello sun18,

It's an interesting problem, for sure! Anyway I worked through it and there are two mistakes in your general approach and perhaps another quibble about your notational semantics. I'll start with notation part.

sun18 said:

Homework Statement


Consider a slender rod satellite in a circular orbit. Show that the position of the center of gravity, ##\rho_G##, can be written in terms of the position of the center of mass, ##\rho_C##, as: ##\rho_G = \rho_C (1+\frac{l^2}{4\rho_C^2})^\frac{1}{4} ##, where ## l## is the length of the rod.

I'm not sure how to best describe this, but the rod is oriented like this, perpendicular to the position vector:
|
|
##\rho_C##--------------------##\rho##-----------------------M
|
|

Homework Equations


I know the total force on the rod acts on the CoG:
## F = \frac{\mu mM}{\rho_G^2} ##

m is the total mass of the bar, M is the mass of the attracting body, ##\rho## is a position vector.
I need to set this force equal to an integral, which I believe goes something like this:

## F = \int_{0}^{m} \frac{\mu M dm}{\rho^2} ##
Not bad so far. But let me comment on the notation.

Newton's universal law of gravitation for point mass particles is often expressed as, using \vec \rho as the displacement vector,
\vec F = G \frac{Mm}{\rho^2} \ \hat \rho

That said, it's often expressed using \mu = GM as
\vec F = \frac{\mu m}{\rho^2} \hat \rho

where \mu is called the "standard gravitational parameter" and has the mass of the gravitating body (often Earth when dealing with things orbiting Earth) folded into it along with Newton's gravitational constant.

In the notation you were using, you've simply replaced G with \mu, but didn't fold M into it. The notation that I've seen most often has M folded into the \mu, along with the G.

Of course, this is all just notational semantics, and isn't important on the grand scheme of things. But for what it's worth, I thought I should mention it.

The Attempt at a Solution


I changed to a length integral (half of the bar multiplied by 2), and have this:
##F = 2 \int_{\rho_C}^{\sqrt{\rho_C^2+l^2/4}} \frac{\mu Mmd \rho }{l \rho^2} ##

using the fact that ##dm = \frac{m}{l} d\rho##

There's one problem right there.

The change in mass does not vary linearly with the change in \rho in this case. It does vary linearly in the case where the rod is parallel to the displacement vector; in that case it would work fine. But here, with the rod perpendicular to the displacement vector, it does not hold true.

Let's define x as the distance from \rho_c to some point along the rod. So you can say easily,
dm = \frac{m}{\ell} dx

You should be able to form a relationship between x and \rho knowing that they form a right triangle. \rho^2 = \rho_c^2 + x^2.

Now for the second problem: you are forgetting about the pesky unit vector \hat \rho.

In the case where the rod was parallel to the position vector, the direction of the force wasn't important; the force on every individual dm all pointed in the same direction, so the unit direction vector didn't really matter. But in this case, where the rod is perpendicular, the each of the individual force directions are different.

You need to break up the \hat \rho into its x- and y-components. It's in the form, \hat \rho = a_x \hat x + a_y \hat y. Due to symmetry you know that the horizontal (x-) components will all ultimately cancel. But the y-components will not. You need to put that into your expression: The way you express the unit vector for a given dm should be should be a function of \rho_c and x. [Edit: or alternatively, express things in polar coordinates involving an angle. Either way will work, as long as you take directionality into account one way or the other.]

Good luck! :wink:

[Edit: Looks like @haruspex and others beat me to the punch. I suppose I should either type faster or write shorter. :smile:]
 
Last edited:
  • #13
Oh right of course, my mistake. I was forgetting that μ absorbs G and M. Luckily it should end up cancelling out anyways.

I had a feeling that the directionality was non-trivial. The first part of the problem had the rod colinear with the displacement vector ρ, so it didn't matter then.

So I had the integral ## \vec{F} = 2 \int_{\rho_C}^{\sqrt{\rho_C^2+l^2/4}} \frac{\mu md\rho}{l\rho^2}\hat{\rho} ##

Using ## \rho^2 = \rho_C^2 + x^2 ##, ## \vec{\rho} = \rho_C \hat{y} + x\hat{x}##, ## \hat{\rho} = \vec{\rho}/\rho ##, and changing the bounds of my integral, I have
## \vec{F} = \int_{-l/2}^{l/2} \frac{\mu mdx}{l\rho^2}\hat{\rho} ## (I didn't replace the ρ's yet for brevity)

I evaluated the integral, and as you said, the ##\hat{x}## components cancel, and I'm left with this integral:
## \frac{\mu m \rho_C \hat{y}}{l} \int_{-l/2}^{l/2} \frac{dx}{(\rho_C^2+x^2)^{\frac{3}{2}}} ##

**EDIT: Yep! This got me the right answer. It's a tricky integral to do, I wonder if it's easier to do in polar coordinates?

Thanks again for all the responses, I appreciate it!
 
Last edited:
  • Like
Likes collinsmark
Back
Top