XII Congreso Chileno de Sismología e Ingeniería Sísmica ACHISINA 2019 Valdivia, 03 al 05 de Abril 2019
deslizamientos (cuando grandes masas de terreno se han desplazado hacia las partes inferiores de las
laderas) sino meramente las etapas iniciales en que el suelo se comporta como un medio continuo.
En este estudio se utiliza el programa Plaxis en su versión bidimensional (2D) [Plaxis 2018]. Este
programa ha sido elegido por sus prestaciones y por ser de amplio uso en ingeniería del terreno, tanto
en investigación como en aplicaciones profesionales.
3 Cálculos efectuados en este trabajo
Tal como se discute en la Introducción, el objetivo de este trabajo es analizar la exactitud del método
simplificado de Newmark comparando sus resultados con los proporcionados por procedimientos
supuestamente más exactos, como la simulación del comportamiento dinámico del suelo mediante
algoritmos basados en la mecánica del medio continuo. En este sentido, se selecciona una serie de
casos de estudio que consisten en taludes de distintas configuraciones geométricas y con suelos de
diferentes características (apartado 5); así mismo se selecciona un conjunto de acelerogramas sísmicos
(apartado 4). En los casos mencionados, se llevan a cabo cálculos simplificados usando el método de
Newmark y cálculos de mayor exactitud empleando Plaxis. En relación a estos últimos análisis, se
efectúan dos tipos de operaciones: determinación de las condiciones de seguridad al deslizamiento en
condiciones estáticas (fase 1) y cálculos dinámicos no lineales para obtener la evolución temporal de la
posición de cada punto del talud (fase 2). En los análisis con Plaxis el comportamiento no lineal del
terreno se representa con un modelo elasto-plástico (perfecto, es decir, sin endurecimiento) de Mohr-
Coulomb; éste se utiliza también como criterio de fallo (deslizamiento). Todos estos cálculos se
efectúan en dos dimensiones.
4 Acelerogramas considerados
Buscando una adecuada representatividad de las conclusiones obtenidas, los acelerogramas que
generan los deslizamientos han sido seleccionados intentado cubrir las situaciones más comunes. Con
esta idea, se consideran registros con diferentes características: para suelo blando y rígido, severos y
muy severos, breves y de mayor duración, y con y sin presencia de pulsos de velocidad; se discute a
continuación la manera en que se tratan cada una de estas cuatro características. El tipo de suelo
representa el contenido en frecuencias; se consideran (según la velocidad promedio de la ondas de
corte, vs,30) suelos B (rígido) y C (flexible) de acuerdo con la normativa europea [EN-1998, 2005]
(equivalen a C y D según la normativa americana). La severidad se relaciona con la magnitud de
momento del terremoto (Mw) y, sobre todo, con la aceleración máxima (PGA) y la intensidad de Arias
(IA) [Arias 1970] del acelerograma. La duración del acelerograma se cuantifica con el valor eficaz
[Trifunac, Brady 1975]. La relevancia de los pulsos de velocidad (efecto de directividad o de proximidad
de falla) se expresa mediante los índices PI [Baker 2007] y Ep [Zhai et al. 2013]. Dado que para cada una
de estas características se consideran dos acelerogramas, el número total de registros es 16 (2 × 2 × 2
× 2); éstos se describen en la Tabla 1 (registros de terremotos con magnitud de momento entre 6 y 7)
y en la Tabla 2 (registros de terremotos con magnitud de momento superior a 7).
5
XII Congreso Chileno de Sismología e Ingeniería Sísmica ACHISINA 2019 Valdivia, 03 al 05 de Abril 2019
Tabla 1 - Acelerogramas correspondientes a terremotos fuertes (6 d Mw d 7)
vs,30 IA Duración
No. Terremoto Año Mw Estación PGA [g] PI Ep Pulso
[m/s] [m/s] eficaz [s]
1 Morgan Hill 1984 6.20 Coyote Lake 561 (B) 0.713 2.889 4.08 (breve) 0.870 0.520 P
2 Parkfield 1966 6.19 Cholame (12) 409 (B) 0.060 0.056 28.99 (largo) 0.012 0.170 NP
3 Northridge 1994 6.69 Jensen 373 (B) 0.411 2.570 12.53 (largo) 0.999 0.470 P
4 Morgan Hill 1984 6.21 Gilroy Array 663 (B) 0.120 0.382 7.29 (breve) 0.063 0.24 NP
Imperial
5 1979 6.53 El Centro 264 (C) 0.317 0.867 8.21 (breve) 1.00 0.650 P
Valley
8.46
6 El Alamo 1956 6.8 El Centro 213 (C) 0.340 0.052 36.74 (largo) 0.080 NP
× 10-5
Northern
7 1954 6.50 Ferndale 230 (C) 0.163 0.546 17.33 (largo) 0.990 0.40 P
California
8 Parkfield 1966 6.19 Cholame (5) 290 (C) 0.444 0.863 6.50 (breve) 0.006 0.20 NP
Tabla 2 - Acelerogramas correspondientes a terremotos muy fuertes (7 d Mw)
vs,30 PGA IA Duración
No. Terremoto Año Mw Estación PI Ep Pulso
[m/s] [g] [m/s] eficaz [s]
9 Kocaeli (Turquía) 1999 7.51 Gebze 792 (B) 0.261 0.549 7.55 (breve) 0.999 0.820 P
10 Chi-Chi (Taiwan) 1999 7.62 CHY006 438 (B) 0.359 1.497 26.66 (largo) 0.0032 0.160 NP
11 Chi-Chi (Taiwan) 1999 7.62 TCU036 478 (B) 0.137 0.734 23.40 (largo) 0.999 0.560 P
12 Chi-Chi (Taiwan) 1999 7.62 CHY028 543 (B) 0.637 5.343 8.65 (breve) 0.016 0.120 NP
13 Kocaeli (Turkey) 1999 7.51 Yarimca 297 (C) 0.227 1.330 15.10 (largo) 0.980 0.590 P
14 Chi-Chi (Taiwan) 1999 7.62 CHY025 277 (C) 0.162 1.222 33.77 (largo) 0.055 0.20 NP
15 Chi-Chi (Taiwan) 1999 7.62 TCU065 306 (C) 0.790 7.667 28.55 (largo) 0.943 0.340 P
16 Duzce (Turquía) 1999 7.14 Bolu 294 (C) 0.739 3.723 8.55 (breve) 0.004 0.140 NP
En los acelerogramas representados en la Tabla 1 y la Tabla 2 el período de muestreo es 't = 0.005 s.
La velocidad de las ondas de corte representada en la Tabla 1 y la Tabla 2 se utiliza para determinar el
módulo de deformación del suelo mediante la conocida expresión ݒୱ = ට ܩ/ ɏ en que G es el módulo
de elasticidad transversal y U es la densidad (masa por unidad de volumen). En estos cálculos se supone
un coeficiente de Poisson Q = 0.30 y una densidad U = 2000 kg/m3.
5 Casos de estudio
La Fig. 3 describe el dominio bidimensional en que se lleva a cabo el análisis. Se considera un terreno
de características uniformes y situado sobre una base horizontal; la ladera tiene 20 m de altura con
pendiente uniforme representada por un ángulo T. Se supone que las zonas planas situadas en las
partes superior e inferior de la ladera tienen extensión indefinida (en otras palabras, puede
despreciarse la influencia de las irregularidades existentes); a efectos de garantizar la precisión y
estabilidad numéricas de los cálculos efectuados con Plaxis, se modelan 60 m de terreno a ambos lados
de la ladera. En todos los casos analizados, r = 0, es decir, no se considera presencia de agua.
