Análisis Pseudo tridimensional de una estructura porticada de tres pisos en la ciudad de Chanchamayo
Pseudo three-dimensional analysis of a three-story porticoed structure in the city of Chanchamayo
Análise pseudotridimensional de uma estrutura porticada de três andares na cidade de Chanchamayo
Correspondencia: mrupay@uniscjsa.edu.pe
Ciencias Técnica y Aplicadas
Artículo de Investigación
* Recibido: 23 de octubre de 2023 *Aceptado: 22 de noviembre de 2023 * Publicado: 23 de diciembre de 2023
- Universidad Nacional Intercultural Juan Santos Atahualpa, Chanchamayo, Perú.
- Universidad Nacional Intercultural Juan Santos Atahualpa, Chanchamayo, Perú.
- Universidad Nacional Intercultural Juan Santos Atahualpa, Chanchamayo, Perú.
- Universidad Nacional Intercultural Juan Santos Atahualpa, Chanchamayo, Perú.
Resumen
El presente artículo científico consiste en calcular las deformaciones y desplazamientos por cada entrepiso, hallar las fuerzas en cada pórtico, encontrar las cortantes en cada uno de los pórticos y comprobar los resultados mediante el software estructural Etabs. La metodología aplicada es de carácter cuantitativo, el proceso de solución y los resultados fueron descritos por fórmulas matemáticas, además se agregó los apuntes de clases de la asignatura de Análisis Estructural II brindada en la Universidad Nacional Intercultural de la Selva Central Juan Santos Atahualpa, donde se pudo extraer pasos para lograr desarrollar el ejercicio planteado por los integrantes que desarrollan este artículo, se considera también los criterios y enseñanzas brindadas por el ingeniero a cargo de la asignatura. Para concluir, se tiene planeado plasmar un esquema del desarrollo del ejercicio, donde se va visualizar el método utilizado para el desarrollo del artículo de investigación. Los resultados se van a plasmar a detalle para que pueda ser interpretados por los lectores de este artículo científico. En conclusión, se observará los resultados que se obtuvo del ejercicio aplicativo, hechos manualmente y comparados con el software Etabs para tener una comparación sobre lo resuelto manualmente.
Palabras Clave: Estructura; Etabs; Pórticos.
Abstract
This scientific article consists of calculating the deformations and displacements for each mezzanine, finding the forces in each frame, finding the shears in each of the frames and checking the results using the Etabs structural software. The applied methodology is quantitative in nature, the solution process and the results were described by mathematical formulas, in addition the class notes of the Structural Analysis II subject offered at the National Intercultural University of the Central Jungle Juan Santos Atahualpa were added, where It was possible to extract steps to develop the exercise proposed by the members who develop this article, the criteria and teachings provided by the engineer in charge of the subject are also considered. To conclude, it is planned to outline a diagram of the development of the exercise, where the method used to develop the research article will be visualized. The results will be expressed in detail so that they can be interpreted by the readers of this scientific article. In conclusion, the results obtained from the application exercise will be observed, done manually and compared with the Etabs software to have a comparison of what was solved manually.
Keywords: Structure; Etabs; Porticos.
Resumo
Este artigo científico consiste em calcular as deformações e deslocamentos de cada mezanino, encontrar as forças em cada pórtico, encontrar os cisalhamentos em cada um dos pórticos e verificar os resultados utilizando o software estrutural Etabs. A metodologia aplicada é de natureza quantitativa, o processo de solução e os resultados foram descritos por fórmulas matemáticas, além disso foram acrescentadas as notas de aula da disciplina Análise Estrutural II oferecida na Universidade Nacional Intercultural da Selva Central Juan Santos Atahualpa, onde foi possível extrair passos para desenvolver o exercício proposto pelos membros que desenvolvem este artigo, também são considerados os critérios e ensinamentos fornecidos pelo engenheiro responsável pela disciplina. Para finalizar, pretende-se traçar um diagrama do desenvolvimento do exercício, onde será visualizado o método utilizado para desenvolver o artigo de pesquisa. Os resultados serão expressos detalhadamente para que possam ser interpretados pelos leitores deste artigo científico. Concluindo, serão observados os resultados obtidos no exercício de aplicação, feito manualmente e comparados com o software Etabs para se ter uma comparação do que foi resolvido manualmente.
Palavras-chave: Estrutura; Etabs; Pórticos.
Introducción
El Perú es un país altamente sísmico debido a que está ubicado dentro del Cinturón de Fuego del Pacífico. En esta región ocurren constantemente sismos, es parte de la naturaleza, y nada ni nadie los puede evitar (Instituto Geofísico del Perú, 2022).
Es por eso que el Perú cuenta con una lista larga de sismos, el último gran sismo ocurrió en el 2007 y las zonas más afectadas fueron las provincias de Ica, Pisco y Chincha en el departamento de Ica y la provincia de Cañete en el departamento de Lima (Centro de Estudios y Prevención de Desastres, 2016). Las consecuencias de este sismo dejaron 595 muertos, 1 800 de heridos, 76 000 viviendas totalmente destruidas y cientos de miles de damnificados.
Nace la necesidad de evaluar las posibles deformaciones, desplazamientos y entre otros puntos, de una estructura antes de su construcción; debido a que en el sector de la ingeniería civil se centra también en la seguridad ciudadana y se busca minimizar los daños ocasionados por fenómenos naturales en lo más mínimo.
El objetivo de este estudio es calcular los puntos anteriormente mencionados, y esto con el método de cortes del análisis pseudo tridimensional de una edificación de 3 niveles.
Material y métodos
Se tomó como referencia una edificación de 3 niveles ubicado en la ciudad de San Ramón - Chanchamayo
Se evaluó los planos que presentaba la edificación donde obtuvimos los dimensionamientos de las columnas, altura de los entre pisos y una vista en planta y perfil de la edificación como se observa en la figura 1, luego en función al plano en planta de cada nivel se procedió a plantear tres coordenadas por nivel en su centro de masa: dos desplazamientos traslacionales y un giro.
Después con la ayuda del Excel se procedió armar nuestra matriz de compatibilidad de la edificación en función a los pórticos planteados y en seguida se armó la matriz de rigidez lateral de cada pórtico por el método de corte y se procedió a ensamblar la matriz de la estructura.
Por último se procedió a calcular el vector de desplazamiento con métodos matriciales y determinas las fuerzas cortantes para cada pórtico para ser evaluados con el software ETABS.
1) Sistemas de coordenadas
Según (Basualdo, 2016) los análisis sísmicos pueden realizarse con dos fuerzas horizontales y un momento de torsión en planta por cada piso, en un punto que generalmente es el centro de masa del piso.
2) Matriz de compatibilidad
Según (Marcos, 2023) para cada estructura plana es necesario definir su ubicación en la planta del edificio y el sentido asumido positivo para su sistema local y con ello se analiza para cada pórtico donde se considera los desplazamientos unitarios y el giro en función al sentido asumido.
3) Matriz de rigidez de cada pórtico
Según (Dovas Ingenieria, 2021) Para armar la matriz de rigidez de cada pórtico primero se evaluará los grados de libertad laterales (traslacionales) para cada piso, sin giros.
En función a ello se arma la matriz de rigidez para un pórtico de 3 niveles como podemos ver en la siguiente ecuación.
Donde:
:
Matriz de rididez de cada portico
:
Rigidez del entrepiso 1,2 y 3
4) Ensamblaje de la matriz de rigidez
Según (Salinas Basualdo, s.f.) “En un pórtico de varios pisos, la matriz de rigidez total es una operación repetitiva de ensambles de matrices de elementos, sean vigas, columnas, muros o arriostres”.
Donde:
:
Matriz de rigidez de la edificación
:
Transpuesta de la matriz de compatibilidad de cada pórtico
:
Matriz de rigidez de cada pórtico
:
Matriz de compatibilidad de cada pórtico
5) Vector de deformación
Según (Razura, 2020) se determina con métodos matriciales, ya que teniendo la matriz de rigidez y el vector de fuerzas sísmicas se hallan los desplazamientos.
Donde:
:
Vector de deformación de la edificación
:
Inversa de la matriz de rigidez de la edificación
:
Carga de la edificación en el centro de masa
6) Fuerzas internas
Según (SALVADOR LLANO, 1999) Las fuerzas internas son las fuerzas que mantienen unidas las partículas que forman un objeto, es decir, refleja la interacción entre ellos. Si un cuerpo rígido está estructuralmente compuesto de varios partes, las fuerzas que mantienen unidos los componentes tambien se definen como fuerzas internas.
Donde:
:
Fuerzas internas de cada pórtico
:
Matriz de rigidez de cada pórtico
:
Matriz de compatibilidad de cada pórtico
:
Vector de deformación de la edificación
Análisis de la estructura porticada
Según (GONZÁLEZ CÁRCELES, 1990) “Primero se adopta el esquema estructural, se decide la disposición de los soportes, entre otros”. La estructura está sometida a fuerzas horizontales , las cuales se analizan con el modelo de corte las vigas tienen una rigidez infinita y las columnas son axialmente indeformables.
Figura 1. Estructura de la edificación - Vista en planta y elevación
Sistema de coordenadas en cada nivel
Primer nivel, segundo nivel y tercer nivel:
Figura 2. Sistema de coordenadas por cada nivel de la edificación
Centro de gravedad del primer y segundo nivel:
Tabla 1 Valores de centro de gravedad del primer, segundo y tercer piso
Pisos |
C.G. (1° y 2° piso) |
C.G. (3° piso) |
X |
3.825 |
3.825 |
Y |
5.145 |
7.375 |
Fuente: Elaboración propia.
Rigidez lateral del pórtico por el método de corte
Según se determina usando el modelo tipo corte, considerando las vigas con EI = ∞.
Propiedades de sección de los elementos:
Definimos las propiedades de cada columna tanto en los ejes “X” y “Y”:
Tabla 2 Propiedades de las columnas en los ejes "X" y "Y"
Sección de columna (a ; b) |
X |
Y |
C-1: 25 x 25 |
707.587 |
707.587 |
C-2: 30 x 30 |
1467.252 |
1467.252 |
C-3: 30 x 33 |
1952.912 |
1613.977 |
C-4: 30 x 25 |
849.104 |
1222.71 |
Fuente: Elaboración propia
Altura de entre pisos:
Rigidez lateral de cada pórtico:
- Rigidez lateral - Pórtico 01:
- Rigidez lateral - Pórtico 02:
- Rigidez lateral - Pórtico 03:
- Rigidez lateral - Pórtico 04:
- Rigidez lateral - Pórtico 05:
- Rigidez lateral - Pórtico 06:
- Rigidez lateral - Pórtico 07:
- Rigidez lateral - Pórtico 08:
Matriz de compatibilidad de deformaciones
Tabla 3 Matriz de compatibilidad actuante en cada pórtico.
D1 |
D2 |
D3 |
D4 |
D5 |
D6 |
D7 |
D8 |
D9 |
PÓRTICOS |
|
1 |
0 |
0 |
0 |
0 |
0 |
7.375 |
0 |
0 |
1 |
P1 |
0 |
1 |
0 |
0 |
0 |
0 |
0 |
7.375 |
0 |
2 |
|
1 |
0 |
0 |
0 |
0 |
0 |
2.915 |
0 |
0 |
1 |
P2 |
0 |
1 |
0 |
0 |
0 |
0 |
0 |
2.915 |
0 |
2 |
|
0 |
0 |
1 |
0 |
0 |
0 |
0 |
0 |
5.145 |
3 |
|
1 |
0 |
0 |
0 |
0 |
0 |
-0.53 |
0 |
0 |
1 |
P3 |
0 |
1 |
0 |
0 |
0 |
0 |
0 |
-0.53 |
0 |
2 |
|
0 |
0 |
1 |
0 |
0 |
0 |
0 |
0 |
1.705 |
3 |
|
1 |
0 |
0 |
0 |
0 |
0 |
-4.03 |
0 |
0 |
1 |
P4 |
0 |
1 |
0 |
0 |
0 |
0 |
0 |
-4.03 |
0 |
2 |
|
0 |
0 |
1 |
0 |
0 |
0 |
0 |
0 |
-1.8 |
3 |
|
1 |
0 |
0 |
0 |
0 |
0 |
-7.38 |
0 |
0 |
1 |
P5 |
0 |
1 |
0 |
0 |
0 |
0 |
0 |
-7.38 |
0 |
2 |
|
0 |
0 |
1 |
0 |
0 |
0 |
0 |
0 |
-5.15 |
3 |
|
0 |
0 |
0 |
1 |
0 |
0 |
-3.83 |
0 |
0 |
1 |
P6 |
0 |
0 |
0 |
0 |
1 |
0 |
0 |
-3.83 |
0 |
2 |
|
0 |
0 |
0 |
0 |
0 |
1 |
0 |
0 |
-3.83 |
3 |
|
0 |
0 |
0 |
1 |
0 |
0 |
-0.03 |
0 |
0 |
1 |
P7 |
0 |
0 |
0 |
0 |
1 |
0 |
0 |
-0.03 |
0 |
2 |
|
0 |
0 |
0 |
0 |
0 |
1 |
0 |
0 |
-0.03 |
3 |
|
0 |
0 |
0 |
1 |
0 |
0 |
3.825 |
0 |
0 |
1 |
P8 |
0 |
0 |
0 |
0 |
1 |
0 |
0 |
3.825 |
0 |
2 |
|
0 |
0 |
0 |
0 |
0 |
1 |
0 |
0 |
3.825 |
3 |
Fuente: Elaboración propia.
Ensamblaje de la matriz de la estructura:
Aplicamos la transpuesta de cada pórtico para
reemplazar en la siguiente fórmula para el ensamblaje de la matriz de la
estructura.
Tabla 4 Matriz de rigidez de la estructura
K edif = |
13611.57 |
-9635.12 |
0.00 |
0.00 |
0.00 |
0.00 |
-4205.23 |
2976.72 |
0.00 |
-9635.12 |
18490.06 |
-8854.94 |
0.00 |
0.00 |
0.00 |
2976.72 |
-15458.36 |
-7264.88 |
|
0.00 |
-8854.94 |
8854.94 |
0.00 |
0.00 |
0.00 |
0.00 |
12481.64 |
7264.88 |
|
0.00 |
0.00 |
0.00 |
12469.92 |
-8826.98 |
0.00 |
-111.41 |
78.86 |
0.00 |
|
0.00 |
0.00 |
0.00 |
-8826.98 |
18101.20 |
9274.22 |
78861.00 |
-161.72 |
82.86 |
|
0.00 |
0.00 |
0.00 |
0.00 |
-9274.22 |
9274.22 |
0.00 |
82.86 |
-82.86 |
|
-4205.23 |
2976.72 |
0.00 |
-111.41 |
78.86 |
0.00 |
378545.53 |
-267958.11 |
0.00 |
|
2976.72 |
-15458.36 |
12481.64 |
78.86 |
-161.72 |
82.86 |
-267958.10 |
480506.20 |
-184714.03 |
|
0.00 |
-7264.88 |
7264.88 |
0.00 |
82.86 |
-82.86 |
0.00 |
-184714.03 |
200914.72 |
Fuente: Elaboración propia
Vector de carga
Qedif |
12 |
18 |
22 |
4 |
6 |
8 |
0 |
0 |
0 |
T |
Vector de desplazamiento “D”
Se determina mediante la siguiente fórmula:
D efic = |
0.01298494 |
0.0170941 |
|
0.02062408 |
|
0.00493841 |
|
0.00652323 |
|
0.00738502 |
|
0.00029793 |
|
-0.0004349 |
|
0.00052712 |
Desplazamientos, fuerzas laterales, y fuerzas cortantes para cada pórtico
Aplicamos las siguientes fórmulas para hallar los desplazamientos, fuerzas laterales y fuerzas cortantes que actúan en cada pórtico:
Tabla 5
Fórmulas para hallar en cada pórtico
Fórmulas |
|
Desplazamientos (m) |
|
Fuerzas laterales |
|
Fuerzas cortantes |
|
|
Fuente: Elaboración propia
Sumamos los valores de las fuerzas cortantes actuantes en los cinco primeros pórticos y los tres últimos pórticos como se muestra, obtenemos:
Resultados
En
el análisis se obtuvo como resultado a comparar con el software Etabs tenemos
los desplazamientos o mejor dicho el vector de desplazamientos que fue
calculado con la siguiente formula: :
Tabla 6
Valores que actúan en cada pórtico
Pórticos |
Desplazamientos D |
Fuerzas Laterales q |
Fuerzas Cortantes V |
|
P-01 |
1 |
0.011 |
1.634 |
5.375 |
2 |
0.014 |
3.741 |
||
P-02 |
1 |
0.012 |
3.167 |
12.27 |
2 |
0.016 |
3.726 |
||
3 |
0.018 |
5.377 |
||
P-03 |
1 |
0.013 |
3.049 |
13.308 |
2 |
0.017 |
4.064 |
||
3 |
0.02 |
6.195 |
||
P-04 |
1 |
0.014 |
2.928 |
14.364 |
2 |
0.019 |
4.409 |
||
3 |
0.022 |
7.027 |
||
P-05 |
1 |
0.015 |
1.223 |
6.685 |
2 |
0.006 |
2.06 |
||
3 |
0.023 |
3.401 |
||
P-06 |
1 |
0.006 |
1.134 |
7.115 |
2 |
0.008 |
2.362 |
||
3 |
0.009 |
3.619 |
||
P-07 |
1 |
0.005 |
1.429 |
6.439 |
2 |
0.007 |
2.146 |
||
3 |
0.007 |
2.864 |
||
P-08 |
1 |
0.004 |
1.438 |
4.447 |
2 |
0.005 |
1.492 |
||
3 |
0.005 |
1.517 |
Fuente: Elaboración propia
Tabla 7
Vector de desplazamientos del cálculo manual
Pisos |
Desplazamiento en X (m) |
Desplazamiento en Y (m) |
Desplazamiento en Z (m) |
Piso 1 |
0.01298494 |
0.00493841 |
-0.00029793 |
Piso 2 |
0.0170941 |
0.00652323 |
-0.0004349 |
Piso 3 |
0.02062408 |
0.00738502 |
-0.00052712 |
Fuente: Elaboración propia.
Tabla 8
Vector de desplazamientos calculo en ETABS
Pisos |
Desplazamiento en X (m) |
Desplazamiento en Y (m) |
Desplazamiento en Z (m) |
Piso 1 |
0.014296 |
0.004553 |
-0.000125 |
Piso 2 |
0.018896 |
0.006037 |
-0.000206 |
Piso 3 |
0.022179 |
0.006957 |
-0.000239 |
Fuente: Elaboración propia.
Discusión
Como se pudo observar en los resultados obtenidos tanto en el cálculo manual y calculo por ETABS existe una ligera diferencia de desplazamientos en los 3 pisos especialmente en el eje Z donde la diferencia en más notoria.
Los desplazamientos en el eje “X” y “Y” de los tres pisos tanto en el cálculo manual y ETABS existe mucha similitud lo cual podemos interpretar que por ese lado el cálculo realizado por nosotros es correcto y que la estructura tiene una buena rigidez lateral, pero en el eje Z de los tres pisos tanto en el cálculo manual y ETABS existe una ligera diferencia de resultados lo cual podemos interpretar que el cálculo manual va demandar mayor rigidez lateral que el cálculo echo en Etabs.
Después de obtener los resultados del cálculo manual y ETABS pudimos observar que en los ejes “X” y “Y” existe mucha similitud con lo cual podemos interpretar que por ese lado el cálculo realizado por nosotros es correcto, pero en el eje “Z” tanto en el cálculo manual y ETABS existe una ligera diferencia de resultados lo cual podemos interpretar que el cálculo manual va demandar mayor rigidez lateral.
Conclusión
- A partir de un proceso de análisis y evaluación manualmente, se obtuvieron resultados satisfactorios sobre la estructura, como viene a ser las deformaciones que presenta, sus desplazamientos, fuerzas y cortes.
- Se obtuvieron satisfactoriamente los resultados de deformaciones que la estructura presenta entre cada piso, observando así que la mayor deformación la presenta por giro con un valor de 0.019 m.
- Al analizar los resultados, observamos los desplazamientos por cada pórtico, viendo de esta manera que el mayor desplazamiento lo presenta el pórtico 2.
- Los resultados obtenidos sobre las fuerzas en los pórticos nos muestran que el pórtico 2 presenta mayor valor que los demás pórticos.
- En la obtención de los resultados de las cortantes, se observa que al sumar las cortantes en el eje “x” dan el mismo resultado que la fuerza externa horizontal de valor 52, lo que significa que los resultados de las cortantes son correctos.
- La estructura que fue modelado en el software de Etabs muestra casi los mismos resultados de deformaciones que los resultados obtenidos manualmente, la variación es muy mínima debido a las decimales que se consideró para el proceso de resolución manual.
Contribución de los autores
MJRV Contribución, Revisión Final del Trabajo (Asesor)
JSLY Contribución, Revisión Final del Trabajo (Asesor)
CBGH Contribución, Redacción de articulo, Materiales y Métodos
JPCM Contribución, Recolección de datos, Modelado en ETABS
JEFP Contribución, Recolección de datos, Solución en Mathcad
IGCM Contribución, Discusiones y Conclusiones
HJTL Contribución, Redacción de articulo, Resumen e Introducción
Conflictos de interés
Los autores declaran que no existe conflicto de interés
Referencias
Basualdo, R. S. (2016). CURSO DE SEGUNDA ESPECIALIZACION EN INGENIERIA SISMORRESISTENTE.
Centro de Estudios y Prevención de Desastres. (7 de Octubre de 2016). Predes. Obtenido de Predes: https://predes.org.pe/15-de-agosto-del-2007-el-terremoto-que-develo-la-falta-de-prevencion-en-el-peru/#:~:text=Las%20zonas%20m%C3%A1s%20afectadas%20fueron,en%20el%20departamento%20de%20Lima.
Dovas Ingenieria. (2021). MÉTODO DE RIGIDEZ EN PÓRTICO. Obtenido de Dovasingenieria.site: http://www.dovasingenieria.site/2020/07/2-metodo-de-rigidez-en-portico.html
GONZÁLEZ CÁRCELES, J. (1990). ANÁLISIS DEL PROCESO DE DISEÑO DE ESTRUCTURAS PORTICADAS. MADRID: ESCUELA TÉCNICA SUPERIOR DE ARQUITECTURA.
Instituto Geofísico del Perú. (3 de Junio de 2022). gob.pe. Obtenido de gob.pe: https://www.gob.pe/institucion/igp/noticias/615872-igp-informa-peru-pais-altamente-sismico
Marcos, R. V. (2023). Apuntes de clase Analisis Estructural II.
Razura, G. (2020). Teoría acerca del Análisis Tridimensional de Edificios bajo Cargas por Sismo mediante el Método Simplificado. Obtenido de Academia.edu: https://www.academia.edu/38123384/Teor%C3%ADa_acerca_del_An%C3%A1lisis_Tridimensional_de_Edificios_bajo_Cargas_por_Sismo_mediante_el_M%C3%A9todo_Simplificado
Salinas Basualdo, R. (15 de 3 de s.f.). FUNDAMENTOS DEL ANALISIS DINAMICO DE ESTRUCTURAS. Lima: PRINCIPIOS COMPUTACIONALES EN INGENIERIA. Obtenido de Marcelopardo.com: https://marcelopardo.com/ensamblaje-matriz-rigidez-portico/
SALVADOR LLANO, R. (1999). ESTATICA APLICADA. Mendoza, Mendoza, Argentina.
© 2023 por los autores. Este artículo es de acceso abierto y distribuido según los términos y condiciones de la licencia Creative Commons Atribución-NoComercial-CompartirIgual 4.0 Internacional (CC BY-NC-SA 4.0)
(https://creativecommons.org/licenses/by-nc-sa/4.0/).