PRACTICA 3A - Métodos Matriciales - Tutorial

Teoría de Estructuras. Practicas sobre el Tema 4: Métodos matriciales.
Práctica 3A - Pórtico a
dos aguas con cargas nodales.
PRÁCTICA 3A
CÁLCULO MATRICIAL DE UN PÓRTICO A DOS

AGUAS BIEMPOTRADO CON CARGAS NODALES
Esta práctica consiste en el cálculo de un pórtico biempotrado con cargas

nodales mediante el método matricial de la rigidez en Excel. Este programa
permite observar cómo se desarrollan todos los pasos necesarios hasta
llegar a obtener la distribución de esfuerzos en las barras del pórtico.
En la siguiente figura se representan las dimensiones principales del pórtico,

que es a dos aguas, simétrico, de 25 m de luz, 5 m de altura de pilares y
10% de inclinación de cubierta, con la que la altura de la cumbrera es 6,25
m.
6,25
5m m
25 m
Figura A1.1: Pórtico biempotrado a 2 aguas de nudos rígidos.
En la siguiente figura se representan las cargas nodales que actúan sobre el

pórtico:
20000 N
40000 N
30000 N
10000 Nm
50000 N
Figura A1.2: Cargas aplicadas en los nudos.
Departamento de Ingeniería Mecánica, Energética y de Materiales. 1

Universidad Pública de Navarra. Campus de Tudela.
Teoría de Estructuras. Practicas sobre el Tema 4: Métodos matriciales. Práctica 3A - Pórtico a
El método de la rigidez requiere realizar un predimensionado de los perfiles
que componen la estructura, en este caso, el pórtico. En la tabla siguiente se
indican los perfiles seleccionados como predimensionamiento:
NUDO NUDO AREA M. DE INERCIA

BARRA PERFIL
INICIAL FINAL (m2) (m4)
1 1 2 HEB 280 0,01314 1,927E-4
2 2 3 IPE 450 0,00988 3,374E-4
3 3 4 IPE 450 0,00988 3,374E-4
4 4 5 HEB 280 0,01314 1,927E-4
Tabla A1.1: Tabla de asignación de perfiles.
En la siguiente figura se representa de nuevo el pórtico del ejemplo,

mostrándose los ejes locales de cada barra y los ejes globales de la
3
estructura. También se indica la numeración de nudos y barras propuesta.
2 3
2 4
1 4
Figura A1.3: Ejes globales, ejes locales, numeración de nudos y barras del
pórtico biempotrado a 2 aguas de nudos rígidos.
Con la información anterior se realiza la siguiente tabla, en la cual se

1
determinan 5 que
las longitudes de las barras del pórtico, así como el ángulo
forman los ejes globales con los ejes locales de cada barra.
LONGITUD ÁNGULO
BARRA
Tabla A1.2: (m) (grados) Datos
geométricos de 1 5 90 las barras del
pórtico.
2 12,5623 5,711
En la siguiente 3 12,5623 354,289 figura se
recuerda la representación
de los ejes 4 5 270 globales de la
estructura y locales de cada barra:

Figura A1.4: Ejes globales del pórtico XG, YG, ZG y ejes locales de una barra
XL, YL, ZL, en este caso de la barra 2.
En la siguiente figura se acotan los ángulos  de cada barra: a partir del eje
x local de la barra se gira hasta buscar el eje X global del pórtico.
Figura A1.5: Ejes globales, ejes locales, numeración de nudos y barras del
pórtico biempotrado a 2 aguas de nudos rígidos.
Tras exponer la geometría del pórtico, las cargas actuantes en los nudos y el
predimnesionado del mismo se describirá de forma abreviada la realización
de la práctica paso por paso.

1. Abrir el archivo de EXCEL “PRÁCTICA 3A.xls”. Observar que consta de
10 hojas en las que se desarrollan los cálculos necesarios. En las 4 primeras
hojas se realizan los cálculos de las matrices de rigidez en locales y en
globales de cada barra. Se empezará introduciendo los datos geométricos
de las barras.
2. En la hoja “K(L) y K (G) barra 1” introducir los datos geométricos y

mecánicos de la barra 1 en las celdas de fondo verde B3, C3, D3, E3 y G3,
tal como se reproduce en la tabla siguiente. Las unidades utilizadas deben
ser metros y Newton o unidades derivadas.
CÁLCULO DE LAS MATRICES DE RIGIDEZ EN LOCALES Y GLOBALES BARRA 1

E (N/m2) A (m2) L (m) Iy (m4) ALPHA (RAD) ALPHA (GRD)
205939650000 0,01314 5 1,927E-04 1,570796 90
Tabla A1.3: Introducción de datos geométricos y mecánicos de la barra 1.
NOTA: No modificar la celda en rojo. En la celda F3 con fondo rojo se calcula

el ángulo  de la directriz de la barra en radianes. Excel realiza los cálculos
trigonométricos en radianes.
En esta hoja “K(L) y K (G) barra 1” se realiza primero el cálculo de la

matriz de rigidez en coordenadas locales KL, en concreto en las celdas
B6:G11, como se reproduce en la siguiente tabla:
MATRIZ DE RIGIDEZ EN LOCALES DE LA BARRA 1

541209400 0 0 -541209400 0 0
0 3809719 9524297 0 -3809719 9524297
0 9524297 31747656 0 -9524297 15873828
-541209400 0 0 541209400 0 0
0 -3809719 -9524297 0 3809719 -9524297
0 9524297 15873828 0 -9524297 31747656
Tabla A1.4: Resultados del cálculo de la matriz de rigidez en locales de la
barra 1.
Estos cálculos se realizan a partir de los datos geométricos y mecánicos de

la barra mediante la expresión siguiente:

Un poco más abajo, en las celdas B37:G42 se calcula la matriz de rigidez de

esta barra en coordenadas globales KG, previo cálculo de la matriz de
transformación de coordenadas L y su traspuesta LT.
La ecuación matricial para el cálculo de la matriz de rigidez en globales de

una barra es:
La matriz de transformación de coordenadas locales a globales L se

determina a partir de la siguiente expresión en función del ángulo  de la
barra:
La matriz de transformación de coordenadas globales a locales LT se obtiene

con la siguiente expresión en función del ángulo  de la barra:
Una vez recordadas las expresiones anteriores se llega a la matriz de rigidez

en globales de la barra 1:

MATRIZ DE RIGIDEZ EN GLOBALES DE LA BARRA 1
3809719 0 -9524297 -3809719 0 -9524297
0 541209400 0 0 -541209400 0
-9524297 0 31747656 9524297 0 15873828
-3809719 0 9524297 3809719 0 9524297
0 -541209400 0 0 541209400 0
-9524297 0 15873828 9524297 0 31747656
Tabla A1.5: Resultados del cálculo de la matriz de rigidez en globales de la
barra 1.
3. Ir a la hoja “K(L) y K (G) barra 2” e introducir los datos geométricos y

tal como se indica en la siguiente tabla:

205939650000 9,88E-03 12,5623 3,374E-04 0,099676 5,711
Al igual que en el caso de la barra 1 se pueden revisar los resultados del

cálculo de la matriz de rigidez en locales y globales de la barra 2.
4. En la hoja “K(L) y K (G) barra 3” introducir los datos geométricos y

mecánicos de la barra 3 en las celdas de fondo verde B3, C3, D3, E3 y G3:
en la siguiente tabla se reproducen las entradas:

205939650000 9,88E-03 12,5623 3,374E-04 6,18351 354,289
Al terminar de introducir los datos se observarán los resultados del cálculo

de las matrices de rigidez en locales y globales de la barra 3.
5. Ir a la hoja “K(L) y K (G) barra 4” introducir los datos geométricos y

tal como se aprecia en la siguiente tabla:

205939650000 0,01314 5 1,927E-04 4,712389 270
Tras la introducción de los datos geométricos y mecánicos de la barra 4 se

apreciarán los resultados del cálculo de las matrices de rigidez en locales y
en globales de la barra.
6. Ir a la hoja “K global estructura”. En esta hoja no es preciso introducir

datos, solo observar los resultados. Se comprobará que en un primer
momento se copian las matrices de rigidez en globales de las barras 1 a 4
previamente calculadas y que, posteriormente, se realiza el ensamblaje de la
matriz de rigidez global del pórtico a partir de las matrices de rigidez en
globales de las barras. La matriz de rigidez global es de dimensión 15 x 15.
Se puede verificar la simetría de la misma si se desea.
7. Ir a la hoja “ACCIONES EN NUDOS” e introducir ordenadamente las

acciones exteriores que actúan directamente en los nudos del pórtico,
indicadas en la figura siguiente:
20000 N
40000 N
30000 N
10000 Nm
50000 N
Figura A1.2: Cargas aplicadas en los nudos en unidades del SI.
En la tabla siguiente se indican los datos a introducir observando la

correspondencia de las acciones actuantes con los 15 grados de libertad del
pórtico:
VECTOR FUERZAS EXTERIORES (GLOBALES)

NUDO gdl VALOR UNIDAD
1 1 0,00 N
1 2 0,00 N
1 3 0,00 Nm
2 4 30000,00 N
2 5 20000,00 N
2 6 10000,00 Nm
3 7 0,00 N
3 8 -50000,00 N
3 9 0,00 Nm
4 10 40000,00 N
4 11 0,00 N
4 12 0,00 Nm
5 13 0,00 N
5 14 0,00 N
5 15 0,00 Nm
Tabla A1.9: Introducción de las acciones exteriores.
8. Ir a la hoja “DESPLAZAMIENTOS GLOBALES”. En esta hoja no se

necesita introducir datos, sino observar los cálculos que se realizan y los
resultados. En esta hoja se encuentran los siguientes c:
1. Se copia la matriz de rigidez global obtenida en la hoja “K global

estructura”.
2. A partir de la matriz de rigidez global se obtiene la matriz de rigidez
reducida eliminando las filas y columnas asociadas a grados de libertad
restringidos, que en el caso del pórtico son los grados de libertad 1, 2,
3, 13, 14 y 15. La matriz de rigidez reducida es de dimensión 9 x 9.
3. Se copia el vector de acciones sobre los nudos de la hoja “ACCIONES
EN NUDOS”.
4. Se forma el vector de acciones sobre los nudos libres a partir del
vector de acciones sobre todos los nudos. El vector de acciones sobre
los nudos libres es de dimensión 9 x 1.
5. Se calcula la matriz inversa de la matriz de rigidez reducida.
6. Se calcula el vector de desplazamientos en globales de los nudos
libres.
7. Se forma el vector de desplazamientos en globales de todos los nudos
del pórtico.
Todo este proceso realizado en esta hoja se resume en la siguiente ecuación

matricial:
En la siguiente tabla se reproduce el vector de desplazamientos de los nudos

en coordenadas globales:
VECTOR DESPLAZAMIENTOS DE NUDOS LIBRES

NUDO gdl SIMBOLO VALOR UNIDAD VALOR UNIDAD
1 1 1x 0,00000000 m 0,00000 mm
1 2 1y 0,00000000 m 0,00000 mm
1 3 1z 0,00000000 rad 0,00000 rad
2 4 2x 0,01159131 m 11,59131 mm
-
2 5 2y m -0,00099 mm
0,00000099
-
2 6 2z rad -0,00528 rad
0,00527952
3 7 3x 0,01759036 m 17,59036 mm
- -
3 8 3y m mm
0,06234258 62,34258
3 9 3z 0,00167777 rad 0,00168 rad
4 10 4x 0,02357855 m 23,57855 mm

-
4 11 4y m -0,05440 mm
0,00005440
-
4 12 4z rad -0,00144 rad
0,00144452
5 13 5x 0,00000000 m 0,00000 mm
5 14 5y 0,00000000 m 0,00000 mm
5 15 5z 0,00000000 rad 0,00000 rad
Tabla A1.10: Desplazamientos en coordenadas globales de los nudos del
pórtico.
En la siguiente figura se representa la deformada del pórtico obtenida
mediante el programa CESPLA. Se puede comprobar el orden de magnitud
de los desplazamientos; por ejemplo, los nudos 3 y 4 apenas se desplazan
según el eje Y global y es apreciable su desplazamiento según el eje X
global; el nudo 3 se desplaza apreciablemente según los ejes X e Y y el
orden de magnitud del desplazamiento en Y global es 6 veces el valor de se
desplazamiento según el eje X global.
2 4
1 5
Figura A1.6: Deformada del pórtico obtenida mediante el programa CESPLA

con el fin de apreciar los desplazamientos de los nudos libres.
En la siguiente figura se reproduce el pórtico con la deformada sin

representar las acciones exteriores con el fin de conseguir una ligera
amplificación en la imagen.
2 4
1 5
Figura A1.7: Deformada del pórtico obtenida mediante el programa CESPLA

sin representar las acciones exteriores.

9. Ir a la hoja “DESPLAZAMIENTOS LOCALES”. No se necesita introducir
datos, solo observar los resultados. En esta hoja se calculan los
desplazamientos de los nudos barra por barra en coordenadas locales a
partir de los desplazamientos de los nudos en coordenadas globales. Para
ello se utiliza la siguiente ecuación matricial:
Se observará que en esta hoja se encuentran los siguientes pasos:
1. Recopilación de las matrices de transformación de globales en locales

de cada barra.
2. En las celdas de fondo amarillo se forma barra a barra el vector de
desplazamientos en globales de los nudos a partir de los resultados de
la etapa 8.
3. En las celdas de fondo verde se calculan los desplazamientos en
coordenadas locales de los nudos de cada barra.
En la siguiente tabla se encuentran los resultados de los desplazamientos en

locales de los extremos de las 4 barras del pórtico.
DESPLAZAMIENTOS DE NUDOS EN COORDENADAS LOCALES

Nudos barra 1 Nudos barra 2 Nudos barra 3 Nudos barra 4
NUDO DESPLAZ. NUDO DESPLAZ. NUDO DESPLAZ. NUDO DESPLAZ.
1 0,000000 2 0,011534 3 0,023707 4 0,000054
1 0,000000 2 -0,001154 3 -0,060283 4 0,023579
1 0,000000 2 -0,005280 3 0,001678 4 -0,001445
2 -0,000001 3 0,011299 4 0,023467 5 0,000000
2 -0,011591 3 -0,063784 4 0,002292 5 0,000000
2 -0,005280 3 0,001678 4 -0,001445 5 0,000000
Tabla A1.11: Desplazamientos en coordenadas locales de los nudos del
pórtico agrupados barra a barra (Unidades en metros y radianes).
10. Ir a la hoja “ESFUERZOS EXTREMOS”. No se necesita introducir datos,

solo observar los resultados. En esta hoja se calculan los esfuerzos en los
extremos de las barras en coordenadas locales multiplicando la matriz de
rigidez de cada barra en coordenadas locales por el vector de
desplazamientos de los nudos extremos expresados también en coordenadas
locales. Se utiliza la siguiente ecuación matricial:
Esta ecuación, desarrollada, por ejemplo, para la barra 1, es así:

En las celdas amarillas se calculan los resultados sin necesidad de introducir

nuevos datos o copiar datos de otras hojas.
A partir de los esfuerzos en los extremos se pueden determinar las leyes de

esfuerzos a lo largo de las barras, es decir, las leyes de esfuerzos normales,
esfuerzos cortantes y momentos flectores. Para ello, en primer lugar, se
recopilan los datos obtenidos en esta hoja “ESFUERZOS EXTREMOS”:
ESFUERZOS BARRA 1 ESFUERZOS BARRA 2

N1 536 N N1 37993 N
V1 -6129 N V1 16839 N
M1 26614 Nm M1 67258 Nm
N2 -536 N N2 -37993 N
V2 6129 N V2 -16839 N
M1 -57258 Nm M1 144282 Nm

N1 38881 N N1 29464 N
V1 -25722 N V1 76129 N
M1 -144282 Nm M1 178848 Nm
N2 -38881 N N2 -29464 N
V2 25722 N V2 -76129 N
M1 -178848 Nm M1 201796 Nm
Tabla A1.12: Esfuerzos en coordenadas locales calculados en los extremos
de las 4 barras que conforman el pórtico.
NOTA: A efectos de representación:
- Los axiles (compresión con signo negativo) cambian de signo en el nudo

izquierdo o nudo inicial de la barra. En el nudo final no se cambia de signo.
- Los cortantes cambian de signo a la derecha o nudo final de la barra. En el
nudo inicial el valor del cortante no cambia de signo
- Los momentos cambian de signo a la izquierda o nudo inicial de la barra.
En el nudo final no se cambia el signo.

Para recordar los criterios anteriores se puede utilizar la regla IDI

(Izquierda-Derecha-Izquierda) aplicada a los esfuerzos en el siguiente orden
NVM (axil-cortante-flector), que es el orden en que se calculan.
Por lo tanto, las tablas anteriores, a efectos de representación, se convierten

en las siguientes:

N1 -536 N N1 -37993 N
V1 -6129 N V1 16839 N
M1 -26614 Nm M1 -67258 Nm
N2 -536 N N2 -37993 N
V2 -6129 N V2 +16839 N
M1 -57258 Nm M1 144282 Nm

N1 -38881 N N1 -29464 N
V1 -25722 N V1 76129 N
M1 +144282 Nm M1 -178848 Nm
N2 -38881 N N2 -29464 N
V2 -25722 N V2 +76129 N
M1 -178848 Nm M1 201796 Nm
Tabla A1.13: Valores de los esfuerzos en los extremos de las barras según
los criterios de Resistencia de Materiales.
Las razones de esos cambios de signo se pueden deducir a partir de los

siguientes aspectos:
a. Los esfuerzos que se obtienen en el método matricial tienen los signos

de los ejes locales en cada barra.
b. De acuerdo a los criterios generalizados de Resistencia de Materiales,
los axiles positivos tienden a estirar una barra o una determinada
sección de la misma, los esfuerzos cortantes positivos tienden a
deslizar una porción de una barra en el sentido de las agujas del reloj
y los momentos flectores positivos tienden a comprimir las fibras
superiores de una barra o porción de la misma.
c. Se debe tener en cuenta el equilibrio de una parte de la barra al
seccionarla mediante un plano de corte intermedio.
Se explicarán gráficamente los razonamientos anteriores.

a. En la siguiente figura se representan los esfuerzos en los extremos de

las barras en coordenadas locales obtenidos en el cálculo matricial, en
este caso es la barra 1. Los signos que resultan del cálculo se
corresponden con los ejes locales de la barra.
Figura A1.8: Signos positivos de los esfuerzos obtenidos en los extremos de

una barra en el método matricial. Se expresan en las coordenadas locales de
la barra.
En las siguientes figuras se recuerdan los signos positivos de esfuerzos

axiles, esfuerzos cortantes y momentos flectores aplicados a un tramo de
longitud muy pequeña de una barra. Este criterio no está universalmente
aceptado. Por ejemplo, el programa CESPLA, considera los esfuerzos
cortantes positivos al contrario que lo que se acepta en esta asignatura de
Teoría de Estructuras y Construcciones Industriales. Pero esto es solo un
criterio: nadie está equivocado.

N N
+
Figura A1.9: Signo positivo del esfuerzo axil en una parte de longitud
infinitesimal de una barra.
V + V
Figura A1.10: Signo positivo del esfuerzo cortante en una parte de longitud
M + M
Figura A1.11: Signo positivo del momento flector en una parte de longitud
En las figuras siguientes se corta la barra por una sección situada a la

derecha del nudo inicial. Se dibujan los esfuerzos en el nudo inicial con el
signo positivo tal como se obtienen del método matricial. Los esfuerzos en la
sección donde se realiza el corte se representan con el signo que se utiliza
en Resistencia de Materiales. Los valores de estos esfuerzos en la sección x
se deducen considerando el equilibrio de la parte cortada de la barra.
V1
M1 M(x)
N1 N(x)
V(x)
V50
x
Figura A1.12: Signo positivo del esfuerzo axil, esfuerzo cortante y
momento flector en una sección de la barra situada a una longitud x.
De la figura anterior, considerando el equilibrio del tramo de la barra de

longitud x partiendo del nudo inicial, se deducen las leyes de esfuerzos
axiles, esfuerzos cortantes y momentos flectores a lo largo de la barra del
pórtico que solo tiene cargas nodales. Estas leyes se son:
Si x es tan pequeño que tiende a 0 se comprenden los cambios de signos

que se realizan para convertir valores de esfuerzos en los extremos
deducidos a partir del método matricial en valores de esfuerzos en los
extremos con los signos utilizados en Resistencia de Materiales (regla IDI).
En la figura siguiente se adapta la figura anterior cuando x tiende a 0.

V1 M(x)
M1
N1 N(x)
V(x)
momento flector en una sección de la barra muy próxima a la sección inicial.
En la sección izquierda se representan los esfuerzos que se deducen del

método matricial con los signos positivos.
En la sección derecha, sección que representa el corte de la barra, se

dibujan los esfuerzos positivos de acuerdo a los criterios de Resistencia de
Materiales.
Teniendo en cuenta que esa parte de longitud infinitesimal está en equilibrio,

se puede deducir la regla IDI:
 Al esfuerzo axil obtenido aplicando el método matricial hay que

cambiarle el signo.
 Al esfuerzo cortante obtenido en el método matricial no hay que
cambiarle el signo.
 Al momento flector calculado con el método matricial hay que
cambiarle el signo.
Aislando una parte infinitesimal de la barra partiendo del nudo derecho se

puede terminar de comprender la regla IDI.
En la siguiente figura se realiza el mismo análisis gráfico realizado

anteriormente pero aplicado al nudo derecho de la barra. En el extremo
derecho se dibujan los esfuerzos obtenidos aplicando el método matricial y
en una sección muy próxima situada a la izquierda se sitúan los esfuerzos
con sentido positivo según el criterio de Resistencia de Materiales.

M(x) V(x) V2
M2
N(x) N2

momento flector en una sección de la barra muy próxima al nudo final.
En el caso del nudo derecho de la barra, los cambios de signos que se

realizan para convertir valores de esfuerzos deducidos a partir del método
matricial en valores de esfuerzos en los extremos con los signos utilizados
en Resistencia de Materiales (regla IDI) se deducen del equilibrio de la
porción de barra de la figura anterior:
A la derecha solo cambia de signo el cortante.
A continuación se dibujarán los diagramas de esfuerzos deducidos. Las

figuras que se indican a continuación se han obtenido con el programa
CESPLA:
-536 3
2 3
2 4
-37993 -38881
4
1
-29464
1 5
Figura A1.14: Representación de los diagramas de esfuerzos normales.
La siguiente figura completa la anterior indicando los esfuerzos axiles en las

barras:

3
2 3
N -37993 N -38881
2 4
N -29464
N -536
4
1
1 5
Figura A1.15: Representación de los valores de los esfuerzos normales.
En la siguiente figura se representan los diagramas de esfuerzos cortantes:
-25722
-6129
3
2 3
2 4
+16839
+76129
1
4
1 5
Figura A1.16: Representación de los esfuerzos cortantes.
En la siguiente figura se representan los diagramas de momentos flectores:
-178848
-67258
3
2 3 -178848
2 4
-26614
1
+201796
-57258
1 5
+144282
Figura A1.17: Representación de los momentos flectores.
12. Ir a la hoja “REACCIONES”. No se necesita introducir datos, solo

observar los resultados. En esta hoja se utiliza parte de la matriz de rigidez
global para multiplicarla por los desplazamientos de los nudos libres en
coordenadas globales y obtener así las reacciones en coordenadas también
globales, que se reproducen en la siguiente tabla:
R1X 6129 N
R1Y 536 N
M1Z 26614 Nm
R5X -76129 N
R5Y 29464 N
M5Z 201796 Nm
Tabla A1.14: Valores de las reacciones en coordenadas globales.
En las siguientes figuras se representan las reacciones obtenidas mediante el

programa CESPLA.
2
2
1
26614
536
6129 76129
1
Figura A1.18: Representación de las reacciones en el apoyo izquierdo.
3 29464
2 3
4
4
76129
5
201796
Figura A1.19: Representación de las reacciones en el apoyo derecho.

En la hoja “REACCIONES” se utiliza la siguiente expresión, que solamente es
válida cuando, como en este caso, solo hay cargas en los nudos de la
estructura:
Si hubiera cargas en las barras de la estructura, la expresión anterior no

valdría. En esos casos, una posibilidad de obtener las reacciones se basa en
convertir los esfuerzos de los extremos de las barras que se obtienen en
coordenadas locales a coordenadas globales. Estos esfuerzos expresados en
coordenadas globales coinciden con las reacciones en los apoyos. Este
proceso se realiza a continuación:
El vector de esfuerzos en el extremo inicial de la barra 1 en

coordenadas locales es:
La matriz de trasformación de locales a globales de la barra 1 es:
Las reacciones en el apoyo 1 se obtienen multiplicando la matriz L por el

vector de esfuerzos en el extremo inicial de la barra 1:
El vector de esfuerzos en el extremo final de la barra 4 en coordenadas

locales es:
La matriz de trasformación de locales a globales de la barra 4 es:

Las reacciones en el apoyo 5 se obtienen multiplicando la matriz L por el

vector de esfuerzos en el extremo final de la barra 4:
Se comprueba que los valores de las reacciones obtenidas por este método
coinciden con los valores obtenidos mediante el método propuesto en la hoja
Excel.


PRACTICA 3A - Métodos Matriciales - Tutorial

Cargado por

Copyright:

Formatos disponibles

PRACTICA 3A - Métodos Matriciales - Tutorial

Cargado por

Información del documento

Descripción original:

Título original

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

PRACTICA 3A - Métodos Matriciales - Tutorial

Cargado por

Copyright:

Formatos disponibles

Teoría de Estructuras. Practicas sobre el Tema 4: Métodos matriciales.

CÁLCULO MATRICIAL DE UN PÓRTICO A DOS

Esta práctica consiste en el cálculo de un pórtico biempotrado con cargas

En la siguiente figura se representan las dimensiones principales del pórtico,

Figura A1.1: Pórtico biempotrado a 2 aguas de nudos rígidos.

En la siguiente figura se representan las cargas nodales que actúan sobre el

Figura A1.2: Cargas aplicadas en los nudos.

Departamento de Ingeniería Mecánica, Energética y de Materiales. 1

NUDO NUDO AREA M. DE INERCIA

En la siguiente figura se representa de nuevo el pórtico del ejemplo,

Con la información anterior se realiza la siguiente tabla, en la cual se

Departamento de Ingeniería Mecánica, Energética y de Materiales. 2

Departamento de Ingeniería Mecánica, Energética y de Materiales. 3

2. En la hoja “K(L) y K (G) barra 1” introducir los datos geométricos y

CÁLCULO DE LAS MATRICES DE RIGIDEZ EN LOCALES Y GLOBALES BARRA 1

NOTA: No modificar la celda en rojo. En la celda F3 con fondo rojo se calcula

En esta hoja “K(L) y K (G) barra 1” se realiza primero el cálculo de la

MATRIZ DE RIGIDEZ EN LOCALES DE LA BARRA 1

Estos cálculos se realizan a partir de los datos geométricos y mecánicos de

Departamento de Ingeniería Mecánica, Energética y de Materiales. 4

Un poco más abajo, en las celdas B37:G42 se calcula la matriz de rigidez de

La ecuación matricial para el cálculo de la matriz de rigidez en globales de

La matriz de transformación de coordenadas locales a globales L se

La matriz de transformación de coordenadas globales a locales LT se obtiene

Una vez recordadas las expresiones anteriores se llega a la matriz de rigidez

Departamento de Ingeniería Mecánica, Energética y de Materiales. 5

3. Ir a la hoja “K(L) y K (G) barra 2” e introducir los datos geométricos y

CÁLCULO DE LAS MATRICES DE RIGIDEZ EN LOCALES Y GLOBALES BARRA 2

Al igual que en el caso de la barra 1 se pueden revisar los resultados del

4. En la hoja “K(L) y K (G) barra 3” introducir los datos geométricos y

CÁLCULO DE LAS MATRICES DE RIGIDEZ EN LOCALES Y GLOBALES BARRA 3

Al terminar de introducir los datos se observarán los resultados del cálculo

5. Ir a la hoja “K(L) y K (G) barra 4” introducir los datos geométricos y

CÁLCULO DE LAS MATRICES DE RIGIDEZ EN LOCALES Y GLOBALES BARRA 4

Tras la introducción de los datos geométricos y mecánicos de la barra 4 se

6. Ir a la hoja “K global estructura”. En esta hoja no es preciso introducir

7. Ir a la hoja “ACCIONES EN NUDOS” e introducir ordenadamente las

Figura A1.2: Cargas aplicadas en los nudos en unidades del SI.

En la tabla siguiente se indican los datos a introducir observando la

VECTOR FUERZAS EXTERIORES (GLOBALES)

8. Ir a la hoja “DESPLAZAMIENTOS GLOBALES”. En esta hoja no se

1. Se copia la matriz de rigidez global obtenida en la hoja “K global

Todo este proceso realizado en esta hoja se resume en la siguiente ecuación

En la siguiente tabla se reproduce el vector de desplazamientos de los nudos

VECTOR DESPLAZAMIENTOS DE NUDOS LIBRES

Departamento de Ingeniería Mecánica, Energética y de Materiales. 8

Figura A1.6: Deformada del pórtico obtenida mediante el programa CESPLA

En la siguiente figura se reproduce el pórtico con la deformada sin

Figura A1.7: Deformada del pórtico obtenida mediante el programa CESPLA

Departamento de Ingeniería Mecánica, Energética y de Materiales. 9

Se observará que en esta hoja se encuentran los siguientes pasos:

1. Recopilación de las matrices de transformación de globales en locales

En la siguiente tabla se encuentran los resultados de los desplazamientos en

DESPLAZAMIENTOS DE NUDOS EN COORDENADAS LOCALES

10. Ir a la hoja “ESFUERZOS EXTREMOS”. No se necesita introducir datos,

Esta ecuación, desarrollada, por ejemplo, para la barra 1, es así:

Departamento de Ingeniería Mecánica, Energética y de Materiales. 10

En las celdas amarillas se calculan los resultados sin necesidad de introducir