UNIVERSIDAD AUTÓNOMA DE BAJA CALIFORNIA
Facultad de Ingeniería Arquitectura y Diseño
Servicio social profesional
Supervisor de servicio social:
Dr. Manuel Antonio Barraza Guerrero
Presentado por:
Santiago Caballero Antonio Miguel (351165)
Fecha:
16 de diciembre del 2022
Método de rigidez en Octave
UNIVERSIDAD AUTÓNOMA DE BAJA CALIFORNIA
Facultad de Ingeniería Arquitectura y Diseño
Contenido
1 Introducción 3
2 Fundamento teórico 3
2.1 Método de rigidez 3
2.2 Ecuaciones de elástica 7
3 Solución en software Octave 11
3.1 Método de rigidez 12
3.2 Ecuaciones de momentos y cortante 15
3.3 Representación grafica 16
4 Modo de uso del programa 17
5 Comprobación de resultados 19
6 Bibliografía 21
UNIVERSIDAD AUTÓNOMA DE BAJA CALIFORNIA
Facultad de Ingeniería Arquitectura y Diseño
1 Introducción
La resolución de vigas hiperestáticas suele ser aplicados, dado que lleva
procedimientos sumamente largos, que se complican de acuerdo a las
características de la viga, propiedades mecánicas de inercia y módulo de
elasticidad, longitudes, condiciones de apoyos, fuerzas gravitacionales actuantes
entre otros, además de esto generalmente en la materia de Análisis estructural el
tiempo es insuficiente, lo que ocasiones que no se estudien de la manera correcta
todos los diferentes métodos o simplemente se omitan algunos, lo que ocasiona
que perdamos las habilidades y herramientas para la solución de este tipo de
vigas.
Lo que busca este proyecto de servicio social, es dar a conocer acerca del método
de rigidez para dar solución a vigas hiperestáticas, esto mediante una recopilación
de teoría y un programa que nos ayuda a resolver por este método.
Por lo tanto, el siguiente documento es sustento del programa llamados Método de
rigidez de vigas, el cual es un programa elaborado con el software de licencia libre
Octave.
Dicho programa resuelve vigas isostáticas e hiperestáticas, esto lo hace
encontrando las reacciones que cada nodo, llámense fuerza vertical y momentos,
además presenta la forma general que sigue los diagramas de momentos de
fuerza cortante y momentos.
2 Fundamento teórico
2.1 Método de rigidez
Este método se fundamente en los principios de: continuidad, la ley de Hooke y el
de equilibrio. El principio de continuidad se define a través del establecimiento de
la compatibilidad delas deformaciones de una sección de una barra, con las
secciones circundantes a ella. A través de la ley de Hooke se establecen las
relaciones esfuerzo-deformación de los materiales, mientras estos se comportan
en el rango elástico, la relación será lineal. Finalmente integrar los efectos
producidos por las deformaciones de una junta con los efectos considerados en la
UNIVERSIDAD AUTÓNOMA DE BAJA CALIFORNIA
Facultad de Ingeniería Arquitectura y Diseño
ley de Hooke, es una operación que se debe llevar a cabo cumpliendo con los
requerimientos de equilibrio estático.
Los principios mencionados anteriormente conducen a la relación F=kδ. En donde
“F”, es la fuerza en la dirección del desplazamiento “δ’, que se produce debido a la
rigidez “k”. Al considera un sistema de fuerzas asociados a un campo de
desplazamiento, la expresión en su representación matricial tiene la siguiente
forma: {F}={k}{δ}. La expresión anterior es la forma básica del método de rigidez.
Conociendo esto podemos describir el procedimiento general para definir la matriz
de rigidez de un elemento barra, de sección primatica de sección constante,
sometida a esfuerzo de flexión (consideran dos grados de libertad, rotación
positiva en sentido horario y desplazamiento vertical positivo hacia arriba).
1. Identificación de nodos para subdividir el elemento.
2. Definición de coordenadas globales y del elemento.
3. Indeterminación cinemática de los elementos (desplazamientos verticales y
rotación).
UNIVERSIDAD AUTÓNOMA DE BAJA CALIFORNIA
Facultad de Ingeniería Arquitectura y Diseño
4. Matriz de rigidez
a. Desplazamiento en Y
Cuando se impone un desplazamiento positivo DY1, mientras se evitan otro posible
desplazamiento, se crean fuerzas cortante y momentos de flexión resultantes.
Donde
UNIVERSIDAD AUTÓNOMA DE BAJA CALIFORNIA
Facultad de Ingeniería Arquitectura y Diseño
12𝐸𝐼
𝑄𝑌 = 3
𝐿
6𝐸𝐼
𝑄𝑍 = 2
𝐿
b. Rotación en Z.
Cuando se impone un desplazamiento positivo DZ1, mientras se evitan otro posible
desplazamiento, se crean fuerzas cortante y momentos de flexión resultantes.
Donde:
6𝐸𝐼
𝑄𝑦 = 2
𝐿
4𝐸𝐼
𝑄𝑍 = 𝐿
2𝐸𝐼
𝑄𝑍(𝐵) = 𝐿
UNIVERSIDAD AUTÓNOMA DE BAJA CALIFORNIA
Facultad de Ingeniería Arquitectura y Diseño
Por superposición:
[𝑄𝑌1 𝑄𝑍1 𝑄𝑌2 𝑄𝑍2 ] = ⎡⎢𝑌1
12𝐸𝐼 6𝐸𝐼 12𝐸𝐼 6𝐸𝐼 6𝐸𝐼 4𝐸𝐼 6𝐸𝐼 2𝐸𝐼 12𝐸𝐼 6𝐸𝐼
3 𝑍1 2 𝑌2 − 3 𝑍2 2 2 𝐿
− 2 𝐿
− 3 2 −
⎣ 𝐿 𝐿 𝐿 𝐿 𝐿 𝐿 𝐿 𝐿
Los valores de Q y D, son valores de cargas y desplazamientos, conocidos tanto
como desconocidos.
𝑄 = 𝐾𝐷
[𝑄𝑎 𝑄𝑏 ] = [𝐾11 𝐾12 𝐾21 𝐾22 ][𝐷𝑏 𝐷𝑎 ]
Expandiendo:
𝑄𝑎 = 𝐾11𝐷𝑏 + 𝐾12𝐷𝑎
𝑄𝑢 = 𝐾21𝐷𝑏 + 𝐾22𝐷𝑎
Los desplazamientos desconocidos Db, se determinan a partir de la primera de
estas ecuaciones. Si se usan estos valores, pueden calcular las reacciones de
apoyo Qu para la segunda ecuación.
2.2 Ecuaciones de elástica
En el método de elementos finitos es conveniente colocar nudos en los extremos
de la viga, en los apoyos, donde existan cargas puntuales, y en todo sector donde
aparezca un cambio brusco en el tipo de carga, resuelto el ejercicio se tendrán el
desplazamiento vertical y el giro para cada nudo.
El método de doble integración de la ecuación de la elastica se utiliza cuando
terminamos de analizar los desplazamientos en los nudos de una viga resuelta por
el método de elementos finitos o de rigidez.
UNIVERSIDAD AUTÓNOMA DE BAJA CALIFORNIA
Facultad de Ingeniería Arquitectura y Diseño
Mediante este método se podrá encontrar el diagrama de momento flector o
cortante de cada barra entre dos nudos. Se deberá repetir el procedimiento para
cada barra entre dos nudos de la viga.
El tramo de viga debe tener los siguientes datos:
Se debe tener en cuenta que la carga distribuida mostrada representa el caso más
general de cargas distribuidas. A partir de esta configuración se pueden generar
cargas distribuidas de distinto tipo.
Dicho esto, podemos describir el método para un tramo con datos indicados
conocidos:
Lo primero que debemos hacer es conocer la ecuación de la carga distribuida a
partir de las cargas q1 y q2 en los extremos. Para eso nos ayudamos de la teoría
de geometría analítica donde encontramos los coeficientes A y B de la ecuación
de la recta genérica mostrada:
En base a las coordenadas mostradas se forma un sistema de ecuaciones a partir
del cual podemos encontrar los coeficientes A y B de la ecuación dela recta.
UNIVERSIDAD AUTÓNOMA DE BAJA CALIFORNIA
Facultad de Ingeniería Arquitectura y Diseño
Despejando B de la primera ecuación y reemplazando la en la segunda podemos
obtener los valores de A y B cómo se muestra:
En la ecuación mostrada se reemplazaron además los valores de A y B dentro de
la ecuación de la recta.
Partimos de la ecuación de la elástica para encontrar los diagramas de cortante y
de momento flector.
Se sabe de la teoría que al integrar esta ecuación 4 veces podemos encontrar
todos los datos del cortante, momento flector, pendiente y deflexión de La viga:
UNIVERSIDAD AUTÓNOMA DE BAJA CALIFORNIA
Facultad de Ingeniería Arquitectura y Diseño
De las anteriores ecuaciones la tercera derivada de la deflexión representa el
cortante de la viga y la segunda derivada representa el momento flector.
Acto seguido se reemplazan las condiciones de contorno en las ecuaciones
mostradas. En este caso las condiciones de contorno son los desplazamientos y
giros en los extremos de la barra estudiada. En el extremo izquierdo para la
coordenada x1 se tienen las condiciones de contorno de desplazamiento u1 y de
giro u2. Una cosa similar sucede en el extremo derecho como se muestra a
continuación:
Ahora se reemplaza las condiciones de contorno dentro de las ecuaciones
integradas de la elástica (tercera y cuarta ecuación) para obtener el sistema de 4
ecuaciones con 4 incógnitas. Las incógnitas en este caso son las 4 constantes de
integración C1, C2, C3, C4.
Hacemos los términos independientes llevándolos todos al extremo izquierdo de
las ecuaciones:
UNIVERSIDAD AUTÓNOMA DE BAJA CALIFORNIA
Facultad de Ingeniería Arquitectura y Diseño
Ahora finalmente convertimos el sistema de ecuaciones a un arreglo matricial que
puede resolverse por cualquier método numérico:
De las soluciones de este sistema de ecuaciones C1, C2, C3, C4 se pueden
conocer y a las ecuaciones de cortante y momento flector para el tramo estudiado:
Se debe repetir el mismo procedimiento para todos los tramos que existan dentro
de la viga.
Todo quedará mucho más claro con un ejemplo, que presento en el siguiente post
y además en el vídeo mostrado a continuación:
3 Solución en software Octave
Para la resolución en el software, nos basamos en el procedimiento antes descrito
en donde podemos describir tres procedimientos principales, calcular reacciones y
desplazamientos con el método de rigidez, calcular las ecuaciones de momento y
UNIVERSIDAD AUTÓNOMA DE BAJA CALIFORNIA
Facultad de Ingeniería Arquitectura y Diseño
cortante con ayuda de las ecuaciones de la elástica, y la representación gráfica de
los resultados obtenidos. Para ello se hace una solicitud de valores o datos
conocidos al usuario, de número de elementos, longitudes, inercias, módulo de
elasticidad, cargas distribuidas, cargas puntuales, se definen las características de
los apoyos, es decir, si es un empotramiento, rodillo o articulación, además
definimos el nodo al que pertenecen los elementos, la cuales van a ayudar a
determinar la matriz de rigidez de los elementos.
3.1 Método de rigidez
Para dicha resolución con los valores introducidos se calcula la matriz de rigidez
de un elemento (tramos de viga) cuya dimensión será de 4 X 4, el cual se
almacenará de una matriz de matrices y tendrá la capacidad de almacenar las
demás matrices de elementos subsecuentes.
De la misma manera con los datos obtenidos, se creará una matriz de grados de
libertad, el cual será cuadrada y sus dimisiones será tal que corresponda al
número total de grados de libertad de toda la viga, cada nodo por definición puede
tener dos llámese desplazamiento vertical y rotación.
Posterior a eso las matrices de rigidez de cada elemento se almacenará en una
matriz global de acorde al grado de libertad al que correspondan, en caso de
elementos que compartan nodo, los valores de la matriz de rigidez se sumaran.
UNIVERSIDAD AUTÓNOMA DE BAJA CALIFORNIA
Facultad de Ingeniería Arquitectura y Diseño
Ahora se definirá nuestra matriz incógnita, es decir la matriz de desplazamientos,
para ello con ayuda de las restricciones que hace el usuario se organizan tan que
la matriz sea de una columna y que el número de filas se los grados de libertan de
la viga.
Posterior a eso se calcula un vector de fuerzas, en la cuales se encuentra la
reacción, que son incógnitas, las cargas puntuales, y las cargas distribuidas que a
su vez generan carcas puntuales y momentos. Para ello se realizó un algoritmo
que almacenara en un vector de una columna y que el número de filas
corresponda al grado de libertad, dichas fuerzas. Las puntales se almacenan
según la posición asignada por el usuario, la distribuidas entran en un bucle, y
según a la forma que tenga la carga se calculan las reacciones generadas con el
teorema de momentos de empotramiento perfecto, y se suman a la posición en el
grado de libertad correspondientes, la reacciones en los apoyos (incógnitas) se
manejan en una matriz de incógnitas por practicidad, las cuales estarán sumando
a la matriz calculada anteriormente.
UNIVERSIDAD AUTÓNOMA DE BAJA CALIFORNIA
Facultad de Ingeniería Arquitectura y Diseño
En este punto obtendremos la siguiente operación de matrices, en su forma
general cuando el desplazamiento sea conocido es decir igual a cero, la reacción
es incógnita, caso contrario cuando el desplazamiento es desconocido la reacción
es cero:
[𝐹] + [𝑅] = [𝐾][𝐷]
La matriz [𝐷], es de restricciones, entonces cuando el valor de la restricción sea 0,
se eliminara la fila y columna con dicha posición, el cual se realizara con el
siguiente bucle:
UNIVERSIDAD AUTÓNOMA DE BAJA CALIFORNIA
Facultad de Ingeniería Arquitectura y Diseño
Por lo tanto, al eliminar las columnas con restricción de ceros, entonces
eliminamos las reacciones incógnitas:
[𝐹] = [𝐾][𝐷]
Los desplazamientos reales serán tales que:
−1
[𝐾] [𝐹] = [𝐷]
Por lo tanto, para obtener las reacciones volvemos a la forma original de la matriz
y realizamos la siguiente operación con matrices:
[𝑅] = [𝐾][𝐷] − [𝐹]
3.2 Ecuaciones de momentos y cortante
Siguiendo los descrito en el fundamento teórico hacemos lo siguiente con el
lenguaje de programación
UNIVERSIDAD AUTÓNOMA DE BAJA CALIFORNIA
Facultad de Ingeniería Arquitectura y Diseño
Conociendo a la forma general de las ecuaciones podemos calcular los valores y
colocarlos es su posición, esto con ayuda de los datos administrados por el
usuario
3.3 Representación grafica
Para la representación graficas se usan las ecuaciones calculadas anteriormente,
para ellos se tabulan según el rango en la viga a la que pertenecen.
UNIVERSIDAD AUTÓNOMA DE BAJA CALIFORNIA
Facultad de Ingeniería Arquitectura y Diseño
4 Modo de uso del programa
El uso de programa es intuitivo y el programa se dan una pequeña instrucciones
de la forma de utilización, para ello se hace una solicitud de valores o datos
conocidos al usuario, de número de elementos, longitudes, inercias, módulo de
elasticidad, cargas distribuidas, cargas puntuales, se definen las características de
los apoyos, es decir, si es un empotramiento, rodillo o articulación, además
definimos el nodo al que pertenecen los elementos, la cuales van a ayudar a
determinar la matriz de rigidez de los elementos.
UNIVERSIDAD AUTÓNOMA DE BAJA CALIFORNIA
Facultad de Ingeniería Arquitectura y Diseño
Los valores del número de elementos son de acuerdo a la forma de dividir la viga
analizada, para ello es conveniente colocar nudos en los extremos de la viga, en
los apoyos, donde existan cargas puntuales, y en todo sector donde aparezca un
cambio brusco en el tipo de carga, resuelto el ejercicio se tendrán el
desplazamiento vertical y el giro para cada nudo.
El valor de la longitud, módulo de elasticidad, e inercia, se asignan de acuerdo
a los datos de problemas, cada dato será separado por un espacio.
En la matriz de cargas distribuidas, el primer valor es el número de elementos y
los siguientes son las cargas al inicio y al final del elemento.
Para definir cargas puntuales, se coloca de la misma manera que en el ejemplo,
considerando que las cargas puntuales caen en el grado de liberta que
corresponde al desplazamiento vertical.
UNIVERSIDAD AUTÓNOMA DE BAJA CALIFORNIA
Facultad de Ingeniería Arquitectura y Diseño
Las restricciones se asignan según el tipo de apoyo, para apoyos empotrados 0
0, para apoyos articulados y rodillos 0 1, para lugares en donde no exista apoyos 1
1.
Por ultimo definir a que numero de nodo está conectado cada tramo, para ello en
la primera columna se coloca el número de nodos, y posterior a ello de que
numero de nodo y a que numero de nodos van.
Se corre el programa y podemos ver los resultados de la reacción de acuerdo al
grado de libertad, ya sea reacción o momento, además podemos ver la forma
general de los diagramas de momento y de corte.
Para fines didácticos se puede acceder a los valores de las distintas variables
utilizadas en el programa, para ello basta con escribir en la venta de comandos la
variable deseada, respetando la escritura asignada a esta.
5 Comprobación de resultados
Para fines de corroborar la eficacia del programa se realza una demostración de
su funcionamiento de acuerdo a un ejercicio de libro consultado. Para dicha
comprobación nos apoyaremos del ejemplo 15.4 de libro Análisis estructural, 8va
edición, el cual dice los siguiente:
UNIVERSIDAD AUTÓNOMA DE BAJA CALIFORNIA
Facultad de Ingeniería Arquitectura y Diseño
Para lo cual hacemos las sustituciones en nuestro programa obteniendo los
siguientes resultados.
Las cuales corresponden a las reacciones y momentos en cada uno de los nodos,
empezado de izquierda a derecha, de las mima manera se muestra una
representación grafica de la viga y del comportamiento de los diagramas de
cortante y momentos.
UNIVERSIDAD AUTÓNOMA DE BAJA CALIFORNIA
Facultad de Ingeniería Arquitectura y Diseño
UNIVERSIDAD AUTÓNOMA DE BAJA CALIFORNIA
Facultad de Ingeniería Arquitectura y Diseño
6 Bibliografía
Ojeda J. Apuntes de análisis estructural. Universidad autónoma de baja california,
faculta de ingeniería arquitectura y diseño.
Marcelo Pardo. (2022, 18 de septiembre). Escribe tu propio programa de Vigas.
Solicitaciones de Momento y Cortante. [Archivo de video]:
https://www.youtube.com/watch?v=Eih5VA8Fj2s
Marcelo Pardo. (2022, 18 de septiembre). Teoria de pórticos por método de
rigidez. [Archivo de video]:
https://www.youtube.com/watch?v=Su51D3SuHZg
Marcelo Pardo. (2022, 18 de septiembre). Solicitaciones de Momentos y cortantes
en Vigas a partir de deflexiones - Método de Rigidez. [Archivo de video]:
https://www.youtube.com/watch?v=Ky1cYSU76uw&t=188s
Marcelo Pardo. (2022, 18 de septiembre). Pórticos: Ensamblaje de Matriz de
rigidez en Matlab u Octave por método de rigidez. [Archivo de video]:
https://www.youtube.com/watch?v=7l3D9vObquI&t=1248s
Marcelo Pardo. (2022, 18 de septiembre). Un Matlab gratuito - Diagramas de
Momento flector con Octave. [Archivo de video]:
https://www.youtube.com/watch?v=JUMr8XER9h0
Pardo M. (2022, 15 de diciembre). Apuntes de ingeniería estructura. Marcelo
pardo ingeniería. https://marcelopardo.com/
Mates NoéC. (2022, 15 de diciembre). Gráficos en Octave. [Archivo de video]:
https://www.youtube.com/watch?v=PFL8jcpdk-k
Andrés Weitzel [.archivo]. (2022, 15 de diciembre). 11 SubGráficos | Función
subplot() | GNU Octave / MATLAB. [Archivo de video]:
https://www.youtube.com/watch?v=Fsqq-CXT4rs&t=1s
UNIVERSIDAD AUTÓNOMA DE BAJA CALIFORNIA
Facultad de Ingeniería Arquitectura y Diseño
