UNIVERSIDAD NACIONAL DE SAN

CRISTOBAL DE HUAMANGA
FACULTAD DE INGENIERIA DE MINAS
GEOLOGIA Y CIVIL
ESCUELA DE FORMACION PROFESIONAL DE
INGENIERIA CIVIL

MATRIZ DE RIGIDEZ DE UN PORTICO

CURSO:
INTRODUCCION AL METODO DE LOS
ELEMENTOS FINITOS
DOCENTE:
ING. CRISTIAN CASTRO PEREZ
ALUMNO:
VICTOR HUGO DIAZ VIVANCO
CODIGO:
16080538

AYACUCHO - PERU
31 de octubre de 2012
Indice
1. INTRODUCCION 3
1.1. Propositos . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3
1.2. Justificacion e importancia . . . . . . . . . . . . . . . . . . . . . . . 3
1.3. Resumen . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3

2. FUNDAMENTACION TEORICA 4

3. METODOLOGIA 4
3.1. Metodo de rigidez . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4
3.1.1. Matriz de rigidez de un elemento en coordenadas globales . . 6
3.1.2. Matriz de rigidez de un elemento de seccion constante o variable 7
3.2. Definicion del problema . . . . . . . . . . . . . . . . . . . . . . . . . 8
3.2.1. Justificacion del problema . . . . . . . . . . . . . . . . . . . . 8
3.2.2. Formulacion del problema . . . . . . . . . . . . . . . . . . . . 9

4. RESULTADOS 11
4.1. En MATLAB . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 11
4.2. En SAP2000 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 14
4.3. Comparacion . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 17
4.4. Codigos del programa en Matlab . . . . . . . . . . . . . . . . . . . . 19

5. CONCLUSIONES 27

6. REFERENCIAS BIBLIOGRAFICAS 28

7. ANEXOS 28

2
1. INTRODUCCION
1.1. Propositos
El presente trabajo se hizo con el proposito de hacer un analisis matricial con
metodo de la rigidez de una determinada estructura. Para poder hallar la matriz de
cada elemento de la estructura se sigue una serie de pasos como se detalla mas ade-
lante, posteriormente se ensambla la matriz global de toda la estructura y despues
se calculan las fuerzas y momentos que actuan en cada extremo de las barras. Final-
mente estos resultados que obtenemos lo compararemos con un programa comercial
que es el SAP2000, como veremos mas adelante estos resultados deben asemejarse.

Para poder seguir todos los pasos mencionados lneas arriba nos apoyaremos en
el programa MATLAB y para la comparacion de resultados el programa SAP2000.

1.2. Justificacion e importancia

Es de suma importancia analizar estructuras complejas como son porticos en 3
o 2 dimensiones.Y cuando son estructuras complejas se hace mas tedioso este analisis.

La matriz de rigidez es parte del metodo matricial de la rigidez que esta di-
senado para realizar analisis computarizado de cualquier estructura incluyendo a
estructuras estaticamente indeterminadas. El metodo matricial se basa en estimar
los componentes de las relaciones de rigidez para resolver las fuerzas o los desplaza-
mientos mediante un ordenador.

El empleo de la notacion matricial presenta dos ventajas en el calculo de estruc-

turas. Desde el punto de vista teorico , permite utilizar metodos de calculo en forma
compacta, precisa y, al mismo tiempo, completamente general. Esto facilita el trata-
miento de la teora de estructuras como unidad, sin que los principios fundamentales
se vean oscurecidos por operaciones de calculo.

Desde el punto de vista practico, proporciona un sistema apropiado de analisis de

estructuras y determina una base muy conveniente para el desarrollo de programas
de computacion.

1.3. Resumen
En la introduccion, se hace la respectiva justificacion e importancia del metodo
de la rigidez.

En la fundamentacion teorica se explica como evolucionaron los metodos de

analisis de estructuras a lo largo de la historia, posteriormente debido al avance
tecnologico surge en analisis matricial de estructuras.

En la parte de metodologa se explica teoricamente en que consiste el metodo de

rigidez, posteriormente se da una serie de pasos a seguir para poder hallar la matriz
de rigidez de cada elemento que tiene la estructura.

En los resultados se muestra la interfaz grafica del MATLAB y secuencialmente

nos muestra una ventana donde se llena los datos para despues darnos los resultados

3
de las fuerzas y momentos que actuan en las barras en coordenadas globales.

Para finalizar mostramos las conclusiones a que se llegaron, tambien mostramos

las referencias bibliograficas y los anexos.

2. FUNDAMENTACION TEORICA
Antiguamente y debido al limitado desarrollo tecnologico, el analisis de estruc-
turas se hacia con metodos aproximados como son el Metodo de Hardy Croos, pero
la mayora de veces eran solo aplicables a determinados tipos de estructuras.

La principal objecion a los primeros metodos de analisis fue que los mismos con-
ducan a sistemas con un gran numero de ecuaciones lineales, difciles de resolver
manualmente. Con los computadores, capaces de realizar el trabajo numerico, esta
objecion no tiene ahora sentido. Debido al avance tecnologico los calculos con compu-
tadoras son mas faciles y ya no debemos preocuparnos por las operaciones rutinarias.

Es de suma importancia hacer un modelo adecuado y una debida interpretacion

de los resultados.

3. METODOLOGIA
3.1. Metodo de rigidez
Hipotesis: Estructura lineal - Todos los movimientos y esfuerzos son funciones
lineales de las cargas- Pequenas deformaciones (ecuaciones de equilibrio en la estruc-
tura no distorsionada).

Para estudiar una estructura por el metodo de la rigidez, al igual que en cual-
quier otro problema elastico, disponemos de tres conjuntos de ecuaciones que deben
cumplirse.

Ecuaciones de compatibilidad

Ecuaciones de compatibilidad

Ecuaciones de equilibrio

Las ecuaciones de compatibilidad relacionan las deformaciones de barras con los

desplazamientos nodales. Introduciendo estas relaciones en las ecuaciones constitu-
tivas, relacionamos las fuerzas en los extremos de barras con los desplazamientos
nodales.

Introduciendo estas ultimas relaciones en las ecuaciones de equilibrio se obtiene

un conjunto de ecuaciones de fuerzas nodales en funcion de desplazamientos noda-
les, que pueden ser consideradas como Ecuaciones de Equilibrio de la estructura en
funcion de desplazamientos.

4
 La resolucion de este sistema de ecuaciones nos permite obtener el valor de las
incognitas (desplazamientos nodales), a partir de los cuales se obtienen las solicita-
ciones de las barras de la estructura, as como las reacciones.

Una de las caractersticas mas importantes del metodo de la rigidez es la forma

en que las propiedades elasticas de las piezas, y su orientacion dentro de la estruc-
tura, son introducidas en el calculo antes de que se efectue ninguna consideracion
sobre el equilibrio o la compatibilidad de los nudos.

Esto nos permite establecer relaciones entre las fuerzas de extremo de barras
y los desplazamientos de nudo. Estas relaciones expresadas en forma matricial se
denomina matriz de rigidez de barra.

Al considerar la interrelacion de cada barra con las demas se obtiene un sistema

global de ecuaciones que define el comportamiento de toda la estructura y nos con-
duce a la solucion del problema.

Podemos considerar seis etapas fundamentales en la solucion de un problema:

1. Identificacion estructural.

2. Calculo de la matriz de rigidez de barra y del vector de cargas nodales equi-

valentes.

3. Calculo de la matriz de rigidez global y del vector de cargas global de la

estructura.

4. Introduccion de las condiciones de borde

5. Solucion del sistema de ecuaciones

6. Calculo de solicitaciones en los extremos de barras y reacciones nodales.

Finalmente se utilizo como libro gua y de consulta Analisis Matricial de Estruc-

turas de Roberto Aguiar Falconi. Para poder hallar la matriz de rigidez y las fuerzas
y momentos que actuan en cada barra de la estructuras, que en este caso se hizo de
un portico de 2 tramos se siguieron los siguientes pasos:

Se define la cantidad de grados de libertad, se numera esos grados. Se cuenta

el numero de nodos y barras.

Una vez numerado los grados de libertad se procede hallar el vector de colo-
cacion.

Se halla la matriz de rigidez de cada elemento en coordenadas globales.

Se utiliza el vector de colocacion para poder ensamblar la matriz de rigidez

global.

Se pide como dato el vector de cargas generalizadas que actuan en cada nudo
de la estructura.

Para poder hallar los giros o desplazamientos de la estructuras se multiplica

la inversa de la matriz de rigidez global con el vector de cargas generalizadas.

5
 Para encontrar las deformaciones de las barras en coordenadas globales, se
trabaja con el vector de colocacion.

Se obtienen las fuerzas y momentos en los extremos de las barras en coordena-

das globales, multiplicando la matriz de rigidez del miembro en coordenadas
globales por su vector respectivo de deformaciones hallado en el item anterior.

3.1.1. Matriz de rigidez de un elemento en coordenadas globales

Las coordenadas globales se miden en forma horizontal y vertical, en consecuen-
cia estas ultimas tienen la misma orientacion que las coordenadas de la estructura.

En la figura se tiene un elemento inclinado que forma un angulo con la hori-

zontal, tambien en la figura se muestra el sistema de coordenadas globales.

Figura 1: Coordenadas globales de un elemento

Para simplificar la escritura se denomina C = cos, S = sen. La siguientes

ecuaciones nos dan la matriz de rigidez de un elemento en coordenadas globales.

EA 2 12EI 2
a= C + S (1)
L L3
EA 12EI
 
b= SC (2)
L L3
6EI
c= S (3)
L2
EA 2 12EI 2
d= S + C (4)
L L3
4EI
f= (5)
L
6EI
g= C (6)
L2
Matricialmente se vera as:

6
 a


b d

c g f
k3 =

a b c a


b a g b

d
f
c g 2 c g f

Donde:
L, A, I son las magnitudes geometricas (longitud, area y momento de inercia).
E la constante de elasticidad longitudinal (modulo de Young).

Se utilizo esta matriz para hallar la matriz de rigidez de una columna en

coordenadas globales, que son de seccion constante. No se considero el efecto de
corte.

3.1.2. Matriz de rigidez de un elemento de seccion constante o variable

Las vigas en nuestro problema tienen las mismas coordenadas globales y locales,
en la figura se muestra estas coordenadas (sistema de coordenadas globales).

Figura 2: Coordenadas globales de un elemento

Para hallar la matriz de rigidez en coordenadas globales de las vigas

acarteladas que son de seccion variable se utilizo las siguientes formulas:

ZL ZL  2
(L X)2 dx 1 dx
= f11 = + (7)
L2 EI(x) L GA(x)
0 0

ZL ZL  2
X (L X) dx 1 dx
= f21 = f12 = + (8)
L2 EI(x) L GA(x)
0 0

ZL ZL  2
0 X 2 dx 1 dx
= f22 = + (9)
L2 EI(x) L GA(x)
0 0

ZL
dx
f33 = (10)
EA(x)
0

0
k= (11)
0 2
7

= (12)
0 2

k0 = (13)
0 2
Matricialmente es as:

r

0 t

0 b k
k2 =

r 0 0 r



0 t b 0 t
0 b0 a 0 b0 k 0

Donde:

k+a
b= (14)
L

k0 + a
b0 = (15)
L

b + b0
t= (16)
L

La matriz que acabamos de ver estan es coordenadas globales, operando todas

las integrales mencionadas se obtendra el valor de la matriz de rigidez para cada
elemento (vigas acarteladas).

Donde E es el modulo de elasticidad del material; G el modulo de corte estatico;

I(x) es el momento de inercia de una seccion variable; A(x) es el area de una seccion
transversal y es el coeficiente de forma.

El coeficiente de forma es adimensional que caracteriza la distribucion de las

tensiones tangenciales en la seccion transversal y depende de la forma de la misma.
Para nuestro caso las vigas son de seccion rectangular y el valor de = 1,2.

3.2. Definicion del problema

3.2.1. Justificacion del problema
Este trabajo se desarrollo para poder obtener las fuerzas y momentos en los ex-
tremos de cada una de las barras, estos resultados se compararan con un software
de estructuras como es el SAP.

Comparados ambos resultados, estos no deben alejarse debido a que se utilizo el

metodo de rigidez. Este metodo es muy confiable ya que nos lleva a muy buenos y
exactos resultados.

8
 Figura 3: Grados de libertad del portico

3.2.2. Formulacion del problema

El problema que nos planteamos es el siguiente: la estructura tiene 9 grados de
libertad, 5 barras: 3 columnas y 2 vigas acarteladas como se puede ver la figura 3.

Las columnas son de 45x45cm y tienen una altura de H = 3m. Las vigas tiene
3 tramos como muestra la figura 4.

Figura 4: Viga de 3 tramos (izquierda a derecha)

El area y la inercia de cada tramo vara debido a que la altura h(x) aumenta o
disminuye para cada tramo, las siguientes expresiones nos muestran estas variacion:

(h1 h)
h(x)1 = h1 x
L1
A(x)1 = b(h(x)1 )
b(h(x)1 )3
I(x)1 =
12

h(x)2 = h
9
 A(x)2 = b(h(x)2 )

b(h(x)2 )3
I(x)2 =
12
(L x)(h2 h)
h(x)3 = h2
L3
A(x)3 = b(h(x)3 )

b(h(x)3 )3
I(x)3 =
12

Para hallar la matriz de rigidez global se utiliza el vector de colocacion. El vector

de colocacion de cada barra se muestra a continuacion, los numeros entre parentesis
indican el numero de barra:
h i
V C(1) = 0 0 0 1 2 3

h i
V C(2) = 1 2 3 4 5 6

h i
V C(3) = 0 0 0 4 5 6

h i
V C(4) = 4 5 6 7 8 9

h i
V C(5) = 0 0 0 7 8 9

A continuacion en la figura 5 se muestra las cargas que actuan en la estructura,

tambien se muestra el vector de cargas generalizadas.

Figura 5: Cargas en la estructura

10

100

0

0


0



Q=
400


0

300



0
0

Con todo estos datos se obtiene las fuerzas y momentos en los extremos de las
barras en coordenadas globales.

4. RESULTADOS
4.1. En MATLAB
Se mostraran los resultados que nos dio el programa Matlab, se explicara como
corre el programa y cuando se llena datos.

Entrando al CD se abre la carpeta ELEMENTOS FINITOS, abrimos la carpeta

que se encuentra con el nombre de anexo1 Programa y nos mostrara esto:

11
 Figura 6: Inicio del programa

Figura 7: Modelado de la estructura

12
Figura 8: Llenar todos los datos de la estructura

Figura 9: Resultados que nos muestra el programa

13
 Al observar toda la secuencia del programa, al final nos muestra los resultados
que nos interesa y estas son las fuerzas y momentos que actuan en cada barra del
portico. Estos resultados los compararemos con lo que nos da el SAP2000.

4.2. En SAP2000
Se muestra el modelamiento que se hizo en el programa SAP. Tambien se mues-
tra las deformada, el diagrama de momentos flectores, axiales y cortantes.

Se vera mas adelante una tabla donde se comparara los resultados obtenidos en
el SAP con los obtenidos en el MATLAB.

14
Figura 10: Modelado y deformada de la estructura

Figura 11: Diagrama de Fuerza Axial

15
 Figura 12: Diagrama de Fuerza Cortante

Figura 13: Diagrama de Momentos Flectores

16
4.3. Comparacion
A continuacion se muestra los resultados en dos tablas, teniendo en cuenta las
coordenadas locales como se muestra en la figura 14.

Figura 14: Coordenadas locales

17
 Axial 1 Cort. 1 Mom. 1 Axial 2 Cort. 2 Mom. 2
Elemento 1 43.81 118.55 261.75 43.81 118.55 93.91
Elemento 2 18.55 44.89 93.91 18.55 44.89 83.52
Elemento 3 399.77 159.36 301.75 399.77 159.86 176.33
Elemento 4 177.91 47.65 92.81 177.91 47.65 95.68
Elemento 5 49.11 122.09 270.58 49.11 122.09 95.68

Cuadro 1: Resultados del programa SAP2000

Axial 1 Cort. 1 Mom. 1 Axial 2 Cort. 2 Mom. 2

Elemento 1 -41.99 -119.26 268.63 41.99 119.26 89.15
Elemento 2 -19.25 -41.98 -89.15 19.25 41.99 -78.82
Elemento 3 397.49 -158.33 308.23 -397.49 158.33 166.76
Elemento 4 -177.59 -44.49 -87.95 177.59 44.49 -90.03
Elemento 5 44.49 -122.41 277.20 -44.49 122.41 90.03

Cuadro 2: Resultados del programa Matlab

Al ver estos resultados que se asemejan bastante en los valores obtenidos tanto
en SAP2000 como en el programa Matlab , se concluye que el metodo de la rigidez
es casi exacto y muy confiable.

18
4.4. Codigos del programa en Matlab
% SE RESOLVERA PARA UN CASO PARTICULAR, UNA PORTICO DE 3 TRAMOS
% Y LA VIGA ACARTELADA
% ES UNA VIGA DE 6 GRADOS DE LIBERTAD

%clc, clear all

global h1 h h2 L1 L2 L3 E q1 q5 q7 b BB H

syms x
%datos columnas

A=b^2;
I=b^4/12;

% altura de la seccion constante (centro) viga acart

L=L1+L2+L3; %LONGITUD HORIZONTAL

%ALTURA

n=6; %numero de nudos

nel=5; %numero d barras
gl=9;

CG=zeros(n,3);
CG(3,:)=[1 2 3];
CG(4,:)=[4 5 6]; % MAtriz de Coordenadas generalizadas

VC= [0 0 0 1 2 3; 1 2 3 4 5 6; 0 0 0 4 5 6;4 5 6 7 8 9; 0 0 0 7 8 9];

% vector de colocacion

S1=sin(pi/2);
C1=cos(pi/2);

% DATOS D LA VIGA ACARTELADA

%BB=0.4; %ancho de la viga acartelada

v=0.2;
G=E/(2*(1+v));

for i=1:3

if i==1 %Elemento 1
L=H

k1 =[ E*A*C1^2/L+12*E*I*S1^2/L^3 (E*A/L-12*E*I/L^3)*S1*C1 -6*E*I*S1/L^2

19
-E*A*C1^2/L-12*E*I*S1^2/L^3 -(E*A/L-12*E*I/L^3)*S1*C1 -6*E*I*S1/L^2;
(E*A/L-12*E*I/L^3)*S1*C1 E*A*S1^2/L+12*E*I*C1^2/L^3 6*E*I*C1/L
-(E*A/L-12*E*I/L^3)*S1*C1 -E*A*S1^2/L-12*E*I*C1^2/L^3 6*E*I*C1/L^2;
-6*E*I*S1/L^2 6*E*I*C1/L^2 4*E*I/L 6*E*I*S1/L^2 -6*E*I*C1/L^2 2*E*I/L;
-E*A*C1^2/L-12*E*I*S1^2/L^3 -(E*A/L-12*E*I/L^3)*S1*C1 6*E*I*S1/L^2
E*A*C1^2/L+12*E*I*S1^2/L^3 (E*A/L-12*E*I/L^3)*S1*C1 6*E*I*S1/L^2;
-(E*A/L-12*E*I/L^3)*S1*C1 -E*A*S1^2/L-12*E*I*C1^2/L^3 -6*E*I*C1/L^2
(E*A/L-12*E*I/L^3)*S1*C1 E*A*S1^2/L+12*E*I*C1^2/L^3 -6*E*I*C1/L^2 ;
-6*E*I*S1/L^2 6*E*I*C1/L^2 2*E*I/L 6*E*I*S1/L^2 -6*E*I*C1/L^2 4*E*I/L];
K1 =[ E*A*C1^2/L+12*E*I*S1^2/L^3 (E*A/L-12*E*I/L^3)*S1*C1 -6*E*I*S1/L^2
-E*A*C1^2/L-12*E*I*S1^2/L^3 -(E*A/L-12*E*I/L^3)*S1*C1 -6*E*I*S1/L^2;
(E*A/L-12*E*I/L^3)*S1*C1 E*A*S1^2/L+12*E*I*C1^2/L^3 6*E*I*C1/L
-(E*A/L-12*E*I/L^3)*S1*C1 -E*A*S1^2/L-12*E*I*C1^2/L^3 6*E*I*C1/L^2;
-6*E*I*S1/L^2 6*E*I*C1/L^2 4*E*I/L 6*E*I*S1/L^2 -6*E*I*C1/L^2 2*E*I/L;
-E*A*C1^2/L-12*E*I*S1^2/L^3 -(E*A/L-12*E*I/L^3)*S1*C1 6*E*I*S1/L^2
E*A*C1^2/L+12*E*I*S1^2/L^3 (E*A/L-12*E*I/L^3)*S1*C1 6*E*I*S1/L^2;
-(E*A/L-12*E*I/L^3)*S1*C1 -E*A*S1^2/L-12*E*I*C1^2/L^3 -6*E*I*C1/L^2
(E*A/L-12*E*I/L^3)*S1*C1 E*A*S1^2/L+12*E*I*C1^2/L^3 -6*E*I*C1/L^2 ;
-6*E*I*S1/L^2 6*E*I*C1/L^2 2*E*I/L 6*E*I*S1/L^2 -6*E*I*C1/L^2 4*E*I/L]

k3=k1;
K3=K1

k5=k1;
K5=K1

end

if i==2 %Elemento 2 y 4

L=L1+L2+L3;

hx1=h1-x*(h1-h)/L1;
Ax1=BB*hx1;
Ix1=BB*hx1^3/12;

hx2=h;
Ax2=BB*hx2;
Ix2=BB*hx2^3/12;

hx3=h2-(L-x)*(h2-h)/L3;
Ax3=BB*hx3;
Ix3=BB*hx3^3/12;

% Hallamos los elemntos de la matriz de rigidez, que esta

en funcion de integrales

a= int((L-x)^2/(L^2*E*Ix1)+1.2/(L^2*G*Ax1),0,L1)+int((L-x)^2/(L^2*E*Ix2)
+1.2/(L^2*G*Ax2),L1,L1+L2)+int((L-x)^2/(L^2*E*Ix3)+1.2/(L^2*G*Ax3)

20
,L1+L2,L1+L2+L3);
e=int(x*(L-x)/(L^4*E*Ix1)+1.2/(L^2*G*Ax1),0,L1)+int(x*(L-x)/(L^2*E*Ix2)
+1.2/(L^2*G*Ax2),L1,L1+L2)+int(x*(L-x)/(L^2*E*Ix3)+1.2/(L^2*G*Ax3),
L1+L2,L1+L2+L3);
ap=int(x^2/(L^2*E*Ix1)+1.2/(L^2*G*Ax1),0,L1)+int(x^2/(L^2*E*Ix2)
+1.2/(L^2*G*Ax2),L1,L1+L2)+int(x^2/(L^2*E*Ix3)+1.2/(L^2*G*Ax3),
L1+L2,L1+L2+L3);
f33=int(1/(E*Ax1),0,L1)+L2/(E*Ax2)+int(1/(E*Ax3),L1+L2,L1+L2+L3);

aa=vpa(a);
ee=vpa(e);
aap=vpa(ap);
ff33=vpa(f33);

%Matriz de rigidez local para la viga axartelada

k=aap/(aa*aap-ee^2);
aaaa=ee/(aa*aap-ee^2);
kp=aa/(aa*aap-ee^2);
r=1/ff33;

bbb=(k+aaaa)/L;
bbbp=(kp+aaaa)/L;
t=(bbb+bbbp)/L;

k2=[r 0 0 -r 0 0; 0 t bbb 0 -t bbbp; 0 bbb k 0 -bbb aaaa; -r 0 0 r 0 0;

0 -t -bbb 0 t -bbbp; 0 bbbp aaaa 0 -bbbp kp];
K2=[r 0 0 -r 0 0; 0 t bbb 0 -t bbbp; 0 bbb k 0 -bbb aaaa; -r 0 0 r 0 0;
0 -t -bbb 0 t -bbbp; 0 bbbp aaaa 0 -bbbp kp]

k4=k2;
K4=K2

end
end

% Ensamblando la matriz de rigidez

for p=1:5

% Elemento 1
if p==1

for i=1:nel
for j=1:6
if VC(p,j)==0

k1(i,j)=0;
for y=1:5
y=y+1;
k1(y,j)=0;

21
 k1(j,y)=0 ;

end
end
end
end
for i=1:3
k1(1,:)=[];
k1(:,1)=[];
end

g=zeros(9);
for i=1:3
for j=1:3
g(i,j)=k1(i,j);

end
end
k1=g
end

% Elemento 2
if p==2
for i=1:nel
for j=1:6
if VC(p,j)==0

k2(i,j)=0;
for y=1:5
y=y+1;
k2(y,j)=0;
k2(j,y)=0 ;
end
end
end
end
k2;
gg=zeros(9);
for i=1:6
for j=1:6
gg(i,j)=k2(i,j);
end
end
k2=gg
end

% Elemento 3
if p==3
for i=1:nel
for j=1:6
if VC(p,j)==0

22
 k3(i,j)=0;
for y=1:5
y=y+1;
k3(y,j)=0;
k3(j,y)=0
end
end
end
end

ggg=zeros(9);
for i=1:6
for j=1:6
ggg(i,j)=k3(i,j);
end
end
k3=ggg
end

% Elemento 4
if p==4
for i=1:nel
for j=1:6
if VC(p,j)==0

k4(i,j)=0;
for y=1:5
y=y+1;
k4(y,j)=0;
k4(j,y)=0;
end
end
end
k4

gggg=zeros(9);
for i=1:6
for j=1:6
gggg(i+3,j+3)=k4(i,j);
end
end
end
k4=gggg

end

% Elemento 5
if p==5
for i=1:nel
for j=1:6

23
 if VC(p,j)==0

k5(i,j)=0;
for y=1:5
y=y+1;
k5(y,j)=0;
k5(j,y)=0 ;
end
end
end
end
k5
gf=zeros(9);
for i=1:6
for j=1:6
gf(i+3,j+3)=k5(i,j);
end
end
k5=gf
end
end

% Rigidez encoordenada globales, sumamos...

format short
K=k1+k2+k3+k4+k5
%K= (k1)+vpa(k2)+vpa(k3)+vpa(k4)+vpa(k5)

% vector de cargas generalizadas en los nudos 1 5 7

Q=[q1; 0 ; 0; 0; -q5; 0;q7; 0; 0]

% desplazamientos y giros q=K-1*Q

format short
q=inv(K)*Q

H ;
L=L1+L2+L3;

% Hallaremos el vector P de desplazamientos para cada miembro y

% multiplicamos por la matriz de rigidez global de cada elemento
% obteniedo las fuerzas y momentos de los miembros ebn coordenada globales

clear P1 P2 P3 P4 P5
for i=1
for j=1:6
if VC(i,j)~=0
P1(j)=q(VC(i,j));
else
P1(j)=0;

24
 end
end
P1=P1

disp(Fuerza y momentos en los extremos de la barra 1:)

PF1=K1*P1
end

for i=2
for j=1:6
if VC(i,j)~=0
P2(j)=q(VC(i,j));
else
P2(j)=0;
end
end
P2=P2
PF2=K2*P2

for j=1:6

Pf2(j)=PF2(j)
end
disp(Fuerza y momentos en los extremos de la barra 2:)
PF2=Pf2

end

for i=3
for j=1:6
if VC(i,j)~=0
P3(j)=q(VC(i,j));
else
P3(j)=0;
end
end
P3=P3

disp(Fuerza y momentos en los extremos de la barra 3:)

PF3=K3*P3
end

for i=4
for j=1:6
if VC(i,j)~=0
P4(j)=q(VC(i,j));
else
P4(j)=0;
end
end
P4=P4

25
disp(Fuerza y momentos en los extremos de la barra 4:)
PF4=K4*P4
end

for i=5
for j=1:6
if VC(i,j)~=0
P5(j)=q(VC(i,j));
else
P5(j)=0;
end
end
P5=P5
disp(Fuerza y momentos en los extremos de la barra 5:)
PF5=K5*P5
end

26
5. CONCLUSIONES
Al desarrollar el problema en MATLAB existe un inconveniente que es el
tiempo de calculo, pues demora casi 10 minutos, una alternativa para que
corra mas rapido el programa es recurrir a los metodos numericos y utilizar
algun algoritmo para la integracion como pued ser Gauss o Simpson.

Al observar los resultados tanto de MATLAB como de SAP 2000, se concluye

que estos son casi iguales y por ende el metodo de la rigidez es muy confiable
ya que nos lleva a resultados casi exactos.

La pequena diferencia que existe en los resultados se debe a que no se consi-

dero el efecto de corte para las columnas, pero si para las vigas.

Viendo los resultados que nos da ambos programa (cuadro 1 y 2) , la variacion

de los resultados esta en el margen permisible.

Es necesaria como herramienta fundamental la computadora, ya que sin esta hu-

biera sido muy tedioso el calculo de la estructura por el metodo de la rigidez.

Al introducir datos al programa se puede dar como dato de cualquier magnitud

de carga o momento que actua en cualquier grado de libertad de la estructura.

Es necesario hacer un buen modelado de una estructura para que de esta mane-
ra cualquier programa nos de buenos resultados y a estos saber interpretarlos.

27
6. REFERENCIAS BIBLIOGRAFICAS
1. AGUIAR FALCONI, Roberto - Analisis matricial de estructuras

2. TENA COLUNGA, Arturo - Analisis de estructuras con metodos matriciales

3. http : //es.scribd.com/doc/53788480/Breve Historia Del Analisis

Estructural

4. http : //www.elprisma.com/apuntes/curso.asp?id = 13685

5. http : //es.wikipedia.org/wiki/M atrizd er igidez

6. http : //ing.unne.edu.ar/pub/e3c ap4.pdf

7. Manuales de Matlab

7. ANEXOS
Se adjunta en un CD al presente informe, donde esta contenido el informe hecho
integramente en LATEXy el programa que resuelve un portico de 9 grados de
libertad que esta en Matlab.

Se adjunta 3 carpetas: Anexo1 Programa, Anexo2 Informe y Anexo3 ModenSAP2000,

revisando estas carpetas se encontraran el programa hecho en MATLAB, el
informe en LATEXy el modelamiento de la estructura que se hizo en SAP2000.

28