ANLISIS DINMICO DE UNA ESTRUCTURA UTILIZANDO STAAD Autor: Dr. Joaqun E. Torre M.

, Honduras, Septiembre 2011 INTRODUCCIN: Staad.Pro es un programa de ordenador con la capacidad de ejecutar un anlisis de super posicin modal espectral para lo que es necesario especificar o definir el espectro que se desea utilizar, sin embargo, los resultados que se obtienen de este anlisis, el programa los presenta en valores absolutos para todos los esfuerzos de cada uno de los miembros que componen la estructura de manera que es incorrecto hacer combinaciones de carga en forma directa utilizando estos resultados. En el procedimiento descrito en este documento, se presenta la metodologa necesaria para hacer uso correcto de la facilidad que ofrece el programa en un ejemplo seleccionado. Adems de hacer las comparaciones entre los resultados por el anlisis dinmico utilizando la facilidad de STAAD.Pro y el procedimiento descrito en este documento, se comparan tambin los mismos resultados con los obtenidos por un anlisis dinmico por super posicin modal espectral de otro programa desarrollado por el autor. DESCRIPCIN DEL CASO DE ESTUDIO: En la figura 1 se muestra la estructura seleccionada en la que se puede observar que la misma, no es simtrica. Consta de dos intercolumnios en cada direccin de 6.00m y 8.00m en el sentido Z y 7.50m y 9.00m en el sentido X con vigas de 0.60x0.40m de seccin en cada uno de los ejes y columnas de 0.60x0.60m en la interseccin de cada eje. Cuenta con tres niveles con alturas de 4.00m, 3.50m y 3.50m y las losas de entrepiso se han modelado por elementos finitos. El origen del sistema coordenado de la estructura es el eje de la columna de la esquina superior izquierda si se ve la estructura en planta y su centro de masas respecto de este origen es (x,z)=8.18,6.90 para cada nivel tomando en cuenta que adems de su peso propio, las losas de entrepiso reciben una carga 2 permanente uniformemente distribuida en toda su rea de 400Kg/m . Los vnculos a tierra se han modelado como empotramientos perfectos y los elementos tipo placa con los que se representa las losas de entrepiso se han indicado de 0.20m de espesor. Para el anlisis dinmico, que se ha realizado nicamente en la direccin X, se le ha especificado al programa la definicin del espectro indicado en el Cdigo Hondureo de la Construccin (CHOC-2008) en el apartado 1.3.6.5.7 con una aceleracin pico del suelo de 0.25g, considerando la estructura apoyada sobre roca (S=1.00), con categora y factor de importancia I=1.00 y con un coeficiente estructural RW=8. Una vez obtenidos los resultados del anlisis dinmico de la estructura, se procede a encontrar la suma de las cortantes en cada columna en el nudo correspondiente al nivel del diafragma en consideracin para por suma algebraica encontrar la fuerza ssmica asignada a la estructura en ese nivel. Esas fuerzas, se colocan en el mismo modelo, en un nudo creado al efecto con las coordenadas del centro de masas del nivel del diafragma correspondiente. En este caso, es el mismo para todos los niveles en virtud de que cada nivel es exactamente igual en masas y colocacin de las mismas y son los valores antes indicados. Es importante que los nudos creados para representar el

centro de masas se conviertan en nudos maestros y el resto de los nudos del nivel del diafragma correspondientes sean esclavizados a este para simular un diafragma rgido en el plano X-Z.

Figura 1: Modelo de la estructura en anlisis

Figura 2: Modelo sombreado de la estructura en anlisis

Corriendo de nuevo el modelo con estas fuerzas, obtenemos los mismos resultados en carga axial, cortantes y momentos pero con el signo correcto y consecuentemente, aptos para combinarse con otras condiciones primarias de carga.

Para efectos comparativos, tambin se le ha practicado a la misma estructura objeto de anlisis en este documento, un anlisis modal espectral utilizando un programa tipo macro escrito en Visual Basic con interfase en Excel en el que su algoritmo fue escrito realizar un anlisis de superposicin modal con un grado de libertad por piso, pero tomando la rigidez de cada piso de un anlisis esttico realizado en STAAD.Pro para obtener los desplazamientos conocidas las fuerzas que con valores arbitrarios se han aplicado a los centros de masa. As mimo, se ha realizado el anlisis ssmico pseudo esttico como lo prescribe el CHOC-2008 en su artculo 1.3.5.2, con la intencin de hacer las comparaciones respectivas. RESULTADOS OBTENIDOS: La figura 3 muestra el diagrama de momentos de un eje de acuerdo con el anlisis dinmico realizado, donde se puede observar que todos los valores estn en positivo. Igual resultado se observa para las cortantes y las cargas axiales.

Figura 3: Diagrama de momentos de un eje por anlisis dinmico Los perodos de vibracin de la estructura son:
MODO 1 2 3 4 5 6 FRECUENCIA (ciclos/sec) 1.640 1.712 2.219 5.500 5.655 5.730 PERODO (sec) 0.60962 0.58426 0.45072 0.18181 0.17683 0.17451

Y se muestra en rojo el correspondiente al perodo fundamental de la estructura del cual se muestra su silueta para la estructura completa en la figura 4 donde se puede apreciar la torsin como consecuencia de la excentricidad resultante de no ser una estructura perfectamente simtrica.

Figura 4: Vista en alzado y planta del primer modo de vibracin La participacin de masas y la cortante total en la base es:
MODO 1 2 3 4 5 6 X 86.56 0.06 0.25 7.51 2.47 0.11 Y Z 0.00 0.04 0.00 87.23 0.00 0.23 2.10 0.06 4.14 1.63 0.40 8.01 SUMM-X 86.558 86.614 86.865 94.380 96.853 96.966 SUMM-Y 0.001 0.003 0.003 2.106 6.249 6.651 TOTAL CQC SUMM-Z 0.045 87.271 87.500 87.557 89.190 97.196 V= X 37,189.40 2 5.07 146.35 5,413.25 1,781.07 81.78 ---------37,916.19 Kg

Y la suma de las fuerzas ssmicas en cada diafragma que se encontr de este anlisis es la siguiente: Nivel (m) 11.00 7.50 4.00 Fx (Kg) 18,783.90 12,097.20 7,035.09 Fz (Kg) 1,352.12 613.42 500.06

Para un cortante total en la base de 37,916.19Kg (0.0444W) en la direccin X y 2,465.60Kg en la direccin Z. Si la estructura fuera perfectamente simtrica, los valores en Z seran cero. El peso total de la estructura es W=854,117.07Kg Estas fuerzas se aplican a los centros de masa en el nivel correspondiente para un anlisis esttico de los que se obtienen los siguientes resultados: Desplazamientos en el sentido X por piso: Nivel 11.00 7.50 4.00 Anlisis dinmico (mm) 5.934 4.380 2.089 Anlisis esttico (mm) 6.063 4.440 2.104 Diferencia absoluta (mm) 0.129 0.060 0.015 Diferencia relativa (%) 2.17 1.37 0.71

Las diferencias observadas corresponden a errores por la aproximacin en el clculo del los centros de masa en cada nivel, sin embargo, son diferencias aceptables.

Momentos en el mismo eje de la figura 3:

Figura 5: Diagrama de momentos del mismo eje mostrado en la figura 3 Seleccionando una columna al azar para hacer la comparacin numrica, obtenemos los siguientes resultados: Tipo de anlisis Dinmico Esttico Diferencia Dif (%) Carga axial (Kg) +5,058.85 -5,427.52 368.67 7.28 Vx (Kg) +3.648.408 +3,820.009 171.601 4.70 Vz (Kg) +366.500 -386.455 19.955 5.44 Torsin (Kg-m) +88.257 -24.371 63.886 72.39 Momentos My (Kg-m) +970.130 +1,018.232 48.102 4.96 Momentos Mz(Kg-m) +10,685.949 +11,110.265 424.316 3.97

Las diferencias se muestran en valores absolutos y la diferencia porcentual es con respecto a los valores del anlisis dinmico. Ntese el cambio de signo en la carga axial y la torsin de la columna as como la diferencia en la torsin como consecuencia de la correcta aplicacin de la carga ssmica. Similares diferencias relativas se obtienen con las dems columnas y vigas. Las diferencias en los valores de las cortantes y los momentos flectores en cada eje de la columna seleccionada son del orden del 5% los cuales son aceptables en trminos de la magnitud, sin embargo, los resultados del anlisis equivalente presentan los valores correctos en cuanto a la direccin del esfuerzo. Los valores mostrados en la tabla anterior corresponden a los valores en el nudo de inicio de la columna seleccionada y en el nudo del extremo final, los signos de cortantes y momentos son completamente diferentes como se aprecia entre los diagramas de momentos de las figuras 3 y 5. En un programa desarrollado por el autor utilizando Excel como interfase y en Visual Basic un algoritmo que permita la generacin del sistema de equilibrio dinmico para poder encontrar los valores propios y realizar la superposicin modal, para el mismo espectro y la misma aceleracin pico del suelo, encontramos los siguientes resultados: Perodos de vibracin: T = 0.6142 seg T = 0.2078 seg T = 0.1378 seg Los cuales comparan muy bien con los obtenidos por STAAD.Pro. El programa en Excel solo toma en cuenta un grado de libertad por piso, pero la rigidez de cada piso la

calcula de los desplazamientos como producto de la aplicacin de fuerzas con valor arbitrario que se han colocado en el centro de masas para ese efecto. Como la solucin solo toma un grado de libertad por piso (desplazamiento en el sentido X) no hay un grado de libertad por torsin de manera que los efectos de la misma no son apreciables ni inciden en este anlisis y su objetivo es comparar las fuerzas smicas a aplicar sobre el modelo. Tipo de anlisis STAAD.Pro Programa del autor Fuerza ssmica (Kg) por piso +4.00 +7.50 +11.00 7,035.09 12,097.20 18,783.90 9,824.69 14,191.51 17,593.28 Desplazamientos (mm) +4.00 +7.50 +11.00 5.934 4.380 2.089 5.674 4.811 3.196

El cuadro comparativo anterior muestra algunas diferencias que son producto de la aproximacin realizada por el programa en EXCEL, no obstante, estas diferencias no son significativas, por lo que podemos aceptar los resultados obtenidos por STAAD.Pro como correctos ya que adems de haber resuelto el caso con todos los grados de libertad necesarios y representativos, los valores absolutos entran en el rango de aceptacin al compararse bien. Si al mismo modelo le aplicamos el anlisis pseudo esttico de acuerdo con UBC-94 que es el mismo indicado en el CHOC-2008 y para la misma aceleracin pico del suelo (0.25g), se obtiene que el perodo fundamental encontrado por el criterio de Rayleigh es T=0.606s que compara muy bien con el valor de T=0.60962s que se muestra en la tabla con los resultados del anlisis modal, y la cortante en la base es Vx=46,611.12Kg = 0.0546W equivalentes a 1.23 veces la misma cortante por el anlisis dinmico. DISCUSIN: El anlisis modal espectral realizado por STAAD.Pro no debe de ser usado en forma directa ni mucho menos combinarse con otras condiciones de carga debido a que los valores mostrados del primero estn en valor absoluto y su combinacin ser incorrecta, sin embargo, los resultados de un anlisis equivalente si es posible combinarlos. El anlisis modal espectral que realiza el programa, no tiene la herramienta para poder aplicar el 5% de excentricidad accidental que exige el CHOC-2008, de manera que cuando se crean los nudos correspondientes al centro de masas estos debern tener en sus coordenadas, la excentricidad de 5% en cada direccin y en cada sentido para cumplir con este requisito, cosa que no se tom en cuenta en el ejemplo de este documento para permitir la comparacin. La distribucin de fuerzas ssmicas en la estructura es completamente diferente entre el anlisis modal y el anlisis pseudo esttico, como se aprecia en la grfica 1 a continuacin, y esa es la razn por la que la norma tcnica exige este anlisis como investigacin de la respuesta de la estructura para edificios de alturas superiores a 75m o las que presentan irregularidad geomtrica como se define en el CHOC-2008 en el artculo 1.3.4.8.3 Ntese en la grfica 1 que la distribucin por el mtodo pseudo esttico es lineal y en el mtodo dinmico la distribucin no es lineal lo que era de esperar ya que esta ltima es el producto de la superposicin de los modos de vibracin los cuales se han combinado con el criterio de la combinacin cuadrtica completa (CQC). Tambin es importante notar, que las mayores diferencias entre la distribucin de las fuerzas, se aprecian en los niveles altos. El caso antes estudiado, que se ha seleccionado sencillo y pequeo para mejor claridad en el anlisis de los resultados, parece presentar diferencias aceptables entre los mtodos, pero en la medida en que la estructura presenta mayores alturas, la diferencia se vuelve importante.

12.00 N 10.00 8.00 i v 6.00 e 4.00 2.00 l 0.00 0.00 5,000.00

Dinmico Esttico

10,000.00 15,000.00 20,000.00 25,000.00 Fuerza en diafragmas (Kg) Grfico 1: Distribucin comparativa de fuerzas por sismo en los diafragmas CONCLUSIONES: 1. Queda evidencia que la combinacin de los resultados obtenidos por el anlisis modal espectral que realiza el programa STAAD.Pro con otras condiciones primarias de carga no es correcto, pues corresponde a una superposicin de fuerzas o esfuerzos en los que los signos obtenidos del primero no necesariamente son correctos o no son consistentes a la direccin del estmulo aplicado. Se vuelve crtico especialmente en el anlisis de las columnas, pues aquellas que deben de estar en traccin se combinan errneamente en compresin con el momento resultante que tambin puede estar en signo incorrecto con los correspondientes de otras condiciones primarias de carga lo que puede llevar a sobre diseo en algunos casos y a un sub dimensionamiento en otros. 2. Utilizando el procedimiento descrito en este documento, es posible hacer uso del anlisis dinmico aunque eso consista en algn trabajo adicional, sin embargo, el mismo es necesario para poder tener una completa y correcta visin de la respuesta de la estructura en el anlisis realizado y adems asegurar fcilmente, la correcta combinacin o superposicin con otros escenarios de carga. 3. Es necesario observar el artculo 1.3.5.5 del CHOC-2008 asignando en las coordenadas del nudo que representa el centro de masas, la excentricidad accidental correspondiente al 5% de la dimensin del edificio en la direccin perpendicular a la del anlisis en cada direccin y en cada sentido. 4. Es importante tener claro las limitaciones del programa seleccionado para realizar los anlisis as como la metodologa para cubrir o superar estas limitaciones, especialmente cuando las normas tcnicas exigen este tipo de anlisis cuando la estructura presenta ciertas condiciones o caractersticas. BIBLIOGRAFA: 1. Comit Tcnico (2008) Cdigo Hondureo de la Construccin, Normas Tcnicas Colegio de Ingenieros Civiles de Honduras, Tegucigalpa, Honduras 2. Chopra, A. (2007) Dynamics of Structures 3 ed. Pearson-Prentice Hall, New Jersey, U.S.A. 3. International Conference of Building Officials (1994) Uniform Building Code Whittier, California, U.S.A. 4. Research Engineers, Intl. (2007) Technical Reference, STAAD.Pro Bentley Solutions Center, California, U.S.A. 5. Torre J. (2010) Programa para anlisis dinmico con un GDL por piso Tegucigalpa, Honduras 6. Torre J. (2010) Curso de Dinmica estructural Programa de Maestra en Ingeniera estructural, UNITEC, Tegucigalpa, Honduras.