Método Matricial
Método Matricial
Método Matricial
1. Conceptos básicos
Concepto
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)
Existen muchas estructuras que no necesitan más que unos pocos grados de
libertad por nudo para analizarse satisfactoriamente.
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:
https://3.bp.blogspot.com/-
hwbvEsFsF4o/W2z8JpV7J5I/AAAAAAAACXY/Io7xJBIwrYQw6
TCm7thMc7M041crK4WBwCLcBGAs/s000/elementos%2Bcon%2
B1GDL.jpg
3. Análisis de cerchas
Teoría
Ya que, por equilibrio, la fuerza del nudo uno debe ser igual y opuesta a la
fuerza del nudo 2,
(ecuación 1)
(ecuación 3)
Donde
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}
siendo (x1 ; y1 ) las coordenadas iniciales del nudo y (x2 , y2)las coordenadas
finales.
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
Planteamiento
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.
Reacciones en la cercha
https://3.bp.blogspot.com/-
uvYk8C1wQ6I/W4lZOapc0cI/AAAAAAAACao/BxYcu_jo4sY6X77Mi
vMJnu61-JGGGtfaQCLcBGAs/s000/Cercha.jpg
https://2.bp.blogspot.com/--
Nd9M2oiBvk/W46tsLmuT3I/AAAAAAAACb8/xY4xTS6JVrEd92YIvQJ
vkitASV09vE0tgCLcBGAs/s000/cercha5.jpg
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
https://3.bp.blogspot.com/-BQRYM36Pm_s/W-
RkVSd7ThI/AAAAAAAACmc/93M5FpoDETQz2norMsCrkQ5mUWcag3
_CACLcBGAs/s000/desarrollo.jpg
4. Análisis de vigas
Teoría
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:
además,
q1 = valor de carga distribuida en el nudo izquierdo
Solución:
I = b*h³/12 = 0.05*0.10³ / 12 m4
E = 13000000 KPa
U4 = -0.00623077
U5 = -0.009470085
U6 = 0.006358974
Enunciado
• un apoyo
• un momento de fuerza
• un desplazamiento inducido
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.
Así, por ejemplo, para el grado de libertad 3, las cargas que concurren
en este grado de libertad son:
Reacciones
Figura 1
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.
Se tiene, así:
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.
https://1.bp.blogspot.com/-
U7mV_7KF2RI/W9PK0eSsr7I/AAAAAAAACjM/L0IIYECxQEI4
YUk2sFy8j-nhhUdpjalPgCLcBGAs/s000/pizarra%2B1.jpg
https://3.bp.blogspot.com/-
y91rcWcxflo/W9PMfUspFuI/AAAAAAAACjc/-
YO3ItAkS8kyrUPEBqtzyQigMcPnG90CACLcBGAs/s000/pizarra%2
B2.jpg
https://2.bp.blogspot.com/-
LEzopIUmgBM/W_TGq4UyPiI/AAAAAAAACrY/ksZhIxLd3zcwlzjc
FHr9ZXvA7i7NgH56QCLcBGAs/s000/viga%2Bcon%2Bresortes.jpg
5 análisis de pórticos
Introducción
• 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.
(ecuación 1)
, se tiene:
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 pórtico es hiperestático
Solución
• Longitud
• Inclinación
• Módulo elástico
• Inercia
Posproceso
Coordenadas locales
Solicitación axial
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:
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.