MATRIZ DE RIGIDEZ PARA UN RESORTE 1
DEFINICION DE MATRIZ DE RIGIDEZ
HERAS QUINTERO CARLOS FRANCISCO
17171139
INSTITUTO TECNOLOGICO DE MEXICO
CAMPUS CULIACAN
INGENIERIA ASISTIDA POR COMPUTADORA
M.C. DAGOBERTO TOLOSA MATA
MATRIZ DE RIGIDEZ PARA UN RESORTE 2
DEFINICION DE MATRIZ DE RIGIDEZ
Para una estructura de elementos finitos, la matriz de tigidez contiene la informacion
geometrica y del material que indica la resistencia del elemento a la deformacion cuando esta
sometido a una carga. Para elementos finitos aplicado en analisis no estructurales, como analisis
de flujos o transferencia de calor, el termino matriz de rigidez tambien se usa. Esto porque la matriz
representa la resistencia del elemento a cambiar cuando se le aplican fuerzas externas.
Los pasos necesarioa para resolver una estructura según el metodo de rigidez comienzan
con la definicion de la geometria y de las acciones sobre la estructura. Se identifican, a
continuacion, los movimientos incognita teniendo en cuenta las condiciones de compatibilidad y
se pasa a resolver las piezas individuales en funcion de los movimientos de sus extremos. El
analisis continua estableciendo las condiciones de equilibrio de la estructura, es decir, realizando
el ensamblaje de las matrices elementales.
Una vez calculadas las ecuaciones elasticas de todas las piezas que forman la estructura en
un sistema de referencia comun, el siguiente paso en el metodo de rigidez es la construccion de la
matriz global de la rigidez de la estructura. Esta matriz global se obtiene medianto el ensamblaje
de las matrices elementales de rigidez de las piezas.
Dicho ensamblaje se realiza considerando el equilibrio de fuerzas (y momentos, en su caso)
que actuan sobre cada nudo de la estructura.
A modo de ejemplo para explicar el proceso a seguir, considerese la sigueinte estructura.
Formada por 4 barras y 4 nodos con la numeracion que se indica. Los vectores que representan las
fuerzas exteriores (cargas en los nudos) y los entidos de recorrido de las piezas pueden verse en la
figura.
MATRIZ DE RIGIDEZ PARA UN RESORTE 3
Las ecuaciones de equilibrio de fuerzas interiores (esfuerzos en extremos de piezas) y
exteriores que actuan sobre los nodos pueden escribirse, en el sistema global de referencia:
MATRIZ DE RIGIDEZ PARA UN RESORTE 4
Estas ecuaciones pueden escribir se en forma matricial:
O de forma mas compacta: f=Kd
Donde f es el vector global de las fuerzas exteriores (incluidas las reacciones), K es la
matriz global de rigidez de la estructura y d es el vector global de movimientos en los nudos. De
lo anterior se obtienen algunas conclusiones de validez general:
1. La matriz K es cimetrica, al serlo las matrices kaa y kbb elementales y traspuestas una
de la otra las matrices kab y kba de cada pieza.
2. La matriz de rigidez K es singular, y que las ecuaciiones de equilibrio no se ven
afectadas por un movimiento solido rigido de la estructura.
3. El sistema de ecuaciones tiene solucion, a pesar de ser K singular, ya que f 2, f4, d1, y d4
son conocidos, siendo f1, f3, d2 y d4 las incognitas. Los vectores d1 y d3 son conocidos
por ser justamente los movimientos prescritos por las condiciones de apoyo de la
estructura. Estas condiciones de apoyo impiden los movimientos de solido rigido de la
estructura y matematicamente, evitan la singularidad de la matriz K, haciendo el
problema resoluble en d2 y d4.
Al ser la matriz global de rigidez K siempre simetrica, solo es necesario considerar el
almacenamiento y el ensamblaje de los terminos de la diagonal principal y los situados por
encima (o por debajo) de ella, es decir, el triangulo superior (o el inferior) de la matriz de
MATRIZ DE RIGIDEZ PARA UN RESORTE 5
la estructura. Dicho triangulo superior puede ensamblarse siguiendo dos algoritmos
alternativos:
a) Nodo a nodo. Se ensambla el grupo de ecuaciones correspondientes a cada nudo,
considerando estos, uno a uno de forma consecutiva. Esto equivale a ensamblar la
matriz global “fila a fila“ siguiendo las siguientes reglas :
El elemento de la diagonal principal de la final i-énesima es la suma de las matrices
Kaa ó Kbb de todas las piezas que concurren al nodo i y a la matriz Kbb si es el
extremo b el que coincide con el nodo i.
Los demás elementos de la fila i-ésima, correspondientes a las columnas j ≥ i, son
las matrices Kab de las piezas que unen el nudo i con los nudos j, respectivamente.
Si no existe la pieza ij, el elemento a ensamblar es una matriz nula 0. Debido a la
simetría de la matriz, sólo es necesario considerar las piezas cuyo extremo de
numeración menor es el nudo i.
b) Pieza a Pieza. Se ensamblan las submatrices elementales correspondientes a cada pieza,
considerando éstas una a una, de forma consecutiva. Esto equivale a ensamblar la
matriz de rigidez de la barra ij (que une el nudo i con el nudo j, i<j), en las siguientes
posiciones de la matriz global.
El algoritmo b), ensamblaje pieza a pieza, tiene la ventaja de que cada una de las piezas
sólo es considerada una vez; esto permite calcular las matrices elementales una a una, antes de
ensamblarlas, eliminando la necesidad de almacenarlas. Esta ventaja es importante en problemas
de gran tamaño.
Observando el proceso de formación de la matriz de rigidez, se deduce fácilmente que su
forma varía con la numeración elegida para los nudos. Por tanto, aunque en principio se dijo que
la numeración de los nudos podía ser arbitraria, lo cierto es que incide de manera directa en la
eficacia y rapidez del análisis.
MATRIZ DE RIGIDEZ PARA UN RESORTE 6
Se llama ancho de banda de la matriz global (bw) a la mayor diferencia de numeración
existente entre los nudos extremos ij de las piezas de la estructura.
Es obvio que la numeración óptima es la que corresponde al ancho de banda menor. Cuanto
menor sea la diferencia de numeración entre los extremos de una barra, más estrecha será la banda
de la matriz, es decir, más agrupados se encontrarán los términos no nulos y menor será el número
de operaciones necesarias para el cálculo, resultando el análisis más eficiente.
MATRIZ DE RIGIDEZ PARA UN RESORTE 7
MODULO DE ELASTICIDAD.
Si aplicamos una fuerza F a una barra de longitud l0 el material se deforma
longitudinalmente y se alarga l-l0.
La razon de proporcionalidad entre el esfuerzo (fuerza por
unidad de área) y de froamcion unitaria (deformacion por unidad de
longitud) esta dada por la constante E, denominada modulo de young,
que es caracteristico de cada material.
La ley de hooke relaciona la deformación εx de una barra sometida a esfuerzo axial, con
la tension normal generada por dicho esfuerzo σx, mediante la constante E que se denomina
modulo de elasticidad lineal o modulo de young.
La rigidez de un material queda caracterizada por la relacion entre el esfuerzo σx y
deformacion εx.
CURVA ESFUERZO-DEFORMACION MATERIALES DUCTILES
Todo cuerpo al soportar una fuerza aplicada trate de deformarse en el sentido de aplicación
de la fuerza. En el caso del ensayo de traccion, la fuerza se aplica en direccion del eje de ella y por
eso se denomina axial. Aunque el esfuerzo y la deformacion ocurren simultaneamente en el ensayo,
MATRIZ DE RIGIDEZ PARA UN RESORTE 8
los dos conceptos son completamente distintos. La curva usual esfuerzo-deformacion, expresa
tanto el esfuerzo como la deformacion en terminos de las dimensiones originales de la probeta, un
procedimiento muy util cuando se esta interesado en determinar los datos de resistencia y
ductilidad para proposito de diseño en ingenieria.
El diseño de elementos estructurales implica determinar la resistencia y rigidez del material
estructural, estas propiedades se pueden relacionar si se evalua una barra sometida a una fuerza
axial para la cual se registra simultaneamente la fuerza aplicada y el aalrgamiento producido. Estos
valores permiten determinar el esfuerzo y la deformacion que al graficar originan el denominada
diagrama de esfuerzo y deformacion. Los diagramas son similares si se trata del mismo material y
de manera general permite agrupar los materiales dentro de dos categorias on propiedades afines
que se denominan materiales ductiles y materiales fragiles. Los diagramas de materiales ductiles
se caracterizan por ser capaces de resistir grandes deformaciones antes de la rotura, mientras que
los fragiles representan un alargamiento bajo cuando llegan al punto de rotura.
MATRIZ DE RIGIDEZ PARA UN RESORTE 9
En un diagrama se observa un tramo recta inicial hasta un punto denominado limite de
proporcional. este limite tiene gran importancia para la teroria de los solidos estaticos, ya que se
basa en el citado limite. Este limite es el superior para un esfuerzo admisible. Los puntos
importantes del diagrama de esfuerzo deformacion son:
Limite de proporcionalidad: hasta este punto la relacion entre el esfuerzo y la deformacion
es lineal.
Limite de elasticidad: mas alla de este limite el material no recupera su forma original al
ser descargado, quedando con una deformacion permanente.
Punto de cedencia: aparece en el diagrama un considerable alargamiento o cedencia sin el
correspondiente aumento de carga. Este fenomeno no se observa en los materiales fragiles.
Esfuerzo ultimo: maxima ordenada del diagrama esfuerzo-deformacion.
Punto de ruptura: cuando el material falla. Dado que el limite de proporcionalidad,
elasticidad y punto de cedencia están tan cerca se considera para la mayoria de los casos como el
mismo punto.
De manera que el material al llegar a la cedencia deja de tener un comportamiento elastico
y la relacion lineal entre el esfuerzo y la deformacion deja de existir.
ESFUERZO.
Se denomina esfuerzo a la fuerza por unidad de área, la cual se denota con la letra griega
sigma y es un parametro que permite comparar la resistencia de dos materiales ya que establece
una base comun de referencia.
Donde:
P: Fuerza axial.
A: Area de la seccion transversal.
MATRIZ DE RIGIDEZ PARA UN RESORTE 10
Se denomina tension al esfuerzo interno a que esta sometido un cuerpo por la aplicación
de dos fuerzas que actuan en sentido opuesto y tienden a estirarlo.
MATRIZ DE RIGIDEZ PARA UN RESORTE 11
Referencias
https://w3.ual.es/~mnavarro/Tema%206%20%20Elasticidad.pdf
esfuerzo y deformacion. Instituto universitario politecnico. Yohandry rodriguez. 2018.
https://es.slideshare.net/yohandryRodriguez/esfuerzo-y-deformacion-33984306
Fundamentals of finite element analysis. David V. hutton. Temas 2.1
Analisis matricial de estructuras. E. Blanco, M.cervera y B. Suárez.