Scribd
Upload
443 views3 pages

Matlab Moment Inertia

This document discusses calculating moments of inertia for an elliptical surface. It presents the equations needed to calculate the moment of inertia about the centroidal x-axis (Ix') and the parallel x-axis (Ix). It uses MATLAB functions and the quad() integration tool to numerically evaluate the integrals and calculate that Ix' is 4890 cm4, the ellipse area is 241.2861 cm2, and applying the parallel axis theorem, Ix is 66660 cm4.

Uploaded by

Kovačević Darko
Copyright
© Attribution Non-Commercial (BY-NC)
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
443 views3 pages

Matlab Moment Inertia

This document discusses calculating moments of inertia for an elliptical surface. It presents the equations needed to calculate the moment of inertia about the centroidal x-axis (Ix') and the parallel x-axis (Ix). It uses MATLAB functions and the quad() integration tool to numerically evaluate the integrals and calculate that Ix' is 4890 cm4, the ellipse area is 241.2861 cm2, and applying the parallel axis theorem, Ix is 66660 cm4.

Uploaded by

Kovačević Darko
Copyright
© Attribution Non-Commercial (BY-NC)
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
You are on page 1/ 3

11

Calculating Moments of Inertia


Ref: Hibbeler 10.1-10.2, Bedford & Fowler: Statics 8.1-8.2

Calculating moments of inertia requires evaluating integrals. This can be accomplished either symbolically, or using numerical approximations. Mathcads ability to integrate functions to generate numerical results is illustrated here. Example: Moment of Inertia of an Elliptical Surface Determine the moment of inertia of the ellipse illustrated below with respect to a) the centroidal x axis, and b) the x axis. The equation of the ellipse relative to centroidal axes is
x' 2 8
2

y' 2 14 2

=1
y -x x = x'

dy' y' C x'

d = dy = 16 x

In this problem, x and y have units of cm.

Solution The moment of inertia about the centroidal x axis is defined by the equation
I x ' = y' 2 dA
A

where dA is the area of the differential element indicated in the figure above.
dA = 2 x dy'

So, the integral for moment of inertia becomes

I x ' = y' 2 2 x dy'


A

Furthermore, x (or x) can be related to y using the equation of the ellipse. Note: Because of the location of the axes, x = x in this example.
y' 2 x = x ' = 8 2 1 2 14

The equation for the moment of inertia becomes:


I x' =

y' 2 y ' 2 2 8 2 1 2 14 8
8

dy'

To perform this integration we need to place the integrand in an m-file function and call MATLABs quad() function on the m-file. function Ix_integrand = Moment_Of_Inertia_Integrand(y_prime) %Saved as Moment_Of_Inertia_Integrand.m in the MATLAB %search path. x = sqrt(8.^2 .* (1 - y_prime.^2./14^2)); Ix_integrand = y_prime.^2 .* 2 .* x; Ix_prime = quad('Moment_Of_Inertia_Integrand',-8,8) %cm^4

Ix_prime =
4890 The moment of inertia relative to the original x axis can be found using the parallel-axis theorem.
I x = I x' + A d y 2

Where A is the area of the ellipse, and dy is the displacement of the centroidal y axis from the original y axis. The required area can be calculated by integration in the same fashion as before. function A = Area_Integrand(y_prime) %Saved as Area_Integrand.m in the MATLAB search path. x = sqrt(8.^2 .* (1 - y_prime.^2 ./ 14^2)); A = 2 .* x; A = quad('Area_Integrand',-8,8) %cm^2 A= 241.2861 Then the moment of inertia about x can be determined. dy = 16; %cm

I_x = Ix_prime + A .* dy.^2

%cm^4

I_x =
66660 Annotated MATLAB Script Solution function Ix_integrand = Moment_Of_Inertia_Integrand(y_prime) %Saved as Moment_Of_Inertia_Integrand.m in the MATLAB search path. x = sqrt(8.^2 .* (1 - y_prime.^2./14^2)); Ix_integrand = y_prime.^2 .* 2 .* x; function A = Area_Integrand(y_prime) %Saved as Area_Integrand.m in the MATLAB search path. x = sqrt(8.^2 .* (1 - y_prime.^2 ./ 14^2)); A = 2 .* x;

%Calculate the moment of inertia relative to the x' axis. Ix_prime = quad('Moment_Of_Inertia_Integrand',-8,8);%cm^4 fprintf('\nThe moment of inertia relative to the x'' axis is %1.0f cm^4\n', Ix_prime) %Calculate the area of the ellipse. A = quad('Area_Integrand',-8,8);%cm^2 fprintf('The area of the ellipse is %1.2f cm^2\n', A) %Use the parallel-axis theorm to calculate the moment of inertia relative to the x axis dy = 16;%cm I_x = Ix_prime + A .* dy.^2;%cm^4 fprintf('The moment of inertia relative to the x axis is = %1.0f cm^4\n', I_x)

You might also like