UNIVERSIDAD DE CUENCA
ESTRUCTURAS III
MA. BELÉN YUQUILIMA
DATOS
%PORTICO 3D
l1=4; % ancho del pórtico (m)
l2=6; % ancho del pórtico (m)
h=3; % altura del pórtico (m)
%COLUMNA
cb=0.20; % base de columna (m)
ca=0.40; % altura de columna (m)
Ac=cb*ca; % area de columna (m2)
Ic=cb*ca^3/12; %inercia de columna (m4)
%VIGA
vb=0.20; % base de viga (m)
va=0.40; % altura de viga (m)
Av=vb*va; % area de viga (m2)
Iv=vb*va^3/12; %inercia de viga (m4)
%CARACTERISTICAS MECÁNICAS DE MATERIAL
E=21000;% módulo de eslasticidad [MPa]
mu=0.2;%[-]
%DENSIDAD DE HORMIGON
rho=2400/1000000;%[kg/m^3]
TIPO
% # tipo
Types={1 'beam'};%VIGA Y COLUMNA
SECCIÓN
% [# A ky kz Ixx Iyy Izz]
Sections=[1 Ac Inf Inf 0.4^3*0.2/12+0.2^3*0.4/12 0.2*0.4^3/12 0.4*0.2^3/12; %COLUMNA
2 Av Inf Inf 0.4^3*0.2/12+0.2^3*0.4/12 0.2*0.4^3/12 0.4*0.2^3/12]; %VIGA
1
MATERIAL
% # E mu rho
Materials=[1 E mu rho];
NODOS (EJES GENERALES)
% # x y z
Nodes=[1 0 0 0;
2 l1 0 0;
3 0 0 h;
4 l1 0 h];
Nodes=reprow(Nodes,(1:4),1,[l1,0,l2,0]);
%NODOS DE REFERENCIA
Nodes=[Nodes; 9 2 0 3;
10 2 6 3;
11 4 3 3;
12 0 3 3];
%GRAFICA
figure
plotnodes(Nodes);
title('GRAFICO DE NUDOS');
2
ELEMENTOS
% [# Typ Sec Mat Ni Nf Nr]
Elements=[1 1 1 1 1 3 9;
2 1 1 1 2 4 9;
3 1 1 1 5 7 10;
4 1 1 1 6 8 10;
5 1 2 1 3 4 10;
6 1 2 1 7 8 9;
7 1 2 1 3 7 11;
8 1 2 1 4 8 12];
%GRAFICA
figure
plotelem(Nodes,Elements,Types);
title('GRAFICO DE NUDOS Y ELEMENTOS');
3
%GRAFICA NODOS ELEMENTOS
plotnodes(Nodes)
hold on
plotelem(Nodes,Elements,Types)
hold off
4
GRADOS DE LIBERTAD
DOF=getdof(Elements,Types);
% .01 DESPLAZAMIENTO EN X
% .02 DESPLAZAMIENTO EN Y
% .03 DESPLAZAMIENTO EN Z
% .04 ROTACION EN X
% .05 ROTACION EN Y
% .06 ROTACION EN Z
DOF=getdof(Elements,Types);
seldof=[1.00 2.00 5.00 6.00];
DOF=removedof(DOF,seldof);
CALCULO PESO PROPIO
g=[0 0 9.81];
pesopropio=accel(g,Elements,Types,Sections,Materials);
MATRIZ DE RIGIDEZ
k=asmkm(Nodes,Elements,Types,Sections,Materials,DOF);
5
APLICACION DE CARGAS
%CARGAS DISTRIBUIDAS
% # kix kiy kiz kjx kjy kjz
Dloads=[5 0 0 -0.012 0 0 -0.012
6 0 0 -0.012 0 0 -0.012
7 0 0 -0.012 0 0 -0.012
8 0 0 -0.012 0 0 -0.012];
Dl=elemloads(Dloads,Nodes,Elements,Types,DOF);
Dl1=elemloads(pesopropio,Nodes,Elements,Types,DOF);
SOLUCION DEL SISTEMA
%U=inv(k)*T;
U=k\(Dl+Dl1); %MATRIZ DE DEFORMACIONES
F=elemforces(Nodes,Elements,Types,Sections,Materials,DOF,U,Dloads);
%Fuerzas en los elementos
GRAFICA DE CORTANTES
figure
plotforc('sheary',Nodes,Elements,Types,F,Dloads);
title('Cortante en Y [MN/m]');
6
GRAFICA DE NORMALES
figure
plotforc('norm',Nodes,Elements,Types,F,Dloads);
title('Esfuerzo Normal [MN]');
GRAFICA DE MOMENTOS
figure
plotforc('momz',Nodes,Elements,Types,F,Dloads);
title('Momentos (MN*m)')
7
GRAFICA DE DEFORMACIONES
figure
plotdisp(Nodes,Elements,Types,DOF,U);
title('GRÁFICO DE DEFORMACIONES [m]');
8
ZAPATA
qadm= 50;%q ADMISIBLE SUELO KN/M2
r1=1;%RELACION1 TAMANO ZAPATA
r2=1.5;%RELACION2 TAMANO ZAPATA
Pu=((cb*ca*3*24)+36);% KN FUERZA ULTIMA TRANSMITIDA DE LA COLUMNA AL PLINTO
%DATOS COLUMNAS
% cb= 0.20; % base de la columna (m)
%ch= 0.40; % ancho de la columna (m)
%resolucion
ladob= sqrt(Pu/(r1*qadm));%LADO b con r1
ladoa= ladob*r1;%LADO b CON r1
%ploteado de rectangulos%
fprintf('DIMENSIONES DE ZAPATA RELACION 1:1')
DIMENSIONES DE ZAPATA RELACION 1:1
disp('(relacion 1:1) lado a (metros) es:')
9
(relacion 1:1) lado a (metros) es:
ladoa
ladoa = 0.9139
disp('(relacion 1:1) lado b (metros) es:')
(relacion 1:1) lado b (metros) es:
ladob
ladob = 0.9139
a=ladoa;
lh=ladob;
fprintf('AREA TOTAL ZAPATA (m2) = %5f\n',a*lh)
AREA TOTAL ZAPATA (m2) = 0.835200
x=linspace(0,a,500);y=linspace(0,lh,500);
y1=0;y2=lh;x1=0;x2=a;
plot(x,y1,'b.',x,y2,'b.',x1,y,'b.',x2,y,'b.');
10
