Salas, F., Sáez, E., Ovalle, C. y Fernández, J. (2019).

Análisis dinámico de una excavación profunda

contenida mediante pilotes anclados en la grava de Santiago. 25, 76-82

Análisis dinámico de una excavación profunda contenida mediante pilotes

anclados en la grava de Santiago
Dynamic analysis of a deep excavation supported by anchored piles in Santiago gravel

Fecha de entrega: 10 de diciembre 2018

Fecha de aceptación: 15 de marzo 2019

Felipe Salas1, Estéban Sáez1,2, Carlos Ovalle3 y Juan Manuel Fernández4

1
Departamento de Ingeniería Estructural y Geotécnica, Pontificia Universidad Católica de Chile, Avda. Vicuña Mackenna 4860,
7820436, Macul, Chile, flsalas@uc.cl, esaez@ing.puc.cl
2
Centro de Investigación para la Gestión Integrada del Riesgo de Desastres CIGIDEN, Edificio Hernán Briones, 3er Piso. Av. Vicuña
Mackenna 4860, Macul, Chile
3
Research Institute of Mining and Environment RIME UQAT-Polytechnique, Department of Civil, Geological and Mining Engineering,
Polytecnique Montréal, Room A-356, 2900 Édouard Montpetit Blvd, Montreal, QC, H3T 1J4, Canada, carlos.ovalle@polymtl.ca
4
Pilotes Terratest Chile, Alonso de Córdova 5151, oficina 1401, Las Condes, Santiago, Chile, jmfernandez@terratest.cl

El sistema de contención más común en Santiago de Chile Discontinuous anchored piling support is one of the
para excavaciones profundas corresponde a pilas (o pilotes) most frequently used retaining systems in temporary
discontinuos y arriostrados lateralmente mediante anclajes. deep excavations in Santiago, Chile. The system has
Este sistema se ha comportado satisfactoriamente en presented satisfactory behaviour under static and
condiciones tanto estáticas como sísmicas. Los reducidos seismic conditions, with negligible deformations.
desplazamientos laterales observados sugieren que los These observations suggest that the design procedures
procedimientos de diseño y que los códigos vigentes and the local code requirements are conservative
son conservadores y que un mejor entendimiento del and a better understanding of the physical behaviour
comportamiento físico del problema podría ayudar a would lead to optimizing the design methodologies.
optimizar las metodologías de diseño. Este trabajo expone los This paper presents a static and dynamic analysis
resultados de un análisis numérico en elementos finitos para of a monitored pre-excavated pile-supported deep
un caso estático y dinámico, entregando una caracterización excavation. A finite element model was developed and
de incrementos de empujes sísmicos, desplazamientos post- calibrated based on static displacement measurements,
sísmicos permanentes y esfuerzos internos inducidos en los geophysical and laboratory characterization using the
pilotes. Para esto se utilizó el modelo constitutivo HS-Small, parallel gradation technique. The dynamic responses
calibrado en base a resultados de exploraciones geofísicas to earthquake induced lateral pressures, induced
en terreno y de ensayos de laboratorio usando la técnica de internal forces and post-earthquake permanent
gradación paralela. La respuesta dinámica fue estudiada displacement were studied for a suite of Chilean
para una gama de registros chilenos de aceleraciones, ground motions, scaled to be representative of current
ajustados con el fin de ser representativos con los códigos seismic design codes. Based on these results, some
de diseño vigentes en el país. A partir de estos resultados, conclusions and recommendations on designing
se exponen conclusiones y recomendaciones acerca de los procedures are provided.
actuales procedimientos de diseño.
Keywords: numerical modelling, dynamic analysis,
Palabras clave: modelación numérica, análisis dinámico, model HS-Small, internal forces
modelo HS-Small, esfuerzos internos

Introducción
La grava de Santiago es conocida por su alta competencia debajo de los 50 m de profundidad. Estas condiciones
geotécnica y por presentar una cota de nivel freático por promueven el uso de pilas discontinuas ancladas como

76
Salas, F., Sáez, E., Ovalle, C. y Fernández, J. (2019). 25, 76-82

sistema de contención en excavaciones temporales. Por ensayos triaxiales monótonos drenados (CID) a diferentes
otra parte, Chile se caracteriza por ser un país con alta confinamientos, cuyos resultados fueron comparados
actividad sísmica, lo que motiva a estudiar en detalle con ensayos realizados in situ, a gran escala, a 3.5 km
una verificación sísmica que ayude a evitar daños en de la excavación estudiada (Kort et al., 1979). Una vez
estructuras vecinas ante la eventual ocurrencia de un sismo que se verificó la representatividad del material escalado
significativo. en trayectorias triaxiales, se ejecutaron ensayos de corte
cíclico en la muestra homotética para lograr obtener
Usualmente, el diseño sísmico de una excavación temporal
curvas de degradación de rigidez y amortiguamiento del
se realiza utilizando métodos pseudo-estáticos donde los
material. Basados en esta caracterización, se calibró un
corrimientos o desplazamientos horizontales permanentes
modelo bidimensional en elementos finitos utilizando el
(Δeq) inducidos por un sismo, se relacionan implícitamente
software Plaxis®. Estos resultados permitieron establecer
con el coeficiente pseudo-estático (kh) escogido (Richards
conclusiones y recomendaciones para el proceso actual de
y Elms, 1979). La norma chilena NCh433 (2009) divide el
diseño.
país en 3 zonas de acuerdo a la máxima aceleración efectiva
(A0) esperada, siendo este parámetro representativo a un
evento sísmico con 10% de probabilidad de excedencia en
50 años con un período de retorno de 475 años (Riddell,
1995). La ciudad de Santiago pertenece a la zona II, lo que
implica una A0 igual a 0.3g. Las cargas sísmicas inducidas
se incorporan mediante un kh calculado en algunos casos
como una fracción de A0 (Manual de Carreteras, 2008), o
en otros casos dependiendo un valor de Δeq definido según
la importancia de las estructuras vecinas, variando Δeq
entre 5 y 40 mm según el código chileno NCh3206 (2010).
Con el fin de analizar este actual procedimiento de cálculo,
se estudió una excavación (ejecutada y monitoreada por
Pilotes Terratest S.A. Chile) ubicada en la comuna de Las
Condes, Santiago. La excavación alcanzó una profundidad
de 21.2 m y fue contenida mediante pilotes pre-excavados
con máquina y dos líneas de anclajes temporales. Uno de
los pilotes fue instrumentado mediante un inclinómetro
capaz de entregar un perfil de deformaciones horizontales
Figura 1: Caso de estudio: a) fotografía de la excavación, b)
en profundidad en las distintas etapas constructivas de la plano en planta del proyecto, c) sección del pilote instrumentado
excavación (ver Figura 1). y d) perfil de velocidad de corte Vs en profundidad

El sitio se caracterizó geotécnicamente mediante métodos

geofísicos activos y pasivos (MASW y ESPAC), y
Geometría de la excavación, mediciones en
mediante ensayos de laboratorio a muestras remodeladas terreno y caracterización del sitio
sobre material extraído en terreno. Debido a que el El sitio estudiado se caracteriza por presentar gravas
tamaño de las partículas que componen el terreno areno-arcillosas bien graduadas (GW-GC), de compacidad
es demasiado grande para ser ensayado en equipos alta y con partículas de más de 90 mm de tamaño máximo.
convencionales de laboratorio, se utilizó la técnica de El contenido de finos varía entre 2-3% con un índice de
gradación paralela (o curvas homotéticas) procurando plasticidad entre 5 y 30. La napa freática se ubica a 50 m de
conservar la representatividad de la muestra original en profundidad. Esta unidad geológica es de origen fluvial y
términos de resistencia al corte y dilatancia, limitando el se conoce como deposiciones del río Mapocho. Basados en
contenido de finos a 9% (de la Hoz, 2007). Se ejecutaron 2 la exploración geofísica fue posible distinguir 2 estratos: la

77
Salas, F., Sáez, E., Ovalle, C. y Fernández, J. (2019). Análisis dinámico de una excavación profunda
contenida mediante pilotes anclados en la grava de Santiago. 25, 76-82

segunda depositación, desde la cota cero hasta los 11 m de

profundidad, y la primera depositación, donde se alcanzó
Modelo de elementos finitos en Plaxis
Parámetros de modelación del suelo
una velocidad de propagación de onda de corte cercana
Se utilizó el modelo constitutivo Hardening Soil with
a 900 m/s (roca) a los 33 m de profundidad, (ver Figura
Small Strains (HS-Small) para la modelación de ambos
1d). Ambos estratos varían entre sí primordialmente en su
estratos de suelo. Este modelo es una extensión de la
grado de compacidad y cohesión.
formulación matemática del modelo Hardening Soil (HS)
Para la contención de la excavación se utilizaron pilotes que presenta una dependencia no-lineal hiperbólica para la
de hormigón armado, pre-excavados con máquina, de relación esfuerzo-deformación (Duncan y Chang, 1970),
diámetro 880 mm y espaciados a 3.2 m eje a eje. Estos además de incorporar endurecimiento tanto en trayectorias
fueron arriostrados mediante dos líneas de anclajes isotrópicas como de corte, lo que ayuda a predecir de
provisorios ubicados a la cota -4.5 y -13.0 m, la primera manera más realistas las deformaciones inducidas. El
y segunda línea respectivamente (Figura 1c). La primera modelo HS-Small ofrece la posibilidad de incorporar un
línea fue tensada a una carga de bloqueo de 880 kN, valor de rigidez al corte a muy pequeñas deformaciones
mientras que la segunda a 1245 kN. (Gmax) y degradación no-lineal ante ciclos de carga-
El escalamiento de las muestras, según el método de descarga (Hardin y Drnevich, 1972).
las curvas homotéticas, se hizo utilizando como base Mediante la herramienta de Plaxis® Soil-test se simularon
el material ensayado por Kort et al. (1979). Para ello se los ensayos efectuados en laboratorio, buscando el mejor
extrajo muestras de suelo a una profundidad de 10 m. El ajuste entre las curvas de esfuerzo-deformación y curvas
proceso de escalamiento se realizó procurando conservar de degradación reales con las modeladas (Figura 2). Para
el coeficiente de uniformidad (Cu) y de concavidad (Cc) complementar esta calibración de parámetros, se utilizaron
del material base, además de verificar que las densidades los resultados de geofísica cuya curva simulada se expone
relativas (densidad máxima y mínima) se mantuvieran en la Figura 1d. Los parámetros obtenidos se presentan en
constantes (de la Hoz, 2007). El contenido de finos se la Tabla 1.
limitó a un 9%, con el fin de que la matriz fina no dominara
en la resistencia al corte ni dilatancia del material (Rollins
et al., 1998).
Se ejecutaron 2 ensayos triaxiales monótonos drenados
(CID) a las muestras escaladas utilizando los mismos
confinamientos empleados por Kort et al. (1979) en sus
ensayos a gran escala (44 y 83 kPa) de acuerdo a ASTM
D5311 (2013). Se procuró confeccionar las probetas a la
misma densidad in situ del terreno. Una vez logrado un
correcto ajuste entre ambas curvas (material original a
gran escala y material escalado), se ejecutaron ensayos
cíclicos en probetas de 5 × 10 cm (torsión cíclica y
columna resonante) completando la caracterización del
material. Las curvas de degradación obtenidas fueron
comparadas con curvas de ensayos a gran escala y al
mismo confinamiento para gravas (Rollins et al., 1998).
Los resultados en términos de degradación de rigidez son
muy similares a los reportados por otros autores, mientras
Figura 2: (a) y (b) Comparación entre curvas de ensayo
el amortiguamiento observado es un poco superior al de la triaxial CID: material original (Kort et al., 1979), escalado,
literatura si la deformación angular es inferior al 0.003%. y simulaciones hechas en Plaxis según HS-Small. (c) y (d)
Resultados de ensayos de corte cíclico: material escalado, curvas
de referencia (Rollins et al., 1998) y simulaciones según HS-
Small

78
Salas, F., Sáez, E., Ovalle, C. y Fernández, J. (2019). 25, 76-82

Tabla 1: Propiedades para las capas de suelo usando el modelo para la fase dinámica se emplearon elementos Free-field
HS small (o campo libre) en los límites laterales y Compliant base
Estratos de suelo
en el extremo inferior. Esto último con el fin de evitar
Parámetros Unidad Primera Segunda
depositación depositación la reflexión de las ondas sísmicas dentro del modelo. Se
g’ kN/m3 22.5 22.0 utilizaron elementos de 15 nodos cuyo tamaño fue escogido
c’ kN/m 2
35.0 15.0 para garantizar un mínimo de 8 elementos por longitud de
f’ ° 45.0 45.0 onda, considerando un contenido de frecuencia relevante
ψ ° 15.0 10.0 desde 0 hasta 15 Hz y un Vs,30 = 600 m/s.
m - 0.51 0.55
E 50
ref
MN/m 2
275.0 80.0
E oed
ref
MN/m 2
190.0 130.0
Eur MN/m 2
900.0 350.0
Gmax MN/m 2
750.0 400.0
g0.7 - 0.0011 0.0008
υur - 0.2 0.3
K0 nc
- 0.43 0.65

Elementos estructurales
Los pilotes pre-excavados fueron modelados utilizando
elementos tipo Plate elasto-plásticos según diseño a flexo- Figura 3: Geometría de modelo de EF en 2D
compresión, incorporando el espaciamiento entre pilotes
de 3.2 m indirectamente en los parámetros de rigidez axial Selección y calibración de input dinámico
y flexural. Alrededor de los pilotes se definieron interfaces Uno de los aspectos claves para llevar a cabo un análisis
de suelo-estructura cuyo criterio de ruptura obedece dinámico, es la elección y calibración de una solicitación
al criterio de Mohr-Coulomb. La longitud libre de los dinámica representativa del sitio estudiado. Esta elección es
anclajes fue modelada mediante elementos Node-to-Node complicada debido a la alta variabilidad entre los distintos
Anchors donde la carga de tensado fue aplicada conforme eventos sísmicos. Una posible solución es considerar una
a la secuencia constructiva. La longitud de bulbo de los amplia gama de registros de aceleraciones. En este estudio
anclajes fue modelada mediante elementos Embedded- se analizaron registros medidos en los terremotos del
pile row con el fin de lograr una mejor representación del Maule 2010 (Mw = 8.8) y de Valparaíso 1985 (Mw = 8.0),
traspaso de carga entre el bulbo y el suelo. todos medidos en roca (outcropping rock). Estos registros
fueron escalados e impuestos en la base del modelo. El
Modelo escalamiento consistió en amplificar iterativamente el
Se desarrolló un modelo bi-dimensional en deformaciones registro de aceleraciones de tal forma de obtener una A0
planas utilizando el software Plaxis®, como se muestra = 0.3g en superficie. De esta manera se aseguró que el
en la Figura 3. Se consideró una profundidad del modelo input dinámico impuesto fuese representativo con lo que
de 80 m con el fin de representar roca elástica (Elastic establece la norma chilena NCh433 (2009). Para mayores
bedrock) en su base, según el perfil de velocidad de ondas detalles acerca de esta calibración y los distintos espectros
de corte obtenido (Figura 1d). La extensión lateral del de respuesta obtenidos, revisar Salas (2018).
modelo fue de 115 m, definida con el objetivo de satisfacer
condiciones de campo libre, es decir, lo suficientemente Resultados
lejos de la perturbación generada por la excavación. En Análisis estático
la fase estática se utilizaron bordes con restricción del Para el análisis estático se consideraron 8 etapas diferentes,
movimiento en la dirección normal a estos, mientras que comenzando por una etapa de inicialización de esfuerzos

79
Salas, F., Sáez, E., Ovalle, C. y Fernández, J. (2019). Análisis dinámico de una excavación profunda
contenida mediante pilotes anclados en la grava de Santiago. 25, 76-82

empleando el método de cálculo Gravity loading de Plaxis® Análisis dinámico

para garantizar el equilibrio de esfuerzos en el modelo. Se monitorearon los desplazamientos de la primera línea
Luego de esta etapa se activó la sobrecarga equivalente de anclajes Δeq según la NCh3206 (2010) durante todo
representativa de calles, igual a 12 kPa (NCh3206, 2010). el evento sísmico en el modelo de EF, obteniendo un
Como tercera etapa se modeló la ejecución e instalación de diagrama tiempo-historia de dichas deformaciones. La
los pilotes y activación de los elementos de interfaz. Luego Figura 5 muestra este diagrama para los distintos registros
se respetó el procedimiento de ejecución de los anclajes, el empleados en sus 2 componentes (NS y EW). Se observa
cual consiste en excavar 0.5 m por debajo de la cota de cada que se alcanza un máximo valor instantáneo de 79 mm para
anclaje, instalar el cable y tensarlo a la carga de bloqueo el registro medido en Las Tórtolas. Los desplazamientos
estipulada. Se repitió esta secuencia hasta alcanzar el nivel residuales permanentes Δeq varían entre el rango de 10 a
de sello de fundación (21.2 m). 20 mm.
Las Figuras 4a y 4b muestran una comparación entre los
resultados obtenidos mediante el modelo numérico y lo
medido en terreno mediante un inclinómetro en distintas
etapas. Se observa un buen ajuste entre ambas curvas,
alcanzando un valor máximo de deflexión de 15 mm en
la etapa final de la excavación. Se concluye que el modelo
es capaz de reproducir apropiadamente las deformaciones
producidas en los pilotes inducidas por la excavación. La
Figura 4c muestra la distribución de empujes calculada en
la última etapa del análisis estático y es comparado con el
diagrama de empujes de diseño empleado y con perfiles
teóricos de empujes (activo de Rankine y geostático). Se
concluye que en la medida en que el pilote se deforma,
se produce una relajación de esfuerzos disminuyendo
los empujes horizontales que actúan sobre el mismo. Se
observa además una concentración de esfuerzos en las
cotas de anclaje, esto es debido al tensado de los mismos.

Figura 5: Diagramas tiempo-historia de desplazamienos del

primer anclaje: a) registro de Rapel (Maule 2010), b) Las
Tórtolas (Maule 2010) y c) UTFSM (Valparaíso 1985)

De cada uno de los registros utilizados, se registraron 10

perfiles instantáneos del incremento de empujes laterales
inducidos por la solicitación sísmica durante la duración
significativa con intensidad de Arias entre 5-95% (Arias,
1970). Dado que la envolvente de estos 10 perfiles no
representa un estado de carga envolvente actuando sobre
Figura 4: a) Comparación entre desplazamientos horizontales
el pilote, se calculó un promedio de estos 10 perfiles
medidos con inclinómetros y calculados con modelo de EF a una instantáneos los cuales se muestran en la Figura 6a
profundidad intermedia de excavación, b) profundidad final de con línea segmentada. Estos perfiles se comparan con
excavación y c) perfiles de empujes laterales en la etapa final de
la excavación: distribución geoestática, de Rankine, de diseño y recomendaciones de empujes para muros perfectamente
resultados de modelo de EF arriostrados según el código NCh433 (2009), y a su

80
Salas, F., Sáez, E., Ovalle, C. y Fernández, J. (2019). 25, 76-82

vez con el empuje teórico de Mononobe-Okabe M-O de un factor de seguridad determinado (DIN4085, 2011;
(Mononobe y Matsuo, 1929) siguiendo una distribución EAB, 2013). Los diagramas de esfuerzos internos luego de
triangular invertida (Seed y Whitman, 1970) utilizando kh aplicar esta metodología también son expuestos. El valor
= 0.15 y kv = 0.0. Se observa que el promedio de los perfiles máximo positivo obtenido del análisis dinámico alcanza
instantáneos calculados no se ajusta a ninguna de las dos los 937 kNm en el nivel de sello de fundación, mientras
teorías, esto es debido a que la distribución que se propone que el máximo negativo es igual a -875 kNm y se alcanza
en la NCh433 (2009) se utiliza en muros perfectamente entre la segunda línea de anclajes y el nivel de sello de
arriostrados, por ejemplo, muros de subterráneos y no en fundación. Se obtiene un incremento sísmico de un 157%
pilotes que admiten deformación. Por otro lado, la teoría en el diagrama de momento con respecto al perfil estático
de M-O asume un muro gravitacional, infinitamente rígido y de un 53% con respecto al diagrama obtenido según la
y deslizante. Situación muy distinta a la de un pilote NCh433 (2009). A esta profundidad el efecto en el diagrama
arriostrado mediante anclajes. de momento de la sobrecarga escogida para el método de
M-O es despreciable. El diagrama obtenido a partir de un
análisis LEM muestra un máximo valor positivo de 1117
kNm en la segunda línea de los anclajes. Valor que es un
20% mayor al calculado mediante un análisis dinámico.
Con respecto a los diagramas de corte, se observa una
tendencia similar a la de los momentos entre todas las
curvas, obteniendo un máximo valor igual a 1098 kN,
presentando un incremento sísmico de un 49% con respecto
al valor estático.
Se observó que la evolución de carga en los anclajes
durante las solicitaciones sísmicas, no superó un 14% de
incremento con respecto al valor inicial (última fase de
análisis estático).

Figura 6: Comparación entre metodologias descritas. a) Perfiles

Conclusiones
de incremento de empujes sísmicos, b) diagrama de momentos y El presente artículo presenta una investigación cuyo ob-
c) diagrama de corte jetivo fue caracterizar los incrementos de empujes latera-
les inducidos por un evento sísmico, los desplazamientos
Con respecto a los esfuerzos internos de los pilotes, permanentes inducidos, y los esfuerzos internos sísmicos
se calculó una envolvente del diagrama de momento y inducidos. Las principales conclusiones son:
de corte una vez terminada la etapa dinámica como se (a) El método de las curvas homotéticas (o gradación
muestra en las Figuras 6b y 6c. Estos valores resumen paralela) es capaz de mantener la representatividad
la máxima respuesta dinámica observada en los pilotes, del material original en términos de resistencia al
por lo que son estos valores los que deberían ser tomados corte y dilatancia en deformaciones menores a un
en cuenta para el diseño estructural. Estos diagramas 0.3%. Basados en esta metodología se propusieron
son comparados con el diagrama de esfuerzos estático curvas de degradación de rigidez y amortiguamiento
obtenido y con diagramas de esfuerzos obtenidos después para la grava de Santiago.
de aplicar la sobrecarga sísmica de la NCh433 (2009) (b) El análisis dinámico llevado a cabo con registros
y el empuje de M-O previamente descrito. El diseño de previamente calibrados con el fin de ser
los pilotes se realizó mediante el software alemán GGU- representativos con la normativa vigente, muestra
Retain®, el cual utiliza la metodología de equilibrio límite que los desplazamientos laterales permanentes
(LEM) para determinar el largo de los anclajes, en función obtenidos en la cabeza de la primera línea de anclaje

81
Salas, F., Sáez, E., Ovalle, C. y Fernández, J. (2019). Análisis dinámico de una excavación profunda
contenida mediante pilotes anclados en la grava de Santiago. 25, 76-82

se encuentran en el rango de 10 a 20 mm. Hardin, B.O. and Drnevich, V.P. (1972). Shear modulus and
damping in soils: measurement and parameter effects. Journal
(c) Se obtuvo un máximo momento sísmico de un
of the Soil Mechanics and Foundation Division 98(6), 603-624
157% mayor con respecto al valor máximo obtenido
en la etapa estática. Este valor es un 20% menor al Kort, I., Musante, H. and Fahrenkrog, C. (1979). In-situ
calculado mediante una metodología de LEM, actual mechanical properties measurements of gravelly soil used in
procedimiento de diseño. an interaction and foundation model for the Santiago Metro. 6th
Panamerican Conference on Soil Mechanics and Foundation
(d) Se observó que la evolución de carga en los anclajes
Engineering, Lima, Peru, vol. 2, 217-224
durante las solicitaciones sísmicas no superó un 14%
de incremento con respecto al valor inicial. Manual de Carreteras (2008). Instrucciones y criterios de diseño.
Volumen 3. Ministerio de Obras Públicas, Santiago, Chile
(e) Se concluye que el diseño de excavaciones profundas
en áreas de alta sismicidad debe ir respaldado de un Mononobe, N. and Matsuo, H. (1929). On the determination of
análisis dinámico con el fin de conocer con mayor earth pressure during earthquakes. World Engineering Congress,
exactitud las deformaciones sísmicas inducidas. Beijing, China, vol. 9, 177-185
Para el caso de estudio, este tipo de análisis hubiera NCh3206 (2010). Geotecnia- Excavaciones, entibaciones y
permitido reducir las cargas sísmicas sobre las pilas. socalzados- Requisitos. Instituto Nacional de Normalización
INN, Santiago, Chile

Agradecimientos NCh433 (2009). Diseño sísmico de edificios. Instituto Nacional

Esta investigación fue parcialmente financiada por el de Normalización INN, Santiago, Chile
Centro para la Gestión Integrada del Riesgo de Desastres Richards, R.J. and Elms, D. (1979). Seismic behavior of gravity
(CIGIDEN) CONICYT/FONDAP/15110017. El primer retaining walls. Journal of the Geotechnical Engineering
autor también fue parcialmente financiado por Pilotest Division 105(4), 449-464
Terratest S.A.
Riddell, R. (1995). Inelastic design spectra accounting for soil
conditions. Earthquake Engineering & Structural Dynamics
Referencias 24(11), 1491-1510
Arias, A. (1970). A measure of earthquake intensity. Hansej R.J.
Rollins, K.M., Evans. M.D., Diehl, N.B. and Daily, W.D. (1998).
editor. Seismic design for nuclear power plants. Cambridge,
Shear modulus and damping relationships for gravels. Journal of
USA
Geotechnical and Geoenvironmental Engineering 124(5), 396-
ASTM D5311 (2013). Standard test method for load controlled 405
cyclic triaxial strength of soil. ASTM International, West
Salas, F. (2018). Monitoring and dynamic analysis of a
Conshohocken, PA, USA
discontinuous pile-supported deep excavation in Santiago gravel.
de la Hoz, K. (2007). Estimación de los parámetros de resistencia Tesis MSc, Pontificia Universdiad Católica de Chile, Santiago.
al corte en suelos granulares gruesos. Tesis MSc, Universidad
Seed, H.B. and Whitman, R.V. (1970). Design of earth retaining
de Chile, Santiago
structures for dynamic loads. ASCE Specialty Conference on
DIN4085 (2011). Subsoil- calculation of earth pressure. Lateral Stresses in the Ground and Design of Earth Retaining
Deutsches Institut für Normung, Berlin, Germany Structures, vol. 1, Cornell University, Ithaca, NY, USA, 103–147
Duncan, J.M. and Chang, C. (1970). Nonlinear analysis of stress
and strain in soils. Journal of the Soil Mechanics and Foundation
Division 96(5), 1629-1653

EAB (2013). Recommendations on excavations. German

Geotechnical Society, Wiley, Germany

82