Análisis Estructural de una Armadura por medio de Matlab y Robot Structural

Analysis

Morales Castillo, Rodrigo Andrés

Aguilar Alzamora, Rosmery
Abanto Torres, Edwar Franklin
Santé Pisconti, Bruno Alonso

Departamento: Lima
Institución: Universidad Peruana de Ciencias Aplicadas
Fecha: martes 9 de junio de 2019
Análisis Estructural de una Armadura por medio de Matlab y
Robot Structural Analysis

Morales Castillo Rodrigo Andrés

Aguilar Alzamora Rosmery
Abanto Torres, Edwar Franklin
Santé Pisconti, Bruno Alonso
Departamento: Lima
Institución: Universidad Peruana de Ciencias Aplicadas
Fecha: martes 9 de junio de 2020

Resumen

Para el presente análisis estructural se trabajó tanto con el Software Matlab, así como
Dynamo, un programa dentro de Revit, y Robot Structural Analysis Professional, programas
comerciales de cálculo y diseño estructural, aplicados en el campo de la ingeniería en general y en
este caso en particular en la rama estructural. Matlab es un sistema de cómputo numérico que
ofrece un entorno de desarrollo integrado con un lenguaje de programación propio; es más,
comprende el análisis numérico, el cálculo matricial, el procesamiento de señales y los gráficos; por
otro lado Robot Structural Analysis Professional es un software especializado en el cálculo,
simulación y diseño de estructuras de ingeniería civil y arquitectura que permite realizar un análisis
exhaustivo del comportamiento estructural de una edificación o de cualquiera de sus partes,
utilizando el método de los elementos finitos (FEM) y varios tipos de algoritmos no lineales y
Dynamo es un programa de scripts visuales que permite crear algoritmos personalizados a fin de
procesar datos y generar geometrías. En otras palabras, Dynamo nos permite programar mediante la
utilización de elementos visuales (nodos), en vez de usar las típicas líneas de código

Detail such as name of institution ii Detail such as date

Contents

Resumen ii

Introducción al Elemento a Estudiar 1

Marco Teórico 2
Matriz de Rígidez 3
Metodología 5
Programación en Matlab 3
Programación Generacional en Dynamo
Analisís estructural en Robot

Resultados
Matriz de rigidez en Matlab
Líneas de Influencia por Sobrecarga de Uso
Comparativa del Caso 1
Comparativa del Caso 2
Comparativa del Caso 3

Anexo 1: Código en Matlab

Anexo 2: Guide de Matlab

Anexo 3: Programación generacional

Referencias Bibliográficas

Universidad Peruana de Ciencias Aplicadas iii martes 9 de junio de 2020

Introducción al Elemento a estudiar

Normas a Utilizar
En este trabajo será necesario el uso de la norma E.020 para establecer la sobrecarga de uso
en la losa y las combinaciones de cargas.

Descripción de la Estructura
Se pide analizar el eje E-E del pórtico de un restaurante. La estructura soportará una carga
sísmica de 75 toneladas, además de acabados y un muro de tabiquería ubicado encima del pórtico.
Además, se tendrá consideración de la sobrecarga por uso, en este caso restaurante, según la norma
E.020. Las medidas del pórtico se expondrán en la tabla 1.

Table 1. Tabla de datos para los grupos, el presente trabajo es realizado por el grupo 8. Fuente: Ing.
Elvis Jony Mamani Yana.

Imagen 1. Elevación del pórtico. Fuente: Ing. Elvis Jony Mamani Yana.

Universidad Peruana de Ciencias Aplicadas 1 28 de junio de

2020
Imagen 2. Vista en planta de la losa. Fuente: Ing. Elvis Jony Mamani Yana.

Imagen 3. Extracto de la tabla 1 de la norma E.020, la sobrecarga de uso a usar para el pórtico es 400 kg/m2. Fuente: Ing. Elvis Jony
Mamani Yana.

Resumen de datos del grupo 8

 L1(m): 5.50 m
 L3(m): 4.50 m
 L4(m): 5.00 m
 L5(m): 4.50 m
 L6(m): 6.00 m
 h1(m): 3.50 m
 Peso de tabiquería: 6.75 kg/cm
 Peso de acabado: 120 kg/m2
 Peso de losa: 350 kg/m2
 Sobrecarga de uso según la Norma E.020: 400 kg/m2
 Sección de columna: 30x30 cm
 Sección de viga: 30x60 cm

Objetivo
El objetivo del trabajo planteado es familiarizar al alumno en la resolución de problemas
estructurales de pórticos implementado Softwares (Matlab y Robot).

Marco Teórico

Matriz de Rigidez

Metodología

Modelo Matemático
Para realizar el modelo matemático se toma en cuenta las longitudes entre ejes de los
elementos. En la imagen 4 se muestra el modelo matemático a trabajar en los programas.

Imagen 3. Modelo matemático del eje E-E del pórtico. Fuente: Propia, AutoCAD.
 Imagen 4. Sentido de los elementos para las matrices. Fuente: Propia, AutoCAD.

Metrado de cargas.
Carga Muerta:
5.50
∗350+ 120
2
+6.75=19.675 kg/cm …(ec . 1)
100

Carga Viva:
400
∗5.50
100
=11 kg /cm… (ec .2)
2

Carga Sismo:

75000 kg …(ec .3)

Programación en Matlab
A continuación, se pasará a explicar a detalle el proceso de programación que hemos tenido
en cuenta para la resolución del problema planteado. Se adjuntan imágenes para que el lector tenga
una mejor comprensión y sea más didáctico.

a. Empezamos creando un guide capaz de mostrar correctamente los resultados obtenidos (Ver
anexo 2).

b. Dentro del script, procedemos a la recolección de datos del guide.

 Imagen 5. Código para la recolección de datos y cálculo de inercias. Fuente: Propia, Matlab.

c. Procedemos a colocar las coordenadas y los grados de libertad en una matriz.

Imagen 6. Código para la introducción de coordenadas de los extremos de los elementos. Fuente: Propia, Matlab.

d. Realizamos la topología y determinamos la longitud de los elementos.

Imagen 7. Código para establecer la topología y las longitudes. Fuente: Propia, Matlab.

e. Realizamos el arreglo para determinar la matriz de rigidez y su ensamble.

Imagen 8. Determinación de la matriz de desplazamientos. Fuente: Propia, Matlab.

f. Procedemos a realizar nuestra matriz de fuerzas. Para ello, metramos las cargas según la
sección a analizar, introduciendo los valores dentro de la matriz de fuerzas según el GDL
en la que actúan.

g. Procedemos a calcular la matriz de desplazamientos.

Imagen 11. Determinación de la matriz de desplazamientos. Fuente: Propia, Matlab.

h. Procedemos a calcular las fuerzas axiales y ver si los elementos se encuentran en tracción o
compresión.

i. Realizamos el código para graficar la armadura deformada con la armadura original.

Para la obtención de todos los resultados provenientes del código programado en Matlab se
procede a ejecutar el programa, seleccionando el tipo de sección a analizar y la ubicación de la
carga de uso, finalmente apretamos el botón calcular y obtendremos los resultados, así como una
gráfica comparativa entre la armadura inicial y la armadura deformada.
Imagen 12. Resultados obtenidos por el código cuando se analiza una sección interna de la armadura y no se considera
la carga de uso. Fuente: Propia, Matlab.

Análisis Estructural en Robot Structural Analysis

Para realizar el análisis en Robot, trazamos vigas y columnas según el modelo matemático
que se muestra en la imagen 3. La imagen 13 muestra el resultado del trazado. Lo siguiente que
vamos a hacer es crear las secciones de vigas y columnas y las asignamos a los elementos.

Imagen 14. Vista frontal del dibujo del pórtico sin asignar secciones. Fuente: Propia, Robot Structural Analysis.

Imagen 14 y 15. Creación de secciones de columna y viga. Fuente: Propia, Robot Structural Analysis.
Imagen 16. Elevación del pórtico con sus secciones asignadas. Fuente: Propia, Robot Structural Analysis.

Lo que sigue es crear los 3 casos de carga: carga muerta, viva y sísmica. Luego colocamos
las cargas según el metrado (ver ecuaciones 1, 2 y 3) y su ubicación.

Imagen 16. Carga Muerta aplicada al pórtico. Fuente: Propia, Robot Structural Analysis.

Imagen 17. Carga Viva aplicada al pórtico. Fuente: Propia, Robot Structural Analysis.

Imagen 17. Carga Viva aplicada al pórtico. Fuente: Propia, Robot Structural Analysis.

Una vez colocadas las cargas procedemos a realizar la combinación de cargas y

finalmente determinar la envolvente.
Imagen 17. Cuadro de combinaciones de cargas. Fuente: Propia, Robot Structural Analysis.

Resultados

Matriz de Rigidez en Matlab

Cuando oprimimos el botón “Calcular” dentro del programa realizado en Matlab, se
ejecutará el código para el análisis. Una tabla nos mostrará la matriz de rigidez de la sección con
dimensión de 24x24, ya que fueron 24 grados que se ingresaron. Los 3 últimos grados están
restringidos. La matriz resultante es simétrica y presenta una diagonal positiva, características
propias de la matriz de rigidez. La imagen 18 muestra la matriz del elemento.

Diagramas en ROBOT
Para obtener los resultados se aprieta el botón “calculate” y se selecciona los diagramas
que se desean mostrar. Las siguientes imágenes muestran los diagramas por casos.
Caso 1: Carga Muerta

Imagen 19. Diagrama de fuerzas cortantes en kg resultado de las cargas muertas. Fuente: Propia, Robot Structural Analysis.

Imagen 19. Diagrama de momento flector en tn/n resultado de las cargas vivas. Fuente: Propia, Robot Structural Analysis.
 Caso 2: Carga Viva

Imagen 19. Diagrama de fuerzas cortantes en kg resultado de las cargas vivas. Fuente: Propia, Robot Structural Analysis.

Imagen 19. Diagrama de momento flector en tn/n resultado de las cargas muertas. Fuente: Propia, Robot Structural Analysis.

Caso 3: Carga Sísmica

Imagen 19. Diagrama de fuerzas cortantes en kg resultado de las cargas sísmica. Fuente: Propia, Robot Structural Analysis.
Imagen 19. Diagrama de momento flector en tn/n resultado de las cargas sísmica. Fuente: Propia, Robot Structural Analysis.

De la misma forma calculamos los diagramas para las combinaciones.

Combinación 1: 1.4CM+1.7CV

Imagen 19. Diagrama de fuerzas cortantes en kg resultado de combinación 1. Fuente: Propia, Robot Structural Analysis.

Imagen 19. Diagrama de momento flector en tn/n resultado de combinación 1. Fuente: Propia, Robot Structural Analysis.

Combinación 2: 1.25(CM+CV) +CS

Imagen 19. Diagrama de fuerzas cortantes en kg resultado de combinación 1. Fuente: Propia, Robot Structural Analysis.

Imagen 19. Diagrama de momento flector en tn/n resultado de combinación 1. Fuente: Propia, Robot Structural Analysis.

Combinación 3: 1.25(CM+CV)-CS

Imagen 19. Diagrama de fuerzas cortantes en kg resultado de combinación 1. Fuente: Propia, Robot Structural Analysis.
Imagen 19. Diagrama de momento flector en tn/n resultado de combinación 1. Fuente: Propia, Robot Structural Analysis.

Combinación 4: 0.9CM+CS

Imagen 19. Diagrama de fuerzas cortantes en kg resultado de combinación 1. Fuente: Propia, Robot Structural Analysis.

Imagen 19. Diagrama de momento flector en tn/n resultado de combinación 1. Fuente: Propia, Robot Structural Analysis.

Combinación 4: 0.9CM-CS
Imagen 19. Diagrama de fuerzas cortantes en kg resultado de combinación 1. Fuente: Propia, Robot Structural Analysis.

Imagen 19. Diagrama de momento flector en tn/n resultado de combinación 1. Fuente: Propia, Robot Structural Analysis.

ENVOLVENTE

Imagen 19. Diagrama de fuerzas cortantes en kg resultado de combinación 1. Fuente: Propia, Robot Structural Analysis.
Imagen 19. Diagrama de momento flector en tn/m resultado de combinación 1. Fuente: Propia, Robot Structural Analysis.

Anexo 1: Código en Matl

Anexo 2: Guide de Matlab

Imagen A2-1. Guide realizado para mostrar y ejecutar el código. El guide presenta un axes para la proyección de las
gráficas de las armaduras, 2 popupmenu para escoger la sección y la ubicación de la carga de uso, 4 tablas para mostrar
la matriz de rigidez, los desplazamientos, las fuerzas axiales y el estado de las barras al ser sometidas a estas fuerzas; y
3 botones: Calcular ejecuta el código, Nuevo limpia el guide, dejándolo listo para otro análisis; y Salir sirve para cerrar
el programa. Fuente: Propia, Matlab.
References

1. Análisis y Diseño de Estructuras con AUTODESK ROBOT STRUCTURAL ANALYSIS 2016 , 3 th ed. by
Carlos Villaroel, https://drive.google.com/file/d/1PTEC-
xZlGOITuwtTpado_EkpP2_hpNL5/view?fbclid=IwAR3d1Hijt0k5AHCqVTTqxHEvaVlw0ygs
pCh-5dR8FpYtzKyYjc0lDCrLGm4(accessed 01 June 2020).

3.
Marc E. Herniter, “Programing in Matlab,” (Thomson learning , 2001)

4.
Peregrina Estevez, “Matemáticas en Ingeniería con Matlab,” (Universidade de Santiago de
Compostela , 2000)