Programas Diseñados para Calcular Momentos de Inercia de Diferentes Figuras-Carla Toapantaleón

Descargar como docx, pdf o txt
Descargar como docx, pdf o txt
Está en la página 1de 9

ESCUELA POLITÉCNICA NACIONAL

I. DATOS INFORMATIVOS:

Nombre : Carla Toapanta León


Facultad: Ingeniería Civil y Ambiental
Carrera : Ingeniería Civil
Fecha de Entrega: 2015/07/15
Tema : Programas para calcular momentos de inercia de
diferentes figuras simétricas, variando los valores de sus dimensiones.

Programación

Deber N°

EJERCICIO N° 1

clc
clear all
close all
disp('Ingrese 1 para seleccionar un círculo o cualquier otro número para
seleccionar un rectángulo: ')
a =input('Ingresar a=');
if a==1
r= input ('Ingresar el radio del círculo: ');
Ix= (pi*r^4)/4
Iy= Ix
else
disp('Ingrese 1 para calcular la inercia del rectángulo con respecto
al centroide u otro número para calcular con respecto al vértice inferior
izquierdo')
f=input('Ingresar f=');
b=input('Ingresar base b=');
h=input('Ingresar altura h=');
if f==1
Ixcg=(1/12)*(b^3)*h
Iycg=(1/12)*(h^3)*b
else
Ix=(1/3)*(b^3)*h
Iy=(1/3)*(h^3)*b
disp('Fin del programa.’)
end
end
EJERCICIO N° 2

%. Realizar un programa que permita calcular el momento de inercia de la


figura con forma de T, siendo esta simétrica.

clc
clear all
close all
disp('"Inercia para calcular una figura con forma de T, siendo esta
simétrica"')
a= input('ingrese valor m: ');
b= input('ingrese valor n: ');
c= input('ingrese valor o: ');
d= input('ingrese valor p: ');
A1=m*n;
A2=o*p;
A=A1+A2;
x1=m/2;
x2=m/2;
y1=p+n/2;
y2=p/2;
Ax=A1*x1+A2*x2;
Ay=A1*y1+A2*y2;
Xc=Ax/A;
Yc=Ay/A;
Ix1=1/12*m*n^3;
Ix2=1/12*o*p^3;
Iy1=1/12*n*m^3;
Iy2=1/12*p*o^3;
Ix=Ix1+Ix2;
Iy=Iy1+Iy2;
dx1=(y1-Yc);
dx2=y2-Yc;
dy1=x1-Xc;
dy2=x2-Xc;
Adx=A1*dx1^2+A2*dx2^2;
Ady=A1*dy1^2+A2*dy2^2;
disp('El Momento de Inercia con respecto al eje x es:')
Ixxcg=Ix+Adx
disp('El Momento de Inercia con respecto al eje y es:')
Iyycg=Iy+Ady
disp('Fin del programa.')
EJERCICIO N° 3

%. Realizar un programa que permita calcular el momento de inercia de la


figura con forma de I, siendo esta simétrica.

clc
clear all
close all
disp('"Inercia para calcular una figura con forma de I, siendo esta
simétrica"')
i=input('Ingrese el valor i: ');
j=input('Ingrese el valor j: ');
k=input('Ingrese el valor k: ');
l=input('Ingrese el valor l: ');
A1=i*k;
A2=(2*k+l)*j;
At=4*A1+A2;
x1=i/2;
x2=1.5*i+j;
x3=i+j/2;
y1=1.5*k+l;
y2=k/2;
y3=k+l/2;
Ax=2*A1*(x1+x2)+A2*x3;
Ay=2*A1*(y1+y2)+A2*y3;
xc=Ax/At;
yc=Ay/At;
Ix1=(i*k^3)/12;
Ix2=(j*(2*k+l)^3)/12;
Iy1=(k*i^3)/12;
Iy2=((2*k+l)*j^3)/12;
Ix01=Ix1+A1*y1^2;
Ix02=Ix1+A1*y2^2;
Ix03=Ix2+A2*y3^2;
Iy01=Iy1+A1*x1^2;
Iy02=Iy1+A1*x2^2;
Iy03=Iy2+A2*x3^2;
Ixxcg=2*Ix01+2*Ix02+Ix03;
Iyycg=2*Iy01+2*Iy02+Iy03;
disp('El Momento de Inercia con respecto al eje x es:')
Ixc=Ixxcg-At*yc^2
disp('El Momento de Inercia con respecto al eje y es:')
Iyc=Iyycg-At*xc^2
disp('Fin del programa.')
EJERCICIO N° 4

%. Realizar un programa que permita calcular el momento de inercia de la


figura con forma de C(SIMPLE, siendo esta simétrica.

clc
clear all
close all
disp('"Inercia para calcular una figura con forma de C(SIMPLE), siendo
esta simétrica"')
r=input('ingrese el valor r: ')
s=input('ingrese el valor s: ')
t=input('ingrese el valor t: ')
u=input('ingrese el valor u: ')
A1=r*s
A2=u*t
A3=r*s
x1=r/2
x2=u/2
x3=r/2
y1=s+t+s/2
y2=s+t/2
y3=s/2
Ax=A1*x1+A2*x2+A3*x3
Ay=A1*y1+A2*y2+A3*y3
A=A1+A2+A3
Xc=Ax/A
Yc=Ay/A
Ix1=(s^3*r)/12
Ix2=(t^3*u)/12
Ix3=(s^3*r)/12
Iy1=(s*r^3)/12
Iy2=(u^3*t)/12
Iy3=(s*r^3)/12
Ixx1=Ix1+A1*y1^2
Ixx2=Ix2+A2*y2^2
Ixx3=Ix3+A3*y3^2
Iyy1=Iy1+A1*x1^2
Iyy2=Iy2+A2*x2^2
Iyy3=Iy3+A3*x3^2
Ixxt=Ixx1+Ixx2+Ixx3
Iyyt=Iyy1+Iyy2+Iyy3
disp('El Momento de Inercia con respecto al eje x es:')

Ixcg=Ixxt-A*(Yc)^2
disp('El Momento de Inercia con respecto al eje y es:')
Iycg=Iyyt-A*(Xc)^2
disp('Fin del programa.')
EJERCICIO N° 5

%. Realizar un programa que permita calcular el momento de inercia de la


figura con forma de L, siendo esta simétrica.

clc
clear all
close all
disp('"Inercia para calcular una figura con forma de L, siendo esta
simétrica"')
w=input('ingrese el valor de w=');
x=input('ingrese el valor de x=');
y=input('ingrese el valor de y=');
z=input('ingrese el valor de z=');
b1=y;
h1=z;
b2=w;
h2=x;
A1=b1*h1;
A2=b2*h2;
At=A1+A2;
x1=b1/2;
y1=h1/2;
x2=b2/2;
y2=h1+(h2/2);
Ax1=A1*x1;
Ax2=a2*x2;
Ay1=a1*y1;
Ay2=a2*y2;
Axt=Ax1+Ax2;
Ayt=Ay1+Ay2;
xc=Axt/At;
yc=Ayt/At;
Ix1=(b1*(h1^3))/12;
Ix2=(b2*(h2^3))/12;
Iy1=((b1^3)*h1)/12;
Iy2=((b2^3)*h2)/12;
Ixxt=Ix1+Ix2;
Iyyt=Iy1+Iy2;
dx1=y1-yc;
dx2=y2-yc;
dy1=x1-xc;
dy2=x2-xc;
adx12=a1*(dx1^2);
adx22=a2*(dx2^2);
ady12=a1*(dy1^2);
ady22=a2*(dy2^2);
adx2=adx12+adx22;
ady2=ady12+ady22;
disp('El Momento de Inercia con respecto al eje x es:')
Ixcg=Ixxt+adx2
disp('El Momento de Inercia con respecto al eje y es:')
Iycg=Iyyt+ady2
disp('Fin del programa.')
EJERCICIO N° 6

%. Realizar un programa que permita calcular el momento de inercia


respecto a un eje de simetría ubicado en la parte inferior izquierda de
la figura con forma de Z invertida su dirección, siendo esta simétrica
clc
clear all
close all
disp('"Inercia para calcular una figura con forma de Z invertida su
dirección, siendo esta simétrica"')
a=input('ingrese el valor a: ')
b=input('ingrese el valor b: ')
c=input('ingrese el valor c: ')
d=input('ingrese el valor d: ')
e=input('ingrese el valor e: ')
f=input('ingrese el valor f: ')
A1=c*(e+d);
A2=(a+b)*d;
A3=b*(f+d);
A4=a*d;
A5=c*(d+e);
x1=1/2*c;
x2=c+1/2*(a+b);
x3=c+a+1/2*b;
x4=c+a+b+1/2*a;
x5=c+a+b+a+1/2*c;
y1=1/2*(d+e);
y2=1/2*d;
y3=d+1/2*(f+d);
y4=d+f+1/2*d;
y5=d+f;
Ax=A1*x1+A2*x2+A3*x3+A4*x4+A5*x5;
Ay=A1*y1+A2*y2+A3*y3+A4*y4+A5*y5;
A=A1+A2+A3+A4+A5;
Xc=Ax/A;
Yc=Ay/A;
Ix1=1/12*c*(e+d)^3;
Ix2=1/12*(a+b)*d^3;
Ix3=1/12*b*(f+d)^3;
Ix4=1/12*a*d^3;
Ix5=1/12*c*(d+e)^3;
Iy1=1/12*c^3*(e+d);
Iy2=1/12*(a+b)^3*d;
Iy3=1/12*b^3*(f+d);
Iy4=1/12*a^3*d;
Iy5=1/12*c^3*(d+e);
Ixxt=Ix1+Ix2+Ix3+Ix4+Ix5;
Iyyt=Iy1+Iy2+Iy3+Iy4+Iy5;
dx1=y1-Yc;
dx2=y2-Yc;
dx3=y3-Yc;
dx4=y4-Yc;
dx5=y5-Yc;
dy1=x1-Xc;
dy2=x2-Xc;
dy3=x3-Xc;
dy4=x4-Xc;
dy5=x5-Xc;
Adx12=A1*(dx1)^2;
Adx22=A2*(dx2)^2;
Adx32=A3*(dx3)^2;
Adx42=A4*(dx4)^2;
Adx52=A5*(dx5)^2;
Ady12=A1*(dy1)^2;
Ady22=A2*(dy2)^2;
Ady32=A3*(dy3)^2;
Ady42=A4*(dy4)^2;
Ady52=A5*(dy5)^2;
Adxt2=Adx12+Adx22+Adx32+Adx42+Adx52;
Adyt2=Ady12+Ady22+Ady32+Ady42+Ady52;
disp('El Momento de Inercia con respecto al eje x es:')
Ixxcg=Ixxt+Adxt2
disp('El Momento de Inercia con respecto al eje y es:')
Iyycg=Iyyt+Adyt2
disp('Fin del programa.')
EJERCICIO N° 7

%. Realizar un programa que permita calcular el momento de inercia de la


figura con forma de C (COMPUESTA) definida en los costados de la misma,
siendo esta simétrica.

clc
clear all
close all
disp('"Inercia para calcular una figura con forma de C (COMPUESTA)
definida en los costados de la misma, siendo esta simétrica"')
disp('condiciones q<g<(h,p))y((h=p)o(hno=p));para que cumpla la figura C
descrita anteriormente')
g=input('ingrese la variable g:');
h=input('ingrese la variable h:');
p=input('ingrese la variable p:');
q=input('ingrese la variable q:');
A1=p*h;
A2=(p-2*q)*(h-2*q);
A3=q*(p-2*g);
At=A1-A2-A3;
x1=h/2;
x2=((h-2*q)/2)+q;
x3=h-(q/2);
y1=p/2;
y2=((p-2*q)/2)+q;
y3=((p-2*g)/2)+g;
Ax1=A1*x1;
Ax2=A2*x2;
Ax3=A3*x3;
Ay1=A1*y1;
Ay2=A2*y2;
Ay3=A3*y3;
Axt=Ax1-Ax2-Ax3;
Ayt=Ay1-Ay2-Ay3;
Ix1cg=(h*(p^3))/12;
Ix2cg=((h-2*q)*((p-2*q)^3))/12;
Ix3cg=(q*((p-2*g)^3))/12;
Iy1cg=((h^3)*p)/12;
Iy2cg=(((h-2*q)^3)*(p-2*q))/12;
Iy3cg=((q^3)*(p-2*g))/12;
Ay21=A1*(y1^2);
Ay22=A2*(y2^2);
Ay23=A3*(y3^2);
Ax21=A1*(x1^2);
Ax22=A2*(x2^2);
Ax23=A3*(x3^2);
Ixx1o=Ix1cg+Ay21;
Ixx2o=Ix2cg+Ay22;
Ixx3o=Ix3cg+Ay23;
Iyy1o=Iy1cg+Ax21;
Iyy2o=Iy2cg+Ax22;
Iyy3o=Iy3cg+Ax23;
Ixxt=Ixx1o-Ixx2o-Ixx3o;
Iyyt=Iyy1o-Iyy2o-Iyy3o;
Xc=Axt/At;
Yc=Ayt/At;
disp('El Momento de Inercia con respecto al eje x es:')
Ixxcg=Ixxt-(At*(Yc^2))
disp('El Momento de Inercia con respecto al eje y es:')
Iyycg=Iyyt-(At*(Xc^2))
disp('Fin del programa.')

ANEXOS
 Archivo adjunto con las figuras acotadas.(PDF de cada uno).

También podría gustarte