• Support PF! Buy your school textbooks, materials and every day products via PF Here!

Inertia tensor of cone around its apex

3
0
Thread moved from the technical forums, so no Homework Template is shown
Im trying to calculate the principals moments of inertia (Ixx Iyy Izz) for the inertia tensor by triple integration using cylindrical coordinates in MATLAB.

inertiatensor.png


Code:
% Symbolic variables
syms r z theta R h M; % R (Radius) h(height) M(Mass)  

% Ixx
unox = int((z^2+(r*sin(theta))^2)*r,z,r,h); % First integration
dosx = int(unox,r,0,R); % Second Integration
tresx = int(dosx,theta,0,2*pi); % Third integration

Ix = tresx / ((pi*R^2*h)/3); % Division by the volume of a cone (Mass since densitiy is 1)
Ix = Ix*M  % Setting up in terms of mass M

% Iyy
unoy = int((z^2+(r*cos(theta))^2)*r,z,r,h);
dosy = int(unoy,r,0,R);
tresy = int(dosy,theta,0,2*pi);

Iy = tresy / ((pi*R^2*h)/3);
Iy = Iy*M

% Izz
unoz = int(r^2*r,z,r,h);
dosz = int(unoz,r,0,R);
tresz = int(dosz,theta,0,2*pi);

Iz = tresz / ((pi*R^2*h)/3);
Iz = Iz*M
I plug in test values in both the formulas from the image and the code and they seem correct except from Iz in which the height (h) still appears and terms like 3h-2r instead of h appear, Am I missing something in the integrations? Thank you in advance.
 

Attachments

Orodruin

Staff Emeritus
Science Advisor
Homework Helper
Insights Author
Gold Member
2018 Award
15,944
5,917
Did you try doing the integrals by hand? (They are not difficult.)

Please show your work in determining the integration boundaries.
 
3
0
Did you try doing the integrals by hand? (They are not difficult.)

Please show your work in determining the integration boundaries.
In cylindrical coordinates for a solid cone the only non-trivial boundarie would be on the first integration, since its a solid of revolution, theta goes grom 0 to 2*pi and r goes from 0 to R (radius of cone).
z would go from z = r (x^2 + y^2) to z = h (height of cone)
 

Orodruin

Staff Emeritus
Science Advisor
Homework Helper
Insights Author
Gold Member
2018 Award
15,944
5,917
z would go from z = r (x^2 + y^2) to z = h (height of cone)
This is incorrect and you have not really provided an argument for why it should be the case. I suggest you take a second look at it.
 
3
0
This is incorrect and you have not really provided an argument for why it should be the case. I suggest you take a second look at it.
Thank you, I found the error, the inferior limit of the integral for z isnt r, its (H/R) * r, The textbook I was reading gives an example for a cone of equal heigth and radius so it ignored the slope in the integral.
 

Want to reply to this thread?

"Inertia tensor of cone around its apex" You must log in or register to reply here.

Related Threads for: Inertia tensor of cone around its apex

Replies
3
Views
11K
Replies
0
Views
2K
  • Posted
Replies
1
Views
5K
  • Posted
Replies
2
Views
4K
  • Posted
Replies
1
Views
3K
Replies
1
Views
982
  • Posted
Replies
8
Views
4K
Replies
6
Views
8K

Physics Forums Values

We Value Quality
• Topics based on mainstream science
• Proper English grammar and spelling
We Value Civility
• Positive and compassionate attitudes
• Patience while debating
We Value Productivity
• Disciplined to remain on-topic
• Recognition of own weaknesses
• Solo and co-op problem solving
Top