Source: https://www.scribd.com/doc/95088788/Libro-de-Dinamica-Estructural
Timestamp: 2016-12-02 23:39:10+00:00

Document:
Libro de Dinamica Estructural
BrowseInterestsBiography & MemoirBusiness & LeadershipFiction & LiteraturePolitics & EconomyHealth & WellnessSociety & CultureHappiness & Self-HelpMystery, Thriller & CrimeHistoryYoung AdultBrowse byBooksAudiobooksArticlesSheet MusicBrowse allUploadSign inJoinDepartamento de Ingeniería CivilFacultad Regional Paraná – Universidad Tecnológica Nacional ANALISIS DE ESTRUCTURAS BAJO ACCIONES DINÁMICAS Arturo M. Cassano Editorial de la Universidad Tecnológica Nacional - edUTecNe
ISBN 978-987-25360-9-1
Introducción - Bibliografía 1 INTRODUCCIÓN Esta publicación pretende brindar en una forma rápida y sencilla los conceptos básicos de la dinámica de estructuras aplicada a las construcciones civiles, enfocada desde un punto de vista numérico. Su organización en capítulos se basa en la secuencia de unidades temá-
ticas del programa de la asignatura homónima que se dicta en la carrera de Ingeniería Civil de la Facultad Regional Paraná de la Universidad Tecnológica Nacional. Todos los temas se desarrollan como síntesis de libros y publicaciones de otros autores (ver Bibliografía) solo que organizados según la conveniencia que mejor se ajusta al curso impartido. El Capítulo 1 hace una muy breve reseña a conceptos sobre sismología y su relación con la ingeniería estructural, caracterización, registros, etc. El Capítulo 2 brinda una aproximación a los conceptos básicos de la dinámica estructural, la importancia de la masa, la relación entre velocidad de carga y de reacción de una estruc-
tura. Avanza sobre los modelos estructurales dinámicos, grados de libertad y la discretiza-
ción espacial y temporal. El Capítulo 3 está dedicado a la caracterización y cálculo de osciladores de un grado de libertad dinámico. Son tratadas las vibraciones libres y forzadas, éstas últimas con cargas armónicas y arbitrarias. Se presenta la resolución numérica mediante integración directa. El Capítulo 4 está dedicado a las estructuras con múltiples grados de libertad dinámicos. Se estudian las vibraciones libres y el cálculo de modos y frecuencias propias, también se presenta la resolución de los sistemas dinámicos mediante diversos métodos como son: descomposición y superposición modal, integración directa y respuesta máxima mediante espectros de respuesta. En el Capítulo 5 se realiza una aproximación al tratamiento dinámico de los efectos del viento sobre las construcciones teniendo como base el reglamento CIRSOC 102. Por último en el Capítulo 6 son tratados los efectos sísmicos y el análisis estructural, desde el punto de la dinámica estructural, enfatizando en la definición de la acción sísmica y la comparación de las fórmulas vistas en capítulos anteriores con los reglamentos vigentes. Un agradecimiento especial a Marlene Jaurena, que con paciencia y dedi-
cación realizó la transcripción de los manuscritos de este trabajo. Introducción - Bibliografía 2 TABLA DE CONTENIDOS Capítulo 1 - NOCIONES BASICAS DE SISMOLOGIA 1-1 Causas que generan los terremotos o sismos 1-1 Los sismos desde el punto de vista de la ingeniería y su caracterización 1-2 Esteva y Rosenblueth: 1-2 Donovan: 1-3 Esteva y Villaverde: 1-3 Esteva: 1-3 Registro de ondas sísmicas. Parámetros utilizados y mapas de riesgo sísmico 1-3 Capítulo 2 - CONCEPTOS BASICOS DE DINAMICA ESTRUCTURAL 2-1 Definición de la acción dinámica 2-1 Acciones y fuerzas dinámicas 2-2 Importancia de la masa en el problema dinámico 2-3 Velocidad de reacción de una estructura 2-5 Modelos dinámicos característicos 2-6 Métodos de modelización dinámica 2-8 Discretización espacial de las estructuras 2-9 Método de las masas concentradas 2-10 Ecuaciones de movimiento 2-13 Principio de Hamilton 2-13 Principio de los trabajos virtuales 2-13 Principio de D’Alembert 2-13 Formulación de la ecuación de movimiento para un sistema de 1 GLD 2-13 Formulación de las ecuaciones de movimiento para modelos con múltiples GLD 2-15 Capítulo 3 - RESPUESTA DE UN OSCILADOR SIMPLE 3-1 Ecuación de movimiento y equilibrio dinámico 3-1 Fórmula de Geiger 3-2 Características dinámicas con amortiguamiento 3-3 Excitación periódica 3-7 Excitación armónica 3-8 Excitación arbitraria. Integral de Duhamel 3-10 Respuesta a un impulso elemental 3-10 Factor de amplificación dinámica 3-12 Espectros sísmicos de respuesta 3-12 Integración numérica de la ecuación de movimiento 3-16 Método de Newmark (1 GLD) 3-16 Desarrollo y forma operativa 3-18 Introducción - Bibliografía 3 Capítulo 4 - RESPUESTA DINÁMICA DE UNA ESTRUCTURA CON MÚLTI-
PLES GRADOS DE LIBERTAD 4-1 Ecuaciones de movimiento y equilibrio dinámico 4-1 Vibraciones libres 4-1 Características dinámicas 4-1 Normalización de los modos 4-3 Obtención de los grados de libertad dinámicos 4-4 Condensación estática de la matriz de rigidez 4-4 Matriz de amortiguamiento 4-6 Matrices de amortiguamiento ortogonales 4-8 Determinación práctica de modos y frecuencias 4-9 Método de Stodola-Vianello 4-9 Resolución de las ecuaciones de movimiento en estructuras con múltiples GLD 4-12 Descomposición y superposición modal 4-12 Integración directa de las ecuaciones de movimiento 4-14 Respuesta máxima utilizando espectros de respuesta 4-15 Capítulo 5 - ANÁLISIS DE CONSTRUCCIONES CON EFECTOS DINÁMI-
COS DE VIENTO 5-1 Acciones paralelas a la dirección del viento 5-1 Acciones perpendiculares a la dirección del viento 5-3 Criterios de confort en edificios que oscilan 5-5 Capítulo 6 - ANÁLISIS DE CONSTRUCCIONES CON EFECTOS SISMICOS 6-1 Definición numérica de la acción sísmica 6-1 Definición mediante espectros de respuesta 6-2 Definición mediante acelerogramas 6-4 Métodos de análisis según INPRES – CIRSOC 103 6-6 Procedimientos con fuerzas estáticas equivalentes 6-6 Métodos dinámicos 6-6 Métodos dinámicos INPRES – CIRSOC 103 6-6 Introducción - Bibliografía 4 BIBLIOGRAFÍA: [1] Estructuras Sometidas a Acciones Sísmicas. Cálculo por Ordenador – A.H. Barbat, J.M. Canet – 2da. Edición, CIMNE, Barcelona, 1994. [2] Diseño Sismorresistente de Edificios – L.M. Bozzo, A. H. Barbat, Editorial Reverté, Barcelona, 2000. [3] Dinámica Estructural – J. Massa, C. Prato, Publicación del Departamento de Estructuras de la Facultad de Cs. Exactas Físicas y Naturales de la Universidad Nacional de Córdoba, 1986. [4] Finite Element Procedures – K. J. Bathe, Prentice-Hall, 1996. [5] Finite Element Modeling in Engineering Practice – C. Spyrakos, Algor Publishing Division, Pitsburg, 1996. [6] Linear and Nonlinear Finite Element Analysis in Engineering Practice – C. Spyrakos, J. Raftoyiannis, Algor Publishing Division, Pitsburg, 1997. [7] Diseño Sísmico de Edificios – E. Bazán, R. Meli, LIMUSA, 2001. Capítulo 1 – Nociones Básicas de Sismología 1-1 CAPÍTULO 1 - NOCIONES BASICAS DE SISMOLOGIA Causas que generan los terremotos o sismos Los terremotos pueden definirse como movimientos de la corteza terrestre, con amplitudes y frecuencias dependientes del tiempo. Las causas que los generan son variadas: Terremotos de colapso: son los originados en cavidades subterráneas por el colapso de las mismas, son de baja intensidad. Terremotos de origen volcánico: la explosión de gases durante las erupciones volcánicas puede producir terremotos que, en general, tienen una intensidad pequeña y afectan a su-
perficies limitadas. Terremotos tectónicos: están causados por la rotura brusca de las capas rocosas a lo largo de superficies de fractura (fallas), son los más fuertes y más frecuentes. Terremotos causados por explosiones: las explosiones producidas por el hombre son capa-
ces de generar vibraciones del terreno, con una intensidad tal que pueda causar movimien-
tos en las estructuras. En general, el movimiento de la corteza se produce por un choque o movimiento brusco ocurrido a una cierta profundidad bajo la superficie terrestre en un punto teórico denomi-
nado foco o hipocentro, a su proyección sobre la superficie terrestre se le denomina epi-
centro. estructura
Fig. 1.1 – Definiciones geométricas de un sismo Capítulo 1 – Nociones Básicas de Sismología 1-2 Los sismos desde el punto de vista de la ingeniería y su caracteri-
zación Los terremotos más importantes son los tectónicos, pues son los que traen consecuencias más desastrosas en las estructuras que afectan, debido a esto, son los que se tienen en cuenta para la elaboración de normas para la contracción de estructuras sismoresistentes. La intensidad sísmica es una medida de los efectos de los terremotos en el entorno y en particular sobre las estructuras. Existen diferentes escalas de intensidades que describen, para cada valor que esta tome, los efectos que produce el terremoto. Una de las más difundidas es la escala de Mercalli Modificada. Algunos de los efectos sobre las estructuras en orden creciente de intensidad son: 1. fisuración de las estructuras de madera 2. agrietamiento de las estructuras débiles de mampostería 3. agrietamiento de las estructuras ordinarias de mampostería 4. colapso parcial de estructuras ordinarias de mampostería; daño en estructuras bien ejecutadas de mampostería no diseñadas para resistir fuerzas sísmicas 5. colapso de estructuras ordinarias de mampostería; las estructuras con diseño anti-
sísmico son seriamente dañadas; daños en cimientos; grietas en el terreno la mayoría de las estructuras son destruidas junto con sus cimientos, daños importantes en presas y diques, grandes deslizamientos del terreno destrucción casi total, grandes masas de rocas desplazadas, etc. Un sismo se caracteriza por su intensidad (parámetro subjetivo) y por su magnitud (pa-
rámetro objetivo). La escala objetiva más popular es la de Ritcher, en la que la magnitud M mide la energía del terremoto en el foco y es el logaritmo decimal de la amplitud del movimiento sísmico, medido en micrones a 100[km] del epicentro, por un sismógrafo Wood-Anderson están-
dar. La magnitud M está relacionada con la energía del terremoto, en ergios, por la expre-
sión: M E 5 , 1 8 , 11 log + = Se han establecido varias relaciones empíricas entre la intensidad I
y la magnitud M, enumeramos algunas a continuación: Esteva y Rosenblueth: R M I
log 46 , 2 45 , 1 16 , 8 − + = R: distancia focal en [km] También se ha relacionado la magnitud M con los valores máximos de las características cinemáticas del movimiento, estas relaciones se han establecido estadísticamente. Capítulo 1 – Nociones Básicas de Sismología 1-3 Donovan: ( )
: aceleración máxima del terreno en [cm/s
] R: distancia focal en [km] Esteva y Villaverde: ( )
: aceleración máxima en [cm/s
: velocidad máxima en [cm/s] R: distancia focal en [km] Esteva: ( )
2 8 , 0
+ = R e a
17 , 0 15
e R e v a
: velocidad máxima en [cm/s] R: distancia focal en [km] Registro de ondas sísmicas. Parámetros utilizados y mapas de riesgo sísmico Los terremotos son fenómenos debidos a la brusca liberación de la energía de deformación acumulada durante largos periodos de tiempo en la zona superficial de la tierra. Los sis-
mos producen ondas de varios tipos, que se propagan desde su foco en todas las direccio-
nes a través de la tierra. Estas ondas son registradas mediante aparatos denominados sis-
mógrafos, diseñados para medir la aceleración, la velocidad o el desplazamiento del mo-
vimiento sísmico. Estos parámetros son relativos, ya que los valores obtenidos están afec-
tados por las características del instrumento registrador y por las condiciones de ruido am-
biental en el lugar de registro. Los mapas de riesgo sísmico representan una síntesis de todos los datos sismológicos y geológicos de un país. Estos mapas se utilizan para determinar el nivel de protección que se debe alcanzar en las estructuras en cada zona de riesgo. Diversos aspectos brindan la subdivisión en zonas, pero los fundamentales son: Capítulo 1 – Nociones Básicas de Sismología 1-4 Estudios geológicos y geotécnicos: proporcionan datos de composición y características dinámicas de las rocas y capas de suelo que componen la corteza terrestre. Estudios sismológicos: sintetizan los parámetros que caracterizan la sismicidad de la zona: 1. ubicación de fallas 2. registro de los terremotos que ocurren en la zona 3. mapas de epicentros 4. datos históricos 5. periodos de retorno (intervalo medio de tiempo en que se espera ocurran dos sismos de igual o mayor intensidad) 6. datos del mecanismo focal 7. correlación de la sismicidad de la zona analizada con la de la macrozona en la que se encuentra Estudios de Ingeniería y Sismología: 1. análisis del efecto que han producido sobre las estructuras y las personas los terremotos ocurridos en el pasado 2. “predicción” estadística de las características más probables de la acción sísmica que se produzca en la zona Es importante destacar que la geología local de la zona puede modificar la propagación de las ondas sísmicas. Las ondas se reflejan y se refractan cuando en su recorrido aparece una discontinuidad, por ejemplo una variación de las características mecánicas del terreno, ello produce cambios en la velocidad. En general, el cálculo y la cuantificación de las acciones sísmicas en la estructuras se rea-
liza en función de protocolos, secuencias y definiciones de acciones dadas por normas y reglamentos. En los capítulos siguientes se ofrecen algunas aplicaciones de este tipo. Capítulo 2 – Conceptos Básicos de Dinámica Estructural 2-1 CAPÍTULO 2 - CONCEPTOS BASICOS DE DINAMICA ES-
TRUCTURAL En un sentido amplio, un sistema dinámico es aquel cuyas variables experimentan varia-
ciones en el tiempo y, si se conocen las influencias externas que actúan sobre el sistema, podrá predecirse el comportamiento de este. influencias externas
conocidas estas acciones
externas, permiten
"predecir" el comportamiento
de las variables temporales
variables con variaciones
En nuestro curso, los sistemas a estudiar serán sistemas estructurales, las variaciones en el tiempo serán vibraciones producidas por cargas dinámicas. permiten evaluar el
la estructura frente
a acciones dinámicas
que gobiernan el
Definición de la acción dinámica Una acción tiene carácter dinámico cuando su variación con el tiempo es rápida y da ori-
gen a fuerzas de inercia comparables en magnitud con las fuerzas estáticas. Algunas fuen-
tes importantes de vibraciones estructurales son: Capítulo 2 – Conceptos Básicos de Dinámica Estructural 2-2 - sismos - viento - olas y corrientes de agua - explosiones e impactos - cargas móviles (vehículos, personas, etc.) La definición de estas cargas externas puede distinguirse entre: determinista y no determi-
nista, ésta última denominada también estocástica o aleatoria. determinista: cuando su variación temporal es perfectamente conocida no determinista: cuando alguno o todos sus parámetros son definidos estadísticamente En nuestro curso trabajaremos con cargas definidas en forma DETERMINISTA. Respuesta dinámica → cualquier magnitud que pueda caracterizar el efecto de una carga dinámica sobre la estructura Una carga definida determinísticamente da origen a una respuesta, también determinista. Fig. 2.1- Definición de la respuesta dinámica: para un punto considerado se calculan: deformaciones, aceleraciones, tensiones, etc. Acciones y fuerzas dinámicas Las acciones dinámicas definidas utilizando representaciones deterministas, son funciones del tiempo cuyo valor en cada instante ES CONOCIDO. Este tipo de representación es apropiado para evaluar el comportamiento de una estructura A POSTERIORI del acontecimiento que dio lugar a dicha acción. Por ejemplo, evaluar el comportamiento de un edificio nuevo ante el terremoto ocurrido en México en 1986 (del que se poseen registros). El diseño de una estructura NO PUEDE encararse en base a ac-
ciones deterministas, pues nada nos asegura que la acción estudiada volverá a repetirse. Capítulo 2 – Conceptos Básicos de Dinámica Estructural 2-3 F(t)
M(t2)
M(t3)
F(t3)
este esquema temporal
de carga debe ser
perfectamente conocido
Fig. 2.2 - Acción y respuesta determinista No considerada como
"carga" sino como una
propiedad intrínseca de
da origen a fuerzas
de "inercia" comparables
con las estáticas
la comparación se
basa en: - velocidad de
- periodo propio
Fig. 2.3 - Acción dinámica y propiedades de la estructura Importancia de la masa en el problema dinámico Aunque la carga varíe con el tiempo, la respuesta de una estructura varía radicalmente se-
gún la masa que vibra con ella. Ante una misma función de carga, una estructura SIN MASA y una CON MASA responden de la siguiente manera: Capítulo 2 – Conceptos Básicos de Dinámica Estructural 2-4 b) Estructura CON MASA ⇒ CON INERCIA
a) Estructura SIN MASA ⇒ SIN INERCIA
Se desarrolla energía
cinética, que modifica
la respuesta y deja vibraciones
rigidez: K
F(t) m/2
masa: m = 0
K x(t) = F(t)
La respuesta seguirá
exactamente la forma
t2 > t1 !!
x0 = F0/K
( ) ( ) ( ) t F t x k t x m = + & &
Fig. 2.4 - Importancia de la masa en la respuesta Capítulo 2 – Conceptos Básicos de Dinámica Estructural 2-5 Velocidad de reacción de una estructura Ante una acción exterior, distintas estructuras reaccionarán de formas diferentes. Esta res-
puesta está íntimamente relacionada con las formas o modos de vibrar y sus correspon-
dientes frecuencias o periodos propios. En el caso de un oscilador de 1 grado de libertad, este periodo propio se obtiene fácilmente. No así para estructuras de múltiples GLD. Como veremos en los capítulos siguientes, los periodos y formas de vibrar dependen de las características geométricas y de materiales (rigidez) y de la inercia que la estructura opone al movimiento (masa). En general si tD >> T ⇒ no es necesario un análisis dinámico
si tD ≅ T ⇒ PROBLEMA DINAMICO
el periodo propio permanece
prácticamente constante
ESTRUCTURA SIN
la amplitud decrece en cada ciclo
Fig. 2.5 - Velocidad de reacción T vs. t
Capítulo 2 – Conceptos Básicos de Dinámica Estructural 2-6 k
a masa constante
k1 > k2
a rigidez constante
; k ; m
T T →
Fig. 2.6 - Velocidad de reacción; varios T
Modelos dinámicos característicos Desde el punto de vista del cálculo numérico, obtener la respuesta dinámica de una estruc-
tura, es el resultado de "filtrar" la señal de excitación a través de la misma estructura y ob-
tener las variaciones de las magnitudes de análisis (desplazamientos, velocidades, acelera-
ciones, momentos, tensiones, etc.) respecto del tiempo. La obtención de la respuesta requiere, previamente, la definición del movimiento del te-
rreno (en caso sísmico) tanto como de las características estructurales del mismo y de la estructura propiamente dicha. El análisis es practicado, no a la propia estructura sino a un modelo mecánico de la misma. La definición del modelo depende del tipo de estructura analizado y pretende brindar una serie de relaciones entre acciones y respuesta que descri-
ban un modelo matemático del problema. Este modelo matemático puede ser resuelto mediante diversas técnicas. En nuestro caso haremos hincapié en los métodos numéricos de análisis. Según la certeza con que fueron formulados los modelos y procedimientos o algoritmos de cálculo durante el análisis, será la precisión de la respuesta obtenida. Capítulo 2 – Conceptos Básicos de Dinámica Estructural 2-7 a3
respuesta en un punto
Fig. 2.7 - “Filtrado” de una señal sísmica Se brindan, a continuación, algunas definiciones típicas del análisis estructural dinámico de una estructura: Grados de libertad (GL) Se definen como grados de libertad (GL) a los puntos de la estructura en los cuales se identifica algún desplazamiento y permiten brindar una deformada de la estructura. Grados de libertad dinámicos (GLD) Son los grados de libertad que tienen asociada masa y para los cuales puede conocerse las vibraciones o movimientos a lo largo del tiempo. Capítulo 2 – Conceptos Básicos de Dinámica Estructural 2-8 menos exacta
RESPUESTA (desplaz. en el piso 1)
mas exacta
[ ] [ ] [ ]a M x K x M − = + & &
Fig. 2.8- Modelización de una estructura Métodos de modelización dinámica Pueden distinguirse modelos dinámicos exactos y modelos dinámicos discretos. En general, para la primera clase, solo pueden resolverse casos muy sencillos y con poca aplicación practica, por lo que a lo largo del curso profundizaremos en modelos discretos. Para estos métodos modelos discretos, se debe tener en cuenta que la subdivisión en domi-
nios finitos es tanto espacial (discretización estructural) como temporal (solución para ins-
tantes de tiempo determinados). Capítulo 2 – Conceptos Básicos de Dinámica Estructural 2-9 nº FINITO DE INSTANTES
se calcula la respuesta
(para cada GLD)
nº FINITO DE PUNTOS
ESPACIALES (GLD)
en CADA PUNTO nº infinito de
DE LA ESTRUCTURA puntos espaciales
en CADA INSTANTE nº imfinito de
DE TIEMPO puntos temporales
Fig. 2.9 - Modelos dinámicos Discretización espacial de las estructuras Fundamentalmente, la diferencia con lo visto en otros cursos de análisis estructural (estáti-
co) radica en que en dinámica estructural, cuando hablamos de discretizar espacialmente, nos referimos a los GLD. Un modelo dinámico exacto (con infinitos GLD) acarrearía más inconvenientes en la reso-
lución matemática que beneficios en su precisión. Además, en estructuras de edificios y en la mayoría de las estructuras civiles, las masas se encuentran más o menos concentradas en lugares conocidos. Es por esto que nuestro principal método de modelización dinámica será el de las MASAS CONCENTRADAS. No obstante, existen otros, como ser: - método de los DESPLAZAMIENTOS GENERALIZADOS - método de los ELEMENTOS FINITOS Capítulo 2 – Conceptos Básicos de Dinámica Estructural 2-10 Método de las masas concentradas nº total de componentes
de desplazamiento según
los cuales vibran las masas
específicas pueden conocerse
mediante los procedimientos
del análisis estático
se calcula la "deformada"
del modelo en cada instante
- MODELOS DE 1 GLD
MULTIPLES GLD
nº de GLD
Modelos con 1 GLD: k
Fig. 2.10 - Modelos con un solo grado de libertad. (a) modelo conservativo; (b) modelo con amortiguamiento; (c) modelo sísmico. x
Fig. 2.11 - Estructuras modelizadas como un sistema de un solo grado de libertad. (a) pórtico; (b) el mismo pórtico con la masa concentrada al nivel de la viga; (c) modelo di-
námico. Capítulo 2 – Conceptos Básicos de Dinámica Estructural 2-11 Modelos con múltiples GLD: xr
Fig. 2.12 -Modelos con varios grados de libertad. (a) modelo conservativo; (b) modelo con amortiguamiento; (c) modelo sísmico. Fig. 2.13 - Estructura con dos grados de libertad: Pórtico de dos pisos y su modelo diná-
mico. a(t)
Fig. 2.14 - Estructura con masa distribuida (antena) y su modelo dinámico discreto con n grados de libertad. Capítulo 2 – Conceptos Básicos de Dinámica Estructural 2-12 Fig. 2.15 - Modelo dinámico de un pórtico de cortante y pórtico espacial modelizado co-
mo un sistema completo (10 grados de libertad) y simplificado ( dos grados de libertad). Fig. 2.16 - Modelo dinámico con grados de libertas de rotación. Capítulo 2 – Conceptos Básicos de Dinámica Estructural 2-13 Ecuaciones de movimiento Las ecuaciones de movimiento son las expresiones matemáticas que gobiernan la respues-
ta dinámica de las estructuras. Pueden obtenerse a partir de cualquiera de los principios de la mecánica clásica: Principio de Hamilton ( ) 0 ;
d dt E dt E E π π (2.1) La primera expresión se denomina funcional de Hamilton; donde E
es la energía poten-
cial, E
es la energía cinética y E
la correspondiente a fuerzas no conservativas. La segun-
da expresión permite establecer el equilibrio a través de una variación funcional nula. Principio de los trabajos virtuales Se trabaja en forma similar a lo visto en análisis estático pero incluyendo las fuerzas de inercia y disipativas. e i
w w δ δ = (2.2) Principio de D’Alembert Proporciona el método más directo para obtener las ecuaciones de movimiento de un sis-
tema dinámico. Puede formularse como sigue: “un sistema dinámico esta en equilibrio cuando todas las fuerzas que actúan en el mismo, incluidas las de inercia y disipativas, cumplen las ecua-
ciones de equilibrio estático en cada instante de tiempo. Formulación de la ecuación de movimiento para un sistema de 1 GLD Tomando el sistema de la figura 2-10, podemos distinguir dos casos: Fig. 2.17- (a) fuerza aplicada; (b) modelo sísmico Para el modelo (a), aplicando el principio de D’Alembert, tendríamos: Capítulo 2 – Conceptos Básicos de Dinámica Estructural 2-14 F(t) Fi(t) Fa(t) Fe(t)
F(t) = Fi(t) + Fa(t) + Fe(t) equilibrio
( ) t x& &
( ) t x
inercia amort. elástica
Fig. 2.18 - Equilibrio de fuerzas para 1 GLD Al aplicar una fuerza exterior F(t), se genera aceleración, velocidad y desplazamiento para un cierto instante “t”; a causa de esto se producen fuerzas: i- de inercia ( ) ( ) t x m t F
= (2.3) ii- de amortiguamiento ( ) ( ) t x c t F
= (2.4) iii- elásticas ( ) ( ) t x k t F
= (2.5) equilibrio en el instante “t” ( ) ( ) ( ) ( ) 0 = − − − t F t F t F t F
(2.6) ( ) ( ) ( ) ( ) t F t F t F t F
= + + (2.7) F x k x c x m = + +
(2.8) Se omite (por simplicidad de notación) la dependencia del tiempo, pero de aquí en adelan-
te ésta se encontrará implícita en toda variable temporal. ( ) ( ) ( )
( ) ( ) . ; ;
; ; : ej
etc F t F a t a
Capítulo 2 – Conceptos Básicos de Dinámica Estructural 2-15 La ecuación (2.8) es la de movimiento correspondiente a 1 GLD con carga exterior y amortiguamiento. Para el modelo (b) de la figura 2-17, el planteo es similar, solo que no tiene fuerza exterior aplicada y la fuerza de inercia se ve afectada por la aceleración total de la masa: [ ] x a m F
+ = (2.9) entonces, la ecuación de movimiento queda: [ ] 0 = − − + − x k x c x a m
(2.10) a m x k x c x m − = + +
(2.11) La (2.11) es la ecuación de movimiento para 1 GLD con aceleración de apoyo (sísmico) y amortiguamiento. Un caso general sería la inclusión de aceleración de apoyo y fuerza exte-
rior: a m F x k x c x m − = + +
(2.12) Formulación de las ecuaciones de movimiento para modelos con múltiples GLD El modelo con varios grados de libertad más sencillo es el de edificios de cortante (fig. 2-
12). Está basado en que las plantas son infinitamente rígidas y en que los únicos movi-
mientos posibles de éstas son los desplazamientos horizontales. Aplicando el principio de D’Alembert en cada GLD (uno por piso) se obtiene: 0 = − − −
er ar ir r
F F F F (2.13) para el piso (r). Planteando el equilibrio para todos los GLD, nos queda un sistema de ecuaciones (vecto-
rial) 0 = − − −
F F F F (2.14) Capítulo 2 – Conceptos Básicos de Dinámica Estructural 2-16 #Nota: en adelante, los vectores y matrices serán representados(en general) con minús-
culas y mayúsculas en negrita, respectivamente. [ ]
f f f ,... ,
= F vector de fuerzas externas { } a
J x M F + − = & & vector de fuerzas de inercia x C F & =
vector de fuerzas disipativas x K F =
vector de fuerzas elásticas entonces, el sistema (2.14) puede escribirse: { } a J M F x K x C x M − = + + & & & (2.15) M → matriz de masas C → matriz de amortiguamiento K → matriz de rigidez J
= [1,1,…1] → vector con todos sus elementos igual a uno Si bien la (2.15) fue deducida para un modelo de edificio cortante, es una expresión AB-
SOLUTAMENTE GENERAL, inclusive para modelos de elementos finitos, y solo varían las formas de M, C y K. Para un modelo sísmico, la (2.15) se reduce a: { } a J M x K x C x M − = + + & & & (2.16) Y en caso general de pórticos 3D, o modelos de elementos finitos, suele sustituirse x por D para indicar que cada GLD puede sufrir desplazamientos en 3 direcciones y respectivos giros. { } a J M F D K D C D M − = + +
(2.17) Nota: en este último caso, J llevará unos en las componentes a las cuales se quiera aplicar el acelerograma a, por ejemplo componente x: J = [1,0,0,0,0,0,1,0,0,0,0,0,…..,1,0,0,0,0,0,…..,1,0,0,0,0,0] etc. Capítulo 3 – Respuesta de un Oscilador Simple 3-1 CAPÍTULO 3 - RESPUESTA DE UN OSCILADOR SIMPLE Ecuación de movimiento y equilibrio dinámico Las características dinámicas de un oscilador de 1 GLD pueden estudiarse mediante un modelo no amortiguado con vibraciones libres, cuya ecuación de movimiento es 0 = + x k x m
(3.1) Fig. 3.1 - Modelo de 1 GLD, no amortiguado. La vibración, del modelo de la fig. 3.1, es inducida por algunas condiciones iniciales, sean desplazamiento, velocidad o aceleración en el instante t = 0. Luego, durante las vibraciones no recibe ningún tipo de perturbación. Dividiendo (3.1) por m y usando la notación: m
ω ; m
= ω (3.2) se obtiene: 0
= + x x ω
(3.3) ω es la “pulsación” o “frecuencia circular” o simplemente frecuencia del modelo estudia-
do. Viene expresada en radianes por segundo (1/s). La frecuencia cíclica viene dada por π
= f (3.4) y se expresa en ciclos por segundo o hertz. Capítulo 3 – Respuesta de un Oscilador Simple 3-2 Finalmente, otra característica es el período natural f
= (3.5) ω
= T (3.6) La solución general de la (3.1) o (3.3) puede escribirse: ( ) Ψ + = t A x ω sen (3.7) donde A es la amplitud del movimiento y ψ el ángulo de fase. Los valores de A y ψ se cal-
culan a partir de las condiciones iniciales del problema, por ejemplo para ( )
0 x x = ; ( )
= Ψ Fórmula de Geiger Sustituyendo g
m = G: peso de m g: aceleración de la gra-
vedad k G
⋅ = = ω SG
= : desplazamiento estático producido por el peso G en la dirección del grado de libertad entonces: SG
⋅ = ω (3.9) Capítulo 3 – Respuesta de un Oscilador Simple 3-3 SG
T ⋅ =
(3.10) Utilizando unidades de S.I. la (3.10) queda: SG
X T 00 , 2 = (3.11) con X
expresado en metros para un peso G en Newton. Características dinámicas con amortiguamiento El amortiguamiento puede definirse estudiando las vibraciones libres del modelo de la fi-
gura 3-2: Fig. 3.2 - Modelo de 1 GLD con amortiguamiento (vibraciones libres) Si se toma la ecuación (2.12) sin cargas ni aceleraciones de apoyo (vibraciones libres) y se divide por m se obtiene: 0 2
= + + x x x ω β & & & (3.12) m
= β 2 (3.13) la solución de (3.12) está dada en la forma: rt
e x = (3.14) que proporciona la ecuación característica: Capítulo 3 – Respuesta de un Oscilador Simple 3-4 0 2
= + + ω β r r (3.15) Ya que (3.12) es una ecuación diferencial de segundo orden, lineal, homogénea a coefi-
cientes constantes. Las soluciones de (3.15) son 0
= − ± − = ω β β r (3.16) Según sea el radicando β
se encuentran tres tipos de amortiguamiento: β
> 0 → SUPERCRITICO: la estructura NO VIBRA β
= 0 → CRÍTICO: caso límite β = ω → m
c ω 2 − = β
< 0 → SUBCRITICO: la estructura VIBRA con amplitud decreciente Este es el caso más frecuente en ingeniería civil, por lo que enfatizaremos su estudio. Para este caso (subcrítico), la cantidad (β
) es negativa, lo que hace que (3.16) tenga raíces complejas: 0 1
= − ± − = ξ ω β i r (3.17) con 1 − = i Llamando frecuencia de vibración amortiguada a: 2
1 ξ ω ω − =
(3.18) se obtiene: v
i r ω β ± − =
(3.19) v
i r ω ξω ± − =
(3.20) En las ecuaciones anteriores aparece la magnitud m
= = (3.21) Capítulo 3 – Respuesta de un Oscilador Simple 3-5 conocida como fracción de amortiguamiento crítico (en estructuras corrientes 0.02<ξ<0.06), también: ω
ξ = o ω
= Volviendo a la resolución de la ecuación de movimiento (3.12), escribimos la solución ge-
neral en la forma: t r t r
+ = (3.22) sustituyendo r
por la expresión 3.20, se obtiene ( ) Ψ + =
t e A x
sen (3.23) Las constantes A y ψ se obtienen de las condiciones iniciales del problema. Fig. 3.3 - Vibraciones libres amortiguadas / 1 GLD La evaluación del amortiguamiento en una estructura es un problema esencial en la diná-
mica estructural. El origen de las fuerzas de amortiguamiento se debe a diferentes causas: - Rozamiento entre superficies de deslizamiento, en donde la fuerza de amortiguamiento es proporcional a la fuerza normal y al coeficiente de rozamiento (hipótesis de Cou-
lomb) - Amortiguamiento debido a fuerzas aero o hidrodinámicas - Debido a fricción interna del material de la estructura Generalmente, en el cálculo dinámico de estructuras, se utiliza un modelo de gran simpli-
cidad que caracteriza el amortiguamiento de toda la estructura. Este modelo denominado de amortiguamiento viscoso se debe a Kelvin-Voigt y es proporcional a la velocidad. x c F
& − = Capítulo 3 – Respuesta de un Oscilador Simple 3-6 Los reglamentos brindan los coeficientes de amortiguamiento para cada tipo de estructura, pero puede obtenerse en forma experimental y con un método relativamente simple: Determinación práctica de ξ Para amortiguamientos bajos (del orden del 10% de crítico) la relación entre dos picos su-
cesivos puede aproximarse: ( )
( ) Ψ +
t sen e A
(3.24) pero v n n
; con v v
T ω π 2 = ∴ v
Ψ + +
t sen e
(3.25) ( )
( ) Ψ + +
t sen e e
ω ω π ξ 2
(3.26) tomando logaritmo natural: ( )
máx ω ω π ξ 2
(3.27) # notar que para amortiguamientos del orden de Capítulo 3 – Respuesta de un Oscilador Simple 3-7 2
1 1 , 0 ξ ω ω ξ − = → =
ω ω ω ω ≈ → ⋅ =
995 , 0 Entonces ( )
(3.28) Para el caso de lecturas separadas por N ciclos: ( ) ( ) [ ]
N n x n x
= (3.29) Excitación periódica En la figura 3-4 pueden observarse diversas funciones de carga. De éstas, nos interesan por ahora, las periódicas y más particularmente las excitaciones armónicas ya que mediante series de Fourier cualquier excitación periódica puede llevarse a una suma de armónicas simples. Fig. 3.4 - Tipos de cargas dinámicas. (a) armónica; (b) periódicas; (c) cuasi periódicas; (d), (e) fuerzas impulsivas; (f) carga dinámica general; (g) aceleración sísmica del terre-
no. Capítulo 3 – Respuesta de un Oscilador Simple 3-8 Excitación armónica Si la carga es de tipo: ( ) t P P Ω = sen
Entonces, la ecuación de movimiento será: ( ) t sen P x k x c x m Ω = + +
& & & (3.30) Ω: frecuencia de la excitación la (3.30) puede escribirse también: ( ) t sen
x x x Ω = + +
2 ω ξω & & & (3.31) La solución general de esta ecuación viene dada por p h g
x x x + = (3.32) t r t r
+ = (3.22) Es la solución de la ecuación diferencial homogénea. t i t t i t
e e c e e c x
ω ξω ω ξω − − −
(3.33) Utilizando matemática para números complejos, esta última ecuación puede escribirse: ( ) t c t sen c e x
(3.34) x
en la (3.32) es la solución particular y lleva la forma t B t A x
Ω + Ω = cos sen (3.35) Capítulo 3 – Respuesta de un Oscilador Simple 3-9 Derivando y reemplazando en (3.30) se obtienen las constantes A y B: Denotando ω
= (3.36) ( ) ( )
ξγ γ
A (3.37) ( ) ( )
B (3.38) Condiciones iniciales Para el caso en que ( )
0 x x & & = Pueden calcularse las partes correspondientes a la solución particular (para t = 0) ( )
(3.39) ( )
& (3.40) Basados en éstas y en (3.32) y (3.34) podemos plantear las siguientes ecuaciones: ( )
x c x x + = = (3.41) ( )
p v p g
x c x x x x & & & + + + − = = ω ξω ξω (3.42) y despejarse c
’ y c
’: [ ]
x x x x c ξω ξω
− − + = & & (3.43) 0 0 2
x x c − = (3.44) Capítulo 3 – Respuesta de un Oscilador Simple 3-10 Forma de operar: 1- Se calculan 0 p
x y 0 p
x& con las 3.39 y 3.40 2- Se computan c
’ con las 3.43 y 3.44 3- Se reemplaza todo en 3.32 y se tiene la respuesta en desplazamientos 4- Para obtener velocidades y aceleraciones se deriva la 3.32 Nota: la parte de la solución correspondiente a la ecuación homogénea incluye el coefi-
ciente e
-ξωt
que es una función decreciente con el tiempo, por lo que se la denomina transi-
toria y desaparece antes o después según sea el valor de ξ. La estructura sigue vibrando con una frecuencia prácticamente igual a la de la excitación Ω. Esta parte de la solución se denomina respuesta en régimen. Excitación arbitraria. Integral de Duhamel Respuesta a un impulso elemental Dada una carga impulsiva como la de la fig. 3-4(d), para condiciones iniciales nulas, pue-
de calcularse la velocidad y posición al finalizar el impulso y a partir de ese instante se trata de un problema de vibraciones libres. P(τ)
α tg =
Fig. 3.5 - Respuesta a un impulso τ d x x x
& & & & + = ; ( ) m P x τ = & & entonces: ( )
= & (a) Capítulo 3 – Respuesta de un Oscilador Simple 3-11 por otra parte ( )
τ τ d x d x x x
& & & + + = (b) ( )
τ d x x
& & = (b) Si dτ es muy pequeño, entonces (b) es despreciable frente a (a), es decir, luego de actuar un impulso la masa queda con velocidad pero prácticamente no se desplaza. La velocidad dada por (a) es la condición inicial para las vibraciones libres siguientes. Reemplazando ésta en la (3.23), para un tiempo infinitesimal, se obtiene: ( )
τ ξω
sen (3.45) Si el sistema es lineal, pueden superponerse las respuestas de una sucesión de impulsos infinitesimales hasta el tiempo genérico “t”. La solución dada por esta superposición es conocida como Integral de Duhamel ( ) ( )
( ) [ ] τ τ ω τ
(3.46) En el caso de excitación sísmica, se transforma en: ( ) ( )
d t sen e a t x
(3.47) Para el caso en que el problema tenga condiciones iniciales 0
≠ x& y 0
≠ x , la solución se obtiene con el mismo procedimiento dado para las cargas armónicas, solo que la solución particular vendrá dada por (3.46) o (3.47). Soluciones a la integral de Duhamel: (a) Para ciertas funciones de P(τ) se encuentran tabuladas las primitivas de (3.46) y (3.47), así como el máximo factor dinámico (FAD) y el instante de máxima respuesta. Capítulo 3 – Respuesta de un Oscilador Simple 3-12 (b) Para funciones aproximadas por “trozos rectos”, también se pueden encontrar solucio-
nes “exactas” a dicha integral (punto 5.3.4 / pág. 154 / ref. [1]) (c) En el caso mas general se puede “integrar numéricamente” utilizando alguno de los métodos conocidos (Trapecios, Simpson, etc.) Factor de amplificación dinámica Es el coeficiente entre el máximo desplazamiento dinámico y el que produciría la carga en forma estática: ( )
= (3.48) Espectros sísmicos de respuesta Puede definirse el espectro d respuesta (para un acelerograma dado) como los máximos valores de la respuesta de un sistema, expresado en función del periodo propio. Esto es muy útil para el diseño de estructuras donde solo los máximos son necesarios. Para un sistema de 1 GLD: ( ) t a m x k x c x m − = + + & & & (2.11) La solución a esta ecuación viene dada por la ec. (3.47) y el máximo valor absoluto será S
. Reemplazando ( )
= Puede obtenerse, para 1 GLD, la expresión exacta del espectro de desplazamientos. Deri-
vando una y dos veces se obtienen las expresiones para velocidad y aceleración: ( )
t x S & = ; ( )
t x S & & = (Ver ref. [1], cap. 5) Capítulo 3 – Respuesta de un Oscilador Simple 3-13 Fig. 3.6 - Espectro sísmico de respuesta para desplazamientos. S
No obstante pueden hacerse ciertas hipótesis que simplifican las expresiones obtenidas. Estas nuevas ecuaciones definen los “Seudo espectros de respuesta”. Hipótesis simplificativas: i- en estructuras civiles, 0,002 < ξ < 0,2, por lo que puede reemplazarse ω
por ω ii- la función coseno que aparece en el espectro de velocidades (al derivar 3.47) puede reemplazarse, sin que ello implique grandes variaciones en los valores máximos, por una tipo seno entonces: ( ) ( )
d t sen e a s τ τ ω τ
, (3.49) ( ) ( )
d t sen e a s τ τ ω τ ξ ω
, (3.50) Capítulo 3 – Respuesta de un Oscilador Simple 3-14 ( ) ( )
d t sen e a s τ τ ω τ ω ξ ω
, (3.51) definen los seudo espectros para desplazamientos, velocidades y aceleraciones. Es importante notar que: d v
S s ω = (3.52) d a
ω = (3.53) Por lo que un espectro completo puede representarse en un gráfico tipo logarítmico como el de la figura 3-7. Fig. 3.7 - Comparación entre los espectros suavizados de Newmark y Hall, los espectros de diseño del UBC y los espectros elásticos para los registros de EL Centro y de la presa de Pacoima. La normativa argentina INPRES-CIRSOC 103 establece la acción sísmica como una fun-
ción del seudo espectro de respuesta de aceleraciones S
. Es por eso que solo aparecen grá-
ficos del tipo de la fig. 3-9, que varían según el tipo de suelo y la zona sísmica considera-
da. Capítulo 3 – Respuesta de un Oscilador Simple 3-15 Fig. 3.8 - Espectro de diseño suavizado de Seed e Idriss (1982). Fig. 3.9 - Seudo espectro y espectro de aceleraciones; (a) corresponde al acelerograma (1). (b) corresponde al acelerograma (2). Capítulo 3 – Respuesta de un Oscilador Simple 3-16 Integración numérica de la ecuación de movimiento Se trata de resolver en forma discreta la ecuación de movimiento para 1 GLD dada por (2.8) y (2.11) utilizando ecuaciones en diferencias que permitan obtener x , x& y x& & en un instante t+∆t en función, únicamente, de los valores en el instante t. Existen varios métodos para plantear las ecuaciones en diferencias antes mencionadas, pe-
ro solo nos dedicaremos en profundidad al de Newmark, dado que es uno de los más di-
fundidos en software de análisis dinámico mediante FEM. Método de Newmark (1 GLD) F x k x c x m = + + & & & (3.54) Fig. - 3.10 con ( ) ( ) [ ] t a m t f F − + = x& &
x& &
ti ti+1 t
t t t − = ∆
haciendo el cambio de variable i
t t − = τ ; t ∆ ≤ τ puede suponerse que el valor de la aceleración de respuesta en un instante τ se expresa como: ( ) ( )( )
x x f x x & & & & & & & & − + =
τ τ (3.55) f(τ) =
0 para τ = 0
1 para τ = ∆t
integrando (3.55) se obtiene la expresión de la velocidad ( ) ( ) ( )
d f x x x x x
& & & & & & & & (3.56) Capítulo 3 – Respuesta de un Oscilador Simple 3-17 ( ) ( )
d f g (3.57) ( )
τ τ γ
d f t (3.58) La expresión (3.56) queda: ( ) ( ) ( ) τ τ τ g x x x x x
& & & & & & & & − + + =
(3.59) que para τ=∆t resulta ( ) t x x t x x x
∆ − + ∆ + =
γ & & & & & & & &
(3.60) ( ) [ ] t x x x x
∆ + − + =
1 & & & & & & γ γ (3.61) Para calcular los desplazamientos, se integra (3.59) y se obtiene ( ) ( ) ( )
d g x x x x x x
& & & & & & & (3.62) particularizando para τ=∆t y llamando ( )
τ τ β
d g t (3.63) se obtiene la relación final en diferencias, propuesta por Newmark 2
& & & & & β β τ (3.64) Las ecuaciones (3.61) y (3.64) juntamente con la ecuación de movimiento (3.54) permiten obtener aceleraciones, velocidades y desplazamientos en un instante t+∆t con solo conocer estas magnitudes en el instante t. Capítulo 3 – Respuesta de un Oscilador Simple 3-18 Los parámetros γ y β surgen de un análisis de estabilidad numérica que escapa a este curso (ver. ref. [1], [2] y [4]) pero pueden tomarse: 25 , 0
Desarrollo y forma operativa La ecuación de movimiento (3.54) valuada en t=t
toma la forma: 1
= + + i
F x k x c x m & & & (3.65) Reemplazando en ésta los valores en diferencias dados por (3.61) y (3.64), agrupando los valores que permanecen constantes en cantidades denominadas por a
; se obtiene el algoritmo de integración paso a paso. A) cálculos iniciales 1- determinar las propiedades de la estructura: m, c, k 2- inicializar 0
x& y 0
x& & utilizando, de ser necesario, la ecuación (3.54) 3- seleccionar el paso de tiempo y los parámetros γ y β ∆t ; γ ; β ∆t ≤ T/10 ; γ=0,5 ; β=0,25 4- calcular las constantes: 2
a ( ) γ − ∆ = 1
t a ; t a ∆ = γ
5- formar la “rigidez efectiva” c a m a k k
+ + = B) para cada paso de integración 1- formar el término de “carga efectiva” en t+∆t ( ) ( )
x a x a x a c x a x a x a m F r & & & & & &
2- resolver el desplazamiento en t+∆t Capítulo 3 – Respuesta de un Oscilador Simple 3-19 1 1 1
= ⋅ r x k
3- calcular aceleraciones y velocidades en t+∆t ( )
x a x a x x a x & & & & &
1 7 6 1 + +
x a x a x x & & & & & & Capítulo 4 – Respuesta Dinámica de una Estructura con Múltiples Grados de Libertad 4-1 CAPÍTULO 4 - RESPUESTA DINÁMICA DE UNA ESTRUCTURA CON MÚLTIPLES GRADOS DE LIBERTAD Ecuaciones de movimiento y equilibrio dinámico Como se vio en el Capítulo 2 un modelo estructural dinámico esta dado por - discretización espacial (por ej. masas concentradas) - discretización temporal (por ej. Newmark) En el Capítulo 3 se estudio el modelo discreto más sencillo (1 GLD) mientras que en la UT2 se plantearon (sin resolver) las ecuaciones de equilibrio dinámico para sistemas de múltiples GLD, utilizando el principio de D’Alembert. En general, el sistema de ecuaciones diferenciales es del tipo ( ) t F D K D C D M = + +
(4.1) En particular, para el caso sísmico: { } a J M D K D C D M − = + +
(4.2) Recordar que J es un vector con 1 en la posición del GLD en el que actúa la aceleración de apoyo a(t). Vibraciones libres Características dinámicas El sistema que gobierna las vibraciones libres en un sistema de múltiples GLD es 0 = + + D K D C D M
(4.3) Como ya se vio, el cambio de frecuencia propia debido a considerar o no el amortigua-
miento, no es relevante para estructuras civiles corrientes, por lo que seguiremos el análi-
sis con el sistema simplificado: Capítulo 4 – Respuesta Dinámica de una Estructura con Múltiples Grados de Libertad 4-2 0 = + D K D M
(4.4) Una solución a este sistema puede ser de la forma ( ) Ψ + = t sen ω A D (4.5) A: vector que contiene las amplitudes de las vibraciones Ψ: ángulo de fase ( ) ( ) Ψ + − = t sen t ω ω A D
(4.6) reemplazada en (4.4) queda: ( ) ( ) 0
= Ψ + + Ψ + − t sen t sen ω ω ω A K A M { } ( ) 0
= Ψ + − t sen ω ω A M K (4.7) para que haya vibraciones, ω≠0, por lo que ( ) 0 sen ≠ Ψ + t ω y puede eliminarse de (4.7), quedando finalmente { } 0
= − A M K ω (4.8) Nos interesan las soluciones de A distintas de la trivial. La ecuación (su resolución) 4.8 representa un problema de autovectores y autovalores, en donde: 0
= − M K ω (4.9) Este determinante puede desarrollarse en la forma polinómica: 0 ...
α ω α ω α ω α ω (4.10) obteniéndose la ecuación característica. Capítulo 4 – Respuesta Dinámica de una Estructura con Múltiples Grados de Libertad 4-3 Cuando K y M son definidas positivas (caso usual en estructuras civiles), de la ecuación característica se obtienen “n” soluciones positivas ω
y en consecuencia “n” valores ω
, siendo “n” el número de GLD de la estructura. Los valores ω
se denominan “frecuencias propias” o “pulsaciones” de la estructura y los “n” periodos propios se calculan 1,...n i
(4.11) Siendo T
el período correspondiente ω
, que es la frecuencia de menor valor, éste se de-
nomina “periodo fundamental” del sistema. Reemplazando cada ω
en la (4.8) se obtiene el correspondiente A
que se denomina vector de “forma” modal o simplemente “modo”. Este vector “modo” contiene la forma que tomara la estructura en cada vibración. Para ca-
da frecuencia vibratoria, la forma modal será diferente. Es importante notar que (4.8) se cumple para cualquier múltiplo de A
, por lo tanto inferimos que no interesa la magnitud de las componentes de A
sino la relación existente entre ellas. Por esto es conveniente “normalizar” los distintos A
. Normalización de los modos 1. Por la máxima componente dado (
A puede dividirse cada elemento por el a
= ϕ ϕϕ ϕ (4.12) (
ϕ ϕϕ ϕ modo i normalizado 2. Según la matriz de masas ( )
M A ϕ ϕϕ ϕ (4.13) Capítulo 4 – Respuesta Dinámica de una Estructura con Múltiples Grados de Libertad 4-4 Con [ ]
M A M A = * (4.14) Esta normalización permite asegurar la condición [ ] 1 =
ϕ ϕϕ ϕ ϕ ϕϕ ϕ M (4.15) Expresión que será de suma importancia más adelante. Obtención de los grados de libertad dinámicos No siempre todos los grados de libertad estáticos que definen el comportamiento de la es-
tructura tienen asociada “masa”, es decir que no son necesarios en el análisis dinámico. Por ejemplo los pórticos de las figuras 2-15 (f) y (h) que son analizados con modelos di-
námicos de menor orden que los respectivos estáticos (GLE). Lo mismo se ve en la figura 2-16. La forma más simple de reducir en número de grados de libertad sin perder precisión es mediante la: Condensación estática de la matriz de rigidez Dada una estructura y su modelo estático puede subdividirse la matriz de rigidez (así como el vector desplazamientos) de manera de separar las ecuaciones que tienen asociada masa (GLD) y las que no la tienen (GLE). modelo dinámico
Fig. 4.1 - Condensación estática Para el ejemplo de la figura 4-1, tenemos que, observando el vector desplazamientos está-
ticos D
podemos determinar que parte queremos eliminar y cuales grados de libertad de-
seamos conservar como GLD. Capítulo 4 – Respuesta Dinámica de una Estructura con Múltiples Grados de Libertad 4-5 (
D: vector con los GLD D
: grados de libertad a eliminar s s s
F D K = (4.1) El sistema de ecuaciones estático puede escribirse (
(4.16) Como F
representa las fuerzas en los grados de libertad a eliminar 0 =
F (4.17) La segunda ecuación (vectorial) de (4.16) puede, entonces, expresarse como 0
D K D K (4.18) de donde D K K D
(4.19) Desarrollando ahora la primera ecuación de (4.16) tendremos { } F D K K K D K = −
(4.20) reordenando { } F D K K K K =
- (4.21) F D K = (4.22) donde K es la matriz de rigidez condensada y vale Capítulo 4 – Respuesta Dinámica de una Estructura con Múltiples Grados de Libertad 4-6 21
= - (4.23) Finalmente el sistema dinámico de la fig. 4-1 estará representado matemáticamente por 0
D M D K
D (4.24) Otra forma de obtener la matriz de rigidez del modelo dinámico es invirtiendo la matriz de flexibilidad de los GLE que tienen asociada masa. Este método, si bien menos formal, puede ser de más sencilla aplicación para procedi-
mientos manuales. Para nuestro ejemplo: 1
Fig. 4-2 - Obtención de K mediante F
Matriz de amortiguamiento Dado el sistema Capítulo 4 – Respuesta Dinámica de una Estructura con Múltiples Grados de Libertad 4-7 ( ) t F D K D C D M = + +
(4.1) solo nos resta conocer (para encarar la resolución) la matriz C, puesto que en el apartado anterior se mostraron dos maneras de obtener K y sabemos (porque usamos un modelo discreto de masas concentradas) que: (
M (4.25) Conviene definir matrices de amortiguamiento ortogonales, pero para esto debemos, pri-
mero incursionar un poco en las condiciones de ortogonalidad de los modos propios. Condiciones de ortogonalidad Las frecuencias y modos propios pueden ordenarse en matrices denominadas espectrales y modales, que son, respectivamente: (
O Ω ΩΩ Ω ; [ ]
ϕ ϕϕ ϕ ϕ ϕϕ ϕ ϕ ϕϕ ϕ K
, = Φ ΦΦ Φ (4.26) Puede demostrarse (ver ref. [1], [2], [3], [4]) que 0 =
ϕ ϕϕ ϕ ϕ ϕϕ ϕ i ≠ j (4.27) Esta condición para los modos puede extenderse a 0 =
ϕ ϕϕ ϕ ϕ ϕϕ ϕ M i ≠ j (4.28) Capítulo 4 – Respuesta Dinámica de una Estructura con Múltiples Grados de Libertad 4-8 0 =
ϕ ϕϕ ϕ ϕ ϕϕ ϕ K i ≠ j (4.29) Si los vectores ϕ ϕϕ ϕ
fueron normalizados según (4.13) las condiciones de ortogonalidad y normalización pueden expresarse como una única condición de ortonormalidad I M = Φ ΦΦ Φ Φ ΦΦ Φ
(4.30) *
K K = Φ ΦΦ Φ Φ ΦΦ Φ (4.31) donde I es la matriz identidad y K* es una matriz diagonal (
k ϕ ϕ K * = (4.32) Matrices de amortiguamiento ortogonales Una de las hipótesis para lograr una representación numérica del amortiguamiento de una estructura está dada por suponer que existe un mecanismo de disipación uniforme de ener-
gía. En tal caso puede desarrollarse una matriz de amortiguamiento que cumpla las condi-
ciones de ortogonalidad respecto a la matriz modal. El “amortiguamiento proporcional” permite definir una matriz que sea proporcional a la de las masas, a la de rigidez o a ambas. K M C
α α + = (4.33) La condición de ortogonalidad: 0 =
ϕ ϕϕ ϕ ϕ ϕϕ ϕ C i ≠ j (4.34) i i i
ξ ω 2 = ϕ ϕϕ ϕ ϕ ϕϕ ϕ C (4.35) Capítulo 4 – Respuesta Dinámica de una Estructura con Múltiples Grados de Libertad 4-9 Los coeficientes α
se calculan a partir de las ecuaciones anteriores (4.35) y (4.33) [ ]
ξ ω α α 2
= + ϕ ϕϕ ϕ ϕ ϕϕ ϕ K M (4.36) i i i
= + ϕ ϕϕ ϕ ϕ ϕϕ ϕ ϕ ϕϕ ϕ ϕ ϕϕ ϕ K M (4.37) i i i
ξ ω ω α α 2
= + (4.38) Recordar que si la normalización de los autovectores ϕ ϕϕ ϕ
fue hecha según (4.13) y (4.14) el producto 2
ω = ϕ ϕϕ ϕ ϕ ϕϕ ϕ K (4.39) Esto se debe a que la (4.8) puede expresarse para A = ϕ ϕϕ ϕ
y premultiplicarse por ϕ ϕϕ ϕ
T [ ] 0
-ω ϕ ϕϕ ϕ ϕ ϕϕ ϕ M K (4.40) 0
ω ϕ ϕϕ ϕ ϕ ϕϕ ϕ ϕ ϕϕ ϕ ϕ ϕϕ ϕ M K (4.41) Determinación práctica de modos y frecuencias Este es, por sí mismo, uno de los problemas más complejos de la dinámica estructural. Generalmente no es necesaria la resolución de todas las frecuencias y sus modos, sino que solo interesan las primeras “q” que representan las posibilidades ciertas de vibrar pues ne-
cesitan menos energía de excitación. Cuando el problema es de pocos GLD es posible resolver el determinante de (4.9) o (4.10) obteniéndose las “n” ω
y luego los “n” ϕ ϕϕ ϕ
, donde “n” es el número de GLD. Para proble-
mas con muchos GLD se utilizan técnicas numéricas (aproximadas) para obtener los pri-
meros “q” pares de (ω
; ϕ ϕϕ ϕ
). En la mayoría de los casos prácticos de ingeniería q << n. En las referencias bibliográficas pueden encontrarse descripciones detalladas de varios métodos. Nosotros desarrollaremos solo uno que es relativamente simple y preciso. Método de Stodola-Vianello Partiendo de la (4.8) y recordando que se cumple para cualquier múltiplo del autovector A podemos escribir Capítulo 4 – Respuesta Dinámica de una Estructura con Múltiples Grados de Libertad 4-10 0
= − A M A K ω (4.42) A M A K
ω = (4.43) Suponiendo un vector inicial A
conocido (cercano a A) 0 1
ω = (4.44) 0 -1 1
A M K A ≈ (4.45) que lleva a la fórmula de recurrencia i -1 1 i
A M K A ≈
(4.46) Nota: una descripción detallada de este método puede encontrarse en [4], parte 11.2 “Vector Iteration Methods”, así como prueba de la convergencia del método. Procedimiento de cálculo 1- Se propone un vector A
inicial. Conviene que los valores correspondan a una defor-
mada suave, no obstante esto solo acelera la convergencia. 2- Se calcula 0 -1 1
A M K A = (4.47) 3- Se normaliza A
= A M A A A
i=1 (4.48) 4- Se mejora la solución i -1 1 i
A M K A =
(4.49) 5- Se normaliza la nueva solución { }
A M A A A (4.50) 6- Se calcula la frecuencia correspondiente 2
ω (4.51) Capítulo 4 – Respuesta Dinámica de una Estructura con Múltiples Grados de Libertad 4-11 7- Se controla la convergencia 4 paso al vuelve se
→  ≤
(4.52) ↓SI 1
A ϕ ϕϕ ϕ ; 1
ω ω (4.53) Se obtuvo el primer autovalor y su correspondiente autovector. Siempre: ∞ → i ; 1
ϕ ϕϕ ϕ →
ω ω →
Obtención de modos y frecuencias superiores El procedimiento es el mismo que para el modo 1, solamente debemos garantizar que el vector de “arranque” sea ortogonal a ϕ ϕϕ ϕ
Se supone que un vector inicial ∑
ϕ ϕϕ ϕ A (4.54) que para hacerlo ortogonal a ϕ ϕϕ ϕ
le “restamos” la componente q
ϕ ϕϕ ϕ q * − = A A (4.55) Para determinar q
, premultiplicamos la (4.54) por [ϕ ϕϕ ϕ
M] y aplicando las propiedades de ortogonalidad nos queda [ ] [ ]
ϕ ϕϕ ϕ ϕ ϕϕ ϕ ϕ ϕϕ ϕ M A M (4.56) 1 1
ϕ ϕϕ ϕ ϕ ϕϕ ϕ
q = (4.57) y en general i
= (4.58) y para los vectores de arranque (ortogonales a los previamente calculados) ∑
q * ϕ ϕϕ ϕ A A (4.59) Capítulo 4 – Respuesta Dinámica de una Estructura con Múltiples Grados de Libertad 4-12 Procedimiento: 1. Se estima A
2. Se calcula q
con (4.58) 3. Se calcula A
con (4.59) 4. Se sigue con el método de Stodola propiamente dicho Resolución de las ecuaciones de movimiento en estructuras con múltiples GLD Descomposición y superposición modal Si bien desarrollaremos el método para el caso sísmico, todo lo expuesto es valido para el caso general. La ecuación que gobierna el comportamiento de una estructura de múltiples GLD es la { } a J M D K D C D M − = + +
(4.2) Las vibraciones libres no amortiguadas se estudian mediante 0 D K D M = +
(4.4) al cual, según lo expuesto en puntos anteriores, le corresponden “n” pares de frecuencias y “modos” que son solución del sistema de ecuaciones algebraicas [ ] 0
= − ϕ ϕϕ ϕ M K ω (4.8) Recordemos que los autovectores ϕ
son ortogonales respecto a las matrices de masa y ri-
gidez. Es por esto que pueden formar una base completa para el espacio de los desplaza-
mientos estructurales, es decir, es posible escribir ( )
ϕ ϕϕ ϕ D (4.60) Capítulo 4 – Respuesta Dinámica de una Estructura con Múltiples Grados de Libertad 4-13 Donde y
(t) es un escalar función del tiempo a determinar, llamado “respuesta generaliza-
da”. Sustituyendo (4.60) en (4.2) obtenemos ( ) ( ) ( ) a t y t y t y
J M K C M − = + +
ϕ ϕϕ ϕ ϕ ϕϕ ϕ ϕ ϕϕ ϕ & & & (4.61) Premultiplicando (4.61) por un ϕ ϕϕ ϕ
cualquiera, se cumple que: *
M M ϕ ϕϕ ϕ ϕ ϕϕ ϕ ϕ ϕϕ ϕ ϕ ϕϕ ϕ (4.62) *
K K ϕ ϕϕ ϕ ϕ ϕϕ ϕ ϕ ϕϕ ϕ ϕ ϕϕ ϕ (4.63) Y si, como es habitual, se está frente a matrices de amortiguamiento proporcionales y or-
togonales (ver punto 4.4) *
C C ϕ ϕϕ ϕ ϕ ϕϕ ϕ ϕ ϕϕ ϕ ϕ ϕϕ ϕ (4.64) Entonces la (4.61) premultiplicada por un ϕ ϕϕ ϕ
cualquiera queda como una ecuación de 1 GLD ( ) ( ) ( ) ( ) t a t y K t y C t y M
J M ϕ ϕϕ ϕ − = + + * * * & & & (4.65) Recordando que: j j j j
ξ ω C 2 * = = ϕ ϕϕ ϕ ϕ ϕϕ ϕ C (4.35) 2
ω K = = ϕ ϕϕ ϕ ϕ ϕϕ ϕ K (4.39) Y dividiendo ambos miembros por j
M ϕ ϕ M = * Capítulo 4 – Respuesta Dinámica de una Estructura con Múltiples Grados de Libertad 4-14 ( ) ( ) ( )
t y ω t y ω ξ t y
2 & & & (4.66) Al valor j
= (4.67) se lo denomina “coeficiente de participación modal”. Finalmente, el sistema (4.2), al ser “proyectado” según la base formada por la matriz mo-
dal Φ ΦΦ Φ = [ϕ ϕϕ ϕ
,… ϕ ϕϕ ϕ
], queda en “n” ecuaciones diferenciales de 1 GLD, del tipo: ( ) ( ) ( ) ( ) t a Q t y ω t y ω ξ t y
2 & & & (4.68) Esta ecuación puede resolverse con cualquiera de los métodos vistos en la UT3. Una vez obtenidos los “n” y
(t), puede obtenerse la solución estructural mediante (4.60). En general, los modos bajos son los que contienen menos energía elástica de deformación y por ende los que más contribuyen a la respuesta estructural. Usualmente ( ) ( )
ϕ ϕϕ ϕ D ; q < n (4.69) Procedimiento de cálculo 1- Se determinan los “q” primeros modos y frecuencias 2- Se resuelven las “q” ecuaciones (4.68) i=1,…q 3- Se obtiene la historia en el tiempo según (4.69) y sus derivadas 4- Se obtienen las fuerzas elásticas, de inercia y de amortiguamiento para cada instante Integración directa de las ecuaciones de movimiento Como su nombre lo indica, este método no requiere ninguna transformación previa de las ecuaciones de movimiento. Consiste, básicamente en obtener la solución en una cierta can-
tidad (discreta) de pasos de tiempo. Es por esto que también es llamada integración “paso a paso”. Capítulo 4 – Respuesta Dinámica de una Estructura con Múltiples Grados de Libertad 4-15 Como ya adelantamos, existen numerosas variantes y algoritmos para la integración numé-
rica de las ecuaciones de movimiento, pero nuestros desarrollos estarán dados en el méto-
do de Newmark. En esencia es igual al planteo dado en el Capítulo 3, solo que en vez de tratarse de una ecuación de 1 GLD se trata de un sistema de “n” GLD. Procedimiento de cálculo A) Iniciales 1. Ensamblar las matrices M, C y K 2. Inicializar 0
D , 0
3. Seleccionar el paso de tiempo y los parámetros γ y β ∆t ≤ T
/10 ; γ=0,5 ; β=0,25 4. Calcular las constantes a
5. Formar la matriz de rigidez efectiva: C M K K
a a + + = (4.70) B) Para cada paso de tiempo 1. Formar el término de “carga efectiva” en t+∆t (*1) ( )
D D D M F r
(4.71) 2. Resolver el desplazamiento en t+∆t 1 1
r D K (4.72) *1: el término ( ) { } t a i i J M F F − + = +1 se refiere a fuerzas generalizadas 3. Calcular aceleraciones y velocidades en t+∆t ( )
a a a D D D D D
(4.73) 1 7 6 1 + +
a a D D D D
(4.74) # Notar que para cada paso de tiempo se debe resolver el sistema (4.72) lo que puede hacer excesivo el costo de cálculo y almacenamiento. Respuesta máxima utilizando espectros de respuesta Las ecuaciones desacopladas (4.68) pueden resolverse utilizando los espectros de respues-
ta: la máxima aceleración será Capítulo 4 – Respuesta Dinámica de una Estructura con Múltiples Grados de Libertad 4-16 ( ) ( )
S Q t y = & & (4.75) entonces, el máximo desplazamiento es ( ) ( )
= (4.76) Podemos, entonces calcular los máximos desplazamientos (en todos los GLD) para el mo-
do “j”: ( ) ( )
ϕ ϕϕ ϕ ϕ ϕϕ ϕ = = D (4.77) ( )
A D = (4.78) Donde A
es el vector de coeficientes de participación modal correspondientes al modo de vibración “j”. Puesto que el máximo en cada grado de libertad no se produce en el mismo instante, la respuesta total máxima no es la suma de los máximos de cada modo!! ∑
D D (4.79) Una forma muy utilizada (y precisa) de evaluar la respuesta máxima (desplazamientos, velocidades, aceleraciones y esfuerzos) como combinación de los máximos modales es: ( )
R R (4.80) denominada SRSS Capítulo 5 – Análisis de Construcciones con Efecto Dinámico de Viento 5-1 CAPÍTULO 5 - ANÁLISIS DE CONSTRUCCIONES CON EFECTOS DINÁMICOS DE VIENTO flujo
T(z): fuerza aerodinámica de empuje, paralela a la dirección del flujo L(z): fuerza aerodinámica de deriva, perpendicular a la direc-
ción del flujo Generalmente son dominantes las fuerzas de empuje, pero para ciertas estructuras: - bajo amortiguamiento - acusada flexibilidad (poca rigidez) - construcciones livianas (poca masa) las acciones perpendiculares de “deriva” pueden tener efectos significativos que deben ser contemplados en el análisis. CIRSOC 102-1 → procedimientos simplificados Validez: a) construcciones prismáticas o cuasiprismáticas b) primer modo dominante y de forma aproximadamente lineal c) periodo fundamental T
> 1 seg d) amortiguamiento ξ < 0,01 Acciones paralelas a la dirección del viento (procedimiento basado en el “factor de ráfaga”) ( ) z x
( ) z X max
( ) z x
( ) ( ) ( ) z x z x z X
max + = ( ) z x : desplazamiento medio ( ) z x
: desplazamiento fluctuante debido a la turbulencia variable con el tiempo Capítulo 5 – Análisis de Construcciones con Efecto Dinámico de Viento 5-2 Se define como factor de ráfaga ( )
+ =1 entonces ( ) ( ) ( ) z x z G z X max ⋅ = Hipótesis básicas para poder reemplazar la acción dinámica del viento turbulento por un procedimiento estático equivalente. a) Comportamiento elástico lineal de la estructura. b) El modo fundamental de vibración es una función lineal de la altura. c) La contribución de los modos superiores al primero en respuesta se considera despreciable, por lo que G(z) = G = cte. d) La velocidad media del viento es promediada sobre intervalos de una hora. e) La variación de la velocidad media del viento varía según ( )
0706 , 0
: velocidad básica de diseño (m/seg) z
: parámetro que depende de la rugosidad i z
para rugosidad I (referencia) Procedimiento de cálculo I) Presiones La presión dinámica que incluye el efecto de la turbulencia del viento se determina mediante 0 2
' q c c G q
⋅ ⋅ ⋅ = G: factor de ráfaga c
: variación por rugosidad y altura (art. 5.2.4.2, CIRSOC 102) c
: factor por cambio de tiempo en velocidad media (Tabla 3/pág. 20, CIRSOC 102-1) q
: presión dinámica básica (art. 5.2.3 , CIRSOC 102) Capítulo 5 – Análisis de Construcciones con Efecto Dinámico de Viento 5-3 )
K G 234 , 1 1 → se calcula mediante tablas y ábacos que utilizan una serie de parámetros auxiliares En cada nivel se comparará: NO → se adopta q
para los esfuerzos q’
SI → se adopta q’
para los esfuerzos CIRSOC 102 → 0
⋅ ⋅ = II) Aceleraciones ( ) ( ) z K z X
σ ⋅ = (paralela al viento) K → fig. 13, CIRSOC 102-1 ( ) z
σ → valor medio cuadrático de las aceleraciones (pág. 21) III) Verificaciones (paralelas al viento) max
a- Dimensionado estructural b- Verificación de confort con gráficos y tablas en función de max
, y T c- Verificación de deformaciones admisibles 350 500
< < Acciones perpendiculares a la dirección del viento I) Resonancia L(z)
torbellinos de Bèrnard - Karman Capítulo 5 – Análisis de Construcciones con Efecto Dinámico de Viento 5-4 La velocidad crítica del viento que produce el fenómeno de resonancia es T S
= d: ancho de la superficie maestra (puede ser variable) S: nº de Strouhal ¹
0,30 a 0,25 prismas
0,20 cilindros
T: período propio ¹
⊥ → ⊥
// // V
> 25 m/seg, entonces puede prescindirse del cálculo de la resonancia. De lo contrario: Fuerzas de deriva: (a la velocidad crítica) ( ) d
q z L
ξ: fracción del amortiguamiento (
000613 , 0
V V Fuerzas de empuje: (a la velocidad crítica) se admite distribución uniforme ( ) d q G C T z T
cr E z
⋅ ⋅ ⋅ ⋅ = = 8 , 0 C
: coeficiente global de empuje CIRSOC 102 G: factor de ráfaga correspondiente a V
→ con V
// Las fuerzas L(z) y T
obtenidas para la velocidad crítica (correspondientes al período per-
pendicular y paralelo respectivamente) se suman de la siguiente manera: ( ) ( )
T z L z F + = y se debe comparar con las correspondientes obtenidas con q’
dadas para la velocidad de diseño. Capítulo 5 – Análisis de Construcciones con Efecto Dinámico de Viento 5-5 Criterios de confort en edificios que oscilan El confort de los ocupantes de los edificios de gran altura que soportan la acción del viento es un tema de primordial importancia en el diseño. La respuesta humana al movimiento oscilatorio de las construcciones abarca una extensa gama de reacciones, con efectos psicológicos y fisiológicos tales como mareos, ansiedad, molestias visuales o temor, llegando hasta náusea aguda. Por el contrario, otras persones sienten placer por el movimiento y la experiencia poco usual de estar en una construcción que oscila. Adicionalmente, las oscilaciones excesivas producen fisuración de la tabiquería, rotura de vidrios de ventanas y caída de revestimien-
tos, lo cual incide negativamente en el valor de una propiedad y su rentabilidad. Sin embargo, un edificio resultaría demasiado costoso si se construye o equipa de modo que pueda soportar sin movimientos perceptibles una tormenta con vientos huracanados o un fuerte sismo. En consecuencia, los movimientos son casi inevitables y el problema del diseño consiste en mantenerlos dentro de los límites aceptables para no perturbar el confort y el bienestar de los usuarios. Por otra parte, el costo del edificio no debe superar los valores normales de una aceptable economía. Los factores que pueden producir vibraciones en un edificio son numerosos, tales como maquinaria en funcionamiento defectuoso, paso de vehículos pesados por el lugar, impac-
tos en rampas, vientos fuertes, sismos, etc. los cuales pueden variar durante la vida de ser-
vicio de la estructura. En general la aceleración es el parámetro predominante para determinar aproximadamente la naturaleza de la respuesta humana a las vibraciones. Las curvas de la fig. 5.3.1 grafican los límites de confort obtenidos del análisis de un gran número de edificios altos, indican-
do las máximas aceleraciones aceptables para diferentes frecuencias, dependiendo del uso o destino del edificio. Los datos se obtuvieron para las aceleraciones pico de las más fuertes tormentas ocurridas durante un periodo de retorno de mas de 5 años. Capítulo 5 – Análisis de Construcciones con Efecto Dinámico de Viento 5-6 0.40
( ) m f ∆
fig. 5.3.1 GRAFICA DE CONFORT AMPLITUD DE OSCILACION EN FUNCION DEL PERIODO PARA DISTINTOS VALORES DE LA ACELERACION EN PORCENTAJE DE g Capítulo 5 – Análisis de Construcciones con Efecto Dinámico de Viento 5-7 La tabla 5.3.1 a continuación da idea de la magnitud de las percepciones que se obtuvieron en experimentos realizados para diferentes niveles de aceleraciones (según Khan Parmelec y según Chang). TABLA 5.3.1 PERCEPCION HUMANA ACELERACION ü (m/seg
) (según Khan y Parmelec) No perceptible ü ≤ 0,004 g Levemente perceptible 0,004 g < ü ≤ 0,0075 g Perceptible 0,0075 g < ü ≤ 0,02 g Molesta 0,02 g < ü (según Chang) No perceptible ü ≤ 0,005 g Perceptible 0,005 g < ü ≤ 0,015 g Desagradable 0,015 g < ü ≤ 0,05 g Muy desagradable 0,05 g < ü ≤ 0,15 g Intolerable 0,15 g < ü Capítulo 5 – Análisis de Construcciones con Efecto Dinámico de Viento 5-8 La sensación subjetiva y el comportamiento humano afectado por las diferentes acelera-
ciones se indica en la tabla 5.3.2 (según Yamada y Goto). TABLA 5.3.2 NIVEL DE PERCEPCION HUMANA RANGO ACELERACION (m/seg
) EFECTO 1 0,05 La gente no percibe el movimiento. 2 0,05 – 0,10 Las personas sensibles perciben el movimiento. Los objetos colgantes pueden moverse algo. 3 0,10 – 0,25 La mayoría de las personas perciben el movimien-
to. La oscilación puede afectar el trabajo de ofici-
na. La exposición de larga duración puede producir malestar. 4 0,25 – 0,40 El trabajo de oficina se vuelve difícil o casi impo-
sible. Aun se puede caminar. 5 0,40 – 0,50 Se percibe fuertemente el movimiento. Hay difi-
cultad para caminar normalmente. Las personas de pie pueden perder el equilibrio. 6 0,50 – 0,70 No se tolera el movimiento y no se puede caminar. 7 - 8 > 0,85 Los objetos caen y pueden lastimar a las personas. Capítulo 5 – Análisis de Construcciones con Efecto Dinámico de Viento 5-9 0.15
1 PERIODO (seg)
0.1 %g
1.1 %g
1.5 %g
15 %g
44 %g
100 %g
fig. 5.3.2 GRAFICA DE CONFORT INTERACCION DE LAS VARIABLES: PERIODO (T), AMPLITUD (∆f) Y ACELERACION (x’’) Capítulo 5 – Análisis de Construcciones con Efecto Dinámico de Viento 5-10 GRADUACIÓN DEL EFECTO SOBRE LAS PERSONAS ACEPTABILIDAD EN EDIFICIOS ESTADO PREVISIBLE DE LA ESTRUCTURA A INTOLERABLE NO COLAPSO B MUY POCO TOLERABLE NO DAÑOS LOCALES C DEMASIADO PERCEPTIBLE SITUACION LIMITE FORMACION DE GRIETAS D MUY PERCEPTIBLE EN TAREAS INDUSTRIALES PESADAS FORMACION DE FISURAS E PERCEPTIBLE EN PERIODOS BREVES EN VIVIENDAS SIN INFLUENCIA EN ESTRUCTURAS CORRIENTES F ESCASAMENTE PERCEPTIBLE PERIODOS LARGOS EN VIVIENDAS SIN INFLUENCIAS G NO PERCEPTIBLE PERIODOS LARGOS EN VIVIENDAS SIN INFLUENCIAS CUADRO ADJUNTO A GRAFICA DE LA fig. 5.3.2 Capítulo 6 – Análisis de Construcciones con Efectos Sísmicos 6-1 CAPÍTULO 6 - ANÁLISIS DE CONSTRUCCIONES CON EFECTOS SISMICOS Definición numérica de la acción sísmica La definición correcta de la acción sísmica es un problema al que se le debe dar la mayor importancia en un análisis estructural sísmico. La solución a este problema parte, generalmente, de los datos experimentales proporcio-
nados por la sismología. Parámetros de terremotos pasados en una región → Predicción de las características sísmicas de los terremotos que afectarán dicha región Tradicionalmente, la fuerza destructiva de un terremoto ha sido expresada en función de la aceleración máxima del terreno, pero se han observado daños moderados a aceleraciones muy altas. Por lo tanto, se deben tener en cuenta otras características: - intensidad, duración - contenidos de frecuencias - secuencia de choques, etc. Por ejemplo, se han utilizado, además de los valores máximos de la aceleración, velocidad y desplazamiento del terreno, el espectro de amplitudes de Fourier, el espectro de seudove-
locidades, intensidad espectral (Housner), el valor medio cuadrático de las aceleraciones en la fase “fuerte” de un acelerograma, etc. Además, la acción definida estará totalmente relacionada al tipo de análisis estructural que se va a realizar: Análisis lineal → desacoplamiento modal → espectros de respuesta Análisis no lineal → integración directa → acelerogramas Capítulo 6 – Análisis de Construcciones con Efectos Sísmicos 6-2 # También existen definiciones de acciones sísmicas a partir de la teoría de procesos esto-
cásticos, en donde la respuesta estructural en sí se obtiene probabilísticamente. No profundizaremos sobre estos temas, pues escapan a nuestro curso. (ver ref. [1 y 2]) Definición mediante espectros de respuesta Es la forma más usual de definir una acción sísmica, dado que se obtienen descripciones de las características más importantes de la respuesta estructural sin la necesidad de dispo-
ner de una historia en el tiempo de la excitación y la respuesta. Otra ventaja es que un espectro puede “modificarse” en base a las características del lugar de emplazamiento de la estructura, sin necesidad de conocer los detalles de la excitación. Espectro de respuesta → representaciones gráficas de valores aproximados de la respuesta máxima de un sistema de 1 GLD lineal y elástico: ( ) t a x x x − = ⋅ + ⋅ +
2 ω ξω (2.11) ( ) ( )
d t e a a S τ τ ω τ
⋅ − ⋅ ⋅ ⋅ ⋅ − =
(3.49) ( ) ( ) ξ ω ω ξ ω , , , , a S a S
⋅ = (3.52) ( ) ( ) ξ ω ω ξ ω , , , ,
⋅ = (3.53) Puede realizarse un análisis de las variaciones de los valores de las curvas S
, en función de las características de la estructura. Recordar que m
= ω Para frecuencias propias altas (en comparación con la del movimiento del terreno), la (2.11) quedaría (despreciando los dos primeros términos): ( )
→ ∞ → ∴
− ≈ ⋅
terreno) del máxima n aceleració la (copia
a S ; si
rígidas muy s estructura
ω t a x
Capítulo 6 – Análisis de Construcciones con Efectos Sísmicos 6-3 suavizado y promediado de valores Para frecuencias bajas, por el contrario, la (2.11) queda: ( )
terreno) del máximos entos desplazami los (copia
d S 0; si
flexibles muy s estructura
Para frecuencias intermedias, se produce una amplificación del movimiento del terreno en su paso a través del filtro estructural. ( ) ( ) ( ) ; , ; , ; ,
T = = = ξ β ξ β ξ β (6.1) Las condiciones locales del terreno, tales como grosor y propiedades de los estratos que se encuentran entre la roca firme y la estructura modifican los espectros de respuesta. Conjunto de espectros de respuesta registrados para una región Espectros de Diseño ↓ INPRES – CIRSOC 103 Parte 1 / Cap. 7 s
Para ξ = 0,05 (línea llena) ( )
⋅ − + = T ≤ T
1 Capítulo 6 – Análisis de Construcciones con Efectos Sísmicos 6-4 b S
< T |
sísmica zona
4 tabla , ,
f T y T b a
Para casos con 0,005 ≤ ξ < 0,05 (línea a trazos) ( )
a b f a S
⋅ − ⋅ + = T ≤ T
1 b f S
: factor de amplificación por amortiguamiento %) en (expresado
f Definición mediante acelerogramas a) Acelerogramas reales Basar un cálculo sísmico en uno o varios registros disponibles en una zona implica un alto riesgo de definición incorrecta de la acción. b) Acelerogramas sintéticos Tiene grandes ventajas respecto al anterior: - posibilidad de generar señales de corta duración que tengan el mismo efecto en las estructuras que el del terremoto real que se quiere simular ⇒ economía computacional - al estar generados en función de un espectro de diseño, se tiene en cuenta (aprox.) las condiciones locales del suelo La generación de acelerogramas sintéticos requiere procedimientos matemáticos basados en procesos estocásticos y expansiones en series sinusoidales. Capítulo 6 – Análisis de Construcciones con Efectos Sísmicos 6-5 Existen muchos procedimientos y algoritmos que permiten generar acelerogramas sintéti-
cos a partir de un espectro dado. INPRES – CIRSOC 103 (Parte 1 – Cap. 14.3.1) Los acelerogramas deben satisfacer: a) ( )
a t a γ ⋅ ≥ a
: ordenada al origen del espectro correspon-
diente γ
: factor de riesgo, art. 5.2 b) El espectro elástico (para ξ = 0,05) deberá cumplir: C
Area Area ⋅ ≥ γ d
S S γ ⋅ ⋅ ≥ 7 , 0 para todos los puntos Cantidad de acelerogramas a emplear: Grupo A y B nº ≥ 3 acelerogramas Grupo A
nº ≥ 4 acelerogramas Para diseño y verificaciones se promediarán las envolventes de solicitaciones y deforma-
ciones obtenidas para cada acelerograma, pero en dicho promedio no se incluirán valores inferiores al 85% del máximo encontrado. Capítulo 6 – Análisis de Construcciones con Efectos Sísmicos 6-6 Métodos de análisis según INPRES – CIRSOC 103 Procedimientos con fuerzas estáticas equivalentes a) Método estático (art. 14.1) b) Procedimiento aproximado (cap. 16) Métodos dinámicos a) Análisis modal espectral b) Superposición modal paso a paso c) Integración directa paso a paso 6... Métodos dinámicos INPRES – CIRSOC 103 a) Análisis modal espectral • excitación sísmica traslacional • materiales lineales y elásticos (según art. 12.1 / pág. 47) • ( ) g S S
⋅ ⋅ = γ a
S : del espectro correspondiente γ
: factor de riesgo g: aceleración de la gravedad • se podrá considerar capacidad de disipación de energía por deformaciones anelásticas de la estructura ( )
= - las deformaciones calculadas con este criterio deberán ser amplificadas multiplicándolas por la ductilidad global µ - para el caso anterior, es imprescindible que la estructura se comporte en forma uniforme, de manera de garantizar la ausencia de concentración de rótulas plásticas Se debe verificar que el corte (con análisis modal espectral) en cada dirección no resulte inferior al 75% del corte obtenido con el método estático; de no verificarse se incrementa-
ran los efectos del análisis modal de la siguiente manera: Capítulo 6 – Análisis de Construcciones con Efectos Sísmicos 6-7 espectral
= Modos a considerar: q ≥ 3, pero además se deberán incluir todos los modos cuya contribu-
ción al total sea mayor que 5% de la contribución del modo fundamental. Modelo de edificio de cortante r
( ) t a
(6.1) con ¦
M # observar que si en (6.1) se susti-
tuye la masa m
por el peso w
el valor de Q
no cambia si 1 = ⋅ ⋅
ϕ ϕ M ∑ ∑
ϕ ϕ Peso efectivo modal ∑
W → si → 3 2 1
la de total peso
w W (6.4) q ⇒ apto Desplazamientos modales máximos ( ) ( )
i i max i i max i
ϕ ϕ ⋅ ⋅ = ⋅ = x (6.5) Capítulo 6 – Análisis de Construcciones con Efectos Sísmicos 6-8 Desplazamientos modales relativos entre pisos ( ) ( )
i r i r ri
− = ∆ (6.6) Fuerzas sísmicas equivalentes (modales) modo : i
nivel r:
m x F =
(6.7) pero ( )
ri t x x
ϕ (6.8) Como ( ) ( )
i a i max
i S Q t x ⋅ =
i a i ri
ri S Q x ⋅ ⋅ =
(6.9) ( )
r i a i ri
m S Q F ⋅ ⋅ ⋅ =ϕ ( )
Q ⋅ ⋅ ⋅ =ϕ (6.10) fuerza sísmica equivalente en el nivel r para oscilaciones en el modo i en función del peso w
y la fracción de g Vectorialmente: ( )
(6.11) ( )
S Q ⋅ ⋅ ⋅ = ϕ M F (6.11) Capítulo 6 – Análisis de Construcciones con Efectos Sísmicos 6-9 Cortante modal En nivel “r”: ∑
F V (6.12) En la base: ∑
(6.13a) o también: ( )
⋅ = W (6.13b) Comparación con INPRES – CIRSOC 103 m am
V W ⋅ ⋅ =
(Fórmula reglamento) γ
: factor de riesgo R: coeficiente de reducción por disipación i = m ( )
f V W ⋅ ⋅ =
(Fórmula 6.13b, multiplicada por un factor f
) Distribución en altura de V
según INPRES – CIRSOC 103: m n
m: modo k: nivel F
: fuerza sísmica equivalente para el modo m en el nivel k Según la fórmula (6.10) y (6.13a): ( ) ( )
r i ri i
w Q V
ϕ ϕ si queremos emplear la notación ( )
⋅ = ⇒ ⋅ =
ri i ri
ϕ ξ ξ comparándola con (6.10) se deduce que Capítulo 6 – Análisis de Construcciones con Efectos Sísmicos 6-10 103 CIRSOC - INPRES notación 1 1
ξ Momento de vuelco modal (para el nivel r) ( )
h h F M
si r = 0 → M
= momento de vuelco en la fundación Momento de torsión modal C.R. C.M.
V e M ⋅ = corresp. al piso “r” Esfuerzos seccionales modales Para cada subestructura vertical con la parte proporcional a su rigidez que le corresponde de la fuerza sísmica equivalente, se calculan los esfuerzos seccionales para cada modo en estudio. subestructura "P"
nivel "r"
α α α α K K m: nº de subestructuras verticales (se deben adicionar los efectos torsionales si los hubiese) CALCULO DE LA RESPUESTA TOTAL ( )
Introducción - Bibliografía
Esta publicación pretende brindar en una forma rápida y sencilla los conceptos básicos de la dinámica de estructuras aplicada a las construcciones civiles, enfocada desde un punto de vista numérico. Su organización en capítulos se basa en la secuencia de unidades temáticas del programa de la asignatura homónima que se dicta en la carrera de Ingeniería Civil de la Facultad Regional Paraná de la Universidad Tecnológica Nacional. Todos los temas se desarrollan como síntesis de libros y publicaciones de otros autores (ver Bibliografía) solo que organizados según la conveniencia que mejor se ajusta al curso impartido. El Capítulo 1 hace una muy breve reseña a conceptos sobre sismología y su relación con la ingeniería estructural, caracterización, registros, etc. El Capítulo 2 brinda una aproximación a los conceptos básicos de la dinámica estructural, la importancia de la masa, la relación entre velocidad de carga y de reacción de una estructura. Avanza sobre los modelos estructurales dinámicos, grados de libertad y la discretización espacial y temporal. El Capítulo 3 está dedicado a la caracterización y cálculo de osciladores de un grado de libertad dinámico. Son tratadas las vibraciones libres y forzadas, éstas últimas con cargas armónicas y arbitrarias. Se presenta la resolución numérica mediante integración directa. El Capítulo 4 está dedicado a las estructuras con múltiples grados de libertad dinámicos. Se estudian las vibraciones libres y el cálculo de modos y frecuencias propias, también se presenta la resolución de los sistemas dinámicos mediante diversos métodos como son: descomposición y superposición modal, integración directa y respuesta máxima mediante espectros de respuesta. En el Capítulo 5 se realiza una aproximación al tratamiento dinámico de los efectos del viento sobre las construcciones teniendo como base el reglamento CIRSOC 102. Por último en el Capítulo 6 son tratados los efectos sísmicos y el análisis estructural, desde el punto de la dinámica estructural, enfatizando en la definición de la acción sísmica y la comparación de las fórmulas vistas en capítulos anteriores con los reglamentos vigentes.
Un agradecimiento especial a Marlene Jaurena, que con paciencia y dedicación realizó la transcripción de los manuscritos de este trabajo.
Capítulo 1 - NOCIONES BASICAS DE SISMOLOGIA Causas que generan los terremotos o sismos Los sismos desde el punto de vista de la ingeniería y su caracterización Esteva y Rosenblueth: Donovan: Esteva y Villaverde: Esteva: Registro de ondas sísmicas. Parámetros utilizados y mapas de riesgo sísmico Capítulo 2 - CONCEPTOS BASICOS DE DINAMICA ESTRUCTURAL Definición de la acción dinámica Acciones y fuerzas dinámicas Importancia de la masa en el problema dinámico Velocidad de reacción de una estructura Modelos dinámicos característicos Métodos de modelización dinámica Discretización espacial de las estructuras Método de las masas concentradas Ecuaciones de movimiento Principio de Hamilton Principio de los trabajos virtuales Principio de D’Alembert Formulación de la ecuación de movimiento para un sistema de 1 GLD Formulación de las ecuaciones de movimiento para modelos con múltiples GLD Capítulo 3 - RESPUESTA DE UN OSCILADOR SIMPLE Ecuación de movimiento y equilibrio dinámico Fórmula de Geiger Características dinámicas con amortiguamiento Excitación periódica Excitación armónica Excitación arbitraria. Integral de Duhamel Respuesta a un impulso elemental Factor de amplificación dinámica Espectros sísmicos de respuesta Integración numérica de la ecuación de movimiento Método de Newmark (1 GLD) Desarrollo y forma operativa 1-1 1-1 1-2 1-2 1-3 1-3 1-3 1-3 2-1 2-1 2-2 2-3 2-5 2-6 2-8 2-9 2-10 2-13 2-13 2-13 2-13 2-13 2-15 3-1 3-1 3-2 3-3 3-7 3-8 3-10 3-10 3-12 3-12 3-16 3-16 3-18
Capítulo 4 - RESPUESTA DINÁMICA DE UNA ESTRUCTURA CON PLES GRADOS DE LIBERTAD Ecuaciones de movimiento y equilibrio dinámico Vibraciones libres Características dinámicas Normalización de los modos Obtención de los grados de libertad dinámicos Condensación estática de la matriz de rigidez Matriz de amortiguamiento Matrices de amortiguamiento ortogonales Determinación práctica de modos y frecuencias Método de Stodola-Vianello Resolución de las ecuaciones de movimiento en estructuras con múltiples GLD Descomposición y superposición modal Integración directa de las ecuaciones de movimiento Respuesta máxima utilizando espectros de respuesta Capítulo 5 -
MÚLTI4-1 4-1 4-1 4-1 4-3 4-4 4-4 4-6 4-8 4-9 4-9 4-12 4-12 4-14 4-15
ANÁLISIS DE CONSTRUCCIONES CON EFECTOS DINÁMICOS DE VIENTO 5-1 Acciones paralelas a la dirección del viento 5-1 Acciones perpendiculares a la dirección del viento 5-3 Criterios de confort en edificios que oscilan 5-5 Capítulo 6 ANÁLISIS DE CONSTRUCCIONES CON EFECTOS SISMICOS Definición numérica de la acción sísmica Definición mediante espectros de respuesta Definición mediante acelerogramas Métodos de análisis según INPRES – CIRSOC 103 Procedimientos con fuerzas estáticas equivalentes Métodos dinámicos Métodos dinámicos INPRES – CIRSOC 103 6-1 6-1 6-2 6-4 6-6 6-6 6-6 6-6
C. Barbat. [3] Dinámica Estructural – J. Cálculo por Ordenador – A. Bazán. LIMUSA. [2] Diseño Sismorresistente de Edificios – L. Bathe. Algor Publishing Division. Exactas Físicas y Naturales de la Universidad Nacional de Córdoba. J. 2001.Introducción . Bozzo. 1996. A. Algor Publishing Division. [7] Diseño Sísmico de Edificios – E. 1986. Barcelona. Barbat. J. Raftoyiannis. Spyrakos.H. 2000. [4] Finite Element Procedures – K.M. Prato. H. Spyrakos. Edición.Bibliografía
[1] Estructuras Sometidas a Acciones Sísmicas. CIMNE. 1994. Pitsburg. [6] Linear and Nonlinear Finite Element Analysis in Engineering Practice – C. Editorial Reverté. J.M. Pitsburg. Meli. Publicación del Departamento de Estructuras de la Facultad de Cs. R. Prentice-Hall. Massa. Barcelona. 1996. Canet – 2da.
. [5] Finite Element Modeling in Engineering Practice – C. 1997.
1 – Definiciones geométricas de un sismo
estructura epicentro
Fig. 1. con una intensidad tal que pueda causar movimientos en las estructuras. Terremotos causados por explosiones: las explosiones producidas por el hombre son capaces de generar vibraciones del terreno. Terremotos tectónicos: están causados por la rotura brusca de las capas rocosas a lo largo de superficies de fractura (fallas). con amplitudes y frecuencias dependientes del tiempo. son los más fuertes y más frecuentes. En general. tienen una intensidad pequeña y afectan a superficies limitadas.Capítulo 1 – Nociones Básicas de Sismología
NOCIONES BASICAS DE SISMOLOGIA
Causas que generan los terremotos o sismos
Los terremotos pueden definirse como movimientos de la corteza terrestre. son de baja intensidad. a su proyección sobre la superficie terrestre se le denomina epicentro. en general. el movimiento de la corteza se produce por un choque o movimiento brusco ocurrido a una cierta profundidad bajo la superficie terrestre en un punto teórico denominado foco o hipocentro. Terremotos de origen volcánico: la explosión de gases durante las erupciones volcánicas puede producir terremotos que. Las causas que los generan son variadas: Terremotos de colapso: son los originados en cavidades subterráneas por el colapso de las mismas.
por un sismógrafo Wood-Anderson estándar. las estructuras con diseño antisísmico son seriamente dañadas.45 M − 2. Una de las más difundidas es la escala de Mercalli Modificada. en la que la magnitud M mide la energía del terremoto en el foco y es el logaritmo decimal de la amplitud del movimiento sísmico. La intensidad sísmica es una medida de los efectos de los terremotos en el entorno y en particular sobre las estructuras. son los que se tienen en cuenta para la elaboración de normas para la contracción de estructuras sismoresistentes. colapso de estructuras ordinarias de mampostería. agrietamiento de las estructuras ordinarias de mampostería 4. medido en micrones a 100[km] del epicentro. etc. Existen diferentes escalas de intensidades que describen. daño en estructuras bien ejecutadas de mampostería no diseñadas para resistir fuerzas sísmicas 5.5M
Se han establecido varias relaciones empíricas entre la intensidad IMM y la magnitud M.8 + 1. por la expresión:
log E = 11. pues son los que traen consecuencias más desastrosas en las estructuras que afectan. Un sismo se caracteriza por su intensidad (parámetro subjetivo) y por su magnitud (parámetro objetivo). estas relaciones se han establecido estadísticamente. enumeramos algunas a continuación: Esteva y Rosenblueth:
I MM = 8. grandes deslizamientos del terreno destrucción casi total. colapso parcial de estructuras ordinarias de mampostería. La escala objetiva más popular es la de Ritcher. daños importantes en presas y diques.46 log R
R: distancia focal en [km]
También se ha relacionado la magnitud M con los valores máximos de las características cinemáticas del movimiento. fisuración de las estructuras de madera 2. agrietamiento de las estructuras débiles de mampostería 3. grandes masas de rocas desplazadas. daños en cimientos. debido a esto. en ergios. Algunos de los efectos sobre las estructuras en orden creciente de intensidad son: 1.16 + 1. los efectos que produce el terremoto. grietas en el terreno la mayoría de las estructuras son destruidas junto con sus cimientos.
. para cada valor que esta tome. La magnitud M está relacionada con la energía del terremoto.Capítulo 1 – Nociones Básicas de Sismología
Los sismos desde el punto de vista de la ingeniería y su caracterización
Los terremotos más importantes son los tectónicos.
ya que los valores obtenidos están afectados por las características del instrumento registrador y por las condiciones de ruido ambiental en el lugar de registro.7 e 0.25
am: aceleración máxima del terreno en [cm/s2] R: distancia focal en [km]
Esteva y Villaverde:
a m 5. que se propagan desde su foco en todas las direcciones a través de la tierra.8 M = g (R + 40)2 vm = 32 e M
(R + 25)1. Estos parámetros son relativos. diseñados para medir la aceleración. 7
am: aceleración máxima en [cm/s2] vm: velocidad máxima en [cm/s] R: distancia focal en [km]
Registro de ondas sísmicas. Los mapas de riesgo sísmico representan una síntesis de todos los datos sismológicos y geológicos de un país.8 M (R + 25)
v m = 15 e M R + 0. pero los fundamentales son:
. Parámetros utilizados y mapas de riesgo sísmico
Los terremotos son fenómenos debidos a la brusca liberación de la energía de deformación acumulada durante largos periodos de tiempo en la zona superficial de la tierra. la velocidad o el desplazamiento del movimiento sísmico.Capítulo 1 – Nociones Básicas de Sismología
am = 1080 e 0.5 M
(R + 25)1.59 M
−1. Los sismos producen ondas de varios tipos. Estos mapas se utilizan para determinar el nivel de protección que se debe alcanzar en las estructuras en cada zona de riesgo. Estas ondas son registradas mediante aparatos denominados sismógrafos.7
a m = 1230 e 0.17 e 0. Diversos aspectos brindan la subdivisión en zonas.
datos históricos 5. registro de los terremotos que ocurren en la zona 3. por ejemplo una variación de las características mecánicas del terreno. periodos de retorno (intervalo medio de tiempo en que se espera ocurran dos sismos de igual o mayor intensidad) 6. Las ondas se reflejan y se refractan cuando en su recorrido aparece una discontinuidad. datos del mecanismo focal 7. mapas de epicentros 4. En general. “predicción” estadística de las características más probables de la acción sísmica que se produzca en la zona Es importante destacar que la geología local de la zona puede modificar la propagación de las ondas sísmicas. correlación de la sismicidad de la zona analizada con la de la macrozona en la que se encuentra Estudios de Ingeniería y Sismología: 1.
. En los capítulos siguientes se ofrecen algunas aplicaciones de este tipo. ello produce cambios en la velocidad. secuencias y definiciones de acciones dadas por normas y reglamentos. ubicación de fallas 2. análisis del efecto que han producido sobre las estructuras y las personas los terremotos ocurridos en el pasado 2. el cálculo y la cuantificación de las acciones sísmicas en la estructuras se realiza en función de protocolos. Estudios sismológicos: sintetizan los parámetros que caracterizan la sismicidad de la zona: 1.Capítulo 1 – Nociones Básicas de Sismología
Estudios geológicos y geotécnicos: proporcionan datos de composición y características dinámicas de las rocas y capas de suelo que componen la corteza terrestre.
los sistemas a estudiar serán sistemas estructurales. si se conocen las influencias externas que actúan sobre el sistema.
SISTEMAS ESTRUCTURALES vibraciones
ecuaciones diferenciales que gobiernan el comportamiento de las vibraciones
permiten evaluar el comportamiento de la estructura frente a acciones dinámicas
resolución de las ecuaciones diferenciales
Definición de la acción dinámica
Una acción tiene carácter dinámico cuando su variación con el tiempo es rápida y da origen a fuerzas de inercia comparables en magnitud con las fuerzas estáticas. las variaciones en el tiempo serán vibraciones producidas por cargas dinámicas.Capítulo 2 – Conceptos Básicos de Dinámica Estructural
CONCEPTOS TRUCTURAL
DINAMICA ES-
En un sentido amplio. un sistema dinámico es aquel cuyas variables experimentan variaciones en el tiempo y.
SISTEMA DINAMICO variables con variaciones temporales
influencias externas sobre el sistema
conocidas estas acciones externas. podrá predecirse el comportamiento de este. permiten "predecir" el comportamiento de las variables temporales
En nuestro curso. Algunas fuentes importantes de vibraciones estructurales son:
Por ejemplo. 2. personas.
. pues nada nos asegura que la acción estudiada volverá a repetirse.)
La definición de estas cargas externas puede distinguirse entre: determinista y no determinista. tensiones. determinista: cuando su variación temporal es perfectamente conocida no determinista: cuando alguno o todos sus parámetros son definidos estadísticamente En nuestro curso trabajaremos con cargas definidas en forma DETERMINISTA.
Fig. Este tipo de representación es apropiado para evaluar el comportamiento de una estructura A POSTERIORI del acontecimiento que dio lugar a dicha acción. aceleraciones.Definición de la respuesta dinámica: para un punto considerado se calculan: deformaciones.Capítulo 2 – Conceptos Básicos de Dinámica Estructural
sismos viento olas y corrientes de agua explosiones e impactos cargas móviles (vehículos. El diseño de una estructura NO PUEDE encararse en base a acciones deterministas. etc.
Acciones y fuerzas dinámicas
Las acciones dinámicas definidas utilizando representaciones deterministas. ésta última denominada también estocástica o aleatoria. Respuesta dinámica → cualquier magnitud que pueda caracterizar el efecto de una carga dinámica sobre la estructura Una carga definida determinísticamente da origen a una respuesta. evaluar el comportamiento de un edificio nuevo ante el terremoto ocurrido en México en 1986 (del que se poseen registros). son funciones del tiempo cuyo valor en cada instante ES CONOCIDO. también determinista.1. etc.
velocidad de la acción ACCION . Ante una misma función de carga.Acción y respuesta determinista
ESTRUCTURA la comparación se basa en: .Acción dinámica y propiedades de la estructura
Importancia de la masa en el problema dinámico
Aunque la carga varíe con el tiempo. M
masa M rigidez K
No considerada como "carga" sino como una propiedad intrínseca de la estructura
Fig. 2.periodo propio DINAMICA de la estructura a(t) da origen a fuerzas de "inercia" comparables con las estáticas depende de K .2 . una estructura SIN MASA y una CON MASA responden de la siguiente manera:
. 2.Capítulo 2 – Conceptos Básicos de Dinámica Estructural
F(t) t3 t1 t2 M(t1)
F(t2) F(t1) F(t3)
este esquema temporal de carga debe ser perfectamente conocido
Fig. la respuesta de una estructura varía radicalmente según la masa que vibra con ella.3 .
4 .Capítulo 2 – Conceptos Básicos de Dinámica Estructural
a) Estructura SIN MASA ⇒ SIN INERCIA x(t) F(t)
rigidez: K masa: m = 0 x K x(t) = F(t) La respuesta seguirá exactamente la forma de la carga
b) Estructura CON MASA ⇒ CON INERCIA x(t) F(t)
rigidez: K masa: m
t1 x xmáx x0
m &&(t ) + k x(t ) = F (t ) x
vibraciones remanentes
Se desarrolla energía cinética.Importancia de la masa en la respuesta
. 2. que modifica la respuesta y deja vibraciones remanentes
t1 t2 t2 > t1 !!
Como veremos en los capítulos siguientes. 2.
con AMORTIGUAMIENTO la amplitud decrece en cada ciclo
x m k x0 x0
ESTRUCTURA SIN AMORTIGUAMIENTO
vibraciones libres F(t) m k F
ESTRUCTURA CON AMORTIGUAMIENTO el periodo propio permanece prácticamente constante
En general si tD >> T ⇒ no es necesario un análisis dinámico si tD ≅ T ⇒ PROBLEMA DINAMICO
Fig. En el caso de un oscilador de 1 grado de libertad.Capítulo 2 – Conceptos Básicos de Dinámica Estructural
Velocidad de reacción de una estructura
Ante una acción exterior. distintas estructuras reaccionarán de formas diferentes. Esta respuesta está íntimamente relacionada con las formas o modos de vibrar y sus correspondientes frecuencias o periodos propios.5 . No así para estructuras de múltiples GLD.Velocidad de reacción T vs. este periodo propio se obtiene fácilmente. los periodos y formas de vibrar dependen de las características geométricas y de materiales (rigidez) y de la inercia que la estructura opone al movimiento (masa). tD
será la precisión de la respuesta obtenida.Velocidad de reacción. El análisis es practicado. etc k1  k2  Fig. La definición del modelo depende del tipo de estructura analizado y pretende brindar una serie de relaciones entre acciones y respuesta que describan un modelo matemático del problema.) respecto del tiempo.
.k. En nuestro caso haremos hincapié en los métodos numéricos de análisis. Este modelo matemático puede ser resuelto mediante diversas técnicas. 2. velocidades. la definición del movimiento del terreno (en caso sísmico) tanto como de las características estructurales del mismo y de la estructura propiamente dicha. no a la propia estructura sino a un modelo mecánico de la misma.Capítulo 2 – Conceptos Básicos de Dinámica Estructural
x0 m k T1
a masa constante k 1 > k2 a rigidez constante m2 > m1 T2
m. es el resultado de "filtrar" la señal de excitación a través de la misma estructura y obtener las variaciones de las magnitudes de análisis (desplazamientos. momentos. previamente. La obtención de la respuesta requiere. obtener la respuesta dinámica de una estructura. variables → m1  m2   → T1  → T2 . etc.6 . aceleraciones. varios Ti
Modelos dinámicos característicos
Desde el punto de vista del cálculo numérico. tensiones. Según la certeza con que fueron formulados los modelos y procedimientos o algoritmos de cálculo durante el análisis.
Grados de libertad dinámicos (GLD) Son los grados de libertad que tienen asociada masa y para los cuales puede conocerse las vibraciones o movimientos a lo largo del tiempo.
.7 . a continuación.“Filtrado” de una señal sísmica Se brindan. algunas definiciones típicas del análisis estructural dinámico de una estructura: Grados de libertad (GL) Se definen como grados de libertad (GL) a los puntos de la estructura en los cuales se identifica algún desplazamiento y permiten brindar una deformada de la estructura.Capítulo 2 – Conceptos Básicos de Dinámica Estructural
ayed axed
respuesta en un punto en particular
onda refractada a2 a1 onda original a1 t onda sísmica original a3 onda reflejada
discontinuidad en el suelo
axed t respuesta
Para estos métodos modelos discretos.
. en el piso 1)
Fig. solo pueden resolverse casos muy sencillos y con poca aplicación practica.Modelización de una estructura
Métodos de modelización dinámica
Pueden distinguirse modelos dinámicos exactos y modelos dinámicos discretos. se debe tener en cuenta que la subdivisión en dominios finitos es tanto espacial (discretización estructural) como temporal (solución para instantes de tiempo determinados). por lo que a lo largo del curso profundizaremos en modelos discretos. En general.8. 2.Capítulo 2 – Conceptos Básicos de Dinámica Estructural
ESTRUCTURA REAL x2 x1
MODELO DINAMICO m2 k2 m1 k1
[M ] && + [K ] x = −[M ]a x
EXCITACION SISMICA PROCEDIMIENTOS NUMERICOS x1 menos exacta t mas exacta RESPUESTA (desplaz. para la primera clase.
Discretización espacial de las estructuras
Fundamentalmente. la diferencia con lo visto en otros cursos de análisis estructural (estático) radica en que en dinámica estructural. nos referimos a los GLD. cuando hablamos de discretizar espacialmente.9 . como ser: . Es por esto que nuestro principal método de modelización dinámica será el de las MASAS CONCENTRADAS.Capítulo 2 – Conceptos Básicos de Dinámica Estructural
Modelo dinámico EXACTO
respuesta dinámica en CADA PUNTO DE LA ESTRUCTURA en CADA INSTANTE DE TIEMPO
nº infinito de puntos espaciales nº imfinito de puntos temporales
Modelo dinámico DISCRETO nº FINITO DE PUNTOS ESPACIALES (GLD) nº FINITO DE INSTANTES DE TIEMPO donde se calcula la respuesta (para cada GLD) DISCRETIZACION ESPACIAL DISCRETIZACION TEMPORAL
Fig. las masas se encuentran más o menos concentradas en lugares conocidos. 2.método de los ELEMENTOS FINITOS
. Un modelo dinámico exacto (con infinitos GLD) acarrearía más inconvenientes en la resolución matemática que beneficios en su precisión.método de los DESPLAZAMIENTOS GENERALIZADOS . en estructuras de edificios y en la mayoría de las estructuras civiles. No obstante. existen otros. Además.
Capítulo 2 – Conceptos Básicos de Dinámica Estructural
Método de las masas concentradas
nº de GLD del modelo
nº total de componentes de desplazamiento según los cuales vibran las masas concentradas se calcula la "deformada" del modelo en cada instante tensiones y deformaciones específicas pueden conocerse mediante los procedimientos del análisis estático
se distinguen: .MODELOS DE MULTIPLES GLD
Modelos con 1 GLD:
m k (a)
m c a(t) (b) k
Fig.10 . (a) modelo conservativo. 2.
. 2.11 . (b) el mismo pórtico con la masa concentrada al nivel de la viga. (a) pórtico.
m k (a) (b) (c)
Fig.Estructuras modelizadas como un sistema de un solo grado de libertad. (b) modelo con amortiguamiento.MODELOS DE 1 GLD . (c) modelo dinámico. (c) modelo sísmico.Modelos con un solo grado de libertad.
Fig. (c) modelo sísmico. 2. 2.Capítulo 2 – Conceptos Básicos de Dinámica Estructural
Modelos con múltiples GLD:
mn kn xr kr mn kn cn xr kr cr mn kn cn xr kr cr
m2 k2 m1 k1 (a)
m2 k2 c1 m1 k1 c1 (b) a(t)
m2 k2 c1 m1 k1 c1 (c)
Fig. (b) modelo con amortiguamiento.12 -Modelos con varios grados de libertad.Estructura con masa distribuida (antena) y su modelo dinámico discreto con n grados de libertad.
Fig. 2.Estructura con dos grados de libertad: Pórtico de dos pisos y su modelo dinámico.
.14 . (a) modelo conservativo.13 .
Fig.15 .Capítulo 2 – Conceptos Básicos de Dinámica Estructural
Fig. 2.Modelo dinámico con grados de libertas de rotación.16 .Modelo dinámico de un pórtico de cortante y pórtico espacial modelizado como un sistema completo (10 grados de libertad) y simplificado ( dos grados de libertad).
Formulación de la ecuación de movimiento para un sistema de 1 GLD
Tomando el sistema de la figura 2-10.(a) fuerza aplicada.Capítulo 2 – Conceptos Básicos de Dinámica Estructural
Las ecuaciones de movimiento son las expresiones matemáticas que gobiernan la respuesta dinámica de las estructuras. donde Ep es la energía potencial. incluidas las de inercia y disipativas. aplicando el principio de D’Alembert.
Se trabaja en forma similar a lo visto en análisis estático pero incluyendo las fuerzas de inercia y disipativas.
δw i = δ w e
Principio de D’Alembert
(2. Ec es la energía cinética y Ed la correspondiente a fuerzas no conservativas. 2. Puede formularse como sigue: “un sistema dinámico esta en equilibrio cuando todas las fuerzas que actúan en el mismo.17. Pueden obtenerse a partir de cualquiera de los principios de la mecánica clásica:
. cumplen las ecuaciones de equilibrio estático en cada instante de tiempo. tendríamos:
. dπ H = 0 p − Ec ) dt + ∫ Ed dt
La primera expresión se denomina funcional de Hamilton. (b) modelo sísmico
Para el modelo (a). podemos distinguir dos casos:
Fig. La segunda expresión permite establecer el equilibrio a través de una variación funcional nula.2)
Proporciona el método más directo para obtener las ecuaciones de movimiento de un sistema dinámico.
pero de aquí en adelante ésta se encontrará implícita en toda variable temporal.5)
equilibrio en el instante “t”
F (t ) − Fi (t ) − Fa (t ) − Fe (t ) = 0 Fi (t ) + Fa (t ) + Fe (t ) = F (t )
& m && + c x + k x = F x
(2. velocidad y desplazamiento para un cierto instante “t”. etc. x(t ) → x .7) (2. F (t ) → F .6) (2.4)
Fe (t ) = k x(t )
(2.18 .Equilibrio de fuerzas para 1 GLD
Al aplicar una fuerza exterior F(t). se genera aceleración.
Fe(t) elástica
x (t ) &
Fig. &&(t ) → && x x a(t ) → a . 2.Capítulo 2 – Conceptos Básicos de Dinámica Estructural
x &&(t )
Fi(t) inercia
Fa(t) amort.
& Fa (t ) = c x(t )
(2. a causa de esto se producen fuerzas: ide inercia
Fi (t ) = m &&(t ) x
& & ej : x(t ) → x .8)
Se omite (por simplicidad de notación) la dependencia del tiempo.
Para el modelo (b) de la figura 2-17. la ecuación de movimiento queda:
& − m [a + &&] − c x − k x = 0 x & m && + c x + k x = −m a x
(2.8) es la de movimiento correspondiente a 1 GLD con carga exterior y amortiguamiento. solo que no tiene fuerza exterior aplicada y la fuerza de inercia se ve afectada por la aceleración total de la masa:
Fi = m [a + &&] x
(2.10) (2.14)
entonces.11)
La (2. 212).11) es la ecuación de movimiento para 1 GLD con aceleración de apoyo (sísmico) y amortiguamiento. Un caso general sería la inclusión de aceleración de apoyo y fuerza exterior:
& m && + c x + k x = F − m a x
(2.Capítulo 2 – Conceptos Básicos de Dinámica Estructural
La ecuación (2. nos queda un sistema de ecuaciones (vectorial)
F T − Fi − Fa − Fe = 0
(2. el planteo es similar. Planteando el equilibrio para todos los GLD.13)
para el piso (r). Aplicando el principio de D’Alembert en cada GLD (uno por piso) se obtiene:
Fr − Fir − Far − Fer = 0
Formulación de las ecuaciones de movimiento para modelos con múltiples GLD
El modelo con varios grados de libertad más sencillo es el de edificios de cortante (fig. Está basado en que las plantas son infinitamente rígidas y en que los únicos movimientos posibles de éstas son los desplazamientos horizontales.
Fi = − M {&& + J a} x & Fa = C x Fe = K x
vector de fuerzas externas vector de fuerzas de inercia vector de fuerzas disipativas vector de fuerzas elásticas
. la (2. Para un modelo sísmico..….0.0.14) puede escribirse:
& M && + C x + K x = F − M {J a} x
(2. inclusive para modelos de elementos finitos. suele sustituirse x por D para indicar que cada GLD puede sufrir desplazamientos en 3 direcciones y respectivos giros.15)
M → matriz de masas C → matriz de amortiguamiento K → matriz de rigidez JT = [1. f 2 .0.0.0.1.15) fue deducida para un modelo de edificio cortante.Capítulo 2 – Conceptos Básicos de Dinámica Estructural
#Nota: en adelante..0] etc. por ejemplo componente x:
J = [1. C y K. es una expresión ABSOLUTAMENTE GENERAL.
&& & M D + C D + K D = F − M {J a}
(2. los vectores y matrices serán representados(en general) con minúsculas y mayúsculas en negrita.17)
Nota: en este último caso.0.…. y solo varían las formas de M. el sistema (2.1.0. o modelos de elementos finitos.1. respectivamente.0.0.0. J llevará unos en las componentes a las cuales se quiera aplicar el acelerograma a.0.0.0.0.16)
Y en caso general de pórticos 3D..0.0.0.0.…1] → vector con todos sus elementos igual a uno
Si bien la (2...
F T = [ f1 ...15) se reduce a:
& M && + C x + K x = − M {J a} x
(2.0.1.
ω= m m
&& + ω 2 x = 0 x
(3.Modelo de 1 GLD. Luego. Viene expresada en radianes por segundo (1/s).2)
(3.Capítulo 3 – Respuesta de un Oscilador Simple
CAPÍTULO 3 . es inducida por algunas condiciones iniciales. 3.
. La frecuencia cíclica viene dada por
(3. La vibración.RESPUESTA DE UN OSCILADOR SIMPLE Ecuación de movimiento y equilibrio dinámico
Las características dinámicas de un oscilador de 1 GLD pueden estudiarse mediante un modelo no amortiguado con vibraciones libres. Dividiendo (3. velocidad o aceleración en el instante t = 0.1 . sean desplazamiento. cuya ecuación de movimiento es
m && + k x = 0 x
ω es la “pulsación” o “frecuencia circular” o simplemente frecuencia del modelo estudiado.4)
y se expresa en ciclos por segundo o hertz.1) por m y usando la notación:
k k . no amortiguado. durante las vibraciones no recibe ningún tipo de perturbación.1. 3. del modelo de la fig.1)
& & x(0 ) = x0
& x  A= x + 0  ω 
ω x0
x0 &
Fórmula de Geiger Sustituyendo
m= G g
G: peso de m g: aceleración de la gravedad
k 1 = g⋅ G g Gk
G = X SG : desplazamiento estático producido k
por el peso G en la dirección del grado de libertad entonces:
1 X SG
ω= g⋅
. otra característica es el período natural
1 f 2π
(3.5) (3. Los valores de A y ψ se calculan a partir de las condiciones iniciales del problema.1) o (3.Capítulo 3 – Respuesta de un Oscilador Simple
Finalmente. por ejemplo para
x(0 ) = x 0
La solución general de la (3.3) puede escribirse:
x = A sen (ω t + Ψ )
donde A es la amplitud del movimiento y ψ el ángulo de fase.
la solución de (3.Modelo de 1 GLD con amortiguamiento (vibraciones libres) Si se toma la ecuación (2.Capítulo 3 – Respuesta de un Oscilador Simple
2π ⋅ X SG g
(3. 3. la (3.11)
con XSG expresado en metros para un peso G en Newton.2 .10)
Utilizando unidades de S.12) está dada en la forma:
x = e rt
que proporciona la ecuación característica:
.I.10) queda:
T = 2.12) sin cargas ni aceleraciones de apoyo (vibraciones libres) y se divide por m se obtiene:
& && + 2 β x + ω 2 x = 0 x 2β = c m
(3.00 X SG
(3.12) (3.
Características dinámicas con amortiguamiento
El amortiguamiento puede definirse estudiando las vibraciones libres del modelo de la figura 3-2:
la cantidad (β2 .17) con i = − 1
Llamando frecuencia de vibración amortiguada a:
ωv = ω 1 − ξ 2
r1. Para este caso (subcrítico).12) es una ecuación diferencial de segundo orden.ω2 se encuentran tres tipos de amortiguamiento: β2-ω2 > 0 → SUPERCRITICO: la estructura NO VIBRA β2-ω2 = 0 → CRÍTICO: caso límite β = ω → cr = −2ω m β2-ω2 < 0 → SUBCRITICO: la estructura VIBRA con amplitud decreciente Este es el caso más frecuente en ingeniería civil.Capítulo 3 – Respuesta de un Oscilador Simple
r 2 + 2β r + ω 2 = 0
(3. 2 = − β ± iω 1 − ξ 2 = 0
(3. 2 = − β ± iω v r1.15)
Ya que (3. 2 = − β ± β 2 − ω 2 = 0
En las ecuaciones anteriores aparece la magnitud
c 2β m = cr 2ω m
(3.16) tenga raíces complejas:
r1.ω2) es negativa.16)
Según sea el radicando β2 .19) (3.18)
(3. lineal.21)
. lo que hace que (3. 2 = −ξω ± iω v
(3. por lo que enfatizaremos su estudio. homogénea a coeficientes constantes. Las soluciones de (3.15) son
3 . se obtiene
x = A e −ξωt sen (ω v t + Ψ )
sustituyendo r1 y r2 por la expresión 3.Amortiguamiento debido a fuerzas aero o hidrodinámicas .20. escribimos la solución general en la forma:
x = c1e r1t + c2 e r2t
& Fa = −c x
Las constantes A y ψ se obtienen de las condiciones iniciales del problema. El origen de las fuerzas de amortiguamiento se debe a diferentes causas: . Este modelo denominado de amortiguamiento viscoso se debe a Kelvin-Voigt y es proporcional a la velocidad. en donde la fuerza de amortiguamiento es proporcional a la fuerza normal y al coeficiente de rozamiento (hipótesis de Coulomb) . en el cálculo dinámico de estructuras.06).Vibraciones libres amortiguadas / 1 GLD
La evaluación del amortiguamiento en una estructura es un problema esencial en la dinámica estructural.Rozamiento entre superficies de deslizamiento.Debido a fricción interna del material de la estructura Generalmente. también:
c 2mω
Volviendo a la resolución de la ecuación de movimiento (3. se utiliza un modelo de gran simplicidad que caracteriza el amortiguamiento de toda la estructura. 3.Capítulo 3 – Respuesta de un Oscilador Simple
conocida como fracción de amortiguamiento crítico (en estructuras corrientes 0.12).02<ξ<0.
pero t n +1 = t n + Tv ∴
t n +1 = t n + 2π
xmáx (n ) = xmáx (n + 1)
e −ξωtn sen(ω v t n + Ψ ) e −ξωtn e
−ξω 2π
sen(ω v t n + 2π + Ψ )
xmáx (n ) 1 = −ξ 2π ω ω v xmáx (n + 1) e
con Tv = 2π ω v entonces
xmáx (n ) = xmáx (n + 1) e
e −ξωtn sen(ω v t n + Ψ )
 2π −ξω  tn +  ωv     
  2π sen ω v t n + ω v + Ψ   ωv  
Los reglamentos brindan los coeficientes de amortiguamiento para cada tipo de estructura.27)
# notar que para amortiguamientos del orden de
tomando logaritmo natural:
 x (n )  ξ 2π ω ω v ln  máx  = ln e  xmáx (n + 1)
(3. . pero puede obtenerse en forma experimental y con un método relativamente simple: Determinación práctica de ξ Para amortiguamientos bajos (del orden del 10% de crítico) la relación entre dos picos sucesivos puede aproximarse:
xmáx (n ) A e −ξωtn sen(ω v t n + Ψ ) = xmáx (n + 1) A e −ξωtn+1 sen(ω v t n+1 + Ψ )
nos interesan por ahora.Capítulo 3 – Respuesta de un Oscilador Simple
ξ = 0. (e) fuerzas impulsivas.
Fig. las periódicas y más particularmente las excitaciones armónicas ya que mediante series de Fourier cualquier excitación periódica puede llevarse a una suma de armónicas simples. De éstas. (d). (g) aceleración sísmica del terreno.
. (a) armónica. (c) cuasi periódicas.1 → ω v = ω 1 − ξ 2 ω v = ω ⋅ 0.4 .Tipos de cargas dinámicas. (f) carga dinámica general.995 → ω v ≈ ω
 x (n )  ln  máx  = 2πξ x máx (n + 1) 
(3. 3.29)
Excitación periódica
En la figura 3-4 pueden observarse diversas funciones de carga. (b) periódicas.28)
Para el caso de lecturas separadas por N ciclos:
ln[x máx (n ) x máx (n + N )] 2πN
xp en la (3.30)
Ω: frecuencia de la excitación la (3.
xh = c1 e −ξωt e iωvt + c2 e −ξωt e −iωvt
(3.30) puede escribirse también:
P0 sen(Ωt ) m
&& + 2ξω x + ω 2 x = & x
Es la solución de la ecuación diferencial homogénea.31)
La solución general de esta ecuación viene dada por
x g = xh + x p xh = c1 e r1t + c2 e r2t
Excitación armónica Si la carga es de tipo:
P = P0 sen (Ωt )
Entonces. esta última ecuación puede escribirse:
xh = e −ξωt (c1 ' senω v t + c2 ' cos ω v t )
(3.32) es la solución particular y lleva la forma
x p = A sen Ω t + B cos Ω t
(3. la ecuación de movimiento será:
& m && + c x + k x = P0 sen(Ωt ) x
Utilizando matemática para números complejos.32) (3.
c2 ' = x0 − x p0
.30) se obtienen las constantes A y B:
Denotando γ =
(3.37) (3.39)
P0 Ω(1 − γ 2 ) & x p (0 ) = ⋅ k (1 − γ 2 )2 + (2ξγ )2
P0 1−γ 2 A= ⋅ k (1 − γ 2 )2 + (2ξγ )2 B= P0 − 2ξγ ⋅ k (1 − γ 2 )2 + (2ξγ )2
Condiciones iniciales Para el caso en que
Basados en éstas y en (3.34) podemos plantear las siguientes ecuaciones:
x g (0 ) = x 0 = c 2 '+ x p 0
& & & x g (0 ) = x0 = −ξω x0 + ξω x p 0 + ω v c1 '+ x p 0
(3.32) y (3.36) (3.42)
y despejarse c1’ y c2’:
c1 ' =
& + ξω x0 − x p 0 − ξω x p 0
(3.43) (3.Capítulo 3 – Respuesta de un Oscilador Simple
Derivando y reemplazando en (3.
Pueden calcularse las partes correspondientes a la solución particular (para t = 0)
P0 − 2ξγ ⋅ 2 2 k (1 − γ ) + (2ξγ )2
x p (0 ) =
(3.41) (3.
Capítulo 3 – Respuesta de un Oscilador Simple
& 1.39 y 3.
P P(τ)
& x f = tg α
f i τ f dτ F(t)
Excitación arbitraria.Respuesta a un impulso
& & x x f = xi + && dτ
α dτ t
2.43 y 3. Esta parte de la solución se denomina respuesta en régimen. && = P (τ ) m x
P(τ ) dτ m
& xf =
.5 .Se computan c1’ y c2’ con las 3. por lo que se la denomina transitoria y desaparece antes o después según sea el valor de ξ. Integral de Duhamel
Respuesta a un impulso elemental Dada una carga impulsiva como la de la fig. 3.Se calculan x p 0 y x p 0 con las 3. puede calcularse la velocidad y posición al finalizar el impulso y a partir de ese instante se trata de un problema de vibraciones libres. La estructura sigue vibrando con una frecuencia prácticamente igual a la de la excitación Ω.Para obtener velocidades y aceleraciones se deriva la 3.32 y se tiene la respuesta en desplazamientos 4.32 Nota: la parte de la solución correspondiente a la ecuación homogénea incluye el coeficiente e-ξωt que es una función decreciente con el tiempo.44 3.Se reemplaza todo en 3. para condiciones iniciales nulas. 3-4(d).
47).46) o (3. para un tiempo infinitesimal. luego de actuar un impulso la masa queda con velocidad pero prácticamente no se desplaza. entonces (b) es despreciable frente a (a). la solución se obtiene con el mismo procedimiento dado para las cargas armónicas. Reemplazando ésta en la (3.46) y (3.Capítulo 3 – Respuesta de un Oscilador Simple
1 2 & x f = xi + x dτ + &&(dτ ) x 2 1 2 x f = &&(dτ ) x 2
Si dτ es muy pequeño. es decir. se transforma en:
x(t ) = −
a(τ )e ωv ∫ 0
−ξω (t −τ )
sen[ω v (t − τ )]dτ
(3. La solución dada por esta superposición es conocida como Integral de Duhamel
1 t −ξω (t −τ ) x(t ) = sen[ω v (t − τ )]dτ ∫ P(τ )e mω v 0
En el caso de excitación sísmica. así como el máximo factor dinámico (FAD) y el instante de máxima respuesta.
Soluciones a la integral de Duhamel: (a) Para ciertas funciones de P(τ) se encuentran tabuladas las primitivas de (3.47)
& Para el caso en que el problema tenga condiciones iniciales x0 ≠ 0 y x 0 ≠ 0 . La velocidad dada por (a) es la condición inicial para las vibraciones libres siguientes. pueden superponerse las respuestas de una sucesión de impulsos infinitesimales hasta el tiempo genérico “t”.
.47).45)
Si el sistema es lineal. se obtiene:
 P (τ )  dx(t ) = e −ξω (t −τ )  dτ sen (ω v (t − τ ))  mω v 
(3. solo que la solución particular vendrá dada por (3.23).
Factor de amplificación dinámica
Es el coeficiente entre el máximo desplazamiento dinámico y el que produciría la carga en forma estática:
FAD = x(t )máx x est
La solución a esta ecuación viene dada por la ec.3. expresado en función del periodo propio. 154 / ref. Para un sistema de 1 GLD:
& m && + c x + k x = −m a (t ) x
(2. [1]. Esto es muy útil para el diseño de estructuras donde solo los máximos son necesarios. Derivando una y dos veces se obtienen las expresiones para velocidad y aceleración:
& S vr = x(t ) máx
Espectros sísmicos de respuesta
Puede definirse el espectro d respuesta (para un acelerograma dado) como los máximos valores de la respuesta de un sistema.Capítulo 3 – Respuesta de un Oscilador Simple
(b) Para funciones aproximadas por “trozos rectos”. 5)
. Reemplazando S d = x(t ) máx Puede obtenerse. etc.47) y el máximo valor absoluto será Sd. para 1 GLD.
r S a = &&(t ) máx x
(Ver ref. (3.4 / pág. Simpson. la expresión exacta del espectro de desplazamientos. cap. también se pueden encontrar soluciones “exactas” a dicha integral (punto 5. [1]) (c) En el caso mas general se puede “integrar numéricamente” utilizando alguno de los métodos conocidos (Trapecios.
6 .2.002 < ξ < 0.Espectro sísmico de respuesta para desplazamientos.50)
. Hipótesis simplificativas: iiien estructuras civiles. por una tipo seno
sd (ω .ξ ) = − ∫ a(τ )e −ξω (t −τ ) sen[ω (t − τ )]dτ
0 máx
(3. 3. Estas nuevas ecuaciones definen los “Seudo espectros de respuesta”.47) puede reemplazarse.ξ ) = −
a(τ )e ω∫
sen[ω (t − τ )]dτ
(3. sin que ello implique grandes variaciones en los valores máximos. 0.Capítulo 3 – Respuesta de un Oscilador Simple
No obstante pueden hacerse ciertas hipótesis que simplifican las expresiones obtenidas.49)
sv (ω . por lo que puede reemplazarse ωv por ω la función coseno que aparece en el espectro de velocidades (al derivar 3.
3-9.53)
Por lo que un espectro completo puede representarse en un gráfico tipo logarítmico como el de la figura 3-7. los espectros de diseño del UBC y los espectros elásticos para los registros de EL Centro y de la presa de Pacoima.Capítulo 3 – Respuesta de un Oscilador Simple
sa (ω .ξ ) = ω ∫ a(τ )e −ξω (t −τ ) sen[ω (t − τ )]dτ
(3.52) (3. 3.
Fig. Es importante notar que:
sv = ω S d sa = ω 2 S d
(3.Comparación entre los espectros suavizados de Newmark y Hall.
definen los seudo espectros para desplazamientos. velocidades y aceleraciones. Es por eso que solo aparecen gráficos del tipo de la fig. La normativa argentina INPRES-CIRSOC 103 establece la acción sísmica como una función del seudo espectro de respuesta de aceleraciones Sa.7 . que varían según el tipo de suelo y la zona sísmica considerada.
.9 .Espectro de diseño suavizado de Seed e Idriss (1982).
Fig.8 . (a) corresponde al acelerograma (1).Seudo espectro y espectro de aceleraciones. 3. (b) corresponde al acelerograma (2). 3.
55) se obtiene la expresión de la velocidad
& & x x(τ ) = xi + &&iτ + (&&i +1 − &&i )∫ f (τ )dτ x x
(3.11) utilizando ecuaciones en diferencias que permitan obtener x .56)
&&i +1 x
∆t = t i +1 − t i
&&i x
τ = t − t i . pero solo nos dedicaremos en profundidad al de Newmark. x y && en un instante t+∆t en función. únicamente. de los valores en el instante t.55)
integrando (3.
Método de Newmark (1 GLD)
(3.8) y (2.3. dado que es uno de los más difundidos en software de análisis dinámico mediante FEM. τ ≤ ∆t
puede suponerse que el valor de la aceleración de respuesta en un instante τ se expresa como:
ti τ
&&(τ ) = &&i + f (τ )(&&i +1 − &&i ) x x x x
0 para τ = 0 f(τ) = 1 para τ = ∆t
(3. .54)
con F = f (t ) + [− m a(t )]
Fig. Existen varios métodos para plantear las ecuaciones en diferencias antes mencionadas.Capítulo 3 – Respuesta de un Oscilador Simple
Integración numérica de la ecuación de movimiento
Se trata de resolver en forma discreta la ecuación de movimiento para 1 GLD dada por & x (2.
Para calcular los desplazamientos.64)
Las ecuaciones (3.64) juntamente con la ecuación de movimiento (3.Capítulo 3 – Respuesta de un Oscilador Simple
g (τ ) = ∫ f (τ )dτ
(3.60) (3.59)
que para τ=∆t resulta
& & x xi +1 = xi + &&i ∆t + (&&i +1 − &&i )γ ∆t x x & & xi +1 = xi + [(1 − γ ) &&i + γ &&i +1 ]∆t x x
. propuesta por Newmark
 1   & xi +1 = xi + xiτ +  − β  &&i + β &&i +1  ∆t 2 x x   2 
se obtiene la relación final en diferencias.61) y (3.54) permiten obtener aceleraciones.57) (3.56) queda:
& & x x(τ ) = xi + &&iτ + (&&i +1 − &&i )g (τ ) x x
(3. velocidades y desplazamientos en un instante t+∆t con solo conocer estas magnitudes en el instante t.58)
γ ∆t =
∫ f (τ )dτ
La expresión (3.59) y se obtiene
& x(τ ) = xi + xiτ + &&i x
+ (&&i +1 − &&i )∫ g (τ )dτ x x
particularizando para τ=∆t y llamando
β ∆t 2 = ∫ g (τ )dτ
(3. se integra (3.
formar el término de “carga efectiva” en t+∆t
& & ri +1 = F i +1 + m(a0 xi + a2 xi + a3 &&i ) + c(a1 xi + a 4 xi + a5 &&i ) x x ˆ
1 β ∆t
1 −1 .25
γ β ∆t
. agrupando los valores que permanecen constantes en cantidades denominadas por a0. se obtiene el algoritmo de integración paso a paso.resolver el desplazamiento en t+∆t
. 2β
γ −1 . [2] y [4]) pero pueden tomarse:
γ = 0.54) x
 ∆t  γ  − 2  2 β  
a 6 = ∆t (1 − γ )
a 7 = γ ∆t
5. x0 y &&0 utilizando.25
Desarrollo y forma operativa La ecuación de movimiento (3. a7.65)
Reemplazando en ésta los valores en diferencias dados por (3. c. k
& 2. A) cálculos iniciales 1.64). ref.54) valuada en t=ti+1 toma la forma:
& m &&i +1 + c xi +1 + k xi +1 = F i +1 x
(3.formar la “rigidez efectiva”
ˆ k = k + a 0 m + a1 c
B) para cada paso de integración 1. de ser necesario.determinar las propiedades de la estructura: m. γ .…. γ ∆t 2 a1 = a4 =
∆t ≤ T/10 .61) y (3. [1].Capítulo 3 – Respuesta de un Oscilador Simple
Los parámetros γ y β surgen de un análisis de estabilidad numérica que escapa a este curso (ver.calcular las constantes:
a0 = a3 = 1 .5 .seleccionar el paso de tiempo y los parámetros γ y β ∆t .5 β = 0. β 4. β=0.inicializar x 0 . γ=0. la ecuación (3.
calcular aceleraciones y velocidades en t+∆t
&&i +1 = a0 ( x+1i − xi ) − a2 xi − a3 &&i & x x & & xi +1 = xi − a6 &&i + a7 &&i +1 x x
.Capítulo 3 – Respuesta de un Oscilador Simple
ˆ k ⋅ x i +1 = r1+1 ˆ
utilizando el principio de D’Alembert.1)
En particular. no es relevante para estructuras civiles corrientes. el sistema de ecuaciones diferenciales es del tipo
&& & M D + C D + K D = F (t )
(4.discretización temporal (por ej.
Características dinámicas El sistema que gobierna las vibraciones libres en un sistema de múltiples GLD es
&& & M D+CD+ K D=0
(4.Capítulo 4 – Respuesta Dinámica de una Estructura con Múltiples Grados de Libertad
CAPÍTULO 4 . Newmark) En el Capítulo 3 se estudio el modelo discreto más sencillo (1 GLD) mientras que en la UT2 se plantearon (sin resolver) las ecuaciones de equilibrio dinámico para sistemas de múltiples GLD. masas concentradas) . por lo que seguiremos el análisis con el sistema simplificado:
. En general.3)
Como ya se vio.2)
Recordar que J es un vector con 1 en la posición del GLD en el que actúa la aceleración de apoyo a(t).discretización espacial (por ej. para el caso sísmico:
&& & M D + C D + K D = − M {J a}
(4. el cambio de frecuencia propia debido a considerar o no el amortiguamiento.RESPUESTA DINÁMICA DE UNA ESTRUCTURA CON MÚLTIPLES GRADOS DE LIBERTAD Ecuaciones de movimiento y equilibrio dinámico
Como se vio en el Capítulo 2 un modelo estructural dinámico esta dado por .
8 representa un problema de autovectores y autovalores. en donde:
K − ω2M = 0
(4..6)
reemplazada en (4.10)
. quedando finalmente
{K − ω
M }A = 0
Este determinante puede desarrollarse en la forma polinómica:
ω 2 n + α 1 ω 2 n − 2 + α 2 ω 2 n − 4 + .8)
Nos interesan las soluciones de A distintas de la trivial.α n −1 ω 2 + α n = 0
obteniéndose la ecuación característica. ω≠0.Capítulo 4 – Respuesta Dinámica de una Estructura con Múltiples Grados de Libertad
&& M D+ K D=0
Una solución a este sistema puede ser de la forma
D = A sen(ω t + Ψ )
(4.7).5)
A: vector que contiene las amplitudes de las vibraciones Ψ: ángulo de fase
&& D(t ) = −ω 2 A sen(ω t + Ψ )
(4. La ecuación (su resolución) 4. por lo que
sen (ω t + Ψ ) ≠ 0
y puede eliminarse de (4.7)
para que haya vibraciones.
(4.4) queda:
− M ω 2 A sen(ω t + Ψ ) + K A sen(ω t + Ψ ) = 0
M }A sen(ω t + Ψ ) = 0
Según la matriz de masas
ϕ i = Ai (M i *)
(4. Reemplazando cada ωi en la (4.13)
. Por esto es conveniente “normalizar” los distintos Ai. Normalización de los modos 1.Capítulo 4 – Respuesta Dinámica de una Estructura con Múltiples Grados de Libertad
Cuando K y M son definidas positivas (caso usual en estructuras civiles). éste se denomina “periodo fundamental” del sistema. siendo “n” el número de GLD de la estructura. que es la frecuencia de menor valor. Es importante notar que (4. Por la máxima componente
 a i1  M    dado Ai = aimáx    M  ain   
puede dividirse cada elemento por el aimáx
ϕi1  M    ϕ i = 1  modo i normalizado   M  ϕin   
Ai aimáx
(4. la forma modal será diferente.8) se obtiene el correspondiente Ai que se denomina vector de “forma” modal o simplemente “modo”.11)
Siendo T1 el período correspondiente ω1.8) se cumple para cualquier múltiplo de Ai. de la ecuación característica se obtienen “n” soluciones positivas ωi2 y en consecuencia “n” valores ωi..12)
(4. Este vector “modo” contiene la forma que tomara la estructura en cada vibración. por lo tanto inferimos que no interesa la magnitud de las componentes de Ai sino la relación existente entre ellas.. Para cada frecuencia vibratoria. Los valores ωi se denominan “frecuencias propias” o “pulsaciones” de la estructura y los “n” periodos propios se calculan
(4. 4.15)
Obtención de los grados de libertad dinámicos
No siempre todos los grados de libertad estáticos que definen el comportamiento de la estructura tienen asociada “masa”. es decir que no son necesarios en el análisis dinámico.Capítulo 4 – Respuesta Dinámica de una Estructura con Múltiples Grados de Libertad
M i * = Ai [M ]Ai
(4. observando el vector desplazamientos estáticos Ds podemos determinar que parte queremos eliminar y cuales grados de libertad deseamos conservar como GLD. Lo mismo se ve en la figura 2-16.
. La forma más simple de reducir en número de grados de libertad sin perder precisión es mediante la: Condensación estática de la matriz de rigidez Dada una estructura y su modelo estático puede subdividirse la matriz de rigidez (así como el vector desplazamientos) de manera de separar las ecuaciones que tienen asociada masa (GLD) y las que no la tienen (GLE).Condensación estática
Para el ejemplo de la figura 4-1. Por ejemplo los pórticos de las figuras 2-15 (f) y (h) que son analizados con modelos dinámicos de menor orden que los respectivos estáticos (GLE). tenemos que.14)
Esta normalización permite asegurar la condición
ϕ i T [M ]ϕ i = 1
Expresión que será de suma importancia más adelante.
d3 kv k2 d1 kv k1 k1 k2
d4 m2 d2 m1
d4 m d2 m
modelo estático 4x4
modelo dinámico 2x2
KD=F
donde K es la matriz de rigidez condensada y vale
La segunda ecuación (vectorial) de (4.20)
-K 12 K 22 K 21 D = F
(4.16) tendremos
K 11 D − K 12 K 22 K 21 D = F
(4. entonces.19)
Desarrollando ahora la primera ecuación de (4.18)
De = − K 22 K 21 D
El sistema de ecuaciones estático puede escribirse
 K 11 K  21 K 12   D   F  = K 22   De   Fe     
Como Fe representa las fuerzas en los grados de libertad a eliminar
(4.16) puede.21) (4. expresarse como
K 21 D + K 22 De = 0
 D  =   De
D: vector con los GLD De: grados de libertad a eliminar
K s Ds = Fs
puede ser de más sencilla aplicación para procedimientos manuales. Este método. Para nuestro ejemplo:
δ m2 1 δ
δ F =  11 δ 12 K = F −1
δ  δ  
Fig. si bien menos formal.24)
Otra forma de obtener la matriz de rigidez del modelo dinámico es invirtiendo la matriz de flexibilidad de los GLE que tienen asociada masa. 4-1 estará representado matemáticamente por
&& K D+ M D=0 ↓ ↓
2× 2 2×2
d  D =  2 d 4 
(4. 4-2 .Obtención de K mediante F-1
.Capítulo 4 – Respuesta Dinámica de una Estructura con Múltiples Grados de Libertad
K = K 11-K 12 K 22 K 21
Finalmente el sistema dinámico de la fig.
Esta condición para los modos puede extenderse a
ϕi M ϕ j = 0
(4.ϕ Kϕ ] O Ω=  1 2 n   0 ωN   
(4. Φ = [ϕ .Capítulo 4 – Respuesta Dinámica de una Estructura con Múltiples Grados de Libertad
(4. Condiciones de ortogonalidad Las frecuencias y modos propios pueden ordenarse en matrices denominadas espectrales y modales. [3].28)
Conviene definir matrices de amortiguamiento ortogonales. primero incursionar un poco en las condiciones de ortogonalidad de los modos propios. [4]) que
ϕi ϕ j = 0
(4. que son.26)
Puede demostrarse (ver ref.1)
solo nos resta conocer (para encarar la resolución) la matriz C. respectivamente:
0 ω1  . [1]. pero para esto debemos. [2]. puesto que en el apartado anterior se mostraron dos maneras de obtener K y sabemos (porque usamos un modelo discreto de masas concentradas) que:
L 0 m1   O   M = M mi M    O   0 L mn   
Capítulo 4 – Respuesta Dinámica de una Estructura con Múltiples Grados de Libertad
ϕi K ϕ j = 0
Si los vectores ϕi fueron normalizados según (4.13) las condiciones de ortogonalidad y normalización pueden expresarse como una única condición de ortonormalidad
ΦT M Φ = I Φ T K Φ = K*
(4.30) (4.31)
donde I es la matriz identidad y K* es una matriz diagonal
k1 *    O    ; ki * = ϕ i T K ϕ i K *=  ki *   O    k n *  
Matrices de amortiguamiento ortogonales
Una de las hipótesis para lograr una representación numérica del amortiguamiento de una estructura está dada por suponer que existe un mecanismo de disipación uniforme de energía. En tal caso puede desarrollarse una matriz de amortiguamiento que cumpla las condiciones de ortogonalidad respecto a la matriz modal. El “amortiguamiento proporcional” permite definir una matriz que sea proporcional a la de las masas, a la de rigidez o a ambas.
C = α1 M + α 2 K
La condición de ortogonalidad:
ϕi C ϕ j = 0
ϕ i C ϕ i = 2ωi ξ i
Los coeficientes α1 y α2 se calculan a partir de las ecuaciones anteriores (4.35) y (4.33)
ϕ i [α1 M + α 2 K ]ϕ i = 2ωi ξ i
α1ϕ i M ϕ i + α 2ϕ i K ϕ i = 2ωi ξ i
α 1 + α 2ω i2 = 2ω i ξ i
Recordar que si la normalización de los autovectores ϕi fue hecha según (4.13) y (4.14) el producto
ϕ i K ϕ i = ωi
Esto se debe a que la (4.8) puede expresarse para A = ϕi y premultiplicarse por ϕiT
ϕ i K-ωi M ϕ i = 0
ϕ i K ϕ i − ωi ϕ i M ϕ i = 0
Determinación práctica de modos y frecuencias
Este es, por sí mismo, uno de los problemas más complejos de la dinámica estructural. Generalmente no es necesaria la resolución de todas las frecuencias y sus modos, sino que solo interesan las primeras “q” que representan las posibilidades ciertas de vibrar pues necesitan menos energía de excitación. Cuando el problema es de pocos GLD es posible resolver el determinante de (4.9) o (4.10) obteniéndose las “n” ωi y luego los “n” ϕi, donde “n” es el número de GLD. Para problemas con muchos GLD se utilizan técnicas numéricas (aproximadas) para obtener los primeros “q” pares de (ωi; ϕi). En la mayoría de los casos prácticos de ingeniería q << n. En las referencias bibliográficas pueden encontrarse descripciones detalladas de varios métodos. Nosotros desarrollaremos solo uno que es relativamente simple y preciso.
Método de Stodola-Vianello
Partiendo de la (4.8) y recordando que se cumple para cualquier múltiplo del autovector A podemos escribir
K A − ω2M A = 0 K A = ω2M A
Suponiendo un vector inicial A0 conocido (cercano a A)
K A 1 = ω 2 M A0 A1 ≈ K -1 M A0
que lleva a la fórmula de recurrencia
A i + 1 ≈ K -1 M A i
Nota: una descripción detallada de este método puede encontrarse en [4], parte 11.2 “Vector Iteration Methods”, así como prueba de la convergencia del método.
Procedimiento de cálculo 1- Se propone un vector A0 inicial. Conviene que los valores correspondan a una deformada suave, no obstante esto solo acelera la convergencia. 2- Se calcula
A1 = K -1 M A0
3- Se normaliza A1
A1 = A1 {A1T M A1 }
i=1 4- Se mejora la solución
A i + 1 = K -1 M A i
5- Se normaliza la nueva solución
Ai +1 = Ai +1 {Ai +1,T M Ai +1 }
6- Se calcula la frecuencia correspondiente
 Ai +1 T K Ai +1   ω i +1 =   Ai +1 T M Ai +1   
7- Se controla la convergencia
ω i +1 − ω i NO ≤ TOL → se vuelve al paso 4 i +1 ω
ϕ1 = Ai +1 ; ω 1 = ω i +1
Se obtuvo el primer autovalor y su correspondiente autovector. Siempre:
Ai → ϕ1
ω i → ω1
Obtención de modos y frecuencias superiores El procedimiento es el mismo que para el modo 1, solamente debemos garantizar que el vector de “arranque” sea ortogonal a ϕ1 Se supone que un vector inicial
A20 = ∑ q j ϕ j
que para hacerlo ortogonal a ϕi le “restamos” la componente q1 ϕ1
A*20 = A20 − q1 ϕ 1
Para determinar q1, premultiplicamos la (4.54) por [ϕ1T M] y aplicando las propiedades de ortogonalidad nos queda
0 T M A2 = ϕ 1 M
]∑ q ϕ
ϕ1 M A20 q1 = T ϕ1 M ϕ1
ϕ i T M Ai0+1 qi = T ϕi M ϕi
y para los vectores de arranque (ortogonales a los previamente calculados)
A*i0 = Ai0 − ∑ q j ϕ j
todo lo expuesto es valido para el caso general. La ecuación que gobierna el comportamiento de una estructura de múltiples GLD es la
(4. Es por esto que pueden formar una base completa para el espacio de los desplazamientos estructurales.4)
al cual.60)
Las vibraciones libres no amortiguadas se estudian mediante
(4. Se estima Ai0k 2.8)
Recordemos que los autovectores ϕi son ortogonales respecto a las matrices de masa y rigidez. Se sigue con el método de Stodola propiamente dicho
Resolución de las ecuaciones de movimiento en estructuras con múltiples GLD
Descomposición y superposición modal
Si bien desarrollaremos el método para el caso sísmico. es posible escribir
D = ∑ ϕ i yi (t )
Procedimiento: 1. Se calcula qi con (4. según lo expuesto en puntos anteriores. es decir. le corresponden “n” pares de frecuencias y “modos” que son solución del sistema de ecuaciones algebraicas
[K − ω M ]ϕ = 0
(4. Se calcula Ai*0k con (4.58) 3.59) 4.
Donde yi(t) es un escalar función del tiempo a determinar. como es habitual.60) en (4.61)
Premultiplicando (4.65)
ϕ j C ϕ j = C j * = 2ω j ξ j ϕ j K ϕ j = K j * = ωj
Entonces la (4.4)
ϕ C ∑ ϕ i = ϕ T Cϕ j = C j * j
T j i =1
(4. Sustituyendo (4.61) por un ϕiT cualquiera.39)
Y dividiendo ambos miembros por M j * = ϕ j M ϕ j
.62) (4.63)
ϕ T K ∑ ϕ i = ϕ T Kϕ j = K j * j j
Y si.61) premultiplicada por un ϕjT cualquiera queda como una ecuación de 1 GLD
& M j * && j (t ) + C j * y j (t ) + K j * y j (t ) = −ϕ T M J a(t ) y j
(4. se cumple que:
ϕ M ∑ ϕ i = ϕ T Mϕ j = M j * j
(4. se está frente a matrices de amortiguamiento proporcionales y ortogonales (ver punto 4.2) obtenemos
& M ∑ ϕ i &&i (t ) + C ∑ ϕ i yi (t ) + K ∑ ϕ i yi (t ) = − M J a y
(4.35) (4. llamado “respuesta generalizada”.
…q 3.68) i=1.
. los modos bajos son los que contienen menos energía elástica de deformación y por ende los que más contribuyen a la respuesta estructural. básicamente en obtener la solución en una cierta cantidad (discreta) de pasos de tiempo.60). el sistema (4. queda en “n” ecuaciones diferenciales de 1 GLD.68)
Esta ecuación puede resolverse con cualquiera de los métodos vistos en la UT3.67)
Finalmente. del tipo:
&&i (t ) + 2ξ i ωi yi (t ) + ωi2 yi (t ) = −Qi a (t ) & y
(4. Es por esto que también es llamada integración “paso a paso”.69) y sus derivadas 4.Se obtienen las fuerzas elásticas. al ser “proyectado” según la base formada por la matriz modal Φ = [ϕ1.66)
ϕ MJ Q j = jT ϕj Mϕj
se lo denomina “coeficiente de participación modal”.2). q < n
(4.Se resuelven las “q” ecuaciones (4.… ϕn].69)
Procedimiento de cálculo 1.
(4. En general. Usualmente
D(t ) = ∑ ϕ i yi (t ) . puede obtenerse la solución estructural mediante (4. de inercia y de amortiguamiento para cada instante
Integración directa de las ecuaciones de movimiento
Como su nombre lo indica.Se obtiene la historia en el tiempo según (4.Capítulo 4 – Respuesta Dinámica de una Estructura con Múltiples Grados de Libertad
ϕ M J a(t ) && j (t ) + 2ξ j ω j y j (t ) + ω y j (t ) = − j T & y ϕj Mϕj
(4. este método no requiere ninguna transformación previa de las ecuaciones de movimiento.Se determinan los “q” primeros modos y frecuencias 2. Consiste. Una vez obtenidos los “n” yi(t).
Inicializar D0 .73) (4.5
. existen numerosas variantes y algoritmos para la integración numérica de las ecuaciones de movimiento. Calcular las constantes a0. pero nuestros desarrollos estarán dados en el método de Newmark.72)
*1: el término F i +1 = F i + {− M J a(t )} se refiere a fuerzas generalizadas 3. solo que en vez de tratarse de una ecuación de 1 GLD se trata de un sistema de “n” GLD.
Respuesta máxima utilizando espectros de respuesta
Las ecuaciones desacopladas (4.70)
B) Para cada paso de tiempo 1.
γ=0.Capítulo 4 – Respuesta Dinámica de una Estructura con Múltiples Grados de Libertad
Como ya adelantamos. Formar el término de “carga efectiva” en t+∆t (*1)
& && ˆ ri +1 = F i +1 + M (a0 Di + a 2 Di + a3 Di & && + C (a1 Di + a4 Di + a5 Di )
(4.…a7 5. Seleccionar el paso de tiempo y los parámetros γ y β
∆t ≤ Tmín/10
. Calcular aceleraciones y velocidades en t+∆t && & && Di +1 = a0 ( Di +1 − Di ) − a 2 Di − a3 Di
& & && && Di +1 = Di − a6 Di + a7 Di +1
(4.72) lo que puede
hacer excesivo el costo de cálculo y almacenamiento.74)
# Notar que para cada paso de tiempo se debe resolver el sistema (4. Ensamblar las matrices M. D0 y D0 3. En esencia es igual al planteo dado en el Capítulo 3. Procedimiento de cálculo A) Iniciales 1.71)
2. Resolver el desplazamiento en t+∆t
ˆ ˆ K Di +1 = ri +1
(4.68) pueden resolverse utilizando los espectros de respuesta: la máxima aceleración será
. Formar la matriz de rigidez efectiva:
ˆ K = K + a0 M + a1C
(4. C y K & && 2.
j Dmáx = A j
(4. aceleraciones y esfuerzos) como combinación de los máximos modales es:
Rmáx =
2 j máx
(4. la respuesta total máxima no es la suma de los máximos de cada modo!!
i Dmáx ≠ ∑ Dmáx i =1 q
entonces. el máximo desplazamiento es
y i (t ) máx =
(4. Puesto que el máximo en cada grado de libertad no se produce en el mismo instante.80)
denominada SRSS
. entonces calcular los máximos desplazamientos (en todos los GLD) para el modo “j”:
j Dmáx = ϕ j y j (t ) máx = ϕ j
ω2 j
&&i (t ) máx = Qi (S a )i y
Una forma muy utilizada (y precisa) de evaluar la respuesta máxima (desplazamientos. velocidades.76)
Podemos.78)
Donde Aj es el vector de coeficientes de participación modal correspondientes al modo de vibración “j”.
Acciones paralelas a la dirección del viento
(procedimiento basado en el “factor de ráfaga”)
X max ( z ) = x( z ) + xmax ( z )
x( z ) : desplazamiento medio x max ( z ) : desplazamiento fluctuante debido a la turbulencia variable con el tiempo
x max ( z )
. pero para ciertas estructuras: .acusada flexibilidad (poca rigidez) . CIRSOC 102-1 → procedimientos simplificados Validez: a) construcciones prismáticas o cuasiprismáticas b) primer modo dominante y de forma aproximadamente lineal c) periodo fundamental T1 > 1 seg d) amortiguamiento ξ < 0.construcciones livianas (poca masa) las acciones perpendiculares de “deriva” pueden tener efectos significativos que deben ser contemplados en el análisis.bajo amortiguamiento . perpendicular a la dirección del flujo
Generalmente son dominantes las fuerzas de empuje.Capítulo 5 – Análisis de Construcciones con Efecto Dinámico de Viento
T(z): fuerza aerodinámica de empuje. paralela a la dirección del flujo
T(z) L(z)
L(z): fuerza aerodinámica de deriva.
a) Comportamiento elástico lineal de la estructura. c) La contribución de los modos superiores al primero en respuesta se considera despreciable. por lo que G(z) = G = cte. b) El modo fundamental de vibración es una función lineal de la altura.2.1      
V0: velocidad básica de diseño (m/seg) z0.1: z0.2. CIRSOC 102-1) q0: presión dinámica básica (art. 20.Capítulo 5 – Análisis de Construcciones con Efecto Dinámico de Viento
Se define como factor de ráfaga
G(z ) = 1 + x max ( z ) x( z )
X max ( z ) = G ( z ) ⋅ x( z )
Hipótesis básicas para poder reemplazar la acción dinámica del viento turbulento por un procedimiento estático equivalente.i para rugosidad I (referencia) Procedimiento de cálculo I) Presiones La presión dinámica que incluye el efecto de la turbulencia del viento se determina mediante
q' z = G ⋅ c z ⋅ c 2 ⋅ q 0
G: factor de ráfaga cz: variación por rugosidad y altura (art.0706   0 . d) La velocidad media del viento es promediada sobre intervalos de una hora.i     z  0.2. CIRSOC 102)
.i V ( z ) = V0 ⋅  10 ln z  0.i: parámetro que depende de la rugosidad i z0.i  ⋅   z 0. e) La variación de la velocidad media del viento varía según
 z ln z  0.4. CIRSOC 102) c2: factor por cambio de tiempo en velocidad media (Tabla 3/pág. 5. 5.3 .
. (z ) → valor medio cuadrático de las aceleraciones (pág.Capítulo 5 – Análisis de Construcciones con Efecto Dinámico de Viento
G = 1 + 1. X max
q' z o qz
a.Verificación de confort con gráficos y
.Dimensionado estructural b.. X max y T c.. CIRSOC 102-1
Aceleraciones (paralela al viento)
X max ( z ) = K ⋅ σ . ( z )
K → fig. .Verificación de deformaciones admisibles H H < X max < 500 350
Acciones perpendiculares a la dirección del viento
I) Resonancia torbellinos de Bèrnard .
tablas en función de X max .
X max . 21)
Verificaciones (paralelas al viento)
.234 ⋅ K ⋅
B+r → J 
se calcula mediante tablas y ábacos que utilizan una serie de parámetros auxiliares
En cada nivel se comparará: q’z > qz NO → se adopta qz para los esfuerzos SI → se adopta q’z para los esfuerzos
CIRSOC 102 → q z = c z ⋅ c d ⋅ q 0 II)
. . 13.Karman
V h T(z) L(z)
⋅ q cr ⋅
z ⋅d h
Vcr = Vcr ⊥
Fuerzas de empuje: (a la velocidad crítica)
se admite distribución uniforme
T ( z ) = T z = 0.25 a 0.000613 ⋅ Vcr  2  m 
L( z ) =
0.Capítulo 5 – Análisis de Construcciones con Efecto Dinámico de Viento
La velocidad crítica del viento que produce el fenómeno de resonancia es
Vcr = d S ⋅T
d: ancho de la superficie maestra (puede ser variable)
cilindros ≈ 0. De lo contrario: Fuerzas de deriva: (a la velocidad crítica) ξ: fracción del amortiguamiento
2  kN  q cr = 0. entonces puede prescindirse del cálculo de la resonancia.8 ⋅ C E ⋅ G ⋅ q cr ⋅ d
CE: coeficiente global de empuje CIRSOC 102 G: factor de ráfaga correspondiente a Vcr qcr → con Vcr// Las fuerzas L(z) y Tz obtenidas para la velocidad crítica (correspondientes al período perpendicular y paralelo respectivamente) se suman de la siguiente manera:
F ( z ) = L( z ) + T z2
y se debe comparar con las correspondientes obtenidas con q’z o qz dadas para la velocidad de diseño.
.20 S: nº de Strouhal  prismas ≈ 0.30 Vcr // → T // T: período propio  Vcr ⊥ → T⊥
Si Vcr > 25 m/seg.
Las curvas de la fig. En consecuencia. otras persones sienten placer por el movimiento y la experiencia poco usual de estar en una construcción que oscila. Adicionalmente.
. Los factores que pueden producir vibraciones en un edificio son numerosos. La respuesta humana al movimiento oscilatorio de las construcciones abarca una extensa gama de reacciones. lo cual incide negativamente en el valor de una propiedad y su rentabilidad. rotura de vidrios de ventanas y caída de revestimientos.1 grafican los límites de confort obtenidos del análisis de un gran número de edificios altos.3. un edificio resultaría demasiado costoso si se construye o equipa de modo que pueda soportar sin movimientos perceptibles una tormenta con vientos huracanados o un fuerte sismo. los movimientos son casi inevitables y el problema del diseño consiste en mantenerlos dentro de los límites aceptables para no perturbar el confort y el bienestar de los usuarios. tales como maquinaria en funcionamiento defectuoso. paso de vehículos pesados por el lugar. Sin embargo. Los datos se obtuvieron para las aceleraciones pico de las más fuertes tormentas ocurridas durante un periodo de retorno de mas de 5 años. etc. con efectos psicológicos y fisiológicos tales como mareos.Capítulo 5 – Análisis de Construcciones con Efecto Dinámico de Viento
El confort de los ocupantes de los edificios de gran altura que soportan la acción del viento es un tema de primordial importancia en el diseño. ansiedad. el costo del edificio no debe superar los valores normales de una aceptable economía. las oscilaciones excesivas producen fisuración de la tabiquería. indicando las máximas aceleraciones aceptables para diferentes frecuencias. En general la aceleración es el parámetro predominante para determinar aproximadamente la naturaleza de la respuesta humana a las vibraciones. impactos en rampas. Por el contrario. Por otra parte. los cuales pueden variar durante la vida de servicio de la estructura. llegando hasta náusea aguda. vientos fuertes. dependiendo del uso o destino del edificio. sismos. molestias visuales o temor. 5.
T(seg) 1 2 3 4 5 6 7 8 9 10
fig.1 GRAFICA DE CONFORT AMPLITUD DE OSCILACION EN FUNCION DEL PERIODO PARA DISTINTOS VALORES DE LA ACELERACION EN PORCENTAJE DE g
INTOLERABLE 15% de g
0.5 % ESTO de g
LE IB PT e g CE % d R PE 0.Capítulo 5 – Análisis de Construcciones con Efecto Dinámico de Viento
∆f (m )
0. 5.3.10
MUY MO LESTO 5 % de g MO L 1.
005 g 0.15 g < ü
.05 g 0. TABLA 5.02 g < ü (según Chang) ü ≤ 0.3.05 g < ü ≤ 0.Capítulo 5 – Análisis de Construcciones con Efecto Dinámico de Viento
La tabla 5.0075 g < ü ≤ 0.0075 g 0.004 g < ü ≤ 0.3.004 g 0.005 g < ü ≤ 0.02 g 0.015 g 0.1 PERCEPCION HUMANA No perceptible Levemente perceptible Perceptible Molesta No perceptible Perceptible Desagradable Muy desagradable Intolerable ACELERACION ü (m/seg2) (según Khan y Parmelec) ü ≤ 0.1 a continuación da idea de la magnitud de las percepciones que se obtuvieron en experimentos realizados para diferentes niveles de aceleraciones (según Khan Parmelec y según Chang).015 g < ü ≤ 0.15 g 0.
3.2 NIVEL DE PERCEPCION HUMANA RANGO ACELERACION (m/seg2) 1 2 3 0. La exposición de larga duración puede producir malestar.50 – 0.85
.70 > 0.25 – 0. Los objetos colgantes pueden moverse algo. No se tolera el movimiento y no se puede caminar. EFECTO
0.2 (según Yamada y Goto).05 0.40 – 0.10 0. Aun se puede caminar. La oscilación puede afectar el trabajo de oficina.40 0.Capítulo 5 – Análisis de Construcciones con Efecto Dinámico de Viento
La sensación subjetiva y el comportamiento humano afectado por las diferentes aceleraciones se indica en la tabla 5.05 – 0.25 La gente no percibe el movimiento. Las personas sensibles perciben el movimiento. TABLA 5. La mayoría de las personas perciben el movimiento.10 – 0. Hay dificultad para caminar normalmente. Los objetos caen y pueden lastimar a las personas. El trabajo de oficina se vuelve difícil o casi imposible.50
0. Las personas de pie pueden perder el equilibrio.3. Se percibe fuertemente el movimiento.
5.3.5 %g 1.01
100 %g 44 %g 15 %g 1.15 0.1 %g
ACELERACION m/seg²
GRAFICA DE CONFORT
2. AMPLITUD (∆f) Y ACELERACION (x’’)
.1 %g 0.11 0.4 1.5 1
INTERACCION DE LAS VARIABLES: PERIODO (T).Capítulo 5 – Análisis de Construcciones con Efecto Dinámico de Viento
10 4.5 1 0.
5.Capítulo 5 – Análisis de Construcciones con Efecto Dinámico de Viento
GRADUACIÓN DEL EFECTO SOBRE LAS PERSONAS A INTOLERABLE B MUY POCO TOLERABLE C DEMASIADO PERCEPTIBLE D MUY PERCEPTIBLE E PERCEPTIBLE F ESCASAMENTE PERCEPTIBLE G NO PERCEPTIBLE
ACEPTABILIDAD EN EDIFICIOS NO NO
ESTADO PREVISIBLE DE LA ESTRUCTURA COLAPSO DAÑOS LOCALES
SITUACION LIMITE EN TAREAS INDUSTRIALES PESADAS EN PERIODOS BREVES EN VIVIENDAS PERIODOS LARGOS EN VIVIENDAS PERIODOS LARGOS EN VIVIENDAS
FORMACION DE GRIETAS FORMACION DE FISURAS SIN INFLUENCIA EN ESTRUCTURAS CORRIENTES SIN INFLUENCIAS SIN INFLUENCIAS
CUADRO ADJUNTO A GRAFICA DE LA fig.2
de los datos experimentales proporcionados por la sismología. velocidad y desplazamiento del terreno. Parámetros de terremotos pasados en una región → Predicción de las características sísmicas de los terremotos que afectarán dicha región Tradicionalmente.secuencia de choques. Además. pero se han observado daños moderados a aceleraciones muy altas. el valor medio cuadrático de las aceleraciones en la fase “fuerte” de un acelerograma. además de los valores máximos de la aceleración. el espectro de seudovelocidades.Capítulo 6 – Análisis de Construcciones con Efectos Sísmicos
ANÁLISIS DE CONSTRUCCIONES CON EFECTOS SISMICOS
Definición numérica de la acción sísmica
La definición correcta de la acción sísmica es un problema al que se le debe dar la mayor importancia en un análisis estructural sísmico. La solución a este problema parte. Por lo tanto. se han utilizado. etc. duración . la fuerza destructiva de un terremoto ha sido expresada en función de la aceleración máxima del terreno.contenidos de frecuencias . la acción definida estará totalmente relacionada al tipo de análisis estructural que se va a realizar: desacoplamiento modal integración directa espectros de respuesta
Análisis no lineal →
. Por ejemplo. generalmente. el espectro de amplitudes de Fourier.intensidad. intensidad espectral (Housner). se deben tener en cuenta otras características: . etc.
11) quedaría (despreciando los dos primeros términos):
ω 2 ⋅ x ≈ −a (t ) estructuras  ∴si ω → ∞. en función de las características de la estructura. representaciones gráficas de valores → aproximados de la respuesta máxima de un sistema de 1 GLD lineal y elástico:
. S a → a max muy rígidas  (copia la aceleración máxima del terreno)
. [1 y 2]) Definición mediante espectros de respuesta Es la forma más usual de definir una acción sísmica. ξ ) = −
⋅ ∫ a (τ ) ⋅ e −ξω ⋅(t −τ ) ⋅ sen ω ⋅ (t − τ ) ⋅ dτ
S v (a. ξ ) S a (a. ω . la (2. ξ )
(3. ω . . ξ ) = ω 2 ⋅ S d (a.
x + 2ξω ⋅ x + ω 2 ⋅ x = −a (t )
S d (a.Capítulo 6 – Análisis de Construcciones con Efectos Sísmicos
# También existen definiciones de acciones sísmicas a partir de la teoría de procesos estocásticos. ξ ) = ω ⋅ S d (a. Sv. No profundizaremos sobre estos temas.. dado que se obtienen descripciones de las características más importantes de la respuesta estructural sin la necesidad de disponer de una historia en el tiempo de la excitación y la respuesta. ω .53)
Puede realizarse un análisis de las variaciones de los valores de las curvas Sd.11) (3. sin necesidad de conocer los detalles de la excitación. pues escapan a nuestro curso. Otra ventaja es que un espectro puede “modificarse” en base a las características del lugar de emplazamiento de la estructura.52) (3. Recordar que ω =
Para frecuencias propias altas (en comparación con la del movimiento del terreno). (ver ref. ω . ω . en donde la respuesta estructural en sí se obtiene probabilísticamente. Sa.
05 (línea llena)
S a = a s + (b − a s ) ⋅ T T1
T ≤ T1
. β v (ξ .Capítulo 6 – Análisis de Construcciones con Efectos Sísmicos
Para frecuencias bajas. se produce una amplificación del movimiento del terreno en su paso a través del filtro estructural.11) queda:
. la (2. .
β a (ξ . a max v max d max
(6. tales como grosor y propiedades de los estratos que se encuentran entre la roca firme y la estructura modifican los espectros de respuesta. T ) = v .  x ≈ −a (t ) estructuras  si ω → 0. Conjunto de espectros de respuesta registrados para una región ↓ INPRES – CIRSOC 103 Parte 1 / Cap. β d (ξ . por el contrario.1)
Las condiciones locales del terreno. T ) = d . T ) =
Sa S S . 7 suavizado y promediado de valores Espectros de Diseño
Para ξ = 0. S d → d max muy flexibles  (copia los desplazamientos máximos del terreno)  
Para frecuencias intermedias.
) las condiciones locales del suelo
La generación de acelerogramas sintéticos requiere procedimientos matemáticos basados en procesos estocásticos y expansiones en series sinusoidales. b) Acelerogramas sintéticos Tiene grandes ventajas respecto al anterior: posibilidad de generar señales de corta duración que tengan el mismo efecto en las estructuras que el del terremoto real que se quiere simular ⇒ economía computacional al estar generados en función de un espectro de diseño.05 (línea a trazos)
S a = as + ( f A ⋅ b − as ) ⋅ Sa = f A ⋅ b
2 T2    T2  3   S a = 1 + ( f A − 1) ⋅  ⋅ b ⋅    T   T     
T ≤ T1 T1 < T ≤ T2 T2 < T
fA: factor de amplificación por amortiguamiento
ξ (expresado en %)
Definición mediante acelerogramas a) Acelerogramas reales Basar un cálculo sísmico en uno o varios registros disponibles en una zona implica un alto riesgo de definición incorrecta de la acción.005 ≤ ξ < 0. se tiene en cuenta (aprox. b.
.Capítulo 6 – Análisis de Construcciones con Efectos Sísmicos
T1 < T ≤ T2
T  Sa = b ⋅  2  T 
 zona sísmica   a s . T1 y T2 → tabla 4 → f   suelo   
Para casos con 0.
pero en dicho promedio no se incluirán valores inferiores al 85% del máximo encontrado.3.7 ⋅ S aD ⋅ γ d
para todos los puntos Cantidad de acelerogramas a emplear: Grupo A y B Grupo A0 nº ≥ 3 acelerogramas nº ≥ 4 acelerogramas
Para diseño y verificaciones se promediarán las envolventes de solicitaciones y deformaciones obtenidas para cada acelerograma. 5.2 b) El espectro elástico (para ξ = 0.
. art.05
Area C ≥ γ d ⋅ Area D
C S a ≥ 0.05) deberá cumplir:
0.1) Los acelerogramas deben satisfacer: a)
a (t )max
a(t )max ≥ a S ⋅ γ d
aS: ordenada al origen del espectro correspondiente γd: factor de riesgo. INPRES – CIRSOC 103 (Parte 1 – Cap. 14.05
0.Capítulo 6 – Análisis de Construcciones con Efectos Sísmicos
Existen muchos procedimientos y algoritmos que permiten generar acelerogramas sintéticos a partir de un espectro dado.
Métodos dinámicos a) Análisis modal espectral b) Superposición modal paso a paso c) Integración directa paso a paso 6. es imprescindible que la estructura se comporte en forma uniforme. de no verificarse se incrementaran los efectos del análisis modal de la siguiente manera:
.. 14.1) Procedimiento aproximado (cap. 47)
(S a )i
= Sa ⋅ γ d ⋅ g
S a : del espectro correspondiente
γd: factor de riesgo g: aceleración de la gravedad • se podrá considerar capacidad de disipación de energía por deformaciones anelásticas de la estructura
Sa ⋅ γ d g R las deformaciones calculadas con este criterio deberán ser amplificadas multiplicándolas por la ductilidad global µ =
para el caso anterior. Métodos dinámicos INPRES – CIRSOC 103
a) Análisis modal espectral • • • excitación sísmica traslacional materiales lineales y elásticos (según art. de manera de garantizar la ausencia de concentración de rótulas plásticas
Se debe verificar que el corte (con análisis modal espectral) en cada dirección no resulte inferior al 75% del corte obtenido con el método estático. 12.1 / pág..Capítulo 6 – Análisis de Construcciones con Efectos Sísmicos
Métodos de análisis según INPRES – CIRSOC 103
Procedimientos con fuerzas estáticas equivalentes a) b) Método estático (art.
1) se susti.5)
. pero además se deberán incluir todos los modos cuya contribución al total sea mayor que 5% de la contribución del modo fundamental. ϕ i = ϕ ri    M  O    ϕ ni  mr    
# observar que si en (6.75 ⋅ V0 est V
0 modal espectral
⋅ Rmodal
Modos a considerar: q ≥ 3. Modelo de edificio de cortante
m ⋅ϕ 1 r =1 ϕ ⋅ M ⋅ { } ∑ r ri
ϕ ⋅ M ⋅ϕi
∑ mr ⋅ ϕ
x r -1
K  O  con M =  mr   0 
0 ϕ 1i   M        .si ϕ iT ⋅ M ⋅ ϕ i = 1 tuye la masa mr por el peso wr el n n valor de Qi no cambia Qi = ∑ mr ⋅ ϕ ri = ∑ wr ⋅ ϕ ri
Peso efectivo modal
 n w ⋅ϕ  ∑ r ri   1 Wi =  r =n 2 ∑ wr ⋅ ϕ ri
r =1 2
∑ Wi ≥ 0.9 ⋅ ∑ wr
(6.Capítulo 6 – Análisis de Construcciones con Efectos Sísmicos
R *modal =
i =1 peso total de la estructura
q ⇒ apto Desplazamientos modales máximos
(x i )max = ϕ i ⋅ yi (t )max = ϕ i ⋅ Qi ⋅
..Capítulo 6 – Análisis de Construcciones con Efectos Sísmicos
Desplazamientos modales relativos entre pisos
∆ ri = ( x r )i − ( x r −1 )i
FriS = ϕ ri ⋅ Qi ⋅ (S a )i ⋅ mr = ϕ i ⋅ Qi ⋅
⋅ wr
fuerza sísmica equivalente en el nivel r para oscilaciones en el modo i Vectorialmente:
en función del peso wr y la fracción de g
 F1iS     M  ϕT ⋅M ⋅J   FiS =  FriS  = M ⋅ ϕ i ⋅ Ti ⋅ (S a )i ϕi ⋅ M ⋅ϕi  M    S  Fni   
x i (t )max = Qi ⋅ (S a )i
.  x ri  = ϕ ⋅ x i (t )   ri max   max
FiS = M ⋅ ϕ i ⋅ Qi ⋅ (S a )i
. FriS =  x ri  = mr     max r: nivel i: modo
(6.9) (6.11)
.  x ri  = ϕ ⋅ Q ⋅ (S )   ri i a i   max .6)
Fuerzas sísmicas equivalentes (modales)
12) (6.13b. multiplicada por un factor f1)
Distribución en altura de Vm según INPRES – CIRSOC 103:
Fkm = Wk ⋅ φ km
∑Wi ⋅ φim
⋅ Vm
m: modo k: nivel Fkm: fuerza sísmica equivalente para el modo m en el nivel k Según la fórmula (6.13b)
Vi = ∑ FriS
Vi = Wi ⋅
Comparación con INPRES – CIRSOC 103
⋅ S am ⋅ Wm
(Fórmula reglamento)
γd: factor de riesgo R: coeficiente de reducción por disipación i=m
Vi = f1 ⋅
⋅ Wi
(Fórmula 6.13a) (6.13a):
Vi = ∑ ϕ ri ⋅ Qi ⋅ wr ⋅
Qi ⋅ (S a )i n ⋅ ∑ ϕ ri ⋅ wr g r =1 Qi ⋅ (S a )i n ⋅ ∑ ϕ ri ⋅ wr g r =1
si queremos emplear la notación
FriS = ξ ri ⋅ Vi
FriS = ξ ri ⋅
comparándola con (6.10) y (6.10) se deduce que
Cortante modal En nivel “r”: En la base: o también:
S Vri = ∑ F ji j =r
subestructura "P"
α p ⋅ FriS
α 1 + α 2 + Kα P + Kα m = 1
m: nº de subestructuras verticales
(se deben adicionar los efectos torsionales si los hubiese) CALCULO DE LA RESPUESTA TOTAL
∑ (R j )max
 ϕ ⋅w W ⋅φ  ξ ri = n ri r = n k km  notación INPRES .R. se calculan los esfuerzos seccionales para cada modo en estudio.M. C.CIRSOC 103 ∑ ϕ ri ⋅ wr ∑ Wi ⋅ φ im  r =1 i =1 
Momento de vuelco modal (para el nivel r)
v M ri =
∑ F jiS (h j − hr )
si r = 0 → M0iv = momento de vuelco en la fundación Momento de torsión modal
C. al piso “r”
Esfuerzos seccionales modales Para cada subestructura vertical con la parte proporcional a su rigidez que le corresponde de la fuerza sísmica equivalente.
t M ri = er ⋅ Vri
More From This UserInforme de Diagnóstico Situación Actual de La Planta de Asfalto Ciber-2Formas de MetrarLimites de SensibilidadDiferencias que existen entre la elaboración de un expediente técnico que se ejecuta por administración directa y un expediente técnico que se ejecuta por la modalidad de ejecución indirectaSeccion Suelos y Pavimentos Manual de CarreterasManual de Carreteras - Especificaciones Tecnicas Generales para Construcción - EG-2013 - (Versión Revisada - JULIO 2013).pdfManual Maestro ObraResidente.Supervisión de Obras.pdfTARJETA-DOBLADO-BACO.pdfTaller I_Wendy Juy Mori
Libro de Dinamica Estructural by Wendy Juy Mori33 viewsEmbedDownloadRead on Scribd mobile: iPhone, iPad and Android.Copyright: Attribution Non-Commercial (BY-NC)List price: $0.00Download as PDF, TXT or read online from ScribdFlag for inappropriate contentMore informationShow less
RelatedDISEÑO SISMICO DE EDIFICIOS - Bazan & Meliby WILBER CUTIMBO CHOQUE DISEÑO SISMICO DE MUROS DE CONTENCIÓN EN GRAVEDAD Yby Ivan Manosalva Aguilaredificio paorticado diseño sismicoby tigreaxulANALISIS DE DISEÑO SISMICO A BASE DE MARCOS CON RECDF 2004.pdfby perezleonardoMetodo Estático y Dinámicoby Wilian Edson Palomino ArangoPractica y Diseno Sismicoby tafo55shernandez DISEÑO ESTANQUE METALICOby PROYECTOMSCálculo Sísmico de Edificiosby vgarciaherrerosoIntroduccion Al Sap2000 Morrisby RamilArtatesDISEÑO AUTOMATIZADO DE ESTRUCTURAS - COLISEO CERRADO DE LA UNIVERSIDAD NACIONAL DE CAJAMARCAby WILBER CUTIMBO CHOQUE 20070430-Ejm%20Edificio%20Alba%20Confinadaby janmsaCFE Sismo 08by Betin NebuchadnezzarSimilar to estruc_dinamDISEÑO SISMICO DE EDIFICIOS - Bazan & MeliDISEÑO SISMICO DE MUROS DE CONTENCIÓN EN GRAVEDAD Yedificio paorticado diseño sismicoANALISIS DE DISEÑO SISMICO A BASE DE MARCOS CON RECDF 2004.pdfMetodo Estático y DinámicoPractica y Diseno Sismicoshernandez DISEÑO ESTANQUE METALICOCálculo Sísmico de EdificiosIntroduccion Al Sap2000 MorrisDISEÑO AUTOMATIZADO DE ESTRUCTURAS - COLISEO CERRADO DE LA UNIVERSIDAD NACIONAL DE CAJAMARCA20070430-Ejm%20Edificio%20Alba%20ConfinadaCFE Sismo 08Analisis Sismico de Edificiosestruc_dinam

References: resolución 
 resolución 
 Resolución 
 resolución 
 resolución 
 Resolución 
 resolución 
 resolución 
 Resolución 

resolución 
 resolución 
 resolución 
 resolución 

Resolución