Analisis de Estructuras Bajo Acciones Di
Analisis de Estructuras Bajo Acciones Di
Analisis de Estructuras Bajo Acciones Di
ISBN 978-987-25360-9-1
2009
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.
1
TABLA DE CONTENIDOS
2
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
3
BIBLIOGRAFÍA:
4
CAPÍTULO 1 - NOCIONES BASICAS DE SISMOLOGIA
sismógrafo
estructura
epicentro
foco
1-1
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.
1-2
Donovan:
1080 e 0,5 M
am =
(R + 25)1,25
am: aceleración máxima del terreno en [cm/s2]
R: distancia focal en [km]
Esteva y Villaverde:
a m 5,7 e 0,8 M
=
g (R + 40)2
32 e M
vm =
(R + 25)1,7
am: aceleración máxima en [cm/s2]
vm: velocidad máxima en [cm/s]
R: distancia focal en [km]
Esteva:
a m = 1230 e 0,8 M (R + 25)
−2
(
v m = 15 e M R + 0,17 e 0,59 M )−1, 7
1-4
CAPÍTULO 2 - CONCEPTOS BASICOS DE DINAMICA ES-
TRUCTURAL
En nuestro curso, los sistemas a estudiar serán sistemas estructurales, las variaciones en el
tiempo serán vibraciones producidas por cargas dinámicas.
SISTEMAS vibraciones
ESTRUCTURALES
ecuaciones diferenciales
cargas dinámicas que gobiernan el
comportamiento de las
vibraciones
permiten evaluar el
comportamiento de resolución de las
la estructura frente ecuaciones diferenciales
a acciones dinámicas
2-1
- sismos
- viento
- olas y corrientes de agua
- explosiones e impactos
- cargas móviles (vehículos, personas, etc.)
Respuesta dinámica → cualquier magnitud que pueda caracterizar el efecto de una carga
dinámica sobre la estructura
2-2
F(t)
t3 F
ACCION
t1 F(t2)
t2
F(t1)
M(t1) F(t3)
t1 t2 t3 t
M(t2)
este esquema temporal
de carga debe ser
M(t3) perfectamente conocido
RESPUESTA
Fig. 2.2 - Acción y respuesta determinista
ESTRUCTURA la comparación se
basa en: - velocidad de
la acción
ACCION - periodo propio
DINAMICA de la estructura
a(t) da origen a fuerzas
de "inercia" comparables
masa M con las estáticas
depende de
rigidez K K , M
No considerada como
"carga" sino como una
propiedad intrínseca de
la estructura
Fig. 2.3 - Acción dinámica y propiedades de la estructura
2-3
a) Estructura SIN MASA SIN INERCIA
x(t)
F(t)
A A' F
F0
rigidez: K t1 t
masa: m = 0
x
K x(t) = F(t) x0 = F0/K
La respuesta seguirá
exactamente la forma
de la carga
t1 t
x(t)
F(t) m/2
F
m/2
F0
rigidez: K t1 t
masa: m
xmáx
m x(t ) + k x(t ) = F (t ) x
x0 vibraciones
remanentes
Se desarrolla energía
cinética, que modifica t1 t2 t
la respuesta y deja vibraciones t2 > t1 !!
remanentes
2-4
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).
con AMORTIGUAMIENTO
la amplitud decrece en cada ciclo
x ESTRUCTURA SIN
AMORTIGUAMIENTO
m x0
k x0
t
tD t tD t
En general si tD >> T no es necesario un análisis dinámico
si tD ≅ T PROBLEMA DINAMICO
2-5
x0 a masa constante
m x0 k 1 > k2
a rigidez constante
k m2 > m1
t
T1 T2
2 1
m;k;
m1 m2
variables → → T1 → T2 , etc
k1 k2
Fig. 2.6 - Velocidad de reacción; varios Ti
2-6
ayed respuesta en un punto
axed en particular
onda
refractada a2 discontinuidad
a1 a3 en el suelo
onda onda
original reflejada
a1 axed
t t
onda sísmica
respuesta
original
2-7
ESTRUCTURA MODELO
REAL DINAMICO
x2 m2 MODELO
k2 MATEMATICO
x1 m1
k1 [M ] x + [K ] x = −[M ]a
EXCITACION
SISMICA
PROCEDIMIENTOS
NUMERICOS
x1
menos exacta
t
mas exacta
RESPUESTA (desplaz. en el piso 1)
2-8
respuesta dinámica
en CADA PUNTO nº infinito de
Modelo DE LA ESTRUCTURA puntos espaciales
dinámico
EXACTO en CADA INSTANTE nº imfinito de
DE TIEMPO puntos temporales
Modelo
dinámico
DISCRETO
nº FINITO DE PUNTOS DISCRETIZACION
ESPACIALES (GLD) ESPACIAL
2-9
Método de las masas concentradas
nº total de componentes
nº de GLD de desplazamiento según
del modelo los cuales vibran las masas
concentradas
se calcula la "deformada"
se distinguen: del modelo en cada instante
- MODELOS DE 1 GLD
- MODELOS DE tensiones y deformaciones
MULTIPLES GLD específicas pueden conocerse
mediante los procedimientos
del análisis estático
m m m
x
c c
k k k
a(t)
m
x
k
a(t)
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.
2-10
Modelos con múltiples GLD:
mn mn mn
kn kn cn kn cn
xr xr xr
mr mr mr
kr kr cr kr cr
m2 m2 m2
k2 k2 c1 k2 c1
m1 m1 m1
k1 k1 c1 a(t) k1 c1
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.
2-11
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).
2-12
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
t2 t2
πH = (E p − Ec ) dt + Ed dt ; dπ H = 0 (2.1)
t1 t1
2-13
F(t) x (t ) Fi(t) Fa(t) Fe(t)
x (t ) inercia amort. elástica
x(t)
F(t) = Fi(t) + Fa(t) + Fe(t) equilibrio
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
Fi (t ) = m x(t ) (2.3)
ii- de amortiguamiento
Fa (t ) = c x(t ) (2.4)
iii- elásticas
Fe (t ) = k x(t ) (2.5)
F (t ) − Fi (t ) − Fa (t ) − Fe (t ) = 0 (2.6)
Fi (t ) + Fa (t ) + Fe (t ) = F (t ) (2.7)
mx+cx+k x = F (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.
2-14
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:
Fi = m [a + x ] (2.9)
− m [a + x ] − c x − k x = 0 (2.10)
m x + c x + k x = −m a (2.11)
m x + c x + k x = F − ma (2.12)
F T − Fi − Fa − Fe = 0 (2.14)
2-15
#Nota: en adelante, los vectores y matrices serán representados(en general) con minús-
culas y mayúsculas en negrita, respectivamente.
M x + C x + K x = F − M {J a} (2.15)
M → matriz de masas
C → matriz de amortiguamiento
K → matriz de rigidez
JT = [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:
M x + C x + K x = − M {J a} (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.
M D + C D + K D = F − M {J a} (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.
2-16
CAPÍTULO 3 - RESPUESTA DE UN OSCILADOR SIMPLE
mx+ k x = 0 (3.1)
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:
k k
ω2 = ; ω= (3.2)
m m
se obtiene:
x +ω2x = 0 (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)
2π
3-1
Finalmente, otra característica es el período natural
1
T= (3.5)
f
2π
T= (3.6)
ω
x = A sen (ω t + Ψ ) (3.7)
resulta:
2
2 x0
A= x + 0
ω
ω x0
tanΨ =
x0
Fórmula de Geiger
Sustituyendo G G: peso de m
m=
g g: aceleración de la gra-
vedad
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
ω= g⋅ (3.9)
X SG
3-2
2π
T= ⋅ X SG (3.10)
g
T = 2,00 X SG (3.11)
x + 2β x + ω 2 x = 0 (3.12)
c
2β = (3.13)
m
x = e rt (3.14)
3-3
r 2 + 2β r + ω 2 = 0 (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
r1, 2 = − β ± β 2 − ω 2 = 0 (3.16)
r1, 2 = − β ± iω 1 − ξ 2 = 0 (3.17)
con i = − 1
ωv = ω 1 − ξ 2 (3.18)
se obtiene:
r1, 2 = − β ± iω v (3.19)
r1, 2 = −ξω ± iω v (3.20)
c 2β m
ξ= = (3.21)
cr 2ω m
3-4
conocida como fracción de amortiguamiento crítico (en estructuras corrientes
0.02<ξ<0.06), también:
β c
ξ= o ξ=
ω 2mω
Volviendo a la resolución de la ecuación de movimiento (3.12), escribimos la solución ge-
neral en la forma:
3-5
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:
pero t n +1 = t n + Tv ; con Tv = 2π ω v
2π
∴ t n +1 = t n + , entonces
ωv
e −ξωtn e ωv
sen(ω v t n + 2π + Ψ )
xmáx (n ) 1
= −ξ 2π ω ω (3.26)
xmáx (n + 1) e v
xmáx (n )
ln
xmáx (n + 1)
[
= ln eξ 2π ω ω v ] (3.27)
3-6
ξ = 0,1 → ω v = ω 1 − ξ 2
ω v = ω ⋅ 0,995 → ω v ≈ ω
Entonces
x máx (n )
ln = 2πξ (3.28)
x máx (n + 1)
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.
3-7
Excitación armónica
Si la carga es de tipo:
P = P0 sen (Ωt )
Entonces, la ecuación de movimiento será:
m x + c x + k x = P0 sen(Ωt ) (3.30)
Ω: frecuencia de la excitación
P0
x + 2ξω x + ω 2 x = sen(Ωt ) (3.31)
m
x g = xh + x p (3.32)
xh = c1 e r1t + c2 e r2t (3.22)
Utilizando matemática para números complejos, esta última ecuación puede escribirse:
3-8
Derivando y reemplazando en (3.30) se obtienen las constantes A y B:
Ω
Denotando γ = (3.36)
ω
P0 1−γ 2
A= ⋅ (3.37)
k (1 − γ 2 )2 + (2ξγ )2
P0 − 2ξγ
B= ⋅ (3.38)
k (1 − γ 2 )2 + (2ξγ )2
Condiciones iniciales
Para el caso en que x(0 ) = x 0 ; x(0 ) = x0
Pueden calcularse las partes correspondientes a la solución particular (para t = 0)
P0 − 2ξγ
x p (0 ) = ⋅ (3.39)
k (1 − γ ) + (2ξγ )2
2 2
P0 Ω(1 − γ 2 )
x p (0 ) = ⋅ (3.40)
k (1 − γ 2 )2 + (2ξγ )2
x g (0 ) = x 0 = c 2 '+ x p 0 (3.41)
x g (0 ) = x0 = −ξω x0 + ξω x p 0 + ω v c1 '+ x p 0 (3.42)
1
c1 ' =
ωv
[x 0 + ξω x0 − x p 0 − ξω x p 0 ] (3.43)
c2 ' = x0 − x p0 (3.44)
3-9
Forma de operar:
1- Se calculan x p 0 y x p 0 con las 3.39 y 3.40
2- Se computan c1’ y c2’ 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
f α
i f i
τ dτ t dτ t
F(t)
P(τ )
entonces: xf = dτ (a)
m
3-10
por otra parte
1
x(dτ )
2
x f = xi + x dτ + (b)
2
1
x(dτ )
2
xf = (b)
2
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.
P (τ )
dx(t ) = e −ξω (t −τ ) dτ sen (ω v (t − τ )) (3.45)
mω v
1 t
x(t ) = P(τ )e −ξω (t −τ ) sen[ω v (t − τ )]dτ (3.46)
mω v 0
1 t
x(t ) = − a(τ )e −ξω (t −τ ) sen[ω v (t − τ )]dτ (3.47)
ωv 0
3-11
(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.)
La solución a esta ecuación viene dada por la ec. (3.47) y el máximo valor absoluto será
Sd.
Reemplazando S d = x(t ) máx
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:
S vr = x(t ) máx ; S ar = x(t ) máx
(Ver ref. [1], cap. 5)
3-12
Fig. 3.6 - Espectro sísmico de respuesta para desplazamientos. Sd
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 ωv 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:
1 t
sd (ω ,ξ ) = − a(τ )e −ξω (t −τ ) sen[ω (t − τ )]dτ (3.49)
ω0 máx
t
sv (ω ,ξ ) = − a(τ )e −ξω (t −τ ) sen[ω (t − τ )]dτ (3.50)
0 máx
3-13
t
sa (ω ,ξ ) = ω a(τ )e −ξω (t −τ ) sen[ω (t − τ )]dτ (3.51)
0 máx
sv = ω S d (3.52)
sa = ω 2 S d (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 Sa. 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.
3-14
Fig. 3.8 - Espectro de diseño suavizado de Seed e Idriss (1982).
3-15
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.
Fig. - 3.10
x ∆t = t i +1 − t i
xi +1
haciendo el cambio de variable
τ = t − t i ; τ ≤ ∆t
xi
puede suponerse que el valor de
∆t la aceleración de respuesta en
ti ti+1 t un instante τ se expresa como:
τ
x(τ ) = xi + f (τ )( xi +1 − xi ) (3.55)
0 para τ = 0
f(τ) =
1 para τ = ∆t
τ
x(τ ) = xi + xiτ + ( xi +1 − xi ) f (τ )dτ (3.56)
0
3-16
τ
g (τ ) = f (τ )dτ (3.57)
0
∆τ
γ ∆t = f (τ )dτ (3.58)
0
xi +1 = xi + xi ∆t + ( xi +1 − xi )γ ∆t (3.60)
xi +1 = xi + [(1 − γ ) xi + γ xi +1 ]∆t (3.61)
τ
τ2
x(τ ) = xi + xiτ + xi + ( xi +1 − xi ) g (τ )dτ (3.62)
2 0
τ
β ∆t 2 = g (τ )dτ (3.63)
0
1
xi +1 = xi + xiτ + − β xi + β xi +1 ∆t 2 (3.64)
2
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.
3-17
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:
γ = 0,5
β = 0,25
m xi +1 + c xi +1 + k xi +1 = F i +1 (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 a0,…, a7; 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 x 0 , x0 y x0 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:
1 γ 1
a0 = ; a1 = ; a2 =
γ ∆t 2 β ∆t β ∆t
1 γ ∆t γ
a3 = −1 ; a4 = −1 ; a5 = −2
2β β 2 β
a 6 = ∆t (1 − γ ) ; a 7 = γ ∆t
5- formar la “rigidez efectiva”
kˆ = k + a 0 m + a1 c
B) para cada paso de integración
1- formar el término de “carga efectiva” en t+∆t
rˆi +1 = F i +1 + m(a0 xi + a2 xi + a3 xi ) + c(a1 xi + a 4 xi + a5 xi )
2- resolver el desplazamiento en t+∆t
3-18
kˆ ⋅ x i +1 = rˆ1+1
3- calcular aceleraciones y velocidades en t+∆t
xi +1 = a0 ( x+1i − xi ) − a2 xi − a3 xi
xi +1 = xi − a6 xi + a7 xi +1
3-19
CAPÍTULO 4 - RESPUESTA DINÁMICA DE UNA ESTRUCTURA
CON MÚLTIPLES GRADOS DE LIBERTAD
M D + C D + K D = F (t ) (4.1)
M D + C D + K D = − M {J a} (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
4-1
M D+ K D=0 (4.4)
D = A sen(ω t + Ψ ) (4.5)
A: vector que contiene las amplitudes de las vibraciones
Ψ: ángulo de fase
− M ω 2 A sen(ω t + Ψ ) + K A sen(ω t + Ψ ) = 0
{K − ω 2
M }A sen(ω t + Ψ ) = 0 (4.7)
sen (ω t + Ψ ) ≠ 0
{K − ω 2
M }A = 0 (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:
K − ω2M = 0 (4.9)
4-2
Cuando K y M son definidas positivas (caso usual en estructuras civiles), de la ecuación
característica se obtienen “n” soluciones positivas ωi2 y en consecuencia “n” valores ωi,
siendo “n” el número de GLD de la estructura.
Los valores ωi se denominan “frecuencias propias” o “pulsaciones” de la estructura y los
“n” periodos propios se calculan
2π
Ti = i = 1,...n (4.11)
ωi
Siendo T1 el período correspondiente ω1, que es la frecuencia de menor valor, éste se de-
nomina “periodo fundamental” del sistema. Reemplazando cada ωi en la (4.8) se obtiene
el correspondiente Ai 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 Ai, por lo tanto inferimos que no interesa la magnitud
de las componentes de Ai sino la relación existente entre ellas. Por esto es conveniente
“normalizar” los distintos Ai.
ϕi1
Ai
ϕi = (4.12) ϕi = 1 modo i normalizado
aimáx
ϕin
2. Según la matriz de masas
ϕ i = Ai (M i *)
−1
2 (4.13)
4-3
Con
M i * = Ai [M ]Ai (4.14)
T
ϕ i T [M ]ϕ i = 1 (4.15)
d3 d4 d4
kv m2 m
k2 k2
d1 d2 d2
kv m1 m
k1 k1
Para el ejemplo de la figura 4-1, tenemos que, observando el vector desplazamientos está-
ticos Ds podemos determinar que parte queremos eliminar y cuales grados de libertad de-
seamos conservar como GLD.
4-4
D D: vector con los GLD
D s
=
D e
De: grados de libertad a eliminar
K s Ds = Fs (4.1)
K 11 K 12 D F
= (4.16)
K 21 K 22 De Fe
Fe = 0 (4.17)
K 21 D + K 22 De = 0 (4.18)
de donde
−1
De = − K 22 K 21 D (4.19)
{ −1
K 11 D − K 12 K 22 K 21 D = F } (4.20)
reordenando
{K 11
−1
-K 12 K 22 K 21 D = F} (4.21)
KD=F (4.22)
donde K es la matriz de rigidez condensada y vale
4-5
−1
K = K 11-K 12 K 22 K 21 (4.23)
K↓ D + M
↓
D=0
2× 2 2×2
(4.24)
d2
D=
d4
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:
δ
m2
δ
1
m1
δ
1
δ δ
F= 11 21
δ 12
δ 22
δ
K = F −1
Matriz de amortiguamiento
Dado el sistema
4-6
M D + C D + K D = F (t ) (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:
m1 0
M= mi (4.25)
0 mn
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:
1 0
Ω= ; Φ = [ϕ 1 ,ϕ 2 ϕn ] (4.26)
0 N
ϕi ϕ j = 0
T
i≠j (4.27)
ϕi M ϕ j = 0 i≠j
T
(4.28)
4-7
ϕi K ϕ j = 0 i≠j
T
(4.29)
ΦT M Φ = I (4.30)
Φ T K Φ = K* (4.31)
k1 *
T
K *= ki * ; ki * = ϕ i K ϕ i (4.32)
kn *
C = α1 M + α 2 K (4.33)
La condición de ortogonalidad:
ϕi C ϕ j = 0 i≠j
T
(4.34)
ϕi C ϕi = 2 (4.35)
T
i i
4-8
Los coeficientes α1 y α2 se calculan a partir de las ecuaciones anteriores (4.35) y (4.33)
ϕi
T
[ 1 M+ 2 K ]ϕ i = 2 i i (4.36)
ϕ i M ϕ i + 2ϕ i K ϕ i = 2 (4.37)
T T
1 i i
α 1 + α 2ω i2 = 2ω i ξ i (4.38)
Recordar que si la normalización de los autovectores ϕi fue hecha según (4.13) y (4.14) el
producto
2
ϕi K ϕi = (4.39)
T
i
Esto se debe a que la (4.8) puede expresarse para A = ϕi y premultiplicarse por ϕiT
ϕ i K-
T
[ i
2
M ϕi = 0 ] (4.40)
2
ϕi K ϕi − ϕi M ϕi = 0
T T
i
(4.41)
1
4-9
K A − ω2M A = 0 (4.42)
K A = ω2M A (4.43)
K A 1 = ω 2 M A0 (4.44)
A1 ≈ K -1 M A0 (4.45)
A i + 1 ≈ K -1 M A i (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 A0 inicial. Conviene que los valores correspondan a una defor-
mada suave, no obstante esto solo acelera la convergencia.
2- Se calcula
A1 = K -1 M A0 (4.47)
3- Se normaliza A1
A1 = A1 {A1T M A1 }
−1
2
(4.48)
i=1
4- Se mejora la solución
A i + 1 = K -1 M A i (4.49)
5- Se normaliza la nueva solución
Ai +1 = Ai +1 {Ai +1,T M Ai +1 }
−1
2
(4.50)
6- Se calcula la frecuencia correspondiente
1
T 2
i +1 i +1
A KA
ω i +1 = T
(4.51)
i +1 i +1
A MA
4-10
7- Se controla la convergencia
ω i +1 − ω i
≤ TOL →
NO
se vuelve al paso 4 (4.52)
ω i +1
↓SI
ϕ1 = Ai +1 ; ω 1 = ω i +1 (4.53)
Se obtuvo el primer autovalor y su correspondiente autovector.
Siempre: i→∞ ; Ai → ϕ1 ; ω i → ω1
Para determinar q1, premultiplicamos la (4.54) por [ϕ1T M] y aplicando las propiedades de
ortogonalidad nos queda
[ϕ ] [ ]
n
T
1 M A20 = ϕ 1T M qj ϕ j (4.56)
j =1
ϕ1 M A20
T
q1 = T (4.57)
ϕ1 M ϕ1
y en general
ϕ i T M Ai0+1
qi = T (4.58)
ϕi M ϕi
y para los vectores de arranque (ortogonales a los previamente calculados)
i −1
A*i0 = Ai0 − qj ϕ j (4.59)
j =1
4-11
Procedimiento:
1. Se estima Ai0k
2. Se calcula qi con (4.58)
3. Se calcula Ai*0k con (4.59)
4. Se sigue con el método de Stodola propiamente dicho
M D + C D + K D = − M {J a} (4.2)
M D+ K D=0 (4.4)
[K − 2
]
M ϕ =0 (4.8)
Recordemos que los autovectores ϕi 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
n
D= ϕ i yi (t ) (4.60)
i =1
4-12
Donde yi(t) es un escalar función del tiempo a determinar, llamado “respuesta generaliza-
da”.
Sustituyendo (4.60) en (4.2) obtenemos
n n n
M ϕ i yi (t ) + C ϕ i yi (t ) + K ϕ i yi (t ) = − M J a (4.61)
i =1 i =1 i =1
n
ϕ MT
j ϕ i = ϕ Tj Mϕ j = M j * (4.62)
i =1
n
ϕ Tj K ϕ i = ϕ Tj Kϕ j = K j * (4.63)
i =1
n
ϕ C T
j ϕ i = ϕ Tj Cϕ j = C j * (4.64)
i =1
Entonces la (4.61) premultiplicada por un ϕjT cualquiera queda como una ecuación de 1
GLD
M j * y j (t ) + C j * y j (t ) + K j * y j (t ) = −ϕ Tj M J a(t ) (4.65)
Recordando que:
ϕ j Cϕ j = Cj* = 2
T
j j (4.35)
2
ϕ j K ϕ j = K j* =
T
j (4.39)
Y dividiendo ambos miembros por M j * = ϕ j M ϕ j
T
4-13
ϕ M J a(t )
T
y j (t ) + 2 y j (t ) + 2
y j (t ) = − j T (4.66)
ϕj Mϕj
j j j
Al valor
ϕ MJ
T
Q j = jT (4.67)
ϕj Mϕj
Finalmente, el sistema (4.2), al ser “proyectado” según la base formada por la matriz mo-
dal Φ = [ϕ1,… ϕn], queda en “n” ecuaciones diferenciales de 1 GLD, del tipo:
yi (t ) + 2 i i yi (t ) + 2
i yi (t ) = −Qi a (t ) (4.68)
Esta ecuación puede resolverse con cualquiera de los métodos vistos en la UT3. Una vez
obtenidos los “n” yi(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
q
D(t ) = ϕ i yi (t ) ; q < n (4.69)
i =1
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
4-14
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 D0 , D0 y D0
3. Seleccionar el paso de tiempo y los parámetros γ y β
∆t ≤ Tmín/10 ; γ=0,5 ; β=0,25
4. Calcular las constantes a0,…a7
5. Formar la matriz de rigidez efectiva:
Kˆ = K + a0 M + a1C (4.70)
4-15
y i (t ) máx = Qi (S a )i (4.75)
Qi
y i (t ) máx = (S ) (4.76)
ω i2
a i
Podemos, entonces calcular los máximos desplazamientos (en todos los GLD) para el mo-
do “j”:
Qj
j
Dmáx = ϕ j y j (t ) máx = ϕ j 2
(S )
a j (4.77)
j
(S )
= Aj
j a j
Dmáx (4.78)
ω 2
j
q
Dmáx ≠ i
Dmáx (4.79)
i =1
q
Rmáx = (R ) j
máx
2
(4.80)
j =1
denominada SRSS
4-16
CAPÍTULO 5 - ANÁLISIS DE CONSTRUCCIONES CON
EFECTOS DINÁMICOS DE VIENTO
T(z): fuerza aerodinámica de
empuje, paralela a la dirección
flujo T(z) del flujo
L(z) 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 T1 > 1 seg
d) amortiguamiento ξ < 0,01
X (z ) X max ( z ) = x( z ) + xmax ( z )
x (z )
max
x( z ) : desplazamiento medio
5-1
Se define como factor de ráfaga
x max ( z )
G(z ) = 1 +
x( z )
entonces
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.
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
z
ln 0 , 0706
z 0,i z 0 ,i
V ( z ) = V0 ⋅ ⋅
10 z 0,1
ln
z 0,i
V0: velocidad básica de diseño (m/seg)
z0,i: parámetro que depende de la rugosidad i
z0,1: z0,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. 5.2.4.2, CIRSOC 102)
c2: factor por cambio de tiempo en velocidad media
(Tabla 3/pág. 20, CIRSOC 102-1)
q0: presión dinámica básica (art. 5.2.3 , CIRSOC 102)
5-2
B+r se calcula mediante tablas y ábacos que utilizan
G = 1 + 1,234 ⋅ K ⋅ → una serie de parámetros auxiliares
J
En cada nivel se comparará:
NO → se adopta qz para los esfuerzos
q’z > qz
SI → se adopta q’z para los esfuerzos
CIRSOC 102 → q z = c z ⋅ c d ⋅ q 0
II) Aceleraciones
..
X max ( z ) = K ⋅ σ . . ( z ) (paralela al viento)
X
V
h T(z)
L(z) d
5-3
La velocidad crítica del viento que produce el fenómeno de resonancia es
d
Vcr =
S ⋅T
d: ancho de la superficie maestra (puede ser variable)
cilindros ≈ 0,20
S: nº de Strouhal
prismas ≈ 0,25 a 0,30
Vcr // → T //
T: período propio
Vcr ⊥ → T⊥
Si Vcr > 25 m/seg, entonces puede prescindirse del cálculo de la resonancia.
De lo contrario:
Las fuerzas L(z) y Tz obtenidas para la velocidad crítica (correspondientes al período per-
pendicular y paralelo respectivamente) se suman de la siguiente manera:
F ( z ) = L( z ) + T z2
2
y se debe comparar con las correspondientes obtenidas con q’z o qz dadas para la velocidad
de diseño.
5-4
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.
5-5
∆f (m )
0.40
0.30
INTOLERABLE
15% de g
LESTO
5 % de g
0.20
MUY MO
1.5 % ESTO
de g
L
MO
0.10 E
IBL
T g
CEP d e
R %
PE 0.5
NO SE
PERCIBE
T(seg)
1 2 3 4 5 6 7 8 9 10
fig. 5.3.1
GRAFICA DE CONFORT
5-6
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/seg2)
(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 < ü
5-7
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 EFECTO
(m/seg2)
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.
5-8
10 100 %g
4.4 44 %g
15 %g
1.5
1 1.5 %g
0.15 1.1 %g
0.1 %g
0.11
0.01
ACELERACION
m/seg
fig. 5.3.2
GRAFICA DE CONFORT
5-9
GRADUACIÓN ACEPTABILIDAD ESTADO
DEL EFECTO EN PREVISIBLE
SOBRE LAS EDIFICIOS DE LA
PERSONAS ESTRUCTURA
A NO COLAPSO
INTOLERABLE
B NO DAÑOS
MUY POCO LOCALES
TOLERABLE
C SITUACION FORMACION
DEMASIADO LIMITE DE
PERCEPTIBLE GRIETAS
D EN TAREAS FORMACION
MUY INDUSTRIALES DE
PERCEPTIBLE PESADAS FISURAS
E EN PERIODOS SIN INFLUENCIA
PERCEPTIBLE BREVES EN EN ESTRUCTURAS
VIVIENDAS CORRIENTES
F PERIODOS SIN
ESCASAMENTE LARGOS INFLUENCIAS
PERCEPTIBLE EN VIVIENDAS
G PERIODOS SIN
NO LARGOS INFLUENCIAS
PERCEPTIBLE EN VIVIENDAS
5-10
CAPÍTULO 6 - ANÁLISIS DE CONSTRUCCIONES CON
EFECTOS SISMICOS
desacoplamiento espectros de
Análisis lineal → →
modal respuesta
integración
Análisis no lineal → → acelerogramas
directa
6-1
# 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])
.. . (2.11)
x + 2ξω ⋅ x + ω 2 ⋅ x = −a (t )
1 t (3.49)
S d (a, ω , ξ ) = − ⋅ a (τ ) ⋅ e −ξω ⋅(t −τ ) ⋅ sen ω ⋅ (t − τ ) ⋅ dτ
ω 0 max
6-2
Para frecuencias bajas, por el contrario, la (2.11) queda:
..
x ≈ −a (t )
estructuras
si ω → 0; S d → d max
muy flexibles
(copia los desplazamientos máximos del terreno)
Para frecuencias intermedias, se produce una amplificación del movimiento del terreno en
su paso a través del filtro estructural.
Sa S S (6.1)
β a (ξ , T ) = ; β v (ξ , T ) = v ; β d (ξ , T ) = d ;
a max v max d max
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 Espectros
suavizado y promediado
respuesta de
de valores
registrados para Diseño
una región
↓
INPRES – CIRSOC 103
Parte 1 / Cap. 7
Sa
g
b
as
T1 T2
6-3
Sa = b T1 < T ≤ T2
2
3
T
Sa = b ⋅ 2 T2 < T
T
zona sísmica
a s , b, T1 y T2 → tabla 4 → f
suelo
6-4
Existen muchos procedimientos y algoritmos que permiten generar acelerogramas sintéti-
cos a partir de un espectro dado.
a(t )max ≥ a S ⋅ γ d
a (t )max
aS: ordenada al origen del espectro correspon-
diente
γd: factor de riesgo, art. 5.2
Area C Area D
CALCULADO DISEÑO
0,05 T0 0,05 T0
S aC S aD
Area C ≥ γ d ⋅ Area D
S aC ≥ 0,7 ⋅ S aD ⋅ γ d
para todos los puntos
6-5
Métodos de análisis según INPRES – CIRSOC 103
Métodos dinámicos
a) Análisis modal espectral
b) Superposición modal paso a paso
c) Integración directa paso a paso
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:
6-6
0,75 ⋅ V0 est
R *modal = ⋅ Rmodal
espectral V espectral
0 modal
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.
xn
n
mr ⋅ ϕ ri
kn ϕ ⋅ M ⋅ {1}
T
r =1
Qi = i
= (6.1)
ϕ ⋅ M ⋅ϕi
T
i
n
mr ⋅ ϕ 2
ri
i =1
mr 0 ϕ 1i
xr
wr kr
x r -1 con M = mr ; ϕ i = ϕ ri
hr m2
0 mr ϕ ni
m1 x1
k1 a (t )
# observar que si en (6.1) se susti- 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
i =1 i =1
r =1 peso total de la
estructura
q apto
(S a )i (6.5)
(x i )max = ϕ i ⋅ yi (t )max = ϕ i ⋅ Qi ⋅ 2
ω i
6-7
Desplazamientos modales relativos entre pisos
∆ ri = ( x r )i − ( x r −1 )i (6.6)
r: nivel
i: modo
pero
.. .. (6.8)
x ri = ϕ ri ⋅ x i (t )max
max
Como
..
x i (t )max = Qi ⋅ (S a )i
.. (6.9)
x ri = ϕ ri ⋅ Qi ⋅ (S a )i
max
(S a )i (6.10)
FriS = ϕ ri ⋅ Qi ⋅ (S a )i ⋅ mr = ϕ i ⋅ Qi ⋅ ⋅ wr
g
fuerza sísmica equivalente en función del peso wr
en el nivel r para oscilaciones y la fracción de g
en el modo i
Vectorialmente:
F1iS
ϕ iT ⋅ M ⋅ J
FiS = FriS = M ⋅ ϕ i ⋅ ⋅ (S a )i
ϕ iT ⋅ M ⋅ ϕ i (6.11)
FniS
Qi
FiS = M ⋅ ϕ i ⋅ Qi ⋅ (S a )i (6.11)
6-8
Cortante modal
En nivel “r”: n
Vri = F jiS (6.12)
j =r
En la base: n
Vi = FriS (6.13a)
r =1
o también: (S a )i
Vi = Wi ⋅ (6.13b)
g
m: modo
k: nivel
Fkm: fuerza sísmica equivalente para el modo m en el nivel k
Según la fórmula (6.10) y (6.13a):
n (S a )i Qi ⋅ (S a )i n
Vi = ϕ ri ⋅ Qi ⋅ wr ⋅ = ⋅ ϕ ri ⋅ wr
r =1 g g r =1
6-9
ϕ ri ⋅ wr Wk ⋅ φ km
ξ ri = n
= n
notación INPRES - CIRSOC 103
ϕ ri ⋅ wr Wi ⋅ φ im
r =1 i =1
F jiS (h j − hr )
n
M riv =
j = r +1
er
M rit = er ⋅ Vri corresp. al piso “r”
C.R. C.M.
6-10