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.
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.
% 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
17.3 KB Views: 2,543