Análisis Matricial

MÉTODO MATRICIAL Y ELEMENTOS FINITOS

1. Conceptos básicos

Grados de libertad de nudos

Concepto

En el Método de Elementos Finitos, se separa la estructura en pequeños

elementos, lo suficientemente pequeños como para lograr la precisión
adecuada en los resultados buscados a partir de los elementos de estudio.
Aunque en el análisis de pórticos y cerchas estos elementos no necesitan ser
tan pequeños.

Al dividir todo el dominio de la estructura en pequeños elementos, estos

elementos están unidos por nudos. Estos nudos tienen distintas propiedades
según las solicitaciones que se necesiten estudiar. De hecho, podría pensarse
que son los nudos los que condensan las propiedades de la estructura
discretizada. La representación de elementos y nudos está mostrada en la
siguiente figura:
Como se dijo, estos nudos tienen diferentes propiedades según la estructura
que se desee estudiar. Son en general los nudos los que tienen la propiedad
de desplazarse y rotar, y es a los nudos donde llega toda la carga de la
estructura.

Las direcciones en las que cada nudo es capaz de rotar y girar son conocidas
como grados de libertad del nudo. En el caso más general en el espacio
tridimensional un nudo es capaz de rotar en tres direcciones y también de
desplazarse en tres direcciones. Entonces se tendrán un total de seis grados
de libertad para el nudo (tres rotaciones y tres desplazamientos rectilíneos)

Reducción de grados de libertad

Sin embargo, el coste computacional de analizar todas las estructuras con 6

grados de libertad por nudo es demasiado alto.

Existen muchas estructuras que no necesitan más que unos pocos grados de
libertad por nudo para analizarse satisfactoriamente.

Estructuras con 1 GDL por nudo

Por ejemplo, elementos sometidos a compresión simple pueden analizarse
con nudos cuyo desplazamiento sea unidireccional. En el caso de la figura
siguiente se analiza la compresión de la torre producto de su peso propio.
Cada nudo lleva consigo el peso propio de una porción de la estructura.

Otro ejemplo de análisis de nudos con un solo grado de libertad es el análisis

dinámico de estructuras que pueden simplificarse de forma tal que la masa
de cada losa se concentre en un solo punto de análisis, teniendo así una
simplificación significativa de la estructura logrando en muchos casos
resultados aceptables. En este caso el desplazamiento a analizarse es en el
sentido de la aceleración. En la mayoría de los casos la traslación es tan
grande comparada con la rotación de las losas, que la rotación se hace
despreciable. Notarás en la siguiente figura que la estructura tiene un solo
grado de libertad por nudo
Estructuras con 2 GDL por nudo

Cuando se trabaja con dos grados de libertad, existen dos tipos principales
de análisis. El análisis de vigas y el análisis de cerchas en el plano.

Cuando se trabaja con cerchas (figura acá arriba), cada nudo necesita tener
la capacidad de desplazarse tanto en el sentido vertical como en el horizontal.
Los desplazamientos en los apoyos no se dibujaron debido a que, según la
configuración del apoyo, estos desplazamientos están restringidos.
Al trabajar con vigas el desplazamiento horizontal en general no tiene mucha
importancia, ya que las cargas son casi siempre verticales. Debido a esta
configuración de cargas es que el desplazamiento vertical es predominante.
Sin embargo, debido a que una viga trabaja a flexión, la forma de los nudos
de medir la magnitud de esta flexión es mediante una combinación de
desplazamiento vertical y rotación de los nudos. Es por eso que una viga
trabaja con dos grados de libertad por nudo. Estos grados de libertad sin
embargo son diferentes a los de cerchas, y en consecuencia las ecuaciones
que rigen el comportamiento de estos nudos serán también diferentes. La
viga abajo tiene dos grados de libertad al igual que la cercha arriba, pero
estos GDL son diferentes a los de la cercha:

(Se ve que en los apoyos no existe la flecha de desplazamiento vertical. Esto

es debido a que en los apoyos de este tipo el desplazamiento vertical está
restringido).

Estructuras con 3 GDL por nudo

Cuando se analizan pórticos en el plano, se necesitarán tres grados de libertad

por nudo. Dos de traslación y uno de rotación en el plano estudiado, como
se muestra en la figura. En este caso tanto la rotación como ambas
traslaciones son necesarias para formas las ecuaciones de equilibrio de cada
nudo (ver siguiente figura).
Seis Grados de Libertad por nudo

Finalmente, los pórticos en 3d cuentan con 6 grados de libertad por nudo,

que es el máximo número de desplazamientos y giros que puede tener un
nudo en análisis estructural. Cualquier programa de cálculo estructural
trabaja con este número de GDL a menos que al programa se le indique lo
contrario.
2. Elementos con 1 grado de libertad por nudo

Ejemplo de elementos con nudos con 1 grado de libertad

La siguiente imagen muestra un ejemplo completo de la resolución de

un sistema de barras conectadas y con una carga en uno de sus
extremos.

https://3.bp.blogspot.com/-
hwbvEsFsF4o/W2z8JpV7J5I/AAAAAAAACXY/Io7xJBIwrYQw6
TCm7thMc7M041crK4WBwCLcBGAs/s000/elementos%2Bcon%2
B1GDL.jpg
 3. Análisis de cerchas

Teoría

Matriz de rigidez de barra con deformación axial solamente

Matriz de rigidez para elemento de cercha en 1D

En el diseño de cerchas por el método matricial se pueden seguir dos caminos
para la deducción de la matriz de rigidez.
El primer método, más complejo y que necesita una base más sólida de
cálculo diferencial, se basa en métodos energéticos igualando el trabajo de
los desplazamientos virtuales producidos por cargas externas, con el trabajo
de las deformaciones virtuales producido por los esfuerzos internos.
El segundo método, que es más sencillo y más directo, se basa en principios
de resistencia de materiales. Este es el método que se explica a continuación.
Dado un elemento de barra en 1 dirección, como el que sigue:

Este elemento de barra tiene una sección transversal A y un módulo elástico

E.
Cualquiera de los nudos puede desplazarse en la dirección x, y producto de
este desplazamiento deformar la barra.
Si asumimos la barra perfectamente elástica, es válida la ecuación:
 σ = E*ε
y como el esfuerzo σ = F/A = fuerza / área
y como la deformación unitaria ε = du/dx = derivada de los desplazamientos
/ derivada de x
Reemplazando en la ecuación anterior queda:
F = E*A (du/dx)

y como para elementos no infinitesimales, para una carga aplicada en el nudo

1, la derivada de los desplazamientos respecto de x puede interpretarse como
la diferencia del desplazamiento del nudo 1 menos el desplazamiento del
nudo dos divididos entre la longitud de la barra:

F1 = E*A (u1 - u2)/L

Ya que, por equilibrio, la fuerza del nudo uno debe ser igual y opuesta a la
fuerza del nudo 2,

Se debe cumplir que,

F1 = -F2
Por tanto
F2 = EA (u2-u1)/L

Finalmente escribiendo ambas fórmulas de forma matricial, el sistema de dos

ecuaciones con dos incógnitas queda de la siguiente manera:
Este es el sistema de ecuaciones que describe el comportamiento de una barra
sometida a fuerzas en una dirección, donde de manera abreviada se escribe
{F}=[k]{u} donde:

{F} = Vector de cargas

[k] = Matriz de rigidez local para un elemento
{u} = Vector de desplazamientos

Aplicando las condiciones de contorno y cargas externas puede resolverse

cualquier problema de elementos cargados axialmente, con deformación
unidireccional.

Transformación de coordenadas para elementos de barra en 2

direcciones

Matriz de Rigidez para Elementos de Barra en 2 Direcciones

Cuando se trabaja con cerchas con barras con deformación axial solamente,
por lo general se necesitan modelos simples en dos dimensiones que
simulen el comportamiento de una estructura como la siguiente:
En estos casos, las fórmulas de matriz de rigidez del tema anterior necesitan
una transformación de coordenadas para poderlas aplicar en la cercha de la
figura.
Entonces, dado un elemento de cercha genérico, como el siguiente, inclinado
un ángulo β, con los superíndices (L) para denotar coordenadas locales:
Si escribimos los desplazamientos locales en función de los desplazamientos
globales u1, u2, u3 y u4, queda:

que, escrito en forma matricial, se muestra así:

(ecuación 1)

De igual manera, en vez de dibujar la barra mencionada en función de los

desplazamientos la dibujamos en función de las fuerzas axiales aplicadas en
los nodos:
esta vez escribiendo las fuerzas globales en función de las fuerzas locales:

y escribiendo las anteriores ecuaciones en forma matricial, se tiene:

(ecuación 2)
Ahora, la fórmula deducida en el anterior Post de la matriz de rigidez de
elementos en una dirección deformados axialmente, es la siguiente:

(ecuación 3)

Reemplazando la ecuación 1 en la ecuación 3, y este resultado finalmente en

la ecuación 2, tenemos el siguiente resultado:

Donde

y [K] es la matriz de rigidez para elementos deformados axialmente, con

desplazamiento en dos direcciones.
Finalmente, la ecuación a resolver será {F} = [K]{u}, que es la forma
abreviada matricial de escribir el sistema de ecuaciones para resolución de
estructuras con cargas en los nudos.
Ejemplos

Ejemplo 1 - Cercha por método Matricial

Cercha por Método Matricial

El siguiente ejemplo consiste en la resolución de una cercha o armadura en
dos dimensiones mediante el método matricial (un método muy similar o
hasta homologo al método de los elementos finitos en su forma más básica
en 1-D).

1. Matrices individuales de cada elemento

2. Matriz global ensamblada
3. Eliminación de grados de libertad (casillas celestes. Ver
planilla)
4. Resolución de la ecuación [Kr]{Ur}={Rr} que corresponde a
la ecuación sin los grados de libertad restringidos.
5. Multiplicación de la matriz [K]{u} ={R} para encontrar las
reacciones en la cercha una vez encontrados los desplazamientos.

Entonces, la cercha es la siguiente:

Se puede apreciar en el gráfico, la enumeración de los grados de libertad de
los nudos, el nombre de los elementos (enumerados), así como las
propiedades de elasticidad y sección de cada una de las barras (constante).

Los pasos a seguir en la resolución de la cercha son:

1. Encontrar la matriz individual de rigidez de cada elemento

mediante las fórmulas:

2.
3. Ensamblar la matriz de rigidez según los grados de libertad de
los nudos de cada elemento (esto está muy bien detallado en la
planilla excel adjunta).
 4. Una vez ensamblada la matriz de rigidez global [K] se
eliminan los grados de libertad restringidos (en este caso, los 1, 2,
11, 12. Se puede observar que se trata de una cercha hiperestática)
5. Eliminados los grados de libertad, se procede a la resolución
del sistema reducido (en este caso un sistema de ecuaciones de 8
ecuaciones con 8 incógnitas)
6. Si llamamos la matriz de rigidez reducida [Kr], entonces el
sistema que se forma es:
7. Kr]{Ur}={Rr} Donde Ur y Rr son los vectores de
desplazamientos y cargas reducidos respectivamente.
entonces,

{Ur}=[Kr]-1 x {Rr}

8. Resueltos los {Ur} se vuelve al sistema global [K]{U}={R} y

al multiplicar la matriz de rigidez por los desplazamientos ya
encontrados de los grados de libertad de los nudos, se encuentran las
reacciones de los apoyos
https://drive.google.com/file/d/1m0XiZPnRO5-
fzii3eaQLYISXhcIBP6yF/view?usp=sharing

Ejemplo 2 - Cercha por método Matricial + SAP2000 final del

curso

Cercha - Método Matricial + SAP2000

En el siguiente ejemplo, veremos la resolución de una cercha en dos

dimensiones, compuesta por elementos de sección tubular metálica.
Solución

Se resolverá esta cercha mediante métodos matriciales, y se hará la

comparación de resultados con el programa SAP2000.

El procedimiento a seguir es el mismo que el de los dos ejemplos

anteriores:

1. Calcular matrices individuales de cada elemento

2. Ensamblar la matriz global

3. Eliminación de grados de libertad (casillas celestes. Ver planilla)

4. Resolución de la ecuación [Kr]{Ur}={Rr} que corresponde a la

ecuación sin los grados de libertad restringidos.

5. Multiplicación de la matriz [K]{u}={R} para encontrar las

reacciones en la cercha una vez encontrados los desplazamientos.

Sin embargo, en este ejemplo se adiciona el Postproceso correspondiente al

cálculo de solicitaciones y esfuerzos en cada barra. Este postproceso se
detalla a continuación:
Los desplazamientos de los nudos de una barra cualquiera, (digamos la
barra 1) son U1, U2, U3, U4. Entonces, las coordenadas desplazadas de los
nudos inicial y final de esta barra serán:

(x1-desp,y1-desp) = (x1+U1 , y1+U2)

(x2-desp,y2-desp) = (x2+U3 , y2+U4)

siendo (x1 ; y1 ) las coordenadas iniciales del nudo y (x2 , y2)las coordenadas
finales.

Entonces, conociendo las coordenadas desplazadas, puede determinarse la

longitud deformada de la barra mediante la ecuación:

y finalmente mediante las ecuaciones de resistencia de materiales, de

deformación axial:

donde

N = solicitación en la barra.
L = longitud de la barra

E = Módulo de elasticidad

e = deformación unitaria.

https://drive.google.com/file/d/1LFY5T383dCLkts2tzgzEuy5m0z6SLI
Op/view?usp=sharing

https://drive.google.com/file/d/1_K5XVWe45-
ltAmGsokr8Vndq2XgqLqI0/view?usp=sharing

Ejemplo 3 - Cercha Hiperestática. Método Matricial

Ejemplo 1 de cercha en 2D por método Matricial

Planteamiento

El siguiente ejemplo es una cercha o armadura calculada mediante el

método matricial. Este método es conocido por su versatilidad al momento
de calcular armaduras hiperestáticas o cuando se tenga un gran número de
elementos en la estructura.

La estructura que calcularemos es la siguiente:

Solución

Matriz de rigidez para un elemento

La teoría del método matricial se puede encontrar en muchos libros de

análisis estructural, pero básicamente se resume en el ensamblaje de
matrices individuales que contengan toda la información de esfuerzo –
deformación de cada elemento. La matriz de cada elemento que contiene
estas propiedades se conoce como Matriz de rigidez [k], y es la mostrada a
continuación
donde:

• E = módulo elástico de cada elemento

• A = Superficie de la sección transversal de cada elemento

• L = Longitud de cada elemento de barra

• α = ángulo de cada barra medida desde el eje x positivo. El ángulo es

antihorario positivo

Reemplazo de datos para cada elemento

Una vez reemplazados los datos de cada matriz de rigidez, para cada
elemento de la cercha (que se ve en la planilla Excel adjunta).
Posteriormente se ensambla la matriz de rigidez global compuesta de cada
una de las matrices elementales. La dimensión de la matriz de rigidez
global es la misma que la cantidad de grados de libertad que existan en la
estructura. En este caso se tienen 8 grados de libertad, por tanto, la matriz
de rigidez será de 8×8.

Si notas en cada matriz de rigidez de cada elemento, existe la numeración

encima indicando de que grado de libertad a qué grado de libertad va esa
respectiva barra. Estas posiciones de la matriz de rigidez local se deben
trasladar a la matriz de rigidez global. Si existen dos elementos cayendo
dentro de la misma casilla en la matriz de rigidez global, entonces se
suman.

Ensamblaje de la matriz Global

ensamblada la matriz de rigidez global se tiene entonces:

Teniendo la matriz de rigidez global, se conforma el sistema de ecuaciones

siguiente:
Resolución del sistema de ecuaciones – Encontrando Desplazamientos

De este sistema de ecuaciones existen datos conocidos (que son las

condiciones de borde del problema). Los desplazamientos de los grados de
libertad 1, 2, 5 son conocidos e iguales a cero. Por otro lado, la carga en el
grado de libertad 4 es igual a 50 hacia abajo, por tanto, negativa. Las cargas
en los grados de libertad 3, 6, 7, 8 son nulas. Finalmente, las reacciones
serán incógnitas en los grados de libertad restringidos (desplazamientos =
0). Por tanto, queda:

el siguiente paso consiste en anular filas y columnas de los grados de

libertad restringidos (1, 2, 5)

Anuladas estas filas y columnas queda el sistema de ecuaciones siguiente:

Multiplicando la matriz de rigidez invertida por las cargas, quedan los
siguientes desplazamientos:

…que son los desplazamientos que faltaban conocer.

Reacciones en la cercha

Con estos desplazamientos ya se pueden conocer las reacciones del sistema

de cercha. Para esto reemplazamos los desplazamientos en el sistema de
8×8 original, donde ya todos los desplazamientos son conocidos:
Luego de la multiplicación respectiva de la matriz de rigidez [k]*{u}
queda:

que son las reacciones de los grados de libertad restringidos.

Ejemplo 4: Desplazamientos en nudos y reacciones en

cercha

Ejemplo: Desplazamientos en nudos de cercha y reacciones

Ejemplo: Desplazamientos en nudos de cercha y reacciones resueltas por

el método matricial o también llamado método de elementos finitos.

Este es un ejemplo que a diferencia de los anteriores, muestra el

procedimiento paso a paso para la resolución de la cercha.
Esta es una cercha isostática, pero el procedimiento es válido incluso para
cerchas hiperestáticas.

https://3.bp.blogspot.com/-
uvYk8C1wQ6I/W4lZOapc0cI/AAAAAAAACao/BxYcu_jo4sY6X77Mi
vMJnu61-JGGGtfaQCLcBGAs/s000/Cercha.jpg

Ejemplo 5 - Solicitaciones en elementos de cercha. Paso

final

Cercha (5) - Ejemplo de cálculo integral de cercha por método

matricial

Cercha (5) - Ejemplo de cálculo integral de cercha por método

matricial

El siguiente ejercicio es el primero en el que presento el cálculo completo

de cercha, desde los desplazamientos hasta las solicitaciones de la cercha
para cada elemento, por el método matricial en base a elementos finitos.

El cálculo es el mismo siempre:

a) se encuentran las matrices elementales
b) se ensambla la matriz global
c) se reduce el sistema de ecuaciones anulando los grados de libertad
restringidos
d) se resuelve el sistema de ecuaciones reducido encontrando los
desplazamientos
e) se reemplazan los desplazamientos en el sistema de ecuaciones global
d) se resuelve la multiplicación [K]{u}={R} para encontrar las reacciones
R
e) se resuelven las solicitaciones de cada elemento multiplicando [Ke]{Ue}
que son: matriz elemental de rigidez por desplazamientos de los nudos del
elemento. El resultado serán las fuerzas en los nudos.
f) se obtiene la resultante de estas fuerzas en los nudos para obtener la
solicitación en la barra

Para seguir el ejemplo, dejo el análisis completo a continuación:

https://2.bp.blogspot.com/--
Nd9M2oiBvk/W46tsLmuT3I/AAAAAAAACb8/xY4xTS6JVrEd92YIvQJ
vkitASV09vE0tgCLcBGAs/s000/cercha5.jpg

Ejemplo 6 - Cercha Hiperestática por método de rigidez.

5 elementos

ejemplo: Cercha hiperestática mediante método matricial

Se procuró realizar este ejemplo lo más didácticamente posible, cuidando

que el lector entienda paso a paso todo el proceso de cálculo de cerchas por
el método matricial. Este método es de preferencia muy sistemático, por
tanto, el que entienda como resolver este ejercicio entenderá como puede
resolver cualquier cercha en 2 dimensiones, sin importar su grado de
hiperestaticidad.

El ejercicio se plantea como sigue:

https://2.bp.blogspot.com/-EM0KaZQ9VHc/W6-
6Csl_s4I/AAAAAAAACgQ/55ZW9MVyjOIuIxtpOfGB2fiFWTIS_1dgwC
LcBGAs/s000/publicacion.jpg
 Ejemplo 7 - Cercha solicitada producto de su peso
propio

Solicitaciones de cercha producto de su peso propio

El siguiente es un ejercicio de cercha deformada solamente bajo su peso

propio. El ejercicio es sencillo, sin embargo, algo largo ya que el
desarrollo está bien detallado.

Para el ejercicio se pide utilizar el método de rigidez o también llamado

método matricial para la resolución de:

-Reacciones en los apoyos

-Desplazamiento de los nudos

-Solicitaciones en una de las barras

Si se deseara encontrar la solicitación en todas las barras, el procedimiento

es repetitivo.

https://3.bp.blogspot.com/-BQRYM36Pm_s/W-
RkVSd7ThI/AAAAAAAACmc/93M5FpoDETQz2norMsCrkQ5mUWcag3
_CACLcBGAs/s000/desarrollo.jpg

4. Análisis de vigas

Teoría

Deducción de matriz de rigidez de una viga

Demostración de matriz de rigidez de elemento de viga. Método

Matricial
Para conocer la matriz de rigidez de un elemento de viga necesitamos
encontrar la interacción de los desplazamientos de la viga en relación
con las solicitaciones.
Si por ejemplo tenemos la siguiente viga, con las deformaciones u1,
u2, u3, y u4 mostradas:

Nótese que los desplazamientos U2 y U4 son en realidad giros.

Estos desplazamientos generan en los extremos de la viga, ciertas

solicitaciones de cortante y de momento flector.

Para deducir la matriz de rigidez de una viga debemos analizar las

solicitaciones resultantes producto de cada uno de los cu atro
desplazamientos mencionados:
Para resolver caso por caso estas solicitaciones provocadas por cada
uno de los desplazamientos mostrados utilizaremos el método de
integración de la ecuación diferencial de la elástica de la viga.
Por ejemplo para las solicitaciones del desplazamiento U1, se tiene:

Integrada 4 veces esta ecuación, se encuentran entonces las constantes

de integración en base a las siguientes condiciones de contorno:

Resolviendo las constantes de integración:

Y con las constantes encontradas se puede conocer las ecuaciones de
momento flector M(x) = d2u/dx2 y de cortante V(x) = d³u/dx³ en
función de x. Para conocer las solicitaciones en los extremos de la viga
se reemplaza x = 0 y además x = L:
Encontradas las solicitaciones en cada extremo, se grafican los
diagramas V y M

Se repite el proceso para cada uno de los desplazamientos U2, U3, U4:
Desplazamiento U2 (giro del nudo izquierdo):
Desplazamiento U3
Desplazamiento U4 (giro del nudo derecho):
Habiendo terminado la obtención de las solicitaciones de M y V
producto de cada caso de desplazamiento, ahora se suman los 4
desplazamientos para generalizar el caso en que el elemento de viga
esté sometido a una combinación de los 4 desplazamientos:

Este sistema de cuatro ecuaciones con 4 incógnitas puede escribirse

más ordenadamente de forma matricial de la siguiente manera:
La matriz 4x4 mostrada constituye la matriz de rigidez de una viga en
2 dimensiones sometida a desplazamientos verticales y giros en los
nudos. Las fuerzas F1, F2, F3, F4 representan en realidad las
solicitaciones V1, M1, V2, M2 respectivamente. Se debe tener en
cuenta que, en el método matricial, la convención de signos para
desplazamientos y solicitaciones es: Positivo hacia arriba, positivo
hacia la derecha, giro positivo antihorario.

Diagramas de V y M en vigas a partir de los desplazamientos

Diagramas de Momento y Cortante a partir de desplazamientos

Cuando se realiza el análisis de vigas mediante elementos finitos o

método matricial, se modela la estructura y los primeros resultados
obtenidos son los correspondientes a los desplazamientos en los nudos
de la estructura. Una vez encontrados los desplazamientos se procede
a encontrar las reacciones en los apoyos.

Hecho esto el siguiente paso es el de encontrar los diagramas de

momentos flectores y diagrama de solicitación de corte en todos los
elementos de la viga. Para esto se procede con el siguiente análisis.
1, 2, 3, 4 son los grados de libertad de un elemento finito de viga. Sus
correspondientes desplazamientos serán u1, u2, u3, u4.

u1 = desplazamiento vertical en el nudo izquierdo (positivo arriba,

negativo abajo)

u2 = giro en el nudo izquierdo (positivo antihorario y negativo horario)

u3 = desplazamiento vertical en el nudo derecho (positivo arriba,

negativo abajo)

u4 = giro en el nudo derecho (positivo antihorario y negativo horario)

además,
q1 = valor de carga distribuida en el nudo izquierdo

q2 = valor de carga distribuida en el nudo derecho

x1 = coordenada del nudo izquierdo

x2 = coordenada del nudo derecho

EI = módulo elástico x inercia de la sección transversal de la viga.

 Ejemplos
Ejemplo 1 - viga hiperestática de dos tramos

Viga resuelta por método Matricial

Dada la siguiente viga:

Resolver las reacciones de la viga mediante el método Matricial

Solución:

Si se desea aplicar el método matricial de rigidez para resolverla,

debemos primero dividir la viga en tres elementos, efectuando las
divisiones en cada punto de discontinuidad de carga o apoyo, de la
siguiente manera:
Los grados de libertad en los extremos de cada elemento se muestran
en la figura.

La matriz de rigidez de un elemento de viga es la siguiente:

En este caso la sección de la viga es constante. Conviene trabajar en

unidades de KN y m, entonces la inercia, módulo elástico y longitudes
se trabajarán en estas unidades. Se calculan:

I = b*h³/12 = 0.05*0.10³ / 12 m4

E = 13000000 KPa

L = según la longitud de cada elemento

Entonces se reemplazan estos valores en la matriz de rigidez y se

ensamblan los elementos según sus grados de libertad.
La matriz ensamblada es la siguiente:

Se puede ver de la figura, que los grados de libertad restringidos (que

no permiten desplazamiento ni giro) son los correspondientes a U1, U2,
U3, U7 y U8

Además, las cargas equivalentes de las cargas distribuidas son las

siguientes (ver fórmulas de momentos de empotramiento perfecto para
ver del detalle de cómo se obtuvieron los valores mostrados)
Se trasladan luego estos valores al sistema de ecuaciones, junto con las
incógnitas correspondientes a los desplazamientos y a las reacciones
desconocidas.

El sistema de ecuaciones queda así:

Se eliminan las filas y columnas con los grados de libertad restringidos
(1, 2, 3, 7, 8) y se resuelve el sistema de ecuaciones para las incógnitas
U4, U5, U6.

Se tiene entonces el resultado:

U4 = -0.00623077

U5 = -0.009470085

U6 = 0.006358974

Reemplazando estos desplazamientos en el sistema de ecuaciones, se

resuelve R1, R2, R3, R7, R8

De donde, despejando las incógnitas se tiene:

R1 = 4.2750 KN (reacción vertical en el apoyo izquierdo)

R2 = 2.0250 KN-m (momento en el apoyo izquierdo)

R3 = 12.5722 KN (reacción vertical en el apoyo central)

R7 = 5.6528 KN (Reacción vertical en el apoyo derecho)

R8 = -2.1583 KN-m (momento en el apoyo derecho)

El siguiente paso consiste en determinar los diagramas de cortante y
Momento flector en función a las reacciones encontradas.

Ejemplo 2 - Viga de tres tramos

Ejemplo de Viga Hiperestática – Método de Elementos Finitos o

de Rigidez

Enunciado

Dada la viga mostrada, se deben encontrar los desplazamientos

producto de las cargas aplicadas a la viga.
Solución

Paso 1 – Enumeración de nudos y elementos

Antes de enumerar los nudos y los elementos, se debe tener presente

donde colocar los nudos. Estos deben colocarse en todos los puntos
donde exista:

• un apoyo

• un cambio en el tipo de carga

• una carga puntual

• un momento de fuerza

• un desplazamiento inducido

Luego, enumeramos cada grado de libertad, cada nudo y cada elemento.

Nota que, al momento de asignar los grados de libertad a los nudos,
solo se asignan dos por nudo. Uno correspondiente a los
desplazamientos verticales y otro para los giros.

En este caso no se necesitan desplazamientos horizontales pues los

elementos de viga, a diferencia de los pórticos, podemos asumir que no
se deforman axialmente

Paso 2 – Se arman las matrices de rigidez de cada elemento

El método matricial consiste en el desglose de la viga en pequeñas

partes. Cada una de estas partes o tramos de la viga contiene
propiedades que matemáticamente y en notación matricial pueden
escribirse según:
Esta es llamada Matriz de Rigidez de un elemento. Esta matriz de
rigidez relaciona las fuerzas y momentos aplicadas en los nudos con
los desplazamientos y giros de estos mismos nudos mediante la
ecuación:

Gráficamente, sucede lo siguiente:

Entonces, con los datos correspondientes a cada tramo de viga (cada
elemento), se deben encontrar las matrices de rigidez correspondientes.
Existen 4 tramos, por tanto 4 elementos, y 4 matrices de rigidez
elementales:

Por ejemplo, el elemento 1 tiene una longitud L = 4m; EI = rigidez =

3000 KNm². Se reemplazan estos datos en la matriz 4×4 de arriba,
resultando:

Nota: Cuando generes una matriz de rigidez elemental, es buena

práctica anotar los grados de libertad a los que corresponde esta matriz
elemental. Esto te ayudará a ensamblar la matriz de rigidez global más
adelante. Por ejemplo, la segunda matriz de rigidez va del nudo 2 al 3.
Los grados de libertad del nudo 2 y 3 respectivamente son 3,4,5,6 y
estos grados de libertad deberás anotarlos ordenadamente sobre y al
lado de la matriz de rigidez elemental, como ves en números vertes en
la imagen arriba.

Te dejo nuevamente la imagen de grados de libertad acá abajo:

Paso 3 – Ensamblaje de la Matriz de rigidez

Para este paso, imagina que existe una matriz vacía de dimensión
10×10 (la misma cantidad de grados de libertad que existen en toda la
viga) donde iremos colocando las matrices elementales una a una,
según la posición que le corresponda.

Por ejemplo, trasladamos la matriz de rigidez del elemento 1 (de color

azul) hacia la matriz vacía de 10×10.
Luego trasladamos la matriz elemental del segundo elemento (pintada
de magenta). Notarás que los grados de libertad que le corresponden en
la matriz global son 3,4,5,6 y se superponen con la primera matriz. En
caso de tener matrices de rigidez elementales superpuestas, su s
elementos se suman.

Se repite el procedimiento para los otros dos elementos

La ecuación generada [K global]{U global}={F global} con todos los

elementos sumados es:
Esta ecuación aun no puede resolverse pues es singular (infinitas
soluciones). Se encontrará una única solución al problema en el
momento en que se apliquen las condiciones del contorno del problema.
En eso consiste el siguiente paso.

Paso 4 – imposición de condiciones de contorno y cargas

Vamos a volver a dibujar la viga con sus condiciones de contorno:

De esta figura, los grados de libertad restringidos son los respectivos a

los apoyos que impiden que la viga gire o se desplace en ese sector.
Por ejemplo, el empotramiento impide que los grados de libertad 1, 2
se desplacen o giren, por tanto, U1 y U2 son cero.
Lo mismo sucede con los grados de libertad 3 y 5. No pueden
desplazarse verticalmente por tanto son cero.

Luego, donde los desplazamientos son cero, las reacciones son

desconocidas o incógnitas.

Para cargar la viga, se deben transformar las cargas distribuidas en

cargas equivalentes puntuales en los extremos del tramo.

Estas cargas equivalentes salen de los Momentos de empotramiento

perfecto. Se pueden deducir o pueden consultarse en bibliografía.

Para cargas triangulares y uniformes, las cargas y momentos

equivalentes en los extremos son los siguientes:
Entonces, cada una de estas cargas equivalentes debe cargarse en el
vector de fuerzas y momentos en el sistema de ecuaciones lineal en el
grado de libertad correspondiente.

Las cargas puntuales como la carga de 30 KN no necesita transformarse

a carga equivalente pues esta ya es una carga aplicada directamente
sobre un nudo.

Así, por ejemplo, para el grado de libertad 3, las cargas que concurren
en este grado de libertad son:

• La carga equivalente de 28 KN producto de la carga triangular

• La carga equivalente de 50 KN producto de la carga distribuida

• La reacción R3 producto del apoyo.

Estas tres fuerzas se anotan en el vector de fuerza en el grado de

libertad 3.

Lo mismo se repite tanto para momentos como para fuerzas en los

grados de libertad que correspondan, quedando, por tanto:
Terminado este paso, el sistema de ecuaciones [K]{u}={F} queda
como en la siguiente imagen:

Resolución del sistema de ecuaciones

Para resolver el sistema, se debe previamente eliminar del sistema, las
filas y columnas correspondientes a los grados de libertad restringidos.
En este caso filas y columnas 1, 2, 3 ,5.

De esta manera nos quedan como incógnitas solamente los

desplazamientos desconocidos:

Esta eliminación de filas y columnas podría parecer arbitraria, sin

embargo, no la es. Se eliminan las columnas 1, 2, 3, 5 ya que al
multiplicar [K] por {u}, los coeficientes que acompañan a los
desplazamientos iguales a cero, obviamente se eliminan en la
multiplicación.

El sistema de ecuaciones reducido mostrado ahora se resuelve de la

manera clásica, mediante inversión de la matriz de rigidez reducida, o
mediante metodos como Gauss-jordan o cualquier método similar.

Aplicando las fórmulas mostradas a nuestro ejercicio se tiene entonces:

Cuyo resultado será el siguiente:

Los desplazamientos encontrados están en unidades de [m] para los U’s

impares y en radianes para los U’s pares

Reacciones

Si se quieren encontrar las reacciones R1, R2, R3, R5 solo se deben

reemplazar los desplazamientos encontrados justo arriba y
reemplazarlos en el sistema de ecuaciones global de 10×10. De esa
manera se tiene:
Multiplicando los dos miembros de la izquierda, se obtiene un vector
de 10×1 que se iguala con el vector 10×1 de la derecha:

Cada término de la izquierda se iguala con cada término de la derecha.

En las filas donde existen las incógnitas se despejan las respectivas
reacciones. En las filas donde no existan incógnitas simplemente se
verifica la identidad.

Y con esto finaliza la solución de las reacciones de la viga.

Reescribiendo las reacciones en la viga se tiene:
¿Qué sigue?

Ejemplo - Líneas de influencia de viga Hiperestática

Ejemplo: Líneas de influencia - Método matricial

Dada la siguiente viga:

Figura 1

El ejercicio consiste en obtener las líneas de influencia de todos los

apoyos y además la línea de influencia de cortante y de momento flector
en el punto H de la viga.

Solución
Primero se procede a dividir la viga en puntos convenientes donde
consideremos que es necesario que se apoye la carga unitaria. Este paso
puede realizarse a criterio, considerando que mientras más divisiones
tenga la viga, existirá más coste computacional.

Entonces, dividiendo la viga de tal forma que los apoyos y el punto H

coincidan con nudos tenemos:

La enumeración de los elementos es correlativa:

Y la enumeración de grados de libertad de todos los elementos es la

siguiente:

Como las divisiones son equidistantes entre apoyos, existen cinco

grupos de elementos de diferentes longitudes, para los cuales deberán
armarse las respectivas matrices de rigidez.

La fórmula de cada matriz de rigidez va como sigue:

Entonces, se tienen las siguientes matrices para los cinco tipos de
elementos:

Cada una de estas matrices de rigidez viene acompañada de sus

respectivos grados de libertad marcados en rojo. Estos grados de
libertad ayudan a armar la matriz de rigidez global.

Finalmente ensamblada la matriz de rigidez y armado el sistema de

ecuaciones se tiene:
Al lado derecho del sistema de ecuaciones se encuentran los vectores
de carga. Por ser este un caso de cargas móviles, se acomodaron lado
a lado los vectores de carga para cada caso de carga (cada salto de la
carga unitaria al siguiente nudo de izquierda a derecha). El arreglo de
vectores de carga mostrado nos sirve para resolver todos los casos de
carga de manera simultánea.

Luego se eliminan filas y columnas de la matriz de rigidez de los grados

de libertad restringidos (5, 13, 25, 33) pues en estos corresponden a los
desplazamientos verticales en los apoyos.

Se resuelve el sistema mediante cualquier método de r esolución de

sistemas de ecuaciones lineales (por ejemplo, Newton Rapson o Jacobi)
y se obtienen los desplazamientos en un arreglo matricial para cada
estado de carga, de la siguiente manera:
Una vez conocidos los desplazamientos en cada nudo, se pueden
reemplazar en el sistema de ecuaciones armado con todos los grados de
libertad siguiente:

Entonces multiplicando la matriz [K] por {u} se obtienen vectores de

las cargas y reacciones ya resueltos.

En los campos donde no exista una Reacción “R” como incógnita, el

vector de la izquierda será igual al vector equivalente de la derecha.
En los campos donde exista una reacción “R” incógnita, se deberá
despejar esta reacción para conocer el resultado de la reacción final.

Se tiene, así:

En resumen, si por ejemplo se deseara dibujar la Línea de Influencia

de reacciones, se puede generar la siguiente tabla resumen:

Y graficando todas las líneas de influencia en una sola gráfica, se

tiene finalmente:
 Ejemplo - Viga con asentamiento impuesto en apoyo

Esfuerzos en viga hiperestática con Asentamiento en un apoyo

Este ejercicio es algo largo, sin embargo bastante útil para el alumno
que esté cursando la materia de análisis matricial de estructuras, ya
que existen pocos ejemplos explicados en toda su extensión,
especialmente cuando se trata de apoyos que se desplazan y pueden
generar esfuerzos producto del desplazamiento del nudo.

La teoría necesaria para entender este ejercicio es la de análisis

matricial en vigas, cuyas fórmulas y lo más relevante de la teoría lo
puedes encontrar en Deducción de la matriz de rigidez de viga
Luego de encontrados los desplazamientos de cada tramo de viga,
necesitarás conocer cómo realizar el postproceso de datos para
encontrar diagramas de cortante y momento a partir de los gi ros y
desplazamientos verticales en los extremos. Para esto necesitarás leer
el siguiente artículo: Diagramas de M y V a partir de giros y
desplazamientos en extremos de un tramo de viga.

Con todo el arsenal listo, ya podrás entender el siguiente ejemplo de

análisis

https://1.bp.blogspot.com/-
U7mV_7KF2RI/W9PK0eSsr7I/AAAAAAAACjM/L0IIYECxQEI4
YUk2sFy8j-nhhUdpjalPgCLcBGAs/s000/pizarra%2B1.jpg

Ahora viene el post proceso donde a partir de estos desplazamientos

se encuentran las solicitaciones en la viga.

https://3.bp.blogspot.com/-
y91rcWcxflo/W9PMfUspFuI/AAAAAAAACjc/-
YO3ItAkS8kyrUPEBqtzyQigMcPnG90CACLcBGAs/s000/pizarra%2
B2.jpg

Ejemplo - Viga apoyada sobre resortes

Resolución de Viga con apoyos de resortes - Método Matricial

En el siguiente ejemplo se resuelve un caso de viga apoyada en

resortes. Una vez entendido el concepto se podría incluso resolver
vigas de fundación apoyadas sobre lecho elástico.

Con la teoría a mano, se puede resolver el ejercicio cuyo enunciado

dice:

Para resolver el ejercicio...

https://2.bp.blogspot.com/-
LEzopIUmgBM/W_TGq4UyPiI/AAAAAAAACrY/ksZhIxLd3zcwlzjc
FHr9ZXvA7i7NgH56QCLcBGAs/s000/viga%2Bcon%2Bresortes.jpg

5 análisis de pórticos

Matriz de Rigidez de pórtico en 2D

Matriz de Rigidez de Pórtico en 2 dimensiones. Método de rigidez

Introducción

Esta ya se trata de una materia de tercer o cuarto año de la carrera de

ingeniería civil, del área de estructuras

Concretamente, necesitas conocimientos de resolución de cerchas o de

vigas por métodos de Rigidez o también conocidos como métodos
matriciales para entender esta publicación.

Lo que pretendo mostrar en esta publicación en concreto es la

demostración de la matriz de rigidez de un pórtico en coordenadas
globales

Matriz de rigidez de pórtico en coordenadas locales

Los elementos de pórtico en 2D no son otra cosa que una combinación

de dos tipos de elementos:

• a) Elementos con deformación axial

• b) Elementos de viga.
Los elementos de cercha, como se vio anteriormente, para los grados
de libertad mostrados tiene la siguiente matriz de rigidez local.

A la vez, vimos que los elementos de viga tienen la siguiente matriz de

rigidez:

Si combinamos ambos efectos en un elemento de barra que contenga

tanto deformaciones axiales, verticales y giros en los nudos,
sobreponiendo efectos de los dos anteriores, se tiene:

… donde se han re enumerado los grados de libertad desde el 1 hasta

el 6 para combinar ambos efectos.

Como se puede ver en la última imagen, la última matriz consiste en la

matriz de rigidez de un elemento de pórtico en coordenadas locales. Se
combinaron ambos efectos tanto de deformación axial como de
deformación transversal para conformar un elemento cuyos nudos se
trasladen en x,y y roten.

Esta matriz de rigidez multiplicada por los desplazamientos es igual al

vector de fuerzas y reacciones del elemento en coordenadas locales.
Ésta es la conocida fórmula del sistema de ecuaciones a resolver por
elementos finitos:

(ecuación 1)

Matriz de transformación de coordenadas

Para poder rotar este elemento de pórtico a cualquier dirección se

necesitan utilizar matrices de rotación. Entonces deberemos poder
convertir los desplazamientos y fuerzas de coordenadas locales a
coordenadas globales.

Para lograrlo debemos concebir ambos sistemas como se muestra en la

siguiente figura:
si queremos convertir por ejemplo el desplazamiento U1 local al
sistema global debemos expresar U1 local como:

Debemos proceder de la misma manera con los grados de libertad U2,

U4, U5 locales. Los grados de libertad U3 y U6 no necesitan
transformación pues corresponden a giros.

Terminado el proceso de transformación y escribiendo las ecuaciones

en sistema matricial se obtiene la siguiente matriz de transformación:
Donde a esta última matriz de senos y cosenos de 6×6 llamaremos
matriz de transformación [T]

Matriz de rigidez de pórtico en coordenadas globales

Tanto para los desplazamientos como para las fuerzas, podemos

escribir los vectores locales en función de los globales de la siguiente
manera:

Y reemplazando {U local} dentro del sistema de ecuaciones 1 que lo

volvemos a mostrar a continuación,

, se tiene:

Si multiplicamos toda esta última ecuación por la inversa de [T] a

ambos lados, se anula el término [T] del lado derecho de la ecuación,
y debido a que para matrices ortogonales como lo es la matriz [T], la
inversa es igual a la transpuesta, finalmente se tiene el sistema de
ecuaciones final:

donde la matriz de rigidez en coordenadas globales es:

Desglosando esta matriz en todas sus componentes se tiene:

Ésta última es la matriz de rigidez en coordenadas globales para

elementos de pórtico en 2 dimensiones.

Con esta matriz resultante de 6×6 elementos se tienen por tanto todos
los datos físico mecánicos de un solo elemento de pórtico, que deberá
ensamblarse con todos los otros elementos de pórtico de la estructura
para obtener la matriz de rigidez global de toda la estructura y
posteriormente resolver los desplazamientos.

Ejemplo

El siguiente paso para entender el uso de esta matriz es mediante un

ejemplo completo. Se sugiere el siguiente ejemplo con cargas puntuales
en los nudos como primer ejemplo de aprendizaje de cálculo de
pórticos por el método de rigidez:

– Ejemplo de pórtico con cargas puntuales – 1ra parte

Ejercicio de Resolución de Pórtico por método de Rigidez (Parte 1)

Planteamiento del problema

Se plantea resolver el siguiente problema mediante el método de

rigidez o también llamado método matricial.
Las características del problema denotan lo siguiente:

• El pórtico es hiperestático

• Existen elementos de pórtico inclinados

• Solo existen cargas puntuales en el enunciado

• Las secciones transversales son constantes

Solución

Paso 1) Numeración de nudos, grados de libertad y barras

En este proceso es conveniente ir de izquierda a derecha en la
numeración de todos los elementos mencionados.

Para la numeración de grados de libertad debe notarse que cada nudo

consta de tres grados de libertad: Desplazamiento en X,
Desplazamiento en Y, Giro.

Paso 2) Datos de cada barra

Los datos necesarios de cada barra para ingresarlos a la matriz de

rigidez de cada elemento son:

• Longitud
 • Inclinación

• Módulo elástico

• Sección transversal (área)

• Inercia

En ese contexto se deben realizar algunos cálculos previos para

encontrar los datos necesarios. Los aspectos a tomar en cuenta son los
siguientes:

• El ángulo de inclinación debe medirse SIEMPRE desde el eje

positivo X, acomodando el nudo de menor numeración en el eje
de coordenadas.

• La sección transversal de las barras debe obtenerse en m², ya que

las cargas están en S.I.

• La inercia debe obtenerse en m4, por el mismo motivo anterior

• El módulo elástico debe calcularse en KPa, ya que las cargas

están en KN

Con esas consideraciones, se tiene lo siguiente:

Paso 3) Se calcula la matriz de rigidez [Ke] de cada elemento

Para calcular la matriz de rigidez de cada elemento se deberá aplicar la

fórmula ya deducida en la publicación de Teoría de matríz de rigidez
de pórticos. La matriz de rigidez deducida allí, se la muestra a
continuación:
Toma en cuenta de esta matriz, que c=cos(α) y s=seno(α). El resultado
de la multiplicación de estas tres matrices será la matriz de rigidez del
elemento, en coordenadas globales.

Calculando entonces cada uno de los elementos en base a la tabla

resumen de datos de cada barra y la matriz de rigidez de la última
imagen, calculamos:
Un aspecto a tomar en cuenta es la numeración en celeste de cada
matriz de rigidez elemental. Estos números corresponden a los grados
de libertad de cada barra. La numeración de los grados de libertad para
cada matriz elemental debe anotarse desde el menor grado de libertad
de esa barra hasta el mayor. Por ejemplo, la barra 4 arranca en el nudo
3 (nudo de menor numeración) con grados de libertad 7 8 9 y termina
en el nudo 5 (nudo con mayor numeración) con grados de libertad 13
14 15. Estos valores van anotados como índices en la matriz de rigidez
de este elemento.

Paso 4) Ensamblaje de la matriz de rigidez global [Kg]

Para ensamblar las matrices de rigidez elementales dentro de la matriz

de rigidez global se debe tomar en cuenta los siguientes pasos:

• Crear una matriz de ceros con dimensión cuadrada igual al

número de grados de libertad del sistema (en este caso 18×18)

• Trasladar cada elemento de las matrices de rigidez elemental a la

global según los índices en celeste anotados arriba y a la
izquierda. Estos índices corresponden a la posición de estos
elementos en la matriz de rigidez global.

• Si en el proceso anterior, dos valores cayeran dentro de la misma

casilla en la matriz de rigidez global, estos valores deben
sumarse.

Con todas las consideraciones necesarias mencionadas arriba, se

obtiene entonces la matriz de rigidez siguiente:
Paso 5) Sistema de ecuaciones y condiciones de contorno

Ensamblada la matriz de rigidez, se arma el sistema de ecuaciones a

resolver siguiente:
Al sistema de ecuaciones lo acompaña un esquema de la estructura con
los grados de libertad mostrados, ya que a partir de acá se impondrán
las condiciones de contorno al sistema de ecuaciones.

Si dibujamos las incógnitas de equilibrio externo de la estructura y sus

cargas, se tiene lo siguiente:
De la imagen con las fuerzas y reacciones mostradas, podemos
comenzar a imponer las condiciones de contorno al sistema de
ecuaciones. Para eso seguimos los siguientes pasos:

• Donde exista apoyo (restricción de desplazamiento o giro) ese

grado de libertad será igual a cero.

• Donde no exista carga ni restricción (apoyo), la fuerza será igual

a cero
 • Donde exista apoyo, para ese grado de libertad existirá siempre
una reacción

• Donde se conocen los desplazamientos, las fuerzas (o reacciones)

son desconocidas

• Donde se desconozcan los desplazamientos, las fuerzas son

siempre conocidas

Paso 6) Sistema de ecuaciones reducido

Una vez que se hayan impuesto todas las condiciones de contorno en

el sistema de ecuaciones se deben anular las filas y columnas de los
grados de libertad restringidos, o sea los grados de libertad cuyos
desplazamientos sean iguales a cero. La siguiente figura muestra este
proceso
Se conforma entonces un sistema de ecuaciones reducido donde todas
las incógnitas están al lado izquierdo del sistema de ecuaciones y todo
las fuerzas conocidas están al lado derecho, y de esta manera se puede
resolver el sistema de ecuaciones.

Paso 7) Resolución del sistema de ecuaciones

No es recomendable utilizar una resolución invirtiendo la matriz

cuadrada pues esto genera errores numéricos que pueden no ser
convenientes y además el uso de procesador es muy alto.

Es mejor aplicar métodos numéricos especializados para la resolución

de sistemas de ecuaciones como por ejemplo el método de Gauss
Jordan.

Realizando Entonces esta resolución nos quedan los siguientes

resultados de desplazamientos:
Paso 8) Reacciones en la estructura

Una vez conocidos los desplazamientos del paso anterior se pueden

encontrar las reacciones del sistema reemplazando los desplazamientos
conocidos en el sistema de ecuaciones original antes de ser reducido,
como se muestra en la imagen:

Hecho esto se multiplica la matriz cuadrada de 18×18 por el vector de

desplazamientos de 18×1 quedando la igualdad mostrada en la
siguiente imagen. En la imagen a continuación se puede ver una
igualdad de dos vectores de donde sólo nos sirve rescatar las posiciones
1 2 3 16 17 y 18. Las demás posiciones desde la fila 4 hasta la fila 15
Simplemente nos van a dar una identidad que corrobora que el sistema
ha sido bien resuelto.

La estructura entonces queda con las reacciones de la siguiente manera

Paso 9) Postproceso

A partir de este momento queda calcular las solicitaciones internas de

las vigas y columnas de la estructura mediante teoría de la elasticidad,
que no se analizará en esta publicación sino en la siguiente publicación:

– Ejemplo de pórtico con cargas puntuales – 2da parte

SOLICITACIONES DE PÓRTICO CON CARGAS PUNTUALES

– MÉTODO DE RIGIDEZ (2da parte)
Enunciado

Se pide encontrar las solicitaciones de momento flector, cortante y

axial del siguiente pórtico hiperestático:
La numeración de nudos, grados de libertad y barras se muestran a
continuación:
Los resultados de desplazamientos de los nudos producto del análisis
Estático por el método de rigidez son los siguientes:
Pasos previos

Para llegar a los desplazamientos de los nudos mostrados en el

enunciado se siguió un procedimiento de resolución mediante el
método de rigidez.

En este artículo se aplican conocimientos de matrice s de rigidez

elementales y globales que utilizaremos también en este artículo.

Posproceso

Suelo llamar «POSTPROCESO» a todo el cálculo de solicitaciones

posteriores al cálculo de desplazamientos y reacciones. Esto se debe a
que concretamente el método de rigidez para encontrar los resultados
del problema termina en el anterior artículo mencionado. A partir de
aquí aplicamos simplemente conocimientos de resistencia de
materiales y de equilibrio estructural para encontrar las solicitaciones
en cada barra.

Coordenadas locales

Los desplazamientos encontrados en toda la resolución del anterior

artículo corresponden a desplazamientos de los nudos en coordenadas
globales. El primer paso para conocer las solicitaciones de cada barra
es transformar estos desplazamientos globales en desplazamientos
locales alineados con el eje de la barra, de la siguiente manera. Por
ejemplo, para la barra 1:
Si bien las barras 2 y 3 no necesitan tal transformación pues su eje
coincide con el eje de las X y por tanto los desplazamientos en
coordenadas locales coinciden con las globales. Por tanto, aplicaremos
la transformación a ejes locales para las barras 1, 4 y 5.

La transformación de coordenadas globales a coordenadas locales se

ejecuta entonces de la siguiente manera:

Aplicado al elemento 1 (y a todos los otros elementos en general) se

obtiene:
Por ejemplo, para el elemento 1, una vez ejecutada la rotación de los
desplazamientos, el elemento queda como en la siguiente imagen. Se
puede a partir de la siguiente imagen, determinar cuáles de los
desplazamientos generan compresión en la barra, y cuáles de los
desplazamientos generan flexión y cortante.

Habiendo desacoplado los desplazamientos separándolos entre los

efectos que generan deformación axial y los que generan deformación
por cortante y flexión, se puede determinar cuál es la solicitación axial
de la barra y cuáles son las solicitaciones de cortante y flexión.

Solicitación axial

A partir de conceptos de resistencia de materiales, se puede determinar

si la barra está a compresión o a tracción en base a la diferencia de
desplazamientos axiales u1 local y u4 local. Además, como la
diferencia de ambos desplazamientos es la deformación, y dividiendo
ésta entre la longitud obtenemos la deformación unitaria, entonces:

De las fórmulas arriba, el esfuerzo axial es por definición igual a la

deformación unitaria por el módulo elástico, y además el esfuerzo axial
se descompone en Fuerza Axial entre la Sección transversal. Por tanto,
despejando la fuerza Axial P se obtiene la fórmula de solicitación axial
de la barra:

Reemplazando los datos para este elemento, se obtiene la solicitación

para la barra. nota que el signo del resultado automáticamente denota
si la barra está a compresión (signo negativo) o si está a tracción (signo
positivo). Además, las unidades de los valores introducidos a la
fórmula son, u1 y u4 en [m], Área en [m²], E en [KN/m²], Longitud de
la barra en [m].
Repetimos el mismo procedimiento para todas las otras barras. Sin
excepción el proceso es idéntico. Obtenemos entonces el resultado
siguiente para las otras barras, y graficamos los resultados.

Solicitaciones de Momento Flector y Cortante

Para conocer las solicitaciones de Momento y Cortante de las barras

continuamos trabajando con los desplazamientos en coordenadas
locales de cada barra.

De los 6 desplazamientos en coordenadas locales, desechamos los

desplazamientos axiales que se utilizaron en el anterior apartado.
Analizamos entonces los desplazamientos locales u2, u3, u5 y u6.

En este procedimiento mencionado los desplazamientos

perpendiculares a la barra son U1 y U3, que nosotros reemplazaremos
por los desplazamientos locales u2 y u5. De la misma manera los
desplazamientos de la teoría U2 y U4 los reemplazaremos por las
rotaciones locales u3 y u6. Queda entonces:

En el sistema de ecuaciones mostrado existen variables que para el caso

de las barras analizadas en coordenadas locales valen X1=0 y
X2=longitud de la barra. Los coeficientes A y B para este ejercicio son
cero ya que las cargas distribuidas son nulas.

Las constantes C1 C2 C3 y C4 obtenidas corresponden a las constantes

de integración de la ecuación de la elástica que sirve para encontrar los
coeficientes del polinomio de las ecuaciones de cortante y Momento
flector.

En resumen, el cortante y momento flector responden a las siguientes

ecuaciones:

M y V Barra 1
Trabajamos con los desplazamientos locales de la barra 1, donde A=0
y B=0 porque no existen cargas distribuidas. Además, x1=0 y
x2=2.7459; A=0.0029m2; I=0.000010174m4; E=200e6.
Reemplazando todos estos datos se obtienen las constantes:

Estas constantes son los coeficientes de la ecuación de V y M:

M y V para las demás barras

La única diferencia entre las otras barras y la barra analizada acá es que
en las otras barras los desplazamientos del nudo 1 no son cero porque
ya no estan empotrados. Por ejemplo para la barra 3:
Finalmente repitiendo el procedimiento para todas las otras barras y
graficando:
Con esto termina el procedimiento del análisis de solicitaciones del
pórtico.

Algoritmo de cálculo de pórticos por método matricial

Un pantallazo global de los pasos a seguir para calcular de principio a

fin un pórtico en 2 dimensiones.El algoritmo es el siguiente:
Leyenda:

E = Módulo Elástico de cada elemento

I = Inercia en dirección de la flexión de cada elemento

A = sección de cada barra

L = longitud de cada elemento

[K]rotada = Matriz de rigidez de cada elemento rotada a coordenadas

globales

[K]local = Matriz de rigidez de cada elemento en coordenadas locale s

[ß] = Matriz de transformación de coordenadas

[K]global = Matriz ensamblada producto de todas las matrices locales

{u} = Vector de desplazamientos globales

{Q} = Vector de cargas y reacciones

[K]reduc = Matriz de ridigez global con los grados de libertad

restringidos eliminados

{u}local = Vector de desplazamientos en coordenadas locales para cada

elemento, alineando los desplazamientos axiales paralelos a la barra y
de flexión perpendiculares a la barra