Propuesta para La Escogencia Del Coeficiente Dinámico para El Análisis Pseudoestático de Estabil
Propuesta para La Escogencia Del Coeficiente Dinámico para El Análisis Pseudoestático de Estabil
Propuesta para La Escogencia Del Coeficiente Dinámico para El Análisis Pseudoestático de Estabil
RESUMEN
El mtodo pseudoesttico para el diseo ssmico de taludes utiliza una fuerza
horizontal producto del peso de la masa a deslizarse y una aceleracin (reflejo del entorno
ssmico y condiciones geotcnicas del sitio), para simular la fuerza de inercia generada por el
sismo. Esta aceleracin deber ser una fraccin de la aceleracin pico correspondiente a la
solicitacin dinmica con fines de diseo estructural del sitio, ya que el mtodo simula con
una fuerza permanente, constante y en una direccin los efectos de la fuerza de inercia real
producida por el sismo que es transitoria, alternante, y aleatoria en magnitud y direccin.
Luego de analizar y comparar las propuestas existentes en la literatura para la determinacin
del coeficiente dinmico, se recomienda utilizar valores de 0.15 a 0.25 g para dicho
coeficiente, como se muestra en el Cuadro 2. Esta propuesta se presenta para discusin y
posible incorporacin en la prxima revisin del Cdigo de Cimentaciones de Costa Rica.
Estos valores propuestos son funcin de la zona ssmica del pas y del tipo de sitio descritos
en el Cdigo Ssmico de Costa Rica (2002). Se recomienda adicionalmente utilizar un
coeficiente vertical de 0.1 g en el anlisis pseudoesttico de estabilidad y obtener como
resultado al menos un valor de 1.15 de factor de seguridad para garantizar la estabilidad
global del mencionado talud.
ABSTRACT
The pseudostatic method for dynamic slope stability uses a horizontal force, the
product of the weight of the potential sliding mass and an acceleration value (that should
reflect the seismic background and the geotechnical properties of the site), to simulate the
inertia force generated by the earthquake. This acceleration should be a fraction of the peak
ground acceleration that characterizes the seismic input for the site for structural design
purposes, since this method models with a constant and permanent force the effects of the
real inertial force that is transient, alternating and random in magnitude and direction. After
examining and comparing the suggestions and guidelines in the existing technical literature
for the determination of the seismic coefficient, values of 0.15 to 0.25g are recommended for
the mentioned coefficient as shown in Table 2. This proposal is set forth to the committee in
charge of reviewing the next version of the Foundations Code of Costa Rica, for its
discussion and possible incorporation. These proposed values are a function of the seismic
zone and type of site described in the Costa Rican Seismic Code (CSCR 2002). A vertical
coefficient of 0.1 g is additionally recommended to be used in the pseudostatic analysis, as
well as a factor of safety of at least 1.15 to insure global slope stability
INTRODUCCIN
El Cdigo de Cimentaciones de Costa Rica y el Cdigo Ssmico de Costa Rica
vigentes no recomiendan valores especficos para la determinacin del coeficiente ssmico
requerido para la aplicacin del mtodo pseudo esttico de taludes. La mayora de los cdigos
de diseo de EEUU tampoco lo hacen y algunos asumen la participacin de un geotecnista
que determine entre otras cosas, la estabilidad ssmica de los taludes. Sin embargo, la ltima
versin del Eurocdigo 8 (Diseo Sismo Resistente de Estructuras, 2003) incluye una
recomendacin especifica para la escogencia de dicho valor, funcin del marco ssmico y de
las condiciones geotcnicas del sitio. A continuacin se comentan las consideraciones
pertinentes para la aplicacin del mtodo pseudoesttico de estabilidad de taludes y la
Terzaghi (en Seed, 1979) comenta las simplificaciones del mtodo ya que una fuerza
permanente tiene que representar el carcter temporal y multi-direccional del sismo.
Menciona que el factor de seguridad, podra tomar valores iguales o hasta menores que 1.0
sin fallar, o fallar a valores mayores de 1.0, dependiendo de la naturaleza de los materiales
del talud. Tambin menciona que los taludes que exhiben un mejor comportamiento son los
que tienen materiales arcillosos (no sensibles) y materiales granulares densos. Los ms
susceptibles a la inestabilidad son los materiales sueltos no cohesivos. Actualmente, la
mayora de los algoritmos clsicos propuestos para la solucin del mtodo de equilibrio
que las deformaciones se dan en el lapso en que la aceleracin del terreno es mayor a la
aceleracin crtica. Esta aceleracin crtica
se puede determinar a partir de un anlisis
de equilibrio lmite pseudo esttico, el cual
se ejecuta reiteradamente hasta que se
obtenga el valor del coeficiente ssmico que
corresponda a un factor de seguridad de la
unidad. Las deformaciones se calculan
luego utilizando la historia de aceleraciones
de diseo, a partir de donde se integran dos
veces los picos de aceleraciones que
excedan la aceleracin crtica.
Figura 2 a: Mtodo de Newmark, cuerpo libre del bloque deslizante
Para llevar a cabo el mtodo simplificado Makdisi y Seed (1977) se requiere de las
propiedades dinmicas del suelo y la fundacin y de una historia de aceleraciones en roca para
llevar a cabo el anlisis. La historia de aceleraciones puede ser un registro obtenido en un sitio de
roca congruente con el sitio de presa, escalado al nivel de aceleracin pico que dicte los
resultados del estudio de amenaza ssmica para el sitio, o puede ser generado sintticamente, de
acuerdo a las propiedades del movimiento del terreno obtenido de los resultados del estudio de
amenaza ssmica.
Los anlisis de respuesta dinmicos rigurosos, como los propuestos por Zienkiewicz
et al (1987, 1988), Prevost et al. (1985), Gazetas (1987), Finn (1986, 1999) pretenden evaluar
las deformaciones permanentes directamente, a partir de relaciones constitutivas que modelan
de una forma ms realista y efectiva la no linealidad del suelo. Se utilizan los modelos
Massing hiperblicos de comportamiento cclico y la terora cinemtica multi-superficie de
plasticidad para evaluar el comportmiento histertico. A pesar de ser estos mtodos
tericamente ms completos, no tienen la aplicabilidad prctica de los que utilizan una
relacin linear equivalente y una posterior integracin de los deformaciones, como en los
descritos en Makdisi y Seed (1977) Serff et al. (1976).
Las metodologas dinmicas ms utilizadas en el anlisis de deformabilidad de presas en
la prctica incorporan un anlisis de respuesta dinmica de la presa que posteriormente requiere
de un procedimiento adicional para evaluar las deformaciones permanentes ya que las relaciones
constitutivas que utilizan no son acopladas al fenmeno de generacin de presin de poro. Con
este enfoque se puede evaluar las deformaciones en zonas extensas dentro del cuerpo de la presa
y no solamente en una superficie de falla bien definida como en los mtodos derivados del de
Newmark. A continuacin se describe brevemente el procedimiento propuesto por Serff et al
(1976). Estos mtodos no son capaces de evaluar las deformaciones generadas provocadas por el
fenmeno de licuacin, por lo que hay que verificar mediante mtodos apropiados si esta
situacin se va a presentar y si es este el caso, recurrir a metodologas capaces de evaluar este
tipo de grandes deformaciones (para ampliar este tema ver Finn, 1999.)
El primer paso del procedimiento propuesto por Serff et al (1976) es la evaluacin de los
esfuerzos estticos iniciales (obtenidos de un anlisis de elemento finito esttico que modele el
proceso constructivo como una colocacin secuencial de capas y utilice modelos constitutivos
no-lineales estticos para el suelo. El segundo paso es el clculo de los esfuerzos cclicos
inducidos por el sismo (obtenidos de un anlisis de respuesta dinmica de elemento finito,
utilizando relaciones constitutivas dinmicas no lineales como el mtodo linear equivalente).
Luego hay que disear y realizar un programa exhaustivo de ensayos dinmicos de laboratorio
para evaluar las deformaciones cclicas (correspondientes a los esfuerzos cclicos determinados
de la respuesta dinmica) en los materiales de la presa superimpuestas a las condiciones
iniciales. Finalmente, la deformabilidad global del cuerpo total de la presa se estima al integrar y
compatibilizar las deformaciones que se evaluaron para los elementos individuales en el
laboratorio. En este ltimo paso se debe llevar a cabo un anlisis esttico de elemento finito con
fuerzas equivalentes en los nudos que correspondan a los potenciales dinmicos de
deformacin evaluados para cada uno de los elementos finitos del cuerpo de la presa a partir de
los resultados de los ensayos dinmicos de laboratorio. Seed en su publicacin clsica de 1979
utiliz este mtodo para predecir con muy buen grado de precisin las deformaciones de las
presas van Norman al fallar durante el terremoto de San Fernando de 1971, lo que le ha dado
gran credibilidad a esta metodologa en el campo ingenierl.
Existen varias versiones de este procedimiento pero aun las ms simplificadas son
laboriosas, costosas y consumen mucho tiempo. Estos anlisis por lo general se aplican a las
presas de alturas mayores de 30 40 m, donde es pertinente evaluar la deformabilidad de la
misma. El mtodo pseudoesttico es el que se aplica a presas y taludes de menor altura y lo hace
para efectos prcticos, el mtodo actualmente ms utilizado en la prctica para evaluar la
estabilidad ssmica de estas estructuras de suelo (y roca).
Los congresos mundiales de ingeniera geotcnica que tratan sobre el diseo de
estructuras de tierra sujetas a movimientos ssmicos (Finn 1986 y 1999, Gazetas 1991) reflejan
que los intereses de la investigacin actual estn enfocados en los mtodos para la evaluacin
permanente de deformaciones. Los temas de investigacin incluyen principalmente la licuacin,
los efectos tridimensionales, los modelos constitutivos plsticos, la generacin de presin de
poros y su acople a los modelos constitutivos. Otros temas de mucha importancia son los
propuestas para incorporar los avances de investigacin a la prctica actual, o para simplificar los
mtodos ms sofisticados para una utilizacin ms global. El mtodo pseudoesttico y la
determinacin del coeficiente dinmico son temas que se consideran superados en la actualidad,
por lo que son las publicaciones clsicas ms viejas las que mejor abordan el tema.
ESCOGENCIA DEL COEFICIENTE DINAMICO
Seed (1973) agrupa los enfoques para determinar el coeficiente dinmico en tres:
emprico, respuesta viscoelstica y respuesta de cuerpo rgido. El enfoque de cuerpo rgido
asume que el terrapln (o talud) se comporta como un cuerpo rgido y la aceleracin es
constante en toda su extensin y es igual a la aceleracin pico del suelo. Sin embargo, salvo
por presas muy pequeas, se ha determinado que las presas de matreriales sueltos se
deforman y se comportan de diferente manera cuando se sujetan a diferentes movimientos en
la base. El movimiento de un sismo es transitorio, alterna en direccin y es aleatorio en
magnitud. Seed (1973) observa que no existe evidencia para corroborar que estos efectos sean
equivalentes a utilizar una fuerza horizontal permanente cuya magnitud sea el producto de la
aceleracin pico y el peso de la masa potencial a deslizarse. El mtodo viscoelstico para
determinar el coeficiente ssmico se coment en prrafos anteriores pero se volvi obsoleto
en la medida en que se desarrollaron metodologas ms precisas y aplicables a la prctica
ingenieril para evaluar el comportamiento dinmico de masas de suelo.
El enfoque emprico deriva sus recomendaciones a partir de casos observados en la
prctica en los que el mtodo pseudoesttico pudo o no predecir los efectos de un sismo en
presas de tierra o taludes. Las experiencias en California y Japn referentes al
comportamiento de presas en sismos han sido muy ilustrativas. A continuacin se enumeran
criterios empricos basados en observaciones de campo, categorizaciones de las condiciones
particulares de cada caso, correlaciones entre el mtodo pseudo esttico y metodologas
simplificadas o ms rigurosas que evalan las deformaciones permanentes en taludes.
Terzaghi en 1950 (Seed, 1979) recomend los siguientes valores de coeficiente dinmico:
= 0.10 g - para sismos severos, escala IX Rossi Forel de Intensidad
= 0.25 g - para sismos destructivos, escala X Rossi Forel de Intensidad (el sismo de
San Francisco de M=8.25 de 1906 tuvo una IRF = X)
= 0.50 g - para sismos catastrficos
De la experiencia del comportamiento de presas sujetas a sismos moderados a grandes,
y posteriores anlisis pseudo estticos y de estimacin de deformaciones, Seed (1979) sugiri
que los coeficientes ssmicos kh, a utilizarse deberan ser 0.1 g para los efectos del campo
cercano de un sismo de M = 6.5 (aceleracin pico mxima 20% g), y 0.15 g para los
efectos del campo cercano de un sismo de M = 8.25 (aceleracin pico mxima 80% g), y el
factor de seguridad mnimo para ambos casos debera ser mayor que 1.15. Seed advirti que
el uso del mtodo pseudo esttico debera ser aplicado solo para presas pequeas y medianas
(de taludes de menos de 40 m de altura), construidas de acuerdo a tcnicas modernas de
colocacin y compactacin de materiales y constituidas por materiales que no sean
susceptibles a la generacin de presiones de poro altas o a la degradacin de resistencia al ser
sometidas a cargas cclicas (que pierdan, como mximo, el 15% de su resistencia total al ser
sometidas a carga cclica). Para presas grandes o construidas hace mucho tiempo sin cumplir
los estndares de hoy en da, este autor recomend llevar a cabo un anlisis de respuesta
dinmica de la presa para evaluar su estabilidad ante las solicitaciones dinmicas.
Marcuson (1981) concluy despus de analizar las publicaciones referentes al tema de
la Segunda Conferencia Internacional sobre Avances recientes en Ingeniera Geotcnica
Kh, g
Kv, g
FS min = 1.15
IX
0.4
0.40 - 0.50
0.29 - 1.15
SB
0.5 amax
Tipo de Sitio
SC
0.45 amax
+ - 0.5Kh
SD
0.4 amax
Kh, g
Kh, g
Tipo de Presa
Tierra
Enrocamiento
Zona Ssmica
Dbil
Fuerte
0.12 - 0.20
0.15 - 0.25
0.10 - 0.15
0.12 - 0.20
Kh, g
Arco
Gravedad
0.20 - 0.30
0.10 - 0.15
Kh, g
Eurocdigo 8, 2003
VIII
0.2
0.35 - 0.43
0.10 - 0.58
Marcuson, 1981
VII
0.1
0.30 - 0.37
0.03 - 0.28
0.24 - 0.40
0.12 - 0.20
Kh
(0.33-0.5) amax
(amplificacin del sitio se debe incorporarse en amax)
Kh
Kv
0.5 amax S
0.5 Kh
S = 1.0 -1.4
para considerar el tipo de sitio
Kh
0.5 amax, citando a Hynes-Griffin y Franklin
Tambin cita a Seed y Marcuson
Tipo de Sitio
Zona II
Zona III
Zona IV
S1
0.15
0.15
0.20
S2
0.15
0.20
0.20
S3
0.15
0.20
0.25
S4
0.15
0.20
0.25
Para esta propuesta se tom en cuenta la magnitud y forma espectral elstica del
movimiento del terreno del CSCR 2002 (Laporte et al., 2000).
La aceleracin pico efectiva es el parmetro con el que se caracterizan las zonas ssmicas
para condiciones en roca y representa una sacudida en el sitio con un perodo de retorno de
500 aos que equivale aproximadamente a una probabilidad de excedencia del 10% en 50
aos. Respectivamente las Zonas II, III y IV tienen una aceleracin pico efectiva de 0.2 g, 0.3
g y 0.4 g. En la Figura 3 se muestran las tres zonas ssmicas y en el Cuadro 3 se describen los
cuatro tipos de sitio del CSCR 2004. Las formas espectrales elsticas en las que se bas el
CSCR reflejan un decaimiento
menos pronunciado para los tipos
de sitio menos competentes
(Laporte et al., 2000) (dicho de
Zona II
otra manera, la grfica espectral
est desplazada hacia perodos ms
altos, donde puede afectar a los
estratos de suelo de 15 metros o
ms de profundidad).
Zona III
Zona IV
Se recomienda seguir la
recomendacin europea de incluir
una carga pseudoesttica vertical
para los casos en ambas direcciones con valor 0.1 g y obtener un factor de seguridad de al
menos 1.15. Se hace la salvedad que el mtodo pseudoesttico no ser aplicable para masas
que incluyan estratos licuables o con arcillas sensibles.
S2
S3
S4
1.
Ambraseys N. N. (1960): The Seismic Stability of Earth Dams. Proceedings, 2nd World
Conference of Earthquake Engineering, Tokyo, Japn.
2.
3.
4.
Building Seismic Safety Council (1997): NEHRP Recommendations Provisions for Seismic
Regulations for New Buildings and Other Structures, Part 1, Provisions Washington D.C.,
EEUU.
5.
Comisin Permanente del Cdigo Ssmico de Costa Rica: (2002): Cdigo Ssmico de Costa
Rica Editorial Tecnolgica de Costa Rica.
6.
7.
8.
9.
10
10. Gazetas, G. y Dakoulas P. (1991): Aspects of Seismic Analysis and Design of Rockfill
Dams. State of the Art Paper. Memorias del Second International Conference on Recent
Advances in Geotechnical Earthquake Engineering and Soil Dynamicas. Univ. Rolla,
Missouri, pp 1851-1888.
11. Gonzlez de Vallejo (coordinador y coautor) (2002). Ingeniera Geolgica Prentice Hall,
Madrid.
12. International Building Code Council (1998). International Building Code. Virginia, EEUU
13. Japanese National Committee on Large Dams (1988): Captulo Earthquake Resistant Design
Features of Dams pp1-30 en Earthquake Resistant Design for Civil Engineering Stuctures in
in Japan 264 pp, compilado por the Earthquake Engineering Committee of the Japan Society
of Civil Engineers publicacin de la Sociedad Japonesa de Ingenieros Civiles. Tokyo, Japn
14. Kramer, S. L. (1996): Geotechnical Earthquake Engineering. Prentice Hall.
15. Laporte, Mara (1999): Presas de Materiales Suelos Consideraciones de Diseo Material
didctico para el Seminario- Taller para establecer un Cdigo de Edificaciones para la
Repblica Dominicana. Banco Interamericano de Desarrollo. Santo Domingo, Repblica
Dominicana.
16. Laporte, M. Santana, G. y Sauter, F. (2000) : Propuesta para Demanda Ssmica y Espectros
Elsticos de Respuesta. Informe preparado por el Comit 1: Demanda Ssmica para la
Comisin Permanente del Cdigo Ssmico de Costa Rica.
17. Makdisi, F. I. y Seed, H. B. (1977): A Simplified Procedure for Estimating EathquakeInduced Deformations in Dams and Embakments, University of California Report No.
UBC/EERC-77/19.
18. Marcuson, W. F. (1981): Informe del Moderador sobre el tema Earth Dams and Stability of
Slopes under Dynamic Loads. Memorias del First International Conference on Recent
Advances in Geotechnical Earthquake Engineering and Soil Dynamicas. Univ. Rolla,
Missouri.
19. Monahan, P. Leveson, McQuarrie, Bean, Henderson y Sy (1999): Seismic Microzonification
Mapping in Greater Victoria, British Columbia, Canada. Specialty Conference on
Earthquake Engineering and Soil Dynamics III. ASCE, Special Geotechnical Publication N
75.
20. Seed H. B. (1973): Stability of Earth and Rockfill Dams during Earthquakes. Embankment
dam Engineering Casagrande Volume, John Wiley and Sons, Nueva York.
21. Seed H. B. (1979): Considerations in the Earthquake Resistant Design of Earth and Rockfill
Dams Gotechnique, Vol 29 No 3.
22. Serff, N.; Seed, H. B.; Makdisi, F. I.; y Chang C. Y. (1976): Earthquake Induced
Deformation of Earth Dams. Report No. EERC 76-4. Earthquake Engineering Research
Center, University of California, Berkeley.
23. Structual Engineers Association of California (1999): Recommended Lateral Force
Requirements and Commentary, Sacramento, California, EEUU.
24. Wang, J. G. Z. Q. and Law, K. T. (1994) Siting in Earthquake Zones, A.A. Balkema,
Rotterdam, pp115. (Se hace referencia a la Tabla 8.3.2 de la pg 94: Horizontal Seismic
Coefficient.
25. Zienkiewicz O. C. (1987): Finite Element Methods in Analysis and Design of Dams, boletn
No. 30a de la International Comission on Large Dams, ICOLD, Paris.
26. Zienkiewicz O. C. (1988): Dam Design Criteria, the Philosophy of Their Selection, boletn
No. 61 de la International Comission on Large Dams, ICOLD, Paris.
11