06 Manual-sismo-Ayudas de Diseño
06 Manual-sismo-Ayudas de Diseño
06 Manual-sismo-Ayudas de Diseño
Tabla 1
h(m) Vs(m/s) γs
2 177 1.42
2 203 1.42
2 181 1.42
2 211 1.42
2 298 1.42
2 354 1.47
2 311 1.47
2 463 1.47
2 530 1.47
2 463 1.52
2 435 1.52
2 482 1.52
2 456 1.52
2 487 1.52
2 508 1.52
8 456 1.80
8 445 1.80
8 436 1.80
8 430 1.80
8 427 1.80
8 425 1.80
∞ 1375 2.00
Para determinar el periodo del terreno, que está idealizado con una estratigrafía en que existen
estratos de alta velocidad intercalados con estratos de baja velocidad, se hace uso de la ecuación 1.2
de la sección 3.1.4.3, dando como resultado
Ts = 0.58s
H s = ∑ i =1 h i = 78 m
N
4 × 78
Vs = = 537.93 m/s
0.58
Por otra parte, el peso volumétrico del estrato equivalente se determina como el promedio de los
pesos volumétricos de los estratos pesados por el espesor de cada uno de ellos, es decir
SECCIÓN 3.1
∑
N
γ hi
i =1 i
γs = = 1.67
Hs
Así, el contraste de impedancias mecánicas se determina como
γ s Vs 1.67 × 537.93
ps = = = 0.33
γ o Vo 2.0 × 1375
Empleando el programa PRODISIS se tiene que la aceleración máxima del terreno en roca para
Guadalajara es
a or = 97 cm/s2
Fd = 97 / 400.0 = 0.2425
El factor de sitio se obtendrá por interpolación de la tabla 1.2 de la sección 3.1.5.2.1. Para ello, se
hará
f sd = 0.58 × ( 0.2425 )
1/ 2
= 0.2856
f sd
0.00 0.05 0.10 0.20 0.50 1.00 2.00 3.00
ps
1.000 1.00 1.00 1.00 1.00 1.00 1.00 1.00 1.00
0.625 1.00 1.08 1.23 1.12 1.00 1.00 1.00 1.00
0.250 1.00 1.18 1.98 1.60 1.40 1.12 1.00 1.00
0.125 1.00 1.20 2.64 2.01 1.69 1.32 1.00 1.00
0.000 1.00 1.22 4.51 3.17 2.38 1.75 1.19 1.00
Siendo Fi , j = 1.60 , Fi , j+1 = 1.40 , Fi +1, j = 1.12 y Fi +1, j+1 = 1.00 , con x = f sd = 0.2856 , y = p s = 0.33 . Las
distancias relativas en la interpolación se determinarán como
0.2856 − 0.20
dx = = 0.2853
0.50 − 0.20
y
0.33 − 0.25
dy = = 0.2133
0.625 − 0.25
Fs = 1.60 × (1 − 0.2853) × (1 − 0.2133) + 1× 0.2853 × 0.2133 + 1.40 × 0.2853 × (1 − 0.2133) + 1.12 × (1 − 0.2853) × 0.2133
= 1.4454
El factor de respuesta se obtendrá por interpolación de la tabla 1.3 de la sección 3.1.5.2.2. Los nodos
entre los que se realizará la interpolación lineal serán
Ts
0.20 0.50 1.00 2.00 3.00
ps
1.000 3.45 3.13 2.73 2.46 2.40
0.625 3.56 3.24 2.77 2.49 2.43
0.250 4.27 3.45 2.85 2.56 2.49
0.125 4.90 3.70 3.06 2.73 2.63
0.000 6.02 4.81 4.05 3.58 3.40
Siendo Fi , j = 3.45 , Fi , j+1 = 2.85 , Fi +1, j = 3.24 y Fi +1, j+1 = 2.77 , con x = Ts = 0.58 s, y = p s = 0.33 . Las
distancias relativas en la interpolación se determinarán como
0.58 − 0.50
dx = = 0.1600
1.00 − 0.50
y
0.33 − 0.25
dy = = 0.2133
0.625 − 0.25
Fr = 3.45 × (1 − 0.16) × (1 − 0.2133) + 2.77 × 0.16 × 0.2133 + 2.85 × 0.16 × (1 − 0.2133) + 3.24 × (1 − 0.16) × 0.2133
= 3.3186
Para determinar el factor no lineal Fn, primero habrá que determinar Fn′ , que se obtendrá por
interpolación de la tabla 1.4 de la sección 3.1.5.3.1, ya que se trata de un suelo granular. Los nodos
entre los que se realizará la interpolación lineal serán
Fd
0.00 0.10 0.20 0.30 0.40 0.50 0.75 1.00
ps
1.000 1.00 0.97 0.93 0.90 0.86 0.83 0.75 0.71
0.625 1.00 0.95 0.91 0.89 0.85 0.82 0.71 0.68
0.250 1.00 0.93 0.87 0.82 0.77 0.73 0.63 0.56
0.125 1.00 0.92 0.84 0.75 0.67 0.64 0.58 0.53
0.000 1.00 0.90 0.78 0.66 0.58 0.54 0.53 0.50
SECCIÓN 3.1
Siendo Fi , j = 0.87 , Fi , j+1 = 0.82 , Fi +1, j = 0.91 y Fi +1, j+1 = 0.89 , con x = Fd = 0.2425 , y = p s = 0.33 . Las
distancias relativas en la interpolación se determinarán como
0.2425 − 0.20
dx = = 0.4250
0.30 − 0.20
y
0.33 − 0.25
dy = = 0.2133
0.625 − 0.25
Fnl′ = 0.87 × (1 − 0.425) × (1 − 0.2133) + 0.89 × 0.425 × 0.2133 + 0.82 × 0.425 × (1 − 0.2133) + 0.91 × (1 − 0.425) × 0.2133
= 0.86
Para determinar el factor de velocidad Fv, primero habrá que determinar Fv′ , que se obtendrá por
interpolación de la tabla 1.6 de la sección 3.1.5.3.2, ya que se trata de un suelo granular. Los nodos
entre los que se realizará la interpolación lineal serán
Fd
0.00 0.10 0.20 0.30 0.40 0.50 0.75 1.00
ps
1.000 1.00 0.99 0.98 0.97 0.96 0.95 0.95 0.95
0.625 1.00 0.98 0.97 0.93 0.90 0.89 0.89 0.89
0.250 1.00 0.97 0.93 0.86 0.81 0.79 0.79 0.79
0.125 1.00 0.97 0.92 0.85 0.80 0.78 0.78 0.78
0.000 1.00 0.97 0.92 0.85 0.80 0.78 0.78 0.78
Siendo Fi , j = 0.93 , Fi , j+1 = 0.86 , Fi +1, j = 0.97 y Fi +1, j+1 = 0.93 , con x = Fd = 0.2425 , y = p s = 0.33 . Las
distancias relativas en la interpolación se determinarán como
0.2425 − 0.20
dx = = 0.425
0.30 − 0.20
y
0.33 − 0.25
dy = = 0.2133
0.625 − 0.25
Fv′ = 0.93 × (1 − 0.425) × (1 − 0.2133) + 0.93 × 0.425 × 0.2133 + 0.86 × 0.425 × (1 − 0.2133) + 0.97 × (1 − 0.425) × 0.2133
= 0.9115
Con estos factores, los parámetros del espectro de diseño se obtienen con las ecuaciones 1.13-1.18
de la sección 3.1.6. Estos son los siguientes
En la figura 1.1 se ilustra el espectro calculado con los parámetros reportados y haciendo uso de la
ecuación 1.12 de la sección 3.1.6, con β=1, ya que es un espectro de peligro uniforme para 5% de
amortiguamiento estructural, carente de efectos de interacción suelo-estructura.
ESTRUCTURAS TIPO 1
ESTRUCTURAS DE EDIFICIOS
ESTRUCTURAS TIPO 1: ESTRUCTURAS DE EDIFICIOS
MÉTODO ESTÁTICO
Se desea obtener los cortantes directos y de torsión en los marcos que conforman el edificio que se
muestra en la figura 3.1, bajo el sismo de diseño, actuando en la dirección X del edificio. El cálculo se
llevará a cabo empleando el método estático. Se acepta que el edificio tiene una fracción de
amortiguamiento igual al 5% del valor crítico.
SECCIÓN 3.3
W1 = 150 ton
W2 = 150 ton
W3 = 120 ton
Como puede apreciarse en la figura, se trata de un edificio de tres niveles, cuyas dimensiones en
planta y altura son:
b X1,3 = 12.0 m
b Y1,3 = 16.0 m
H = 9.0 m
Además, se supone que la ubicación del centro de masa en todos los niveles se ubica en la siguiente
coordenada (en planta):
C.G. = ( 8.0,7.0 ) m
Espectro de diseño
Para este cálculo se supone que el edificio se ubica en una zona de la República Mexicana, con un
tipo de terreno y estratificación tales que se tienen los siguientes parámetros:
ESTRUCTURAS TIPO 1: ESTRUCTURAS DE EDIFICIOS
a 0r = 80cm s 2
Ts = 0.5 s
ps = 0.4
Q = 4.0
a 0 = 0.1094
c = 0.3692
T´s = 0.5092 s
Ta = 0.1782 s
Tb = 0.6110 s
Tc = 3.5000 s
Fd = 0.2
f sd = 0.2236
Fs = 1.3948
Fr = 3.3740
Fnl = 0.9620
Fv = 0.9820
CONDICIONES DE REGULARIDAD
Para que un edificio sea considerado regular, debe cumplir las 11 condiciones que se indican en el
inciso 3.3.2.1 de la sección de recomendaciones:
b Y1 16.0
= = 1.33 < 2.5
b X1 12.0
4. En planta no se tienen entrantes ni salientes cuya dimensión exceda 20% de la dimensión de
la planta
5. Se acepta que en cada nivel se tiene un sistema de piso rígido y resistente.
6. Se acepta que no se tienen aberturas en los sistemas piso con grandes dimensiones, ni
ocasionan asimetrías significativas, ni difieren en posición de un piso a otro.
7. El peso de cada nivel nunca excede el del nivel inmediato inferior.
8. Las áreas de todos los niveles son iguales.
9. En todos los pisos, todas las columnas están restringidas en las dos direcciones.
10. Inicialmente se acepta que la rigidez y la resistencia al corte de cada entrepiso no excede en
más de 50% a la del entrepiso inmediatamente inferior. Esto se comprueba más adelante,
con respecto a la rigidez y la resistencia dependerá del diseño.
11. En todos los entrepisos, la excentricidad torsional es menor del 10% de su dimensión.
∑W n
cβ
Pn = Wn h n n =1
Ne
Q ' Rρ
∑W h
n =1
n n
R = R 0 = 2.5 (por ser estructura de concreto diseñada para Q = 4.0 ,inciso 3.3.1.1)
De aquí,
cβ
= 0.0888785
Q 'Rρ
ESTRUCTURAS TIPO 1: ESTRUCTURAS DE EDIFICIOS
Las fuerzas sísmicas y cortantes por cada nivel obtenidas de esta manera se resumen en la tabla 3.1:
Wn hn Wnhn Pn Vn
NIVEL
(ton) (m) (ton-m) (ton) (ton)
1 150 3 450 6.91 37.33
2 150 6 900 13.83 30.42
3 120 9 1080 16.59 16.59
Suma 420 2430
Los desplazamientos de cada nivel se pueden obtener de resolver la ecuación matricial P=KX,
empleando las fuerzas sísmicas obtenidas sin considerar el período natural. Dichos desplazamientos
son los siguientes:
SECCIÓN 3.3
⎡ 0.568⎤
X = ⎢⎢1.456 ⎥⎥ cm
⎢⎣ 2.125⎥⎦
1/2
⎡ Ne ⎤
⎢ ∑
2π n =1
Wn X 2n ⎥
Te = ⎢ N ⎥ = 0.78505 s
g⎢ e ⎥
⎢ ∑ Pn X n ⎥
⎣ n =1 ⎦
como Te > Tb , las fuerzas sísmicas considerando el período fundamental se deben calcular con la
ecuación 3.15 de la sección de recomendaciones:
a(β)
Pn = Wn ( α1h n + α 2 h 2n )
Q ' Rρ
donde:
Ne
∑W n
α1 = (1 − 0.5r (1 − q ) ) n =1
Ne
∑W h
n =1
n n
Ne
∑W n
α 2 = 0.75r (1 − q ) n =1
Ne
∑W hn =1
n
2
n
q = ( Tb Te ) = 0.7165
r
obteniendo:
α1 = 0.140255
α 2 = 0.007211
Entonces, el cálculo de las fuerzas sísmicas por nivel considerando el período de vibrar se resume en
la tabla 3.2
ESTRUCTURAS TIPO 1: ESTRUCTURAS DE EDIFICIOS
Wn hn Wnhn Wnhn2 Pn Vn
NIVEL
(ton) (m) (ton-m) (ton-m2) (ton) (ton)
1 150 3 450 1350 5.28 33.294
2 150 6 900 5400 11.97 28.017
3 120 9 1080 9720 16.05 16.051
Suma 420 2430 16470
⎡ 0.514 ⎤
X = ⎢⎢1.333 ⎥⎥ cm
⎢⎣1.962 ⎥⎦
Las fuerzas y, por ende, los cortantes directos en cada marco se obtienen como el producto matricial
de su rigidez por el vector de desplazamientos obtenido:
⎡ 2.099 ⎤ ⎡10.626 ⎤
P1,3 = ⎢⎢3.964 ⎥⎥ ton V1,3 = ⎢⎢ 8.527 ⎥⎥ ton
⎢⎣ 4.563⎥⎦ ⎢⎣ 4.563 ⎥⎦
⎡1.080 ⎤ ⎡12.043⎤
P1,3 = ⎢⎢ 4.038⎥⎥ ton V1,3 = ⎢⎢10.963⎥⎥ ton
⎢⎣ 6.925⎥⎦ ⎢⎣ 6.925 ⎥⎦
Estos cortantes pueden obtenerse distribuyendo el cortante total en cada nivel en forma proporcional
a las rigideces de entrepiso de cada uno de los marcos. Estas rigideces se pueden obtener con algún
método aproximado o dividiendo el cortante entre el desplazamiento relativo de cada nivel. Los
resultados obtenidos deben ser equivalentes, dependiendo de la aproximación del método empleado
para obtener las rigideces de entrepiso. Las rigideces de entrepiso, obtenidas a partir de los
desplazamientos relativos y los cortantes de cada marco son las siguientes:
⎡ 0.514 ⎤
X rel = ⎢⎢ 0.819 ⎥⎥ cm
⎢⎣ 0.630 ⎥⎦
⎡ 20.673⎤ ⎡ 23.431⎤
K E1,3 = ⎢⎢10.417 ⎥⎥ ton cm KE2 = ⎢⎢13.393 ⎥⎥ ton cm
⎢⎣ 7.245 ⎥⎦ ⎢⎣10.994 ⎥⎦
Como puede apreciarse, las rigideces de los entrepisos son menores que las de los pisos inferiores,
verificando la suposición realizada al revisar condiciones de irregularidad. Asimismo, en vista que los
cortantes sísmicos son mayores en los entrepisos inferiores, se supone que la resistencia será mayor
en dichos niveles.
Efectos de torsión
Para evaluar los cortantes debidos a los efectos de torsión, como primer paso es necesario obtener
los momentos torsionantes en cada nivel y posteriormente se evalúa la distribución de cortantes en
los marcos, en función de su ubicación y rigidez. Los momentos de torsión se calculan con la
ecuación 3.19 de la sección de recomendaciones:
Los momentos torsionantes en cada nivel, para las dos excentricidades se resumen en la tabla 3.3.
En esta tabla el signo positivo de los momentos es en el sentido de las manecillas del reloj.
Momento (ton-m)
Nivel en+ en-
1 69.92 -3.329
2 58.84 -2.802
3 33.71 -1.605
En cada nivel, los cortantes en los marcos se pueden obtener a partir de sus rigideces de entrepiso y
de su distancia al centro de torsión, como:
ESTRUCTURAS TIPO 1: ESTRUCTURAS DE EDIFICIOS
∑K
i =1,n
Ei di
Vti = M t n
∑K
i =1,n
Ei d i2
Es importante recalcar que la distribución de cortantes por efectos de torsión se debe obtener
considerando los marcos en ambos sentidos, ya que también contribuyen a la rigidez a torsión del
entrepiso. Para fines del ejemplo, se supuso que las rigideces de entrepiso de los marcos en la
dirección Y (marcos extremos y central) es 25% mayor que las de los marcos en la dirección X. Los
resultados obtenidos para cada nivel se muestran en las tablas 3.4 a 3.6.
Ke d Vt (en+) Vt (en-)
Marco eje Ke d K ed 2
(ton/cm) (m) (ton) (ton)
1 9.62 6.00 57.72 346.32 1.933 -0.092
2 11.49 0.00 0.00 0.00 0.000 0.000
3 9.62 -6.00 -57.72 346.32 -1.933 0.092
A 10.90 8.00 87.20 697.60 2.920 -0.139
B 13.08 0.00 0.00 0.00 0.000 0.000
C 10.90 -8.00 -87.20 697.60 -2.920 0.139
SUMA 2087.84
Ke d Vt (en+) Vt (en-)
Marco eje Ke d K ed 2
(ton/cm) (m) (ton) (ton)
1 4.42 6.00 26.52 159.12 1.536 -0.073
2 5.93 0.00 0.00 0.00 0.000 0.000
3 4.42 -6.00 -26.52 159.12 -1.536 0.073
A 5.45 8.00 43.60 348.80 2.525 -0.120
B 7.27 0.00 0.00 0.00 0.000 0.000
C 5.45 -8.00 -43.60 348.80 -2.525 0.120
SUMA 1015.84
Ke d Vt (en+) Vt (en-)
Marco eje Ke d K ed 2
(ton/cm) (m) (ton) (ton)
1 3.63 6.00 21.78 130.68 0.855 -0.041
2 5.17 0.00 0.00 0.00 0.000 0.000
3 3.63 -6.00 -21.78 130.68 -0.855 0.041
A 4.67 8.00 37.36 298.88 1.466 -0.070
B 6.54 0.00 0.00 0.00 0.000 0.000
C 4.67 -8.00 -37.36 298.88 -1.466 0.070
SUMA 859.12
SECCIÓN 3.3
Finalmente, se puede comentar que la distribución de cortantes en los marcos se puede obtener en
forma directa si se emplea un modelo de análisis tridimensional que incluya el grado de libertad de
torsión, ya que bastará aplicar las fuerzas sísmicas y momentos torsionantes en los centros de
gravedad de cada nivel. Debe recordarse que los cortantes y momentos en cada nivel se acumulan
desde la parte superior. Por ello, debe tenerse especial cuidado al trabajar con edificios irregulares,
ya que los centros de torsión de los distintos niveles no necesariamente son coincidentes en planta.
SECCIÓN 3.6
ESTRUCTURAS TIPO 3
MUROS DE RETENCIÓN
ESTRUCTURA TIPO 3: MUROS DE RETENCIÓN
Calcular la reducción en el cortante basal y en el momento de volteo, por unidad de longitud, que
produce un relleno de 10 m de espesor sobre un muro, cuando se considera la rotación de la base. La
velocidades de propagación de de corte (S) en el relleno es vsre = 250 m/s . Además, su peso
volumétrico es γ=1.6 ton/m3 y su relación de Poisson es νre=0.45. Considérese que el amortiguamiento
del relleno es de 5%, que el muro tiene una masa relativa μ = 0.15 y que el sistema muro-relleno
descansa sobre el sitio rocoso cuya aceleración máxima es a 0r = 0.25g . Considérese que la rigidez
rotacional relativa es dθ=1/2. El relleno tendrá un periodo
4H re 4 × 10.0
Tr = = = 0.16 s
vsre 250
De acuerdo con la sección de recomendaciones 3.1.6.3, para terreno rocoso se tiene que Ta=0.1 s y
Tb=0.6 s, por lo tanto, en este problema interesa conocer el valor de la meseta espectral es decir, la
segunda rama de la ec. 6.8 de la sección de recomendaciones, ya que Ta≤Tre<Tb. Por tratarse de un
sitio rocoso como terreno de desplante los valores de CQ y CM se tomarán de la tabla 6.2 de la misma
sección. Se determinarán dos valores de rigidez relativa: dθ=0 (sin rotación) y dθ=1/2 (con rotación)
ν re = 0.45 μ = 0.15 CQ CM
dθ= 0 1.605 0.973
dθ= 1/2 1.389 0.797
De acuerdo con las ecuaciones 6.6 y 6.7, de la sección de recomendaciones, los valores del cortante
basal y el momento de volteo deberán multiplicarse por los factores γ re H re2 a or y γ re H3re a or ,
respectivamente, ya que ρ =γ/g. De esta forma se llega a que el cortante basal y el momento de
volteo, se calculan como:
ΔQ = 8.64 ton / m
ΔM o = 70.4 ton
31
SECCIÓN 3.7
ESTRUCTURAS TIPO 4
CHIMENEAS, SILOS Y
SIMILARES
CHIMENEAS, SILOS Y SIMILARES
1/ 2
0.0155H 2 ⎛ M'⎞
Te = ⎜ ⎟
( 3D0 − DH ) E c ⎝ M ⎠
E c = 10,000 f c' = 10,000 200 = 141421 kg / cm 2 = 13.87 GPa
0.0155 × ( 60 )
2 1/2
⎛ 1696.97 ⎞
Te = ⎜ ⎟ = 0.8 s
( 3 × 9.25 − 7.5) 13.87 ⎝ 1475.21 ⎠
33
SECCIÓN 3.7
Aceleración Espectral
De acuerdo a las recomendaciones para la construcción del espectro de diseño sísmico contenidas
en este manual y por estar la chimenea en terreno rocoso, el espectro de diseño se caracteriza por
los siguientes valores:
a o = 0.14 c = 0.35 Ta = 0.1 s Tb = 0.6 s Tc = 2.0 s k = 2.0 r = 1.0
En vista de que se trata de una chimenea de concreto reforzado, se puede tomar un factor de
comportamiento sísmico Q=2 para propósitos de reducción de las ordenadas espectrales por
ductilidad.
0.45
⎛ 0.05 ⎞
Te < Tc ⇒ β=⎜ ⎟ = 1.26
⎝ 0.03 ⎠
r 1
⎛T ⎞ ⎛ 0.6 ⎞
Tb < Te < Tc ⇒ a ( β ) = β c ⎜ b ⎟ = (1.26)(0.35) ⎜ ⎟ = 0.33
⎝ Te ⎠ ⎝ 0.8 ⎠
34
CHIMENEAS, SILOS Y SIMILARES
a(β)
Ps = 0.15W (1 + 0.5r − 0.5rq )
QR
r 1.0
⎛ T ⎞ ⎛ 0.6 ⎞
q=⎜ b ⎟ =⎜ ⎟ = 0.75
⎝ Te ⎠ ⎝ 0.8 ⎠
0.33
∴ Ps = 0.15 × 16641.50 × (1 + 0.5 × 1.0 − 0.5 × 1.0 × 0.75 ) × = 370.69 kN
2 × 1.25
Para tomar en cuenta el efecto de la aceleración de rotación en la base, el valor espectral de dicha
aceleración es igual a
θ(Te )H e W
Vθ =
Q'R
r 1.0
β ⎛ Tc ⎞ Te 1.26 ⎛ 2.0 ⎞ 0.8
Q ' = 1 + (Q − 1) ⎜ ⎟ = 1 + (2 − 1) ⎜ ⎟ = 1.5
k ⎝ Te ⎠ Tc 2.0 ⎝ 0.8 ⎠ 2.0
En la tabla 7.2 se presentan los cálculos necesarios para determinar, según la ecuación 7.10, los
momentos de volteo de diseño en diferentes secciones de la chimenea.
35
SECCIÓN 3.7
Dovela
hn h 2n Wn Wnhn Wn h n
2
Pn Pθ (Pn2+Pθ2)1/2 Vn
2 2
(m) (m ) (t) (t-m) (t-m ) (t) (t) (t) (t)
10 57 3,249 1,044.40 59,530.80 3,393,255.60 410.13 191.97 452.84 823.53
9 51 2,601 1,148.35 58,565.85 2,986,858.35 392.85 183.88 433.75 1,257.27
8 45 2,025 1,263.09 56,839.05 2,557,757.25 370.94 173.62 409.56 1,666.84
7 39 1,521 1,389.60 54,194.40 2,113,581.60 343.83 160.94 379.63 2,046.47
6 33 1,089 1,528.85 50,452.05 1,664,917.65 310.93 145.53 343.30 2,389.77
5 27 729 1,680.85 45,382.95 1,225,339.65 271.44 127.05 299.71 2,689.48
4 21 441 1,849.52 38,839.92 815,638.32 225.25 105.43 248.71 2,938.18
3 15 225 2,034.87 30,523.05 457,845.75 171.47 80.26 189.33 3,127.51
2 9 81 2,238.85 20,149.65 181,346.85 109.54 51.27 120.94 3,248.45
1 3 9 2,463.42 7,390.26 22,170.78 38.83 18.18 42.88 3,291.33
SUMAS 16641.50 421867.98 15418711.80
⎧ NS
⎪ ∑ Wn
⎪α1 = (1 − 0.5r (1 − q ) ) Nn =1
⎪ S
⎪ ∑ Wn h n
2 a(β) ⎪
Pn = 0.85Wn ( α1h n + α 2 h n )
n =1
; ⎨
Q 'R ⎪ NS
⎪ ∑ Wn
⎪α 2 = 0.75r (1 − q ) NS
n =1
⎪
⎪⎩ ∑
n =1
Wn h 2n
36
CHIMENEAS, SILOS Y SIMILARES
NS
M = ( 0.75 + 0.25 h n / H
r
n )∑ Vm ( h m − h m −1 )
m = n +1
37
SECCIÓN 3.8
ESTRUCTURAS TIPO 5
TANQUES, DEPÓSITOS Y
SIMILARES
TANQUES, DEPÓSITOS Y SIMILARES
Las figuras 8.1 a 8.4 indican, para algunas relaciones de aspecto, los valores de las funciones
adimensionales C 0 (z) y C1 (z) necesarios para el cálculo de las presiones hidrodinámicas en las
paredes del depósito.
39
SECCIÓN 3.8
40
TANQUES, DEPÓSITOS Y SIMILARES
El recipiente y la estructura de soporte poseen las mismas características en las dos direcciones
ortogonales en que se debe analizar la estructura, de manera que el análisis sísmico del tanque
elevado se reduce solamente a una dirección.
41
SECCIÓN 3.8
Se supone que el tanque elevado será destinado al almacenamiento de agua cuyo peso volumétrico
es 9.81 kN / m3 , por lo que la masa del fluido almacenado es igual a
42
TANQUES, DEPÓSITOS Y SIMILARES
Los modos dominantes de vibración del tanque elevado se pueden determinar al resolver el problema
de valores característicos (K S − ω2n M S ) Zn = 0 , cuyas matrices de masa M S y rigidez K S se
definen, según las ecuaciones 8.30 y 8.31, como sigue:
Resolviendo el problema de valores característicos resultante se encuentra que las frecuencias y los
modos naturales de vibración son:
⎧ 1 ⎫
ω1 = 1.29 rad / s Z1 = ⎨ ⎬
⎩7.457 ⎭
⎧ 1 ⎫
ω2 = 3.29 rad / s Z2 = ⎨ ⎬
⎩−0.214 ⎭
El espectro de diseño para el terreno de cimentación y la localización del mismo se caracteriza por los
siguientes valores
Las características de la estructura de soporte son tales que puede tomarse un factor de
comportamiento sísmico Q = 3 para propósitos de reducción de las ordenadas espectrales por
ductilidad.
De acuerdo con la ecuación 5.1, sólo se justificará tomar en cuenta los efectos de la interacción
suelo-estructura cuando se cumpla que
Te H s
< 2.5
Ts H e
43
SECCIÓN 3.8
Z1T M s H 156,174.42
He = = = 21.04 m
Z1T M s J 7, 422.13
4.87 × 78
= 31.13 > 2.5
0.58 × 21.04
Por lo que se pueden despreciar los efectos de interacción suelo-estructura. Además, en tanques
elevados se puede despreciar la interacción líquido-recipiente, lo que se justifica aún más en
recipientes de concreto.
Los desplazamientos máximos que ocurren en el modo fundamental se determinan con la ecuación
8.32, esto es:
φ1 a ( T1 , β ) g
X1 = Z1
ω12 Q ' ( T1 ) R
Z1T M s J 7, 422.13
φ1 = = = 0.158
T
Z1 M s Z1 46,901.63
0.185
Tc 2.0 ⎛ 0.05 ⎞
T1 > Tc ⇒ λ = 0.45 = 0.45 = 0.185; β=⎜ ⎟ = 1.18
Te 4.87 ⎝ 0.02 ⎠
⎛T ⎞
r
⎡ ⎛ T ⎞ ⎤⎛ T ⎞
2 2
T1 > Tc ⇒ a ( T1 , β ) = β c ⎜ b ⎟ ⎢ k + (1 − k ) ⎜ c ⎟ ⎥ ⎜ c ⎟
⎝ Tc ⎠ ⎢⎣ ⎝ Te ⎠ ⎥⎦ ⎝ Te ⎠
⎛ 0.72 ⎞
0.58
⎡ ⎛ 2.0 ⎞ ⎤ ⎛ 2.0 ⎞
2 2
44
TANQUES, DEPÓSITOS Y SIMILARES
2 2
⎛T ⎞ ⎛ 0.72 ⎞
p = k + (1 − k ) ⎜ b ⎟ = 1.42 + (1 − 1.42 ) ⎜ ⎟ = 1.41
⎝ Te ⎠ ⎝ 4.87 ⎠
βp 1.18 × 1.41
T1 > Tc ⇒ Q ' ( T1 ) = 1 + ( Q − 1) = 1 + ( 3 − 1) = 3.16
k 1.42
T1 > Ta ⇒ R=2
Los desplazamientos máximos que ocurren en el modo superior se determinan con la ecuación 8.33,
esto es:
φ a ( T2 , β ) g
X 2 = 22 Z2
ω2 Q ' ( T2 ) R
Z T2 M s J 1,132.47
φ2 = = = 0.842
Z T2 M s Z 2 1,345.47
0.45
⎛ 0.05 ⎞
T2 < Tc ⇒ λ = 0.45 ; β=⎜ ⎟ = 1.51
⎝ 0.02 ⎠
r 0.58
⎛T ⎞ ⎛ 0.72 ⎞
Tb < T2 < Tc ⇒ a ( T2 , β ) = βc ⎜ b ⎟ = 1.51 × 0.45 ⎜ ⎟ = 0.39
⎝ Te ⎠ ⎝ 1.91 ⎠
r 0.58
β ⎛ Tc ⎞ Te 1.51 ⎛ 2.0 ⎞ 1.91
Tb < T2 < Tc ⇒ Q ' ( T2 ) = 1 + ( Q − 1) ⎜ ⎟ = 1 + ( 3 − 1) ⎜ ⎟ = 3.00
k ⎝ Te ⎠ Tc 1.42 ⎝ 1.91 ⎠ 2.0
T2 > Ta ⇒ R=2
45
SECCIÓN 3.8
Las fuerzas de inercia máximas correspondientes a los modos naturales de vibración del sistema se
obtienen con la ecuación 8.34 como sigue:
Pn = K s X n ; n = 1, 2
Las fuerzas cortantes en la base de la estructura de soporte asociadas a cada uno de los modos
naturales de vibración del sistema se obtienen sumando las fuerzas de inercia del modo
correspondiente. De esta forma se tiene que
Los momentos de volteo en la base de la estructura de soporte asociados a cada uno de los modos
naturales de vibración del sistema se obtienen sumando los momentos flexionantes originados por las
fuerzas de inercia del modo correspondiente. De esta forma se tiene que
392.26 915.64
M1 = 22.14 × × 15 + 22.14 × × (15 + 5.97 ) + 104.12 × (15 + 6.44 )
1,307.9 1,307.9
M1 = 2,656.97 kN − m
392.26 915.64
M 2 = 702.12 × × 15 + 702.12 × × (15 + 5.97 ) − 94.11 × (15 + 6.44 )
1,307.9 1,307.9
M 2 = 11, 448.59 kN − m
Cabe recordar que las fuerzas de inercia P11 = 22.14 kN y P12 = 702.12 kN se deben tanto a la masa
impulsiva como a la masa de la plataforma, razón por la cual se tuvieron que distribuir
proporcionalmente a cada una de estas masas a fin de calcular el momento de volteo en la base de la
estructura de soporte.
46
TANQUES, DEPÓSITOS Y SIMILARES
Para estimar la fuerza cortante y el momento de volteo basales debidos a la combinación de los
modos naturales de vibración del sistema se recurre al criterio de la raíz cuadrada de la suma de los
cuadrados de las respuestas modales. De esta forma se tiene que la fuerza cortante y el momento de
volteo de diseño en la base de la estructura de soporte resulta ser
(126.26 ) + ( 608.01)
2 2
V0 = = 620.98 kN
Por último, la estructura deberá diseñarse de acuerdo con lo especificado en la sección 3.8.7 de
recomendaciones. Cabe asentar que la fuerza cortante y el momento de volteo basales debidos a la
acción del componente vertical son nulos.
47
SECCIÓN 3.10
ESTRUCTURAS TIPO 7
PUENTES
PUENTES
Diseñar las pilas del puente de cuatro claros simplemente apoyados que se muestra en la figura 10.1.
Los claros tienen una longitud de 40 m y la longitud total del puente es de 160 m. La superestructura
consiste en una losa de concreto reforzado de 25 cm de espesor y 12 m de ancho total, apoyada sobre
tres trabes armadas de 2.4 m de peralte, separadas a una distancia de 5 m centro a centro. Las trabes
metálicas descansarán sobre apoyos de neopreno reforzado, con dureza Shore 60, con dimensiones
en planta de 40 x 40 cm, y con espesor igual a: 2.5 cm para los apoyos fijos y 7.3 cm para los apoyos
libres. En los extremos del puente se colocarán estribos que se consideran rígidos en las direcciones
transversal y longitudinal. Las pilas están formadas por cinco columnas de 10.0 m de altura libre,
separadas a una distancia de 2.5 m a ejes y unidas mediante un cabezal de 1.5 x 1.0 m. Se considera
que el amortiguamiento de la estructura es de 5%.
b) Media elevación
49
SECCIÓN 3.10
12000 mm
250
2400 mm
1000
10000 mm
250 250 250 250
600
229
25
51
14 1500
19
623
1000
2400
d) Secciones transversales
SOLUCIÓN:
50
PUENTES
b) Acción sísmica (sección 3.10.3 y sección 3.1). Para la construcción del espectro de respuesta se
tienen que definir los siguientes factores (sección 3.1.5):
Factores de terreno rocoso (sección 3.1.5.1.1).- De acuerdo con las coordenadas de localización de
la estructura, se determina la aceleración máxima del terreno en sitios rocosos. En este caso,
a 0r = 485cm / s2
a 0r 485
Fd = = = 1.213
a ref 400
Factores de comportamiento lineal del suelo (sección 3.1.6.1). El factor de sitio se obtiene a partir de
la aceleración máxima en roca y es igual a:
Fs Fnl a 0r
a0 =
981.0
Para obtener Fs es necesario determinar el parámetro fsd = Ts(Fd)1/2 = 0.5(1.213)1/2 = 0.55, y el contraste
de impedancias mecánicas ps = vs / v0 = 450 / 720 = 0.63. De la tabla 1.2 se obtiene Fs = 1.00. Por otro
lado, el factor no lineal (sección 3.1.5.3) se determina mediante:
T 0.5
Fnl = 1 - (1- Fn′ ) s = 1 - (1-0.68 ) = 0.89
Tc 1.5
En la ecuación anterior se empleo un factor Fn′ = 0.68 , tomado de la tabla 1.4 correspondiente a
suelos de tipo granular por lo tanto, la aceleración máxima del terreno es
La aceleración máxima del espectro se obtiene a partir del factor de respuesta Fr, mediante,
c = Fr a 0
Fr se obtiene con la tabla 1.3 para los valores de impedancia ps = vs / v0 = 0.63, y Ts = 0.5. El factor de
respuesta es igual a 3.24, por lo tanto, la aceleración máxima espectral es:
51
SECCIÓN 3.10
T 0.5
Fv = 1 - (1- Fv′ ) s = 1 - (1-0.89 ) = 0.96
Tc 1.5
donde, Fv´ se obtiene de la tabla 1.5 y es igual a 0.89. La velocidad equivalente es:
T 0.5
Ts′ = s = = 0.52 s
Fv 0.96
y,
Ta = 0.35Ts′ = 0.35 ( 0.52 ) = 0.18s
k = 2-0.5 = 1.5
S (T ) T
a = a e = 0.44 + (1.426 - 0.44 ) e para Te < 0.18s
g 0.18
S (T )
a = a e = 1.426 para 0.18 ≤ Te < 0.625s
g
0.5
S (T ) ⎛ 0.625 ⎞
a = a e = 1.426 ⎜ ⎟ para 0.625 ≤ Te < 3
g ⎝ Te ⎠
0.5 ⎡ ⎛ 3 ⎞ ⎤⎛ 3 ⎞
2 2
Sa (Te ) ⎛ 0.625 ⎞ ⎢
a= = 1.426 ⎜ ⎟ 1.5 + (1 −1.5 ) ⎜ ⎟ ⎥ ⎜ ⎟ para Tc ≥ 3s
g ⎝ 3 ⎠ ⎢ ⎝ Te ⎠ ⎥⎦ ⎝ Te ⎠
⎣
Espectro de desplazamientos de diseño (sección 3.1.6.6)
52
PUENTES
981Te 2 ⎛ T ⎞
Sd ( Te ) = ⎜ 0.44 + (1.426 − 0.44 ) e ⎟ para Te < 0.18s
4π 2 ⎝ 0.18 ⎠
981Te2
Sd ( Te ) = (1.426) para 0.18 ≤ Te < 0.625s
4π 2
0.5
981Te 2 ⎛ 0.625 ⎞
Sd ( Te ) = (1.426) ⎜ ⎟ para 0.625 ≤ Te < 2s
4π 2 ⎝ Te ⎠
0.5 ⎡ 2⎤ 2
981Te 2 ⎛ 0.625 ⎞ ⎛ ⎞
⎢1.5 + (1−1.5 ) ⎜ 2 ⎟
⎛ ⎞
⎥⎜ 2 ⎟
Sd ( Te ) = (1.426) ⎜ ⎟ para Tc ≥ 2s
4π 2 ⎝ 2 ⎠ ⎢ ⎝ Te ⎠ ⎥ ⎝ Te ⎠
⎣ ⎦
Por lo que respecta a la sobrerresistencia se aplicará un factor de 1.5 a las ordenadas del espectro de
diseño, como se especifica en la sección 3.10.3.
90
80
70
Desplazamiento (cm)
60
50
40
30
20
10
0
0.00 1.00 2.00 3.00 4.00 5.00 6.00 7.00 8.00 9.00 10.00
Periodo (s)
Figura 10.2. Espectro de desplazamientos
53
SECCIÓN 3.10
Si se admite que la respuesta del puente se debe fundamentalmente al primer modo de vibrar, se
puede aplicar, en primera instancia, el método de análisis estático. Para ello es necesario verificar
que se cumplen los requisitos establecidos en la sección 3.10.4.2.
En vista de que se cumplen los requisitos indicados, se procederá al análisis estático lineal (AEL). En
primer lugar, se idealiza la estructura como un sistema de un grado de libertad, donde la masa
corresponde a la masa de la superestructura, del cabezal y de la mitad superior de las columnas.
Dada la gran rigidez del cabezal con respecto a la rigidez de las columnas, se puede suponer que
éstas se encuentran empotradas en ambos extremos, razón por la cual la rigidez en dirección
transversal kt puede obtenerse con:
5 (12 ) E Ief
kt = = 76219 kg / cm 2
L3
Para determinar la rigidez se consideró que la pila es de 125cm de diámetro, que la resistencia del
concreto es de 300 kg/cm2, que el módulo de rigidez a flexión E = 265,000 kg/cm2 y que el momento de
inercia efectivo es 0.4 Ig (sección 3.10.4.6.1)
54
PUENTES
Sd = 7.14 cm.
Diámetro φ = 125 cm
Porcentaje de acero longitudinal ρ = 1.5%
Resistencia del concreto fc′ = 300 kg / cm 2
P
Porcentaje de carga axial = 0.04
f c′ Ag
Estribos circulares No. 4 s = 20 cm (en la zona de la articulación plástica)
Recubrimiento del acero r = 5 cm
A partir del análisis momento–curvatura se determinó la gráfica de la figura 10.3 y los valores de
curvatura siguientes:
500
450
400
350
Momento (t-m)
300
250
200
150
100
50
0
0 0.00005 0.0001 0.00015 0.0002 0.00025 0.0003 0.00035 0.0004
ϕ (rad/cm)
55
SECCIÓN 3.10
Por lo tanto, el desplazamiento de fluencia de la pila (ecuación 10.20), considerando que las
columnas se deforman en doble curvatura es:
Demanda máxima de ductilidad (sección 3.10.5.1.5).- Para el nivel de prevención del colapso, los
valores máximos de ductilidad son (tabla 10.4):
Δ pi ≤ Δ p ≤ Δ py
Δ pi
≤ μ Δ ≤1
Δ py
Posteriormente debe revisarse si la pila tiene la suficiente capacidad para alcanzar μΔ = 5. Para ello
se determina la rotación inelástica última de la pila. Dicha rotación se obtiene con (sección
3.10.5.1.4):
56
PUENTES
donde
Lap se obtuvo con la ecuación 10.18
Δ máx 17.8
μΔ = = = 7.6 > 5.0
Δy 2.33
Se concluye que la pila es capaz de desarrollar la demanda de ductilidad impuesta como límite por la
norma.
La fuerza cortante que actúa sobre cada una de las columnas (Vcj) es:
k cj
Vcj = V = 0.2 V = 0.2 k t Δ p = 0.2(76, 219)10.7 = 163,109 kg
5
∑ k ci
i =1
Para fines de la revisión por resistencia, las ordenadas del espectro deben reducirse por efecto de la
sobrerresistencia y la ductilidad. El factor de sobrerresistencia es R = 1.50 (sección 3.10.3.3), y el
factor de ductilidad puede considerarse como mínimo Q = 2.0 para pilas que forman parte de un
marco. Por lo tanto el momento para revisar la capacidad de la pila es:
815.5
M máx = = 271.8 t - m
1.5 x 2
De la gráfica momento – curvatura (figura E2) se observa que la resistencia de la pila es mayor. Para
evitar fallas de tipo frágil, la fuerza cortante debe corresponder al momento de diseño incrementado
por el factor de sobrerresistencia λ = 1.35 (sección 3.10.5.2.2),
57
SECCIÓN 3.10
λ M 1.35 x 271.8
V= = = 73.4 t
L 5
2
Es necesario revisar además la dirección longitudinal del puente y considerar el estado límite de de
servicio.
58
SECCIÓN 3.11
ESTRUCTURAS TIPO 8
TUBERÍAS
TUBERÍAS
La tubería esta enterrada con una cubierta de Hc de 3 ft en arena seca cuyo peso volumétrico es de
130 pcf y un ángulo de fricción interna de 33º. El ángulo de fricción en la interfaz suelo-tubo es de 28º y
el coeficiente de presión de tierras en reposo k0 es 0.5.
La curva esfuerzo deformación para la tubería API-5LX60 esta caracterizada por los siguientes
parámetros: Un modulo elástico inicial de 30000 ksi y coeficientes de Ramberg-Osgood de σ0 = 67.5
ksi, α = 30.26 y r =21.02.
Δx = 6.34 ft
Δy = 13.6 ft
Δz = 0 ft
3) Suponiendo un esfuerzo axial en la falla, σ a ,como 71.1 ksi . La correspondiente fuerza axial, P
es de 5,202 kips.
31
SECCIÓN 3.11
El esfuerzo flexionante ε b , se calcula como 0.58% (solo curvatura lateral) y el esfuerzo axial,
ε a , como 1.46% .Las tres verificaciones que se deben hacer son:
La longitud requerida de la tubería a cada lado de la falla para transferir la fuerza axial a al
suelo por fricción esta dada por la ecuación 11.34 como 1,830 ft . Este es tomado como la
longitud desde la falla al punto de anclaje, se asume para este ejemplo que no existen puntos
de anclaje efectivos dentro de esta distancia sobre cada lado de la falla.
ΔL c = 0.51 ft
ΔL s1 = ΔL s 2 = 2.08 ft
Los pasos para la el segundo ciclo interactivo están dados como sigue:
32
TUBERÍAS
33
SECCIÓN 3.12
ESTRUCTURAS TIPO 9
PRESAS
PRESAS
Considérese una cortina de gravedad que tiene una altura de H = 50 m y una base de B = 40 m. El
material de la cortina tiene las siguientes propiedades: Módulo de elasticidad E = 28 GPa, densidad
de masa ρ = 2400 kg/m3 y amortiguamiento ξ1 = 0.05. Se considera el nivel NAMO igual a H a = 47 m.
La cimentación rocosa tiene un módulo de elasticidad E c = 56 GPa, densidad de masa ρc = 2100
kg/m3, amortiguamiento ηc = 0.02 y velocidad de compresión Vc = 1200 m/s. Se requiere determinar si
la compresibilidad del agua o la flexibilidad de la cimentación son importantes así como el periodo
fundamental modificado y amortiguamiento modificado del sistema. Considérese la velocidad del
sonido del agua C = 1440 m/s.
El periodo fundamental del vaso está dado por la ecuación 12.9 de la sección de recomendaciones.
4H a 4 ⋅ 47
Ta = = = 0.13 s
C 1440
H 50
C1 = 6 = 6 = 7.5
B 40
H 50
T1 = C1 = 7.5 = 0.11 s
E/ρ 28 ⋅109 / 2400
Como
T1 0.11
= = 0.85 < 2 la compresibilidad del agua es importante (ecuación 12.11).
Ta 0.13
Como
2400 ⋅1200 − 1000 ⋅1440
α= = 0.67 < 0.9 la flexibilidad del fondo del vaso es importante
2400 ⋅1200 + 1000 ⋅1440
(Ecuaciones 12.12 ).
Como
E c 56
= = 2 < 4 la flexibilidad de la cimentación es importante (ecuación 12.13).
E 28
H a 47
Con los datos = = 0.94 , E = 28 GPa y α = 0.67 y las tablas 12.4 y 12.5 de comentarios, se
H 50
obtiene R a = 1.268 y ξa = 0.0415 .
101
SECCIÓN 3.12
Ec
Con los datos de = 2 y la tabla 12.6 de Comentarios se obtiene R c = 1.099 y con el dato adicional
E
de ηc = 0.10 y la misma tabla se obtiene ξc = 0.028 .
1 1 1 1
ξ = ξ + ξa + ξc = 0.05 + 0.0415 + 0.028 = 0.099
Ra (Rc ) 3 1
1.268 (1.099 )3
Considérese una cortina de sección homogénea con una altura de H = 46 m. El material de la cortina
tiene las siguientes propiedades: Módulo de corte máximo G max = 175 MPa, densidad de masa
ρ = 2040 kg/m3. Las variaciones del módulo cortante y amortiguamiento se muestran en la fig. E2.1.
La cortina está cimentada sobre terreno firme. Los espectros de diseño se muestran en la fig. E.2.2
para los amortiguamientos de 5, 10, 15 y 20%. La aceleración máxima del terreno es de 0.2g.
Primera iteración.
Paso 1. De la fig. E2.1 se determina G / G max = 0.4 y ξ =13%. EL módulo de corte es igual a
G = 0.4 × G max = 0.4 ×175 =70 MPa.
H 46
T1 = 2.62 = 2.62 =0.656 s
vs 183.5
Paso 4. De la fig. E2.2 y usando los valores del periodo y amortiguamiento T1 y ξ se determina la
aceleración espectral para el primer modo Sa1 = 0.2 × 1.3 = 0.26 g
102
PRESAS
H 46
γ c = 0.20. S = 0.20.
2 a1
0.26 × 9.81 = 0.068%
vs 183.52
Segunda iteración.
Paso 1. De la fig. E2.1 se determina G / G max = 0.36 y ξ =14%. EL módulo de corte es igual a
G = 0.36 × G max = 0.36 ×175 =63 MPa.
H 46
T1 = 2.62 = 2.62 =0.575 seg
vs 174.1
Paso 4. De la fig. E2.2 y usando los valores del periodo y amortiguamiento T1 y ξ se determina la
aceleración espectral para el primer modo Sa1 = 0.2 ×1.22 = 0.244 g.
H 46
γ c = 0.20. S = 0.20.
2 a1
0.244 × 9.81 = 0.07%
vs 174.12
La diferencia entre 0.07 y 0.68% está dentro del error permitido se concluye el proceso iterativo.
Se determina los periodos naturales de vibración para el segundo y tercer periodos de vibración
(Ecuación 12.34 de Recomendaciones)
46
T2 = 1.14 =0.301 segundos
174.1
46
T3 = 0.73 = 0.192 segundos
174.1
103
SECCIÓN 3.12
Usando los valores de T2 , T3 y ξ las aceleraciones espectrales se determinan de la fig. E2.2 como
Sa 2 = 0.2 ×1.22 = 0.32 g y
Sa3 = 0.2 ×1.47 =0.294g
Las aceleraciones máximas para cada modo se calculan con la ecuación 12.35 como
a corona,máx = a1,max
2
+ a 2,max
2
+ a 3,max
2
= 0.3902 + 0.3392 + 0.2522 =0.575g
104
SECCIÓN 3.13
AISLAMIENTO SÍSMICO Y
DISIPACIÓN DE ENERGÍA
AISLAMIENTO SÍSMICO Y DISIPACIÓN DE ENERGÍA
Diseñar el sistema de aislamiento con apoyos de neopreno con corazón de plomo (ALNP) para un
edificio de tres niveles de concreto reforzado destinado a oficinas, cuyo periodo sobre base fija es
Te=0.38 s. La planta y elevación del edificio se muestran en la figura 13.1 y las reacciones de cada
apoyo en la tabla 13.1.
a) Planta
30x55 cm
2.70m
30x55 cm
2.70m
30x55 cm
2.70m
45x45 cm 45x45 cm 45x45 cm 45x45 cm 45x45 cm
b) Elevación
El edificio se localiza en un terreno duro, alejado de una falla activa, con una aceleración máxima del
terreno de acuerdo con su ubicación geográfica de 485 cm/s2, periodo de vibración Ts = 0.5 y velocidad
de ondas de cortante Vs = 710 m/s.
101
SECCIÓN 3.13
SOLUCIÓN
Clasificación de la estructura
En vista de que se trata de una estructura convencional, en la que no existe riesgo de pérdida de un
gran número de vidas o pérdidas económicas cuantiosas, el edificio se clasifica como estructura del
grupo B.
Acción sísmica
Las ecuaciones que definen el espectro de diseño de acuerdo con el presente Manual son:
T
Sa ( Te ) = a 0 + ( β c - a 0 ) e para Te < Ta
Ta
Sa ( Te ) = β c para Ta ≤ Te < Tb
r
⎛T ⎞
Sa ( Te ) = β c ⎜ b ⎟ para Tb ≤ Te < Tc
⎝ Te ⎠
102
AISLAMIENTO SÍSMICO Y DISIPACIÓN DE ENERGÍA
⎛ Tb ⎞ ⎡
r 2⎤ 2
⎛ Tc ⎞ ⎛ ⎞
Sa ( Te ) = β c ⎜ ⎟ ⎢ k + (1 − k ) ⎜ ⎟ ⎥ ⎜ Tc ⎟ para Te > Tc
⎝ Tc ⎠ ⎢⎣ ⎝ Te ⎠ ⎥ ⎝ Te ⎠
⎦
a0 = 485 cm/s2
c = 1537 cm/s2
Ta = 0.18 s
Tb = 0.61 s
Tc = 2.0 s
k = 2.0 (por tratarse de un suelo rocoso)
b = factor de amortiguamiento
Te
Sa ( Te ) = 485 + (1537 - 485 ) para Te < 0.18s
0.18
0.5
⎛ 0.61 ⎞
Sa ( Te ) = 1537 ⎜ ⎟ para 0.61 ≤ Te < 2.0s
⎝ Te ⎠
0.5 ⎡ ⎛ ⎞
2⎤
⎛ ⎞
2
⎛ 0.61 ⎞ ⎢ 2.0 + (1 − 2.0 ) ⎜ 2.0 ⎟ ⎥ ⎜ 2.0 ⎟
Sa ( Te ) = 1537 ⎜ ⎟ para Te > 2.0 s
⎝ 2 ⎠ ⎢ ⎝ Te ⎠ ⎥ ⎝ Te ⎠
⎣ ⎦
Área de apoyos
El área mínima de los apoyos se determina con base en las cargas verticales que deben soportar. En
este caso, la descarga máxima para los apoyos obtenida del análisis del edificio puede agruparse en
tres casos distintos.
Para elegir las dimensiones iniciales de los apoyos, se considera que el esfuerzo admisible a
compresión del hule σmáx está comprendido entre 50 kg/cm2 y 100 kg/cm2 (Priestley, 1996). Por lo tanto,
el área mínima del apoyo considerando σmáx = 60 kg/cm2 es:
103
SECCIÓN 3.13
P 66, 700
Ab = = = 1,112 cm 2
σ max 60
Se propone emplear apoyos de 35 x 35 cm. Como estrategia de diseño se decide colocar apoyos con
las mismas dimensiones en todas las columnas para ahorrar costos de fabricación.
Para determinar las propiedades de rigidez efectiva (ecuación 13.28) y amortiguamiento equivalente
(ecuación 13.33) de los apoyos, es necesario conocer la demanda de ductilidad. Como esta variable
se desconoce al inicio del procedimiento de análisis, se puede proponer un desplazamiento máximo
inicial, mismo que tendrá que ser verificado al final del procedimiento. En caso de que el
desplazamiento propuesto al inicio sea suficientemente cercano al desplazamiento final, el proceso se
da por terminado, en caso contrario se repite el procedimiento utilizando como desplazamiento inicial
el valor obtenido en el ciclo anterior.
Se propone un desplazamiento máximo igual que la mitad de la dimensión del apoyo en la dirección
de análisis, es decir,
B 35
x max = = = 17.5 cm
2 2
104
AISLAMIENTO SÍSMICO Y DISIPACIÓN DE ENERGÍA
Conocido el desplazamiento inicial, se procede a determinar la fuerza de fluencia del apoyo Fy, que se
elige de manera que no se produzcan desplazamientos excesivos en condiciones de servicio,
ocurrencia de sismos moderados, efectos del viento, etc. Para fines del ejemplo se considera que la
fuerza de fluencia debe ser mayor que el 3% del peso total del edificio Fy = 17,342 kg, por lo tanto, el
área del elemento de plomo, considerando que se colocarán 15 apoyos y que el esfuerzo de fluencia
del apoyo es σy = 90 kg/cm2, está dada por (ecuación 13.17):
Fy 17, 342
A pl = = = 12.85 cm 2
σyN 90(15)
4A pl 4(12.9)
d= = = 4.04 cm
π π
Conocida la fuerza de fluencia se propone la ductilidad esperada en el apoyo, de manera que se logre
un trabajo eficiente del disipador. Para el ejemplo se propone una ductilidad μ = 12. El
desplazamiento de fluencia correspondiente es:
x 17.5
x y = max = = 1.46 cm
μ 12
Fy
xy =
GA b
10
Tr
Por lo tanto se propone un espesor total H = 15.3 cm, que corresponde a: Tr = 12.3 cm y 10 placas de
acero de 0.3 cm.
15.3 cm
45 cm
105
SECCIÓN 3.13
Rigidez equivalente
La rigidez equivalente de los aisladores se determina con la ecuación 13.28 que se muestra a
continuación:
1 + α ( μ − 1)
( k )ef = ( k )e
μ
Amortiguamiento equivalente
El amortiguamiento histerético del aislador se obtiene con la ecuación 13.33 del Manual:
⎛ 1⎞
2 (1 − α ) ⎜ 1 − ⎟ 2 (1 − 0.1) ⎛⎜ 1 − ⎞⎟
1
⎝ μ ⎠= ⎝ 12 ⎠ = 0.25
( ζ )his =
π ⎡⎣1 + α ( μ − 1) ⎦⎤ π ⎡⎣1 + 0.1 (12 − 1) ⎤⎦
El amortiguamiento del sistema conjunto edificio – aislador se obtiene con (ecuación 13.35):
Se supone que el edificio arriba del sistema de aislamiento se mueve como cuerpo rígido, de manera
que la suma de la ecuación anterior se reduce al nivel de aislamiento, es decir n=1. Por lo tanto, la
rigidez de la estructura (k)es, se determina con:
4π 2 m 4π 2 ( 589.26 )
( k )es = = = 161,100 kg / cm
T2 0.382
106
AISLAMIENTO SÍSMICO Y DISIPACIÓN DE ENERGÍA
Conocidas las propiedades equivalentes del modelo (rigidez y amortiguamiento), se realizó el análisis
de la estructura para obtener el desplazamiento final del aislador, los desplazamientos relativos y los
elementos mecánicos. Se empleó un modelo tridimensional utilizando el programa SAP2000
considerando elementos barra con propiedades elásticas para representar a los elementos del
edificio, pues se pretende que la estructura se mantenga dentro del intervalo elástico. Los aisladores
se representan por resortes lineales con las propiedades equivalentes que se obtienen en cada etapa
del análisis. El procedimiento seguido para determinar el desplazamiento final consiste de los
siguientes pasos:
1. La rigidez inicial de los resortes se determina con la ductilidad que resulta del valor de
desplazamiento inicial propuesto. Se realiza un análisis estático del edificio sobre resortes
equivalentes, empleando el coeficiente sísmico del espectro reducido por amortiguamiento
(factor b), para un periodo objetivo de la estructura sobre base aislada igual a tres veces el
periodo sobre base fija, es decir, T = 3*0.38 = 1.14 s. Como resultado del análisis se obtiene un
periodo de la estructura aislada Ti+1 = 1.11 s y un desplazamiento de los resortes xi+1 = 15.16
cm. Los resultados de esta primera etapa se muestran en el primer renglón de la tabla 2.
2. Como el desplazamiento inicial (17.5cm) difiere del desplazamiento obtenido del análisis
(15.16cm), se determina nuevamente la ductilidad de los apoyos y la rigidez y
amortiguamiento equivalentes, usando como desplazamiento inicial el valor final de la etapa
anterior. Se realiza el análisis estático con el coeficiente sísmico reducido por amortiguamiento
(0.478) y para un periodo de 1.11 s. El periodo de la estructura es ahora de 1.08s y el
desplazamiento de los resortes se reduce a 13.68 cm. Los resultados de la segunda etapa se
muestran en la tabla 13.2.
4. Para fines del ejemplo se considera que la aproximación es suficiente y los valores del análisis
serán los que se utilicen para fines de diseño.
107
SECCIÓN 3.13
Una vez encontrada la convergencia se decide si los valores de diseño finales son adecuados, con lo
que se acepta el aislador propuesto. En caso contrario se reinicia el análisis con un nuevo aislador de
base.
En la figura 13.4 se presenta el espectro para un amortiguamiento del 5% y el espectro reducido por
el factor de amortiguamiento b=0.471. También se indican los valores de seudoaceleración de diseño
para el edificio sobre base fija (1537cm/s2) y para el edificio con aisladores (547 cm/s2).
Figura 13.4. Espectros de diseño para edificio con base fija y para edificio con aisladores.
108
SECCIÓN 3.14
ESTRUCTURAS TIPO 11
TORRES DE
TELECOMUNICACIÓN
TORRES DE TELECOMUNICACIÓN
1.40
1.20
Modo 1
1.00
Modo 2
0.80 Modo 3
Ft
0.60
0.40
0.20
0.00
0 1 2 3 4 5 6 7 8 9 10
DO/DL
101
SECCIÓN 3.14
1.5
1.0 C1
0.5 C2
0.0 C3
Ci -0.5
-1.0 C4
-1.5 C5
-2.0 C6
-2.5
0.0 0.2 0.4 0.6 0.8 1.0
Altura relativa
0.2
C2
0.1
C3
Ci
C5
0.0
C6
-0.1
0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8
Altura relativa
102
TORRES DE TELECOMUNICACIÓN
1.5
C1
1.0
0.5 C2
0.0 C3
Ci
-0.5 C4
-1.0 C5
-1.5
C6
-2.0
0.0 0.2 0.4 0.6 0.8 1.0
Altura relativa
0.2
0.1 C2
C3
Ci 0.0
C5
-0.1 C6
-0.2
0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9
Altura relativa
103
SECCIÓN 3.14
1.5
C1
1.0
0.5 C2
0.0 C3
Ci
-0.5 C4
-1.0 C5
-1.5
C6
-2.0
0.0 0.2 0.4 0.6 0.8 1.0
Altura relativa
0.2
C2
0.1
C3
Ci
C5
0.0
C6
-0.1
0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8
Altura relativa
104
SECCIÓN 3.15
INTRODUCCIÓN
Los parámetros necesarios para evaluar un túnel en condiciones dinámicas se deben determinar de
acuerdo a lo expuesto en la sección de recomendaciones (sección 3.15); sin embargo, se pueden
emplear los valores presentados en la siguiente tabla para un diseño preliminar, considerando que
éstos corresponden a intervalos probables de valores o valores promedio que pueden presentar un
grado de imprecisión.
Tabla 3.15.1. Parámetros de algunas propiedades de las rocas (adaptada de González, 2002).
Peso Módulo de elasticidad Velocidad de propagación Relación de
Tipo de volumétrico de onda en compresión poisson
dinámico
roca
γ (kN/m3) Ed (MPa) vp (m/s) ν
Arenisca 22.56 – 25.50 4 903.3 – 54 917.2 1 400 – 4 200 0.24 - 0.31
Basalto 26.48 – 28.44 40 207.3 – 85 317.9 4 500 – 6 000 0.19 – 0.38
Caliza 22.56 – 25.50 7 845.3 – 97 085.8 2 500 – 6 000 0.25 – 0.33
Diabasa 28.44 58 839.9 – 96 105.2 5 500 – 7 000 0.28
Dolomía 24.52 – 25.50 21 574.6 – 84 337.2 5 000 – 6 000 0.29 – 0.37
Gneiss 26.48 – 29.42 24 516.6 – 102 969.8 3 100 – 5 500 0.20 – 0.30
Granito 25.50 – 26.48 9 806.70 – 82 375.9 4 500 – 6 000 0.18 – 0.24
Lutita 21.57 – 25.50 9 806.70 – 68 646.6 1 400 – 3 000 0.25 – 0.29
v máx ⋅ Cs
a= (15.1)
3 ⋅ a máx ⋅ R
22
SECCIÓN 3.15
a=
( 0.12 ) ⋅ ( 2500.0 ) = 14.0766
3 ⋅ ( 2.21) ⋅ ( 3.2 )
⎛ 8 ⋅ a 2 − 4 ⋅ a +1 −1 ⎞
φc = arcsin ⎜ ⎟ (15.2)
⎜ 4 ⋅ a − 2 ⎟
⎝ ⎠
⎛ 8 ⋅ (14.0766 )2 − 4 ⋅ (14.0766 ) + 1 − 1 ⎞
φc == arcsin ⎜ ⎟ = arcsin ⎛ 38.1139 ⎞ = 44.57º
⎜ ⎟
⎜ 4 ⋅ (14.0766 ) − 2 ⎟ ⎝ 54.3064 ⎠
⎝ ⎠
Para obtener la deformación unitaria máxima del revestimiento, debida a la onda de corte con un
ángulo de incidencia, se sustituyen los valores en la siguiente ecuación:
v máx a
ε=± ⋅ sin ( φ ) ⋅ cos ( φ ) ± máx2 ⋅ R ⋅ cos3 ( φ ) (15.3)
Cs Cs
0.12 2.21
ε=± ⋅ sin ( 44.57 ) ⋅ cos ( 44.57 ) ± 2 (
3.2 ) ⋅ cos3 ( 44.57 )
2500.0 ( 2500.0 )
Para estimar las deformaciones por ovalización, se determina la distorsión angular y la deformación
unitaria diametral máxima:
v máx
γ máx = (15.4)
Cs
0.12
γ máx = = 48 ×10−6
2500.0
⎛ ΔD ⎞
⎜ ⎟ = 2 ⋅ γ máx ⋅ (1 − ν m ) (15.5)
⎝ D ⎠máx
23
TÚNELES
⎛ ΔD ⎞
⎟ = 2 ( 48 × 10 ) (1 − 0.25 ) = 72 ×10
−6 −6
⎜
⎝ D ⎠ máx
Para determinar la deformación unitaria por ovalización:
v máx ⎡ t 1 E (1 − ν L 2 ) R ⎤
ε ov = ⋅ ⎢3 ⋅ (1 − ν m ) ⋅ + ⋅ m ⋅ ⋅ ⎥ (15.6)
Cs ⎣ R 2 E L (1 + ν m ) t ⎦
⎛ ΔD ⎞
ε ov > ⎜ ⎟ → 0.00008673 > 0.000072
⎝ D ⎠ máx
Solución B. Empleando la gráfica
R 3.2
= = 10.67
t 0.3
E m (1 − ν L 2 ) ⎛ 10.7 ⎞ (1 − 0.22 )
F' = ⋅ =⎜ ⎟⋅ = 0.399 ≈ 0.40
E L (1 − ν m 2 ) ⎝ 27.46 ⎠ (1 − 0.252 )
Empleando la gráfica (figura 15.1), se observa que no hay un valor para la curva F’=0.40, pero con la
ecuación de deformación unitaria normalizada se pueden determinar las líneas que corresponden a
las deformaciones por fuerza normal y por momento, respectivamente.
εov t 1 E m (1 − ν L 2 ) R t 1 R
= 3⋅ + ⋅ ⋅ ⋅ = 3 ⋅ + ⋅ F '⋅
v máx R 2 E L (1 − ν m ) t
2
R 2 t (15.7)
⋅ (1 − ν m )
Cs
La línea descendente corresponde a los valores que aporta la fuerza normal (primer sumando de la
ecuación) y la línea ascendente corresponde a los valores que aporta el momento flexionante
(segundo sumando de la ecuación). Se debe hacer notar que el primer sumando siempre será el
mismo para cualquier condición, por otra parte el segundo sumando depende de F’ pero la pendiente
de esta “curva linearizada” es la misma para cualquier caso.
Tomando en cuenta lo anterior, se obtienen dichas líneas, para el primer sumando se determina que
los valores de la relación, R/t, de 0.3 y 300 corresponden a valores de la deformación unitaria
24
SECCIÓN 3.15
normalizada de 10 y 0.01, respectivamente. Para el segundo sumando, valores de, R/t, igual a 0.01 y
50 corresponden a valores de 0.02 y 10, respectivamente, sin olvidar que F’=0.4. Con estas líneas, se
pueden determinar cada uno de los valores de la deformación inducida por la fuerza normal (igual a
2.1) y el momento flexionante (igual a 0.28).
ε ov
= 2.1 + 0.28 = 2.38
v máx
⋅ (1 − ν m )
Cs
F' =0.67
F' =0.15
F' =0.07
F' =0.5
F' =0.3
F' =0.2
F' =0.1
F' =1
10 F' =0.05
9
8
7
6 F' =0.03
5
4 F' =0.02
3 F' =0.015
Límite
2 F' =0.01
1.5
F' =0.007
Deformación unitaria normalizada (ovalización),
1 F' =0.005
0.9
0.8
0.7
0.6
/C )
0.5
m máx s
0.4
ε /((1-ν )v
0.3
ov
0.2
0.15
0.1
0.09
0.08
0.07
0.06
0.05
0.04
0.03
0.02
0.015
4.8 10.67
0.01
0.1 1 10 100
Relación, R/t
25
TÚNELES
Se verifica fácilmente que con esta geometría y con estas propiedades, la sección transversal no
cumple con la condición de deformación transversal debido a que es mayor a 2; además, si se
requiere estimar el valor de la deformación unitaria, εov:
v máx 0.12
ε ov = ( 2.38 ) ⋅ ⋅ (1 − ν m ) = ( 2.38 ) ⋅ ⋅ (1 − 0.25 ) = 0.0000857
Cs 2500.0
Dicho valor es muy aproximado al que se calculó con las fórmulas (0.00008673).
Por lo tanto, se concluye que para estas condiciones, el túnel no cumple con las condiciones
requeridas bajo solicitud sísmica. Además, se demuestra que las gráficas propuestas sirven para
revisar de una manera fácil y rápida si el túnel cumple o no con las condiciones de deformación por
ovalización.
Solución C: Comparando el valor de R/t con sus valores máximos y mínimos permisibles
⎛R⎞ R ⎛R⎞
⎜ ⎟ ≤ ≤⎜ ⎟ (15.8)
⎝ t ⎠ mín t ⎝ t ⎠máx
donde:
⎛R⎞ 2 − 4 − 6 ⋅ F'
⎜ ⎟ = (15.9)
⎝ t ⎠mín F'
⎛R⎞ 2 + 4 − 6 ⋅ F'
⎜ ⎟ = (15.10)
⎝ t ⎠ máx F'
2
0 < F' < (15.11)
3
Entonces, para este ejemplo se verifica fácilmente que F’=0.4 y cumple con la primera condición.
Finalmente, se obtienen los valores máximos y mínimos de la relación R/t, considerando que el valor
es igual a 10.67 para la geometría del túnel de este ejemplo.
26
SECCIÓN 3.15
Se verifica que no cumple con las deformaciones unitarias permisibles por ovalización, ya que no se
encuentra entre el rango de valores mínimo y máximo. Aunque no se obtienen los valores de la
deformación unitaria, es útil debido a que es muy sencillo evaluar si el túnel va cumplir con la
condición de deformación unitaria permisible.
Resumen de resultados:
Se incluyen el resumen de resultados de otros ejemplos, incluyendo el antes desarrollado:
Revestimiento
Concreto Acero
vmáx (m/s) 0.12 0.12 0.12 0.12
amáx (m/s2) 2.21 2.21 2.21 2.21
Cs (m/s) 2500 2500 2500 2500
EL (GPa) 10.70 10.70 10.70 10.70
Em (GPa) 27.46 27.46 205.94 205.94
Datos
27
TÚNELES
F' =0.5
F' =0.3
F' =0.2
F' =0.1
F' =1
10 F' =0.05
9
8
7
6 F' =0.03
5
4 F' =0.02
3 F' =0.015
Límite
2 F' =0.01
1.5
F' =0.007
Deformación unitaria normalizada (ovalización),
1 F' =0.005
0.9
0.8
0.7
0.6
F' =0.0025
/C )
0.5
m máx s
0.4
ε /((1-ν )v
0.3
ov
0.2
0.15
0.1
0.09
0.08
0.07
0.06
0.05
0.04
0.03
0.02
0.015
0.01
0.1 1 10 100
Relación, R/t
28
SECCIÓN 3.15
F' =0.67
F' =0.15
F' =0.07
F' =0.5
F' =0.3
F' =0.2
F' =0.1
F' =1
10 F' =0.05
9
8
7
6 F' =0.03
5
4 F' =0.02
3 F' =0.015
Límite
2 F' =0.01
1.5
F' =0.007
Deformación unitaria normalizada (ovalización),
1 F' =0.005
0.9
0.8
0.7
0.6
/C )
0.5
m máx s
0.4
ε /((1-ν )v
0.3
ov
0.2
0.15
0.1
0.09
0.08
0.07
0.06
0.05
0.04
0.03
0.02
0.015
0.01
0.1 1 10 100
Relación, R/t
Figura 15.3. Ábaco para ayuda de prediseño, para generar curvas con otros valores
de F’ o los valores que aportan tanto la fuerza normal como el momento flexionante
se hace de acuerdo al ejemplo antes descrito.
29