1UNIVERSIDAD NACIONAL DE PIURA FACULTAD DE INGENIERÍA CIVIL TRABAJO DE INVESTIGACIÓN CIENTÍFICA Realizado por: Ing° Carlos Silva Castillo Tema: MÉTODO GENERAL PARA DETERMINAR LA RIGIDEZ LATERAL DE UN PÓRTICO RECTANGULAR 2 ÍNDICE Resumen 01 Parte I. Discusión del método. Cálculos y deducciones 01 1 Preliminares 01 2 Coeficientes de rigidez 05 3 Matriz de rigidez 06 4 Condensación de la matriz de rigidez 06 5 Casos particulares 09 5.1 Pórtico simétrico con columnas de igual altura 09 5.2 Pórtico simétrico de la misma sección transversal: 6 Estudio de los casos extremos 11 6.1 Valor mínimo 11 6.2 Valor máximo 12 Parte II. Algunas Aplicaciones 13 Parte II. Conclusiones y recomendaciones 16 Parte III. Bibliografía 17 Anexos 18 3 MÉTODO GENERAL PARA DETERMINAR LA RIGIDEZ LATERAL DE UN PÓRTICO RECTANGULAR Resumen Dos fueron las razones que impulsaron la elaboración del presente trabajo. Primero, tiene por finalidad desarrollar un método general que permita determinar la rigidez lateral de un pórtico rectangular que trabaja dentro del rango elástico y que se supone axialmente rígido, que sea del material que fuese (concreto armado, acero, madera, etc.) siempre y cuando el pórtico en cuestión esté hecho de un mismo material. La importancia de conocer o calcular la rigidez lateral, estriba en que prácticamente todos los parámetros dinámicos del pórtico (tal como la frecuencia natural de vibración, por ejemplo) dependen de este valor. Segundo, la información existente en el medio se presenta, casi siempre, en otros idiomas que no sea el español, lo que hace tedioso el estudio para un alumno que, además de la Dinámica de Estructuras, tiene que estudiar el idioma Inglés ya que la información brindada en nuestro idioma materno no siempre está completa o se da de una manera sesgada e incompleta. Para su deducción se ha empleado el método de la condensación estática de la matriz de rigidez, se deducen y estudian varios casos particulares, terminando con algunos ejemplos prácticos donde se ilustran el poder del método obtenido. Parte I: Discusión del Método. Cálculos y deducciones 1 Preliminares: Se hace necesario conocer, para ensamblar la matriz de rigidez, los coeficientes de rigidez para un miembro sea este una viga o columna sometido a diferentes efectos tales como la flexión o el corte. Es posible aplicar los diversos métodos del Análisis Estructural y demostrar que en tales circunstancias se tienen los siguientes resultados: Fig N° 01: Coeficientes de rigidez Usando este resultado elemental, estamos en la capacidad de formar la matriz de rigidez para el pórtico mostrado: 4 Fig N° 02: Pórtico empotrado en sus bases y con características generales Donde: H = Altura de la columna de la izquierda h = Altura de la columna de la derecha L = Longitud de la viga E = Módulo de elasticidad (se supone que todo el pórtico es del mismo material) Ic 1 = Momento de inercia de la columna de la izquierda Ic 2 = Momento de inercia de la columna de la derecha Iv = Momento de inercia de la viga El sistema tiene tres grados de libertad (en el futuro, GDL) tal como de muestra en la figura siguiente: Fig N° 03: Grados de libertad del pórtico suponiéndolo rígido axialmente 5 2 Coeficientes de rigidez: Dado que el sistema tiene 3 GDL, la matriz de rigidez será de 3 × 3 y para obtener cada columna de dicha matriz, por definición de coeficiente de rigidez, se dan desplazamientos o rotaciones unitarias según la dirección del GDL empleándose en cada caso los resultados de la Fig. N° 01. Para determinar la primera columna de la matriz de rigidez, se realiza un desplazamiento unitario a lo largo del GDL x 1 , manteniéndose los otros GDL nulos. O sea 1 2 3 1 y 0 x x x = = = , obteniéndose: Fig N° 04: Coeficientes de rigidez cuando 1 2 3 1 y 0 x x x = = = Para determinar la segunda columna de la matriz de rigidez, se realiza una rotación unitaria a lo largo del GDL x 2 , manteniéndose los otros GDL nulos. O sea 2 1 3 1 y 0 x x x = = = , obteniéndose: Fig N° 05: Coeficientes de rigidez cuando 2 1 3 1 y 0 x x x = = = 6 Para determinar la tercera columna de la matriz de rigidez, se realiza una rotación unitaria a lo largo del GDL x 3 , manteniéndose los otros GDL nulos. O sea 3 1 2 1 y 0 x x x = = = , obteniéndose: Fig N° 06: Coeficientes de rigidez cuando 3 1 2 1 y 0 x x x = = = 3 Matriz de rigidez: Habiendo calculado los coeficientes de rigidez, estamos en condiciones de formar la matriz de rigidez del pórtico, obteniendo al ensamblar: 1 2 1 2 3 3 2 2 1 1 2 2 2 2 6 6 12 6 2 4 6 2 4 Ic Ic EIc EIc E H h H h EIc Ic Iv EIv K E H H L L EIc Ic EIv Iv E h L h L 1 | | + | 1 \ . 1 1 | | = + 1 | \ . 1 1 | | + 1 | \ . ¸ ] (1.1) 4 Condensación de la matriz de rigidez: El término condensación se refiere a la disminución en tamaño de un sistema de ecuaciones por la eliminación de determinados GDL. Físicamente equivale a sustituir todo el pórtico por un sistema cuyo modelo matemático está dado por una masa puntual y un resorte de rigidez k* (condensada de la matriz). Para determinar una expresión general de condensación de la matriz de rigidez para una fuerza lateral, como se supone que el trabajo del pórtico es bajo el rango elástico, es aplicable la Ley de Hooke: j ¦ { ; { ; K F A = (1.2) 7 Donde: j ¦ K = Matriz de rigidez { ; A = Matriz de desplazamientos { ; F = Matriz de fuerzas La matriz columna de fuerzas tiene sólo un elemento no nulo que es el primero+ ya que sólo se está suponiendo un comportamiento lateral del pórtico, esto es a la largo del GDL x 1 . Debido a esto las otras dos componentes son nulas (o sea, las que se ubicarían a lo largo de los GDL x 2 y x 3 ) ya que no hay fuerzas en dichas direcciones. Representando la ecuación (1.2) en su forma particionada, se obtiene: 0 aa ab a ba bb b k k f k k A 1 ¦ ¹ ¦ ¹ = ´ ` ´ ` 1 A ¹ ) ¸ ] ¹ ) (1.3) Desarrollando el producto matricial, obtenemos el sistema: j ¦ { ; j ¦ { ; { ; aa a ab b k k f A + A = (1.4) j ¦ { ; j ¦ { ; { ; 0 ba a bb b k k A + A = (1.5) Despejando { ; b A de la ecuación (1.5): { ; j ¦ j ¦ { ; 1 b bb ba a k k ÷ A = ÷ A (1.6) Reemplazando en (1.4) y factorizando { ; a A : j ¦ j ¦ j ¦ j ¦ , { ; { ; 1 aa ab bb ba a k k k k f ÷ ÷ A = (1.7) Definiendo la expresión entre paréntesis como matriz condensada k* , y arreglando para que tenga el aspecto de la ley de Hooke: j ¦ { ; { ; * a k f A = (1.8) De donde, por definición de matriz condensada: j ¦ j ¦ j ¦ j ¦ j ¦ 1 * aa ab bb ba k k k k k ÷ = ÷ (1.9) 8 Fig N° 07: Interpretación física de la condensación de la matriz de rigidez Condensando la matriz de rigidez del pórtico dada por la ecuación (1.1): 1 1 1 2 1 2 1 2 3 3 2 2 2 2 2 2 6 4 6 6 * 12 6 2 4 Ic Iv EIv EIc E H L L Ic Ic EIc EIc H k E EIc H h H h Ic EIv Iv E h L h L ÷ 1 | | 1 + | 1 1 \ . | | 1 1 = + ÷ 1 | 1 1 \ . ¸ ] | | 1 + 1 | 1 ¸ ] \ . ¸ ] (1.10) Llevando a cabo las operaciones matemáticas y simplificando algebraicamente: , { ; , 3 2 4 4 2 2 2 2 2 2 2 1 1 1 2 1 1 3 3 2 2 1 2 1 2 12 ( ) (3 ) 3 4 3 4 * 3 4 4 E H LIc HIv LIc h Ic Iv HIv LIc hIc H Iv HLIc Iv h hH H h L Ic k h H hHIv LIv hIc HIc L Ic Ic 1 + + + + + + + + ¸ ] = 1 + + + ¸ ] (1.11) Dividiendo a numerador y denominador entre Iv tenemos: , , 2 3 2 4 4 2 2 2 2 1 1 2 1 1 2 1 3 3 2 1 2 1 2 12 ( ) (3 ) 3 4 3 4 * 3 4 4 Ic Ic E H LIc H L h Ic HIv LIc hIc H Iv HLIc h hH H h L Iv Iv k Ic Ic h H hHIv L hIc HIc L Iv ¦ ¹ 1 + + + + + + + + ´ ` 1 ¸ ] ¹ ) = 1 + + + 1 ¸ ] (1.12) Esta expresión es precisamente la que andamos buscando, es aquella que nos da el valor de la rigidez lateral del pórtico en estudio. Observamos su aspecto exterior complicado y de forma poco agradable para su memorización, razón por la cual se prefiere en un ejercicio práctico comenzar por la formación y ensamblaje de la matriz de rigidez y luego llevar a cabo el cálculo de su condensada, antes que usar una fórmula general de muy complicado aspecto. Pero veamos algunos casos particulares en que esta expresión adopta formas más manejables y simples. En lo que sigue, siempre que se hable de rigidez del pórtico se refiere al valor que se obtiene al condensar la matriz de rigidez, por lo simplemente hablaremos de rigidez. O sea * k k = 9 5 Casos particulares: 5.1 Pórtico simétrico con columnas de igual altura: Para el caso en que las dos columnas del pórtico tienen la misma altura, teniendo cada una diferente sección transversal, esto es si = ℎ, reemplazando en (1.6) se obtiene: Fig N° 08: Pórtico simétrico, con columnas de diferente momento de inercia , 2 3 2 4 4 3 2 2 1 1 2 1 1 2 1 6 2 2 1 2 1 2 12 ( ) (3 ) 3 11 3 4 4 Ic Ic E h LIc h L h Ic hIv LIc hIc h Iv h LIc h L Iv Iv k Ic Ic h h Iv hL Ic Ic L Iv 1 | | + + + + + + 1 | \ . ¸ ] = 1 + + + 1 ¸ ] (1.13) Podemos observar que aún conserva su aspecto poco amigable aunque sólo un poco simplificada. 5.2 Pórtico simétrico con columnas de igual altura y de la misma sección transversal: Además de la condición = ℎ, debe cumplirse que 1 2 Ic Ic Ic = = . Imponiendo esta condición en (1.7), obtenemos: Fig N° 09: Pórtico simétrico, con columnas de igual momento de inercia 10 3 12 6 3 2 Ic EIc h L Iv k Ic h h L Iv | | + | \ . = | | + | \ . (1.14) Dividiendo a numerador y denominador entre L y reordenado, se tiene: 3 6 12 3 2 h Ic EIc L Iv k h Ic h L Iv | | + | = | | + \ . (1.15) Para simplificar la ecuación anterior y hacerla más manejable, introducimos los siguientes parámetros Ic Iv = y h j L = , la ecuación anterior queda: 3 12 6 3 2 EIc j k h j | | + = | + \ . (1.16) Llamando 6 3 2 j j + = + , estudiemos la variación de como función de los parámetros y j , se obtiene el gráfico siguiente para diferentes valores del parámetro j : Fig N° 10: Variación del parámetro 11 Analizando este gráfico llegamos a determinar los valores extremos que puede tomar este parámetro , el mínimo valor es 0.5 y el máximo es 2 0.5 2 s s (1.17) Este resultado es muy importante, quiere decir que el valor de la rigidez lateral de un pórtico simétrico rectangular que tiene idénticas secciones transversales en las columnas, está limitado entre dos extremos. De la ecuación (1.16) y de la definición del parámetro : 3 12EIc k h = (1.18) Introduciendo los valores extremos que puede tomar el parámetro y que se resumen en la ecuación (1.17), obtenemos finalmente los valores extremos de la rigidez: 3 3 6 24 EIc EIc k h h s s (1.19) 6 Estudio de los casos extremos: 6.1 Valor mínimo m n 3 6 í EIc k h = (1.20) Matemáticamente hablando, se tiene este valor ante dos posibilidades: a. 0 j = que ocurre, a su vez, cuando la altura del pórtico es nula ( 0 h = ) o despreciable en comparación con la longitud de la viga ( L = · ) b. = · que solamente puede ocurrir si el momento de inercia de la viga es nula o despreciable ( 0 Iv = ). Esto quiere decir, cuando la viga no tiene rigidez alguna. Formalmente hablando, este se da cuando la rigidez de la viga ( EIv ) es cero o despreciable ( 0 EIv = ) Se podría haber llegado más rápidamente a este resultado analizando el pórtico en estudio para el caso en que la viga carece de rigidez ( 0 EIv = ), como se puede observar en la siguiente Fig N° 11. En este caso, la rigidez lateral está dada por la suma de rigideces de columnas en ausencia de viga: 3 3 3 6 columnas EIv EIv k h h = = _ (1.21) 12 Fig N° 11: Caso del pórtico con viga de rigidez despreciable 6.2 Valor máximo m 3 24 áx EIc k h = (1.22) Matemáticamente hablando, se tiene este valor ante dos posibilidades: a. j = · que ocurre, a su vez, cuando la altura del pórtico es infinita ( h = · ) o muy grande en comparación con la longitud de la viga ( 0 L = ) b. 0 = que solamente puede ocurrir si el momento de inercia de la viga es infinita o muy grande en comparación con el de las columnas ( Iv = ·). Esto quiere decir, cuando la viga es completamente rígida. Formalmente hablando, este se da cuando la rigidez de la viga ( EIv ) es muy grande o infinita 1 ( EIv = · ). Es posible que también se halla llegado a este resultado cuando el momento de inercia de las columnas es nula o despreciable ( 0 EIc = ), pero por obvias razones de estabilidad y seguridad estructural es absurdo que se tenga una columna carente de rigidez ¡Pues el pórtico se desploma! Se podría haber llegado más rápidamente a este resultado analizando el pórtico en estudio para el caso en que la viga se puede suponer como completamente rígida ( EIv = · ), como se puede observar en la siguiente Fig N° 12. En este caso, la rigidez lateral está dada por la suma de rigideces de columnas: 3 3 12 24 columnas EIv EIv k h h = = _ (1.23) 1 Esto es lógico por cuanto la rigidez de un elemento no sólo está dada por la forma de la sección transversal sino también por su módulo de elasticidad. Pensemos en dos barras de las mismas dimensiones pero una de plástico y la otra de acero. Tienen el mismo momento de inercia, pero diferente módulo de elasticidad, razón por la cual más rígida será la barra de acero que la de plástico, ya que tiene mayor valor del módulo de elasticidad. 13 Fig N° 12: Caso del pórtico con viga completamente rígida Parte II: Algunas aplicaciones Ejemplo 01: Determinar la rigidez lateral del siguiente pórtico de concreto armado con 2 ' 210 / c f kg cm = las dimensiones de los elementos son Col 1 = 0.40m×0.40m, Col 2 = 0.35m×0.35m y viga = 0.30m×0.60m. La manipulación matemática se detalla en la siguiente hoja de cálculo, usando la Ec. (1.12): 14 Ejemplo 02: Determinar la frecuencia natural de vibración para un pórtico de concreto armado con 2 ' 280 / c f kg f cm = ÷ las dimensiones de los elementos son Col 1 = 0.50m×0.45m, Col 2 = 0.45m×0.40m y viga = 0.30m×0.50m. Se puede suponer una masa distribuida uniformemente en el techo de 25 t/m, despreciar la masa del pórtico. La manipulación matemática se detalla en la siguiente hoja de cálculo, usando la Ec. (1.13): Calculamos la masa como una cantidad puntual: 25 4.20 105 105000 m t kg = × = = Finalmente hallamos la frecuencia natural de vibración: 2 3447.6815 / 0.1812 rad/ 105000 k kg s s m kg = = = 15 Ejemplo 03: Hallar una expresión general que permita determinar la rigidez lateral de un pórtico rectangular simétrico, donde tanto la viga como las columnas tienen el mismo momento de inercia I; además, la longitud de la viga es el doble del valor para la altura de las columnas. Este ejemplo se resuelve aplicando directamente la Ec. (1.16), donde los valores de los parámetros ahí encontrados son 1 Ic I Iv I = = = y 1 2 2 h h j L h = = = . Reemplazando estos resultados en la citada ecuación: 3 3 1 6 1 12 96 2 1 7 3 2 1 2 EI EI k h h | | × + | = = | | × + × \ . Ejemplo 04: Hallar el gráfico tridimensional de la variación del parámetro en función de y j. O sea, graficar la función 6 3 2 j j + = + 16 Parte III: Conclusiones y recomendaciones 1. El conocimiento del parámetro llamado rigidez lateral es fundamental para la determinación de muchos parámetros mecánicos de un pórtico, sin cuyo valor nada podría hacerse para calcular el periodo de vibración, por ejemplo. 2. Se ha desarrollado un método general para determinar la rigidez lateral de un pórtico rectangular axialmente rígido y que trabaja en el rango elástico, empleando fórmulas y métodos algebraicos en vez de métodos matriciales, pero deducidos de éstos, fundamentalmente del método estático de condensación de la matriz de rigidez. 3. Los ejemplos resueltos en este trabajo empleando las fórmulas algebraicas obtenidas para el cálculo de la rigidez lateral de un pórtico rectangular, han sido verificados empleando los softwares Mathematica ® y la hoja de cálculo Excel ® . Los resultados de dichas verificaciones se adjuntan en el anexo correspondiente. 4. Los ejemplos planteados y resueltos, tuvieron (dentro de la precisión de los cálculos y del programa usado) los mismos resultados, quedando verificado que los métodos algebraicos y deducidos en este trabajo son fiables y seguros. 17 Parte IV: Bibliografía 1.0 PAZ, Mario – LEIGH, William. “STRUCTURAL DYNAMICS – Theory and Computation – Updated with SAP2000” KLUGER ACADEMIC PUBLISHERS Printed in the United States of America – 2004 – 844 pp 2.0 ZALKA, Karoly A. “GLOBAL STRUCTURAL ANALYSIS OF BUILDINGS” E & FN SPON London – 2000 – 334 pp 3.0 CHOPRA, Anil K. “DYNAMICS OF STRUCTURES - Theory and Aplications to Earthquake Engineering” PRENTICE HALL Printed in the United States of America – 1995 – 730 pp In[4]:= K 12 G [Ic1 / H 3 Ic2 / h 3 j 6 G Ic1 / H 2 6 G Ic2 / h 2 6 G Ic1 / H 2 4 G (Ic1 / H Iv / L) 2 G Iv / L 6 G Ic2 / h 2 2 G Iv / L 4 G (Ic2 / h Iv / L) Out[4]= 12 G Ic1 H 3 Ic2 h 3 , 6 G Ic1 H 2 , 6 G Ic2 h 2 , 6 G Ic1 H 2 , 4 G Ic1 H Iv L , 2 G Iv L , 6 G Ic2 h 2 , 2 G Iv L , 4 G Ic2 h Iv L In[5]:= G 15000 210 Out[5]= 15000 210 In[14]:= Ic1 40 40 3 / 12 Out[14]= 640 000 3 In[15]:= Ic2 35 35 3 / 12 Out[15]= 1500 625 12 In[16]:= Iv 30 60 3 / 12 Out[16]= 540 000 In[17]:= H 420 Out[17]= 420 In[18]:= h 310 Out[18]= 310 In[19]:= L 515 Out[19]= 515 In[20]:= K Out[20]= 117 152 040 625 10 21 4379 277 , 16000 000 10 21 7 , 112 546 875 210 961 , 16000 000 10 21 7 , 202 000 000 000 10 21 103 , 3240 000 000 210 103 , 112 546 875 210 961 , 3240 000 000 210 103 , 278 162 187 500 210 3193 In[21]:= [12 G [Ic1 / H 3 Ic2 / h 3 jj [ 6 G Ic1 / H 2 6 G Ic2 / h 2 j. Inverse_] 4 G (Ic1 / H Iv / L) 2 G Iv / L 2 G Iv / L 4 G (Ic2 / h Iv / L) ¿_. 6 G Ic1 / H 2 6 G Ic2 / h 2 In[22]:= __ 215 247 461 950 557 846 875 10 21 8645 537 871 461 967 28967 052 932 409 375 105 2 117 626 365 598 122 __ In[24]:= N_ 215 247 461 950 557 846 875 10 21 8645 537 871 461 967 28967 052 932 409 375 105 2 117 626 365 598 122 _ Out[24]= 15396.2 2 Untitled-1 In[4]:= K 12 G [Ic1 / H 3 Ic2 / h 3 j 6 G Ic1 / H 2 6 G Ic2 / h 2 6 G Ic1 / H 2 4 G (Ic1 / H Iv / L) 2 G Iv / L 6 G Ic2 / h 2 2 G Iv / L 4 G (Ic2 / h Iv / L) Out[4]= 12 G Ic1 H 3 Ic2 h 3 , 6 G Ic1 H 2 , 6 G Ic2 h 2 , 6 G Ic1 H 2 , 4 G Ic1 H Iv L , 2 G Iv L , 6 G Ic2 h 2 , 2 G Iv L , 4 G Ic2 h Iv L In[25]:= G 15000 280 Out[25]= 30000 70 In[26]:= Ic1 50 45 3 / 12 Out[26]= 759 375 2 In[27]:= Ic2 45 40 3 / 12 Out[27]= 240 000 In[28]:= Iv 30 50 3 / 12 Out[28]= 312 500 In[29]:= H 325 Out[29]= 325 In[30]:= h 325 Out[30]= 325 In[31]:= L 420 Out[31]= 420 In[32]:= K Out[32]= 14277 600 70 2197 , 109 350 000 70 169 , 69120 000 70 169 , 109 350 000 70 169 , 20882 500 000 10 7 13 , 312 500 000 10 7 , 69120 000 70 169 , 312 500 000 10 7 , 16189 000 000 10 7 13 In[33]:= [12 G [Ic1 / H 3 Ic2 / h 3 jj [ 6 G Ic1 / H 2 6 G Ic2 / h 2 j. Inverse_] 4 G (Ic1 / H Iv / L) 2 G Iv / L 2 G Iv / L 4 G (Ic2 / h Iv / L) ¿_. 6 G Ic1 / H 2 6 G Ic2 / h 2 Out[33]= 791 359 496 773 200 70 188 392 976 291 In[34]:= N_ 791 359 496 773 200 70 188 392 976 291 _ Out[34]= 35144.6 2 Untitled-1
Comments
Report "Método general para determinar la rigidez lateral de un pórtico - sísmica (1)"