Antonio - Soriano - Martinez Cimentaciones en Puertos
Antonio - Soriano - Martinez Cimentaciones en Puertos
Antonio - Soriano - Martinez Cimentaciones en Puertos
TESIS DOCTORAL
Madrid, 2015
DEPARTAMENTO DE INGENIERÍA Y MORFOLOGÍA
DEL TERRENO
Madrid, 2015
UNIVERSIDAD POLITÉCNICA DE MADRID
Presidente D. __________________________________________________
Secretario D. __________________________________________________
Vocal D. ______________________________________________________
Vocal D. ______________________________________________________
Vocal D. ______________________________________________________
Calificación: ________________________
EL PRESIDENTE: VOCAL:
VOCAL: VOCAL:
EL SECRETARIO:
A Nuria, Isabel y C.Matías
AGRADECIMIENTOS
Cuando he pensado en las personas a las que tengo que agradecer la ayu-
da prestada para que esta tesis se haya escrito he visto que, realmente, son
una gran cantidad que, con su sabiduría, paciencia y cariño han conseguido
que haga este trabajo; desde las monjas que empezaron a enseñarme a leer
en la guardería hasta los profesores de los cursos de doctorado incluyendo a
todo el entorno (familia y amigos) que también han sido fundamentales en mi
formación.
No puedo concluir este apartado sin agradecer a las personas más impor-
tantes en mi vida: mi familia. El apoyo, ánimo y cariño que me han mostrado
durante toda mi vida. Gracias a mis padres Antonio (van tres veces agradecido
y me parecen pocas) y Mª Carmen, a mis hermanas Isabel, Bárbara y Lucia, a
mi señora esposa Nuria y a mis hijos Isabel y C. Matías.
Los diques de abrigo verticales son estructuras monolíticas que sirven para
reflejar el oleaje creando, de este modo, una zona abrigada en el lado tierra.
Son estructuras de contención de gravedad, es decir, su peso es el elemento
resistente fundamental.
Se conoce que en los suelos saturados, bajo cargas cíclicas, tiende a pro-
ducirse una elevación de las presiones intersticiales y una reducción de tensio-
nes efectivas así como una degradación del módulo de deformación tangencial
en función del número de ciclos de carga, pudiéndose producir el fenómeno
conocido como licuefacción (arenas) o movilidad cíclica (suelos más finos).
Una vez investigado el estado del arte sobre este tema, se propone un
procedimiento de cálculo en el que, utilizando el programa comercial FLAC, se
establecen las acciones cíclicas sobre un dique de abrigo vertical tipo con dis-
tintas condiciones de apoyo, aplicando, para el cimiento, un modelo de com-
portamiento tipo hiperbólico con generación de presiones intersticiales cuyos
parámetros pueden obtenerse de ensayos de campo y laboratorio.
Loads acting on these structures are dynamic (waves) and consist essen-
tially in a frontal thrust and an uplift pressure on the base, which both vary
over time.
It is known that pore pressures tend to rise in saturated soils under cyclic
loading and, consequently, there is a reduction of effective stress. A degrada-
tion of the shear modulus also occurs depending on the number of load cycles.
All of these effects can bring about the phenomenon known as liquefaction in
sands or cyclic mobility in fine-grained soils.
The aim of the thesis is to explore the possibility that current technology
provides to analyze the stability of vertical breakwaters founded on soft soils,
and to suggest a method to evaluate the dynamic stability in this type of
works.
To this end, a review has been made of procedures currently used to de-
fine the actions in calculations, the main models of dynamic behaviour of satu-
rated soils available and of calculation procedures.
Once the state of the art on this subject has been reviewed, a method of
calculation is proposed that uses the commercial program FLAC and is applied
to a typical vertical breakwater on a range of different foundation conditions.
For the foundation soil, a hyperbolic constitutive model with pore pressure
generation has been employed, whose parameters can be obtained from field
and laboratory tests.
ÍNDICE GENERAL
Págs.
III.1 INTRODUCCIÓN..................................................................... 87
III.1.1 Análisis estático previo .................................................... 87
-I-
Estudio dinámico de cimentaciones en obras portuarias
- II -
ÍNDICE GENERAL
- III -
ÍNDICE DE FIGURAS
ÍNDICE DE FIGURAS
Págs.
-V-
Estudio dinámico de cimentaciones en obras portuarias
- VI -
ÍNDICE DE FIGURAS
- VII -
Estudio dinámico de cimentaciones en obras portuarias
- VIII -
ÍNDICE DE FIGURAS
- IX -
Estudio dinámico de cimentaciones en obras portuarias
Figura nº IV.4. Sección tipo del dique de la nueva bocana del puerto
de Barcelona, el día de la rotura ..................................... 151
Figura nº IV.5.- Planta del dique actual ............................................. 153
Figura nº IV.6.- Perfil geológico-geotécnico por el eje del dique
actual ......................................................................... 154
Figura nº IV.7.- Datos de humedad natural y porcentaje de finos de
los suelos del delta del Llobregat. Zona de la nueva
bocana del puerto de Barcelona ...................................... 155
Figura nº IV.8.- Gráfico de plasticidad de Casagrande .......................... 156
Figura nº IV.9.- Valores de ensayos SPT en función de la
profundidad ................................................................. 157
Figura nº IV.10.- Ángulos de rozamiento cv, medidos en ensayo
triaxial CU ................................................................... 158
Figura nº IV.11.- Resistencia al corte simple sin drenaje ...................... 159
Figura nº IV.12.- Resultado de ensayos de corte simple sin drenaje,
en función de la humedad natural ................................... 159
Figura nº IV.13.- Gráficos de ajuste de la deformabilidad en la
interpretación de los ensayos triaxiales CU según la
teoría del modelo hiperbólico .......................................... 160
Figura nº IV.14.- Resistencia dinámica .............................................. 161
Figura nº IV.15.- Números de módulo en función del ensayo SPT .......... 163
Figura nº IV.16.- Relación de rotura (Rf) en función del ensayo SPT ....... 163
Figura nº IV.17.- Esquema del modelo para ensayos numéricos ............ 164
Figura nº IV.18.- Comparación de resistencia dinámica ........................ 166
Figura nº IV.19.- Relación σ'v/σ'v0 en función del número de ciclos......... 166
Figura nº IV.20.- Esquema de módulo de rigidez transversal y su
reducción con el incremento de la deformación
tangencial .................................................................... 167
Figura nº IV.21.- Gmax en función de SPT corregido. Ohta y Goto
(1976) y Seed et al. (1986) ........................................... 170
Figura nº IV.22.- Gmax en función de SPT corregido. Imai y Tonouchi
(1982) ........................................................................ 171
Figura nº IV.23.- Curvas de reducción de G para IP=0 ......................... 173
Figura nº IV.24.- Curvas de reducción de G para IP = 50% .................. 174
-X-
ÍNDICE DE FIGURAS
- XI -
Estudio dinámico de cimentaciones en obras portuarias
- XII -
ÍNDICE DE FIGURAS
- XIII -
ÍNDICE DE TABLAS
ÍNDICE DE TABLAS
Págs.
- XV -
CAPÍTULO I
INTRODUCCIÓN Y OBJETO
CAPÍTULO I.- INTRODUCCIÓN Y OBJETO
-1-
Estudio dinámico de cimentaciones en obras portuarias
1
altura, sobre la base del dique, del orden de H y la subpresión suele actuar
2
1
a una distancia del centro de la base igual a B.
6
El peso del espaldón hace que el momento del peso propio respecto al cen-
tro de la base sea hacia el lado mar compensando así del todo o parcialmente,
el momento que provoca la subpresión. Si esa compensación fuese exacta se
tendría que la desviación, respecto al centro de la base de apoyo, de las ac-
ciones sobre el cimiento sería aproximadamente:
1 EH
e
2 W
1
Es conveniente que esa excentricidad sea inferior a B para que el cajón
6
no tenga tendencia a despegarse en el lado mar de su apoyo. Esto conduce a:
L E 1 B
e< B <
6 W 3 H
Por otro lado, para asegurar que el cajón no desliza se suele exigir un coe-
ficiente de seguridad F > 1,5 suponiendo que el coeficiente de rozamiento ca-
jón banqueta es = 0,6. Con esto resulta:
E
< 0,4
W
-2-
CAPÍTULO I.- INTRODUCCIÓN Y OBJETO
Durante la retirada de la ola (seno del oleaje) el empuje del agua en el la-
do abrigado es mayor que en el lado mar y el peso efectivo del dique aumenta
al quedar mayor volumen de su cuerpo emergido (subpresión negativa respec-
to al mar en reposo). Esta situación puede ser condicionante cuando el espal-
dón es de gran volumen. En esta exposición previa se supone que la situación
condicionante es la correspondiente al empuje de la cresta del oleaje.
-3-
Estudio dinámico de cimentaciones en obras portuarias
su = · z
-4-
CAPÍTULO I.- INTRODUCCIÓN Y OBJETO
-5-
Estudio dinámico de cimentaciones en obras portuarias
Para soportar el peso del cajón tras el fondeo se requiere una resistencia al
corte sin drenaje:
p 480
su = = 4,7 t/m2
B ( 2) 20 ( 2)
480
su =
(20 t) ( 2)
2s uo ' B *
su,eq = 3 suo
6s uo ' B *
-6-
CAPÍTULO I.- INTRODUCCIÓN
resulta:
0,5 t 5 t 10
su,eq = 0,75 t t
1,5 t 5 4 t 13,3
t (t 10) 480
4 t 13,3 (20 t) ( 2)
-7-
Estudio dinámico de cimentaciones en obras portuarias
20 m
-8-
CAPÍTULO I.- INTRODUCCIÓN
Existen, en consecuencia dos tipo de rotura: una brusca (que también po-
dría denominarse frágil) que se asocia a la licuefacción del suelo y otra con
grandes deformaciones, sería una rotura dúctil, también denominada movili-
dad cíclica.
-9-
Estudio dinámico de cimentaciones en obras portuarias
- 10 -
CAPÍTULO I.- INTRODUCCIÓN
- 11 -
Estudio dinámico de cimentaciones en obras portuarias
Existe una diferencia fenomenológica clara entre los dos tipos de rotura cíclica
encontrados en la práctica. Los casos de licuefacción en los que la presión inters-
ticial crece hasta hacer el ensayo inestable y los casos de “movilidad cíclica” en
los que existe un exceso de deformación aunque la presión intersticial quede es-
tabilizada en un nivel alto pero mantenido, que aumenta poco con los ciclos.
Cuando existe una tensión de corte inicial “o” la amplitud de la onda cícli-
ca, c, que supone el suelo sin cambio de volumen es generalmente menor pe-
ro el valor de max = o + c es mayor que en el caso de tensión tangencial pre-
existente nula, esto es o = 0.
- 12 -
CAPÍTULO I.- INTRODUCCIÓN
Lógicamente las formas de estas curvas se espera que sean específicas pa-
ra cada tipo de suelo.
- 13 -
Estudio dinámico de cimentaciones en obras portuarias
Se quiere en este caso comparar el valor del empuje estático que rompe el
cimiento del dique de manera estática, Eest, con el que rompe el dique cuando
se produce quince veces el empuje dinámico como carga alternativa, E15.
Pero en un caso real, con geometría no tan simple, como la del ensayo de
laboratorio y cuando además existe una tensión de corte previa, la relación
E15/Eest normalmente será distinta.
- 14 -
CAPÍTULO I.- INTRODUCCIÓN
a) CÁLCULO PSEUDOESTÁTICO
b) CÁLCULO DINÁMICO
- 15 -
Estudio dinámico de cimentaciones en obras portuarias
- 16 -
CAPÍTULO II
ESTADO DEL ARTE
CAPÍTULO II.- ESTADO DEL ARTE
- 19 -
Estudio dinámico de cimentaciones en obras portuarias
- 20 -
CAPÍTULO II.- ESTADO DEL ARTE
- 21 -
Estudio dinámico de cimentaciones en obras portuarias
C 3 2vd
vd C1 C 2 vd
C 4 vd
Presumiblemente vd debería ser cero si es cero; esto implica que las
constantes están relacionadas de forma que C1 C2 C4 = C3.
vd
C vd
C1 2
- 22 -
CAPÍTULO II.- ESTADO DEL ARTE
12
10
8
evd
g·C1
6
0
0 20 40 60 80 100 120 140 160 180 200
Número de ciclos
- 23 -
Estudio dinámico de cimentaciones en obras portuarias
C 1 7600 Dr
2,5
1
Dr 15 N1 60
2
esto es,
0,4
En este caso, C2 se calcula como C 2
C1
ds ij dp ij
d ije
2Gmax 3K
- 24 -
CAPÍTULO II.- ESTADO DEL ARTE
II.3, se definen en el espacio de tensiones J-p con una función de forma g()
en p = plano constante ( es el ángulo de Lode). Se propone una superficie de
tipo flat-cap, p-pm = 0 para la simulación de la consolidación y entumecimien-
to. La variable pm es el máximo valor de p experimentado en la historia de
- 25 -
Estudio dinámico de cimentaciones en obras portuarias
p
Gmax G0 p a Ve
pa
donde:
V e
2,973 e 2
1e
1 e p
K pa
pa
R
m
Hr Gmax hr C f
1
R m
- 26 -
CAPÍTULO II.- ESTADO DEL ARTE
dependiente del estado para carga cíclica irregular. Para carga cíclica de ten-
sión controlada, las relaciones tensión de corte - deformación de corte obteni-
das en base a las ecuaciones anteriores conforman las bien conocidas reglas
de Masing. C() es una función de la medida del endurecimien-
to/ablandamiento. C() se incrementa (durante el endurecimiento) o se reduce
(durante el ablandamiento) en función de .
- 27 -
Estudio dinámico de cimentaciones en obras portuarias
- 28 -
CAPÍTULO II.- ESTADO DEL ARTE
- 29 -
Estudio dinámico de cimentaciones en obras portuarias
- 30 -
CAPÍTULO II.- ESTADO DEL ARTE
- 31 -
Estudio dinámico de cimentaciones en obras portuarias
- 32 -
CAPÍTULO II.- ESTADO DEL ARTE
Si bien los modelos basados en la teoría del estado crítico explican el efec-
to del índice de huecos en el comportamiento del suelo, estos modelos pare-
cen haber caído en desuso. El problema con estos modelos es que no reprodu-
cen bien la dilatancia existente en suelos densos que son, precisamente, los
de mayor interés práctico desde el punto de vista de la ingeniería de la cons-
trucción. Estos modelos tampoco presentan buenas predicciones de compor-
tamiento en arenas sueltas ni del comportamiento frente a licuefacción.
- 33 -
Estudio dinámico de cimentaciones en obras portuarias
2. Los suelos se "mueven" hacia esa línea de estado crítico con la de-
formación tangencial (γ).
Ψ = e - ec
- 34 -
CAPÍTULO II.- ESTADO DEL ARTE
Los datos recogidos en un amplio rango de arenas indican que existen infi-
nitas líneas de compresión normal. Esto hace que, para definir el estado del
suelo sea necesario disponer de dos parámetros: ψ y R.
- 35 -
Estudio dinámico de cimentaciones en obras portuarias
En el modelo la teoría del estado crítico, también existe una relación entre
la tensión media y la tensión tangencial, esta relación es de la forma:
q,c = M · m,c
- 36 -
CAPÍTULO II.- ESTADO DEL ARTE
- 37 -
Estudio dinámico de cimentaciones en obras portuarias
- 38 -
CAPÍTULO II.- ESTADO DEL ARTE
El modelo queda definido, según indica su autor, por los siguientes pará-
metros:
- 39 -
Estudio dinámico de cimentaciones en obras portuarias
Otro contraste, de más interés, se realizó con otro ensayo en modelo redu-
cido, en este caso de mayor tamaño. El modelo tiene una altura de 3,765 m y
recibe el impulso de una ola con valor máximo de 3,5 t/ml. La cimentación del
cajón se establece sobre una banqueta de escollera de 0,45 m de espesor que,
a su vez, apoya sobre una capa de arena de 1,40 m. Suponiendo ciertos pa-
rámetros de las relaciones constitutivas de la banqueta y de la arena, el pro-
grama en cuestión reproduce bien el desplazamiento horizontal máximo regis-
trado en el ensayo que es del orden de 0,3 mm.
- 40 -
CAPÍTULO II.- ESTADO DEL ARTE
- 41 -
Estudio dinámico de cimentaciones en obras portuarias
- 42 -
CAPÍTULO II.- ESTADO DEL ARTE
max
d A d A
A A
- 43 -
Estudio dinámico de cimentaciones en obras portuarias
D 2 0,212
y que las tensiones de corte en el punto de equilibrio (=0)
3
también lo son; eq/max = 2 1 . Para que esos datos clave sean otros es
a
(1, 2) =
1
2
2
2
1 s o 1 so
eq
eq
- Tensiones de corte de equilibrio “eq”. Punto de corte del ciclo con el eje
vertical. Valores de para = 0.
- 44 -
CAPÍTULO II.- ESTADO DEL ARTE
eq
eq =
max
so = 2,52 D2/3
La forma del ciclo de histéresis es importante. Pero hay otro aspecto que lo
es aún más cuando se pretende simular la rotura del terreno, que es la forma
de establecer la generación de presiones intersticiales. Esto es el mecanismo
que hace que ’v decrezca y que max crezca de ciclo en ciclo y que eventual-
mente se produzca una rotura.
Este aspecto es tratado con el modelo RG, de una forma tradicional basada
en el concepto del “daño acumulado” establecido hace ya mucho tiempo y con
gran aceptación práctica. El problema general es complejo, pero para una se-
rie de N ciclos iguales de tensión de corte, se puede suponer que tanto max
como la presión intersticial generada, u, vienen dados por las siguientes leyes:
C2 c
c C3
'c
max = C1 o N
'c
u = c (a1 + a2 ln N) + a3 o
Los ensayos utilizados por Rubén Galindo para desarrollar su modelo indi-
can que para ese suelo (limo arcilloso de la zona del Prat del Llobregat, Barce-
- 45 -
Estudio dinámico de cimentaciones en obras portuarias
C1 = 2 · e1,5 C2 = 2,8 C3 = 1 + e
Para otros suelos sería necesario realizar una potente campaña de ensa-
yos para determinar estas leyes de dependencia que podrían ser, como en ese
caso, función exclusiva del índice de poros.
C1 = C3 = 2 C2 = 2,8
a1 = 1 a2 = 0,7 a3 = 1,2
2
2,8 c
max = 2 c N 'c
'c
u = c (1 + 0,7 ln N)
nula, o = 0.
- 46 -
CAPÍTULO II.- ESTADO DEL ARTE
'
dP = d , d vP = (sen cv - ) dP
GP
2
P
G =G P
o 1 R f
f
- 47 -
Estudio dinámico de cimentaciones en obras portuarias
Los valores de GPo y Rf son los valores bien conocidos en las aplicaciones
geotécnicas del modelo hiperbólico.
np
'
GP
o = Kp pa
pa
Las constantes elásticas del modelo en condiciones estáticas son las con-
vencionales del bien conocido “modelo hiperbólico”.
- 48 -
CAPÍTULO II.- ESTADO DEL ARTE
- 49 -
Estudio dinámico de cimentaciones en obras portuarias
II.2.1 Temporales
Los temporales están formados por una serie de olas que varían en su al-
tura, período y dirección. La forma más sencilla de representar este tren de
olas irregular es mediante la altura, período y dirección de propagación que se
puedan considerar representativos de su variación estadística.
Para la zona del mar abierto es razonable aplicar el modelo gaussiano. Sin
embargo, cerca de la costa, donde empiezan a aparecer problemas de refle-
xión, difracción y rotura, este modelo deja de tener precisión. Esto es debido a
que el oleaje comienza a ser asimétrico respecto al nivel medio.
La ROM 1.0-09, Puertos del Estado (2009) indica una posible distribución
estadística de la altura de ola, según la ecuación de Rayleigh, que se considera
a continuación.
La altura de ola de paso por cero se puede definir como la suma aritmética
de la amplitud de cresta y de la amplitud de seno. Para el caso particular de
profundidades indefinidas y se considera que el oleaje es tipo swell la función
de densidad de la altura de ola es la función de Rayleigh, que se indica a con-
tinuación:
2
H
H
2 H e Hrms ,
p H 0
Hrms Hrms
2
- 50 -
CAPÍTULO II.- ESTADO DEL ARTE
2
H
H
Pr 1 e Hrms
, H 0
Hrms
2
Hq
Hq
e Hrms , Hq 0
q 1 F
Hrms
- 51 -
Estudio dinámico de cimentaciones en obras portuarias
mulas empíricas que permiten estimar los diagramas de presiones que actúan
sobre la estructura.
En 1919 Hiroi diseña un primer diagrama de presiones definido por una ley
rectangular desde el pie del paramento vertical hasta una altura de 1,25 veces
la altura de ola de diseño (H) o la altura del paramento vertical (se tomaría el
menor de ellos). El valor de la presión se define como sigue.
Pmax 1,50 w H
Esta fórmula toma como punto de partida la analogía con la presión hidro-
dinámica y está modificada por los resultados de presiones medidas en campo.
- 52 -
CAPÍTULO II.- ESTADO DEL ARTE
- 53 -
Estudio dinámico de cimentaciones en obras portuarias
1
dos (1/25 < h/L < , siendo L la longitud de onda y h la altura desde el pie
2
de la estructura al nivel del agua).
p 2 w h H 0
p1
h H 0
w H
p2
2 h
cosh
L
H2 2 h
0 coth
L L
Para el caso del seno de la onda se supone que la presión es nula a la cota
del nivel del mar en reposo (cota h) y tiene un valor p3 a la cota h+0-H. El
valor de la presión p3 es:
p 3 w H 0
No existe consenso sobre el valor de la ola de diseño que hay que utilizar.
Actualmente se recomienda considerar la ola máxima para el diseño de la
obra.
- 54 -
CAPÍTULO II.- ESTADO DEL ARTE
- 55 -
Estudio dinámico de cimentaciones en obras portuarias
Hb d s
Pm 101 w D ds
L hs D
- 56 -
CAPÍTULO II.- ESTADO DEL ARTE
2
H
w h b
P Hb 2
Rt m
3 2
3
H
w h b
P Hb h 2
Mt m
3 6
Los valores que se obtienen son del orden de 15-18 veces superiores a las
olas no rotas.
- 57 -
Estudio dinámico de cimentaciones en obras portuarias
H2
h0 coth k h
L
2
donde k k
L
h H h0 h p1 h2
Rc
2 2
Mc
h H h0 2 h p1 h3
2 6
1 H
Siendo p1
2 cosh k h
Además, Rundgren (1958) también estima las resultantes para el paso del
seno de la ola:
h2 h h0 H h p1
Rs
2 2
h3 h h0 H h p1
2
Mc
6 6
P 0,70 w H H<d
H
P w H 0,15 0,55 H>d
d
- 58 -
CAPÍTULO II.- ESTADO DEL ARTE
tanh h H
Fmax
1
2
h2 H2 H
k
tanh k h
Fmax
1
2
h2 H2 H
k
tanh k h
Fmax
1
2
h2 1,3 H 0,15 H h H
2
- 59 -
Estudio dinámico de cimentaciones en obras portuarias
Con la idea de proponer una única formulación válida tanto en zona de on-
das estacionarias como en la zona de rotura, Goda (1974) realiza una nueva
propuesta con la intención de emplear una única formulación tanto en la zona
de ondas estacionarias como de olas en rotura.
El método Goda tiene una base teórica que posteriormente fue contrastada
y ajustada con una serie de ensayos de laboratorio y análisis de fallos históri-
cos. La forma geométrica de la ley de empuje puede consultarse en la Figura
nº II.19.
- 60 -
CAPÍTULO II.- ESTADO DEL ARTE
p1 0,5 1 cos 1 2 cos2 Hd
p1
p2
2 h
cosh
L
gT 2 2 h
L tanh
2 L
- 61 -
Estudio dinámico de cimentaciones en obras portuarias
Siendo h la profundidad del fondo del mar respecto al nivel del agua y T
el período del oleaje.
p 3 3 p1
2
1 4 h /L
1 0,6
2 senh4 h / L
h d H 2 2 d
2 min b d ,
3 hb d Hd
donde hb es la profundidad del fondo del mar a una distancia del pie lado
mar de la banqueta de 5 veces la altura de ola significante.
h' 1
3 1 1
h 2 h
cosh
L
* 0,75 1 cos Hd
- 62 -
CAPÍTULO II.- ESTADO DEL ARTE
p u 0,5 1 cos 1 3 H d
p1 Cb2 T
- 63 -
Estudio dinámico de cimentaciones en obras portuarias
Por tanto, la ley de presiones tendría la misma forma que la de Goda pero
modificando el valor de los parámetros. En concreto, las principales modifica-
ciones introducidas serían las siguientes:
- 64 -
CAPÍTULO II.- ESTADO DEL ARTE
Con todo ellos los nuevos parámetros que habría que considerar según Ta-
kahashi serían:
* 0,75 1 cos 1 Hd
p1 0,5 1 cos 1 1 2 * cos2 Hd
p 3 3 p1
p 4 4 p1
p u 0,5 1 cos 1 1 3 H d
i i0 i1
H
i0 para H < 2d y i0 = 2 para H > 2d
d
- 65 -
Estudio dinámico de cimentaciones en obras portuarias
3,134
H
Fh 15 w h2s si
hs
Donde hs es la diferencia de cota entre el pie del cajón lado mar y el nivel
libre del agua y Hsi la altura de ola.
En esta misma línea, en 2010 Cuomo et al. propusieron una fórmula para
la fuerza horizontal impulsiva si la pendiente del lecho marino 1/250:
1
PROVERBS es un proyecto de la Unión Europea creado para desarrollar e implementar las he-
rramientas de probabilidad en el diseño del dique vertical en los Grupos de Trabajo MAST III.
Ver Oumeraci et al. (2000).
- 66 -
CAPÍTULO II.- ESTADO DEL ARTE
hb h s
Fh,imp,1 / 250 w Hmo L hs 1
hs
1 Hmo
hb arctan h
k 0,14 L hs
Esta expresión es válida si Hmo está comprendida entre 0,2 m y 0,7, hs en-
tre 0,5 m y 1,3 m y Tm entre 2 s y 3,7 s.
hb d
Fh,imp(1 / 250 ) C 1r,65 g Hm0 L h s 1
d
2
Fh,qs(1 / 250) 4,8 g Hm 0
Fh
Fh*
w Hb2
- 67 -
Estudio dinámico de cimentaciones en obras portuarias
Fh*
1 ln P Fh*
- 68 -
CAPÍTULO II.- ESTADO DEL ARTE
En las últimas décadas han sido desarrollados modelos numéricos para re-
producir la interacción del oleaje con las estructuras. A continuación se resu-
men algunos de os trabajos más interesantes.
- 69 -
Estudio dinámico de cimentaciones en obras portuarias
Concluida la revisión del estado del arte sobre la acción del oleaje sobre los
diques verticales se pueden resumir las siguientes conclusiones:
- 70 -
CAPÍTULO II.- ESTADO DEL ARTE
- 71 -
Estudio dinámico de cimentaciones en obras portuarias
El primer paso es establecer los criterios para obtener los valores del olea-
je para verificar la seguridad del dique de abrigo mediante un procedimiento
de Nivel I. Existen dos posibilidades tal como se explica a continuación.
- 72 -
CAPÍTULO II.- ESTADO DEL ARTE
Por una parte suponer que las probabilidades asignadas a los modos de fa-
llo geotécnicos son bajas, inferiores al 5%. En este caso el valor de cálculo del
oleaje es igual al valor característico correspondiente a un periodo de retorno
de 50 años afectado por el coeficiente de mayoración de acciones y el coefi-
ciente de compatibilidad. Y para esta hipótesis de cálculo se aceptarían los
coeficientes de seguridad propuestos en la ROM 0.5-05 para los distintos mo-
dos de fallo geotécnicos.
- 73 -
Estudio dinámico de cimentaciones en obras portuarias
h
FHU 0
u w dz
Por una parte existen numerosos cálculos que intentan reproducir median-
te un modelo numérico los empujes producidos por el oleaje sobre el para-
mento vertical del dique. Las referencias se pueden consultar en el apartado
de este estado del arte dedicado al oleaje.
- 74 -
CAPÍTULO II.- ESTADO DEL ARTE
- 75 -
Estudio dinámico de cimentaciones en obras portuarias
También Moradi et al. (2015) han utilizado el programa FLAC con para es-
tudiar el comportamiento de un aluvial compuesto por arena en condiciones
saturadas. Y también la generación de presiones intersticiales es realizada
mediante el modelo Finn (1977).
En este apartado del estado del arte sobre metodologías de cálculo es in-
teresante incluir la referencia de Madrid et al (2013). En este procedimiento
simplificado se estudia en primer lugar el comportamiento dinámico del suelo
mediante ensayos triaxiales cíclicos y ensayos de corte simple. A partir de es-
tos datos se puede conocer el número de ciclos (impactos del oleaje en este
caso) necesarios para alcanzar la rotura del terreno en función de la tensión
tangencial preexistente normalizada y la tensión cíclica normalizada (diagrama
de interacción).
- 76 -
CAPÍTULO II.- ESTADO DEL ARTE
Para cada cálculo se supone una resistencia al corte sin drenaje deducida
del diagrama de interacción obtenido con los ensayos de laboratorio. Y el re-
sultado de cada cálculo es un coeficiente de seguridad.
- 77 -
Estudio dinámico de cimentaciones en obras portuarias
- 78 -
CAPÍTULO II.- ESTADO DEL ARTE
de una longitud suficiente sobre las que actúan cargas contenidas en el plano
de cálculo (sección transversal).
Adicionalmente el programa tiene una opción del cálculo para tensión pla-
na (para materiales elásticos o plásticos con criterio de rotura de Mohr-
Coulomb). En este tipo de modelo, las tensiones normales a la sección trans-
versal (de cálculo) son nulas. Este caso puede darse, por ejemplo, en placas
delgadas cargadas únicamente en su plano.
Por último existe otra opción para realizar los modelos: geometría axilsi-
métrica. En este caso e utilizan coordenadas cilíndricas para definir la malla,
siendo el eje de simetría el correspondiente a x=0, el sentido positivo de x
correspondería a la coordenada radial, la dirección según el eje y corresponde-
ría a la coordenada axial y la dirección del eje z a la coordenada circunferen-
cial. Este modelo de cálculo puede utilizarse para reproducir, por ejemplo,
aquellos ensayos de laboratorio que se realizan sobre probetas cilíndricas.
- 79 -
Estudio dinámico de cimentaciones en obras portuarias
- 80 -
CAPÍTULO II.- ESTADO DEL ARTE
Para realizar cada ciclo alrededor del bucle representado se utiliza un paso
de tiempo determinado. Una cuestión importante es que cada recuadro pre-
sentado en la Figura nº II.25 actualiza todas las variables de la malla a partir
de los valores conocidos que permanecen fijos cuando el control se encuentra
en el recuadro determinado. Por ejemplo, en el recuadro inferior, se toman las
velocidades ya calculadas y, para cada elemento, calcula las nuevas tensiones,
asumiendo que las velocidades se encuentran “congeladas” para la realización
de este cálculo. Se supone que las nuevas tensiones no afectan a las velocida-
des en este recuadro. Esto puede no parecer razonable ya que se sabe que si
cambian las tensiones en alguna zona, esto afectará tanto a esa zona como a
las vecinas que cambiarán sus velocidades. Para minimizar este hecho, el pro-
grama elige un paso de tiempo muy reducido de tal modo que la “información”
de cambio de tensión no pueda propagarse a elementos contiguos en el inter-
- 81 -
Estudio dinámico de cimentaciones en obras portuarias
En general puede decirse que los métodos explícitos son mejores para sis-
temas con un comportamiento “malo” (modelos no lineales, de grandes de-
formaciones, físicamente inestables,…) mientras que no son tan eficientes pa-
ra problemas lineales en pequeñas deformaciones.
- 82 -
CAPÍTULO II.- ESTADO DEL ARTE
- 83 -
CAPÍTULO III
PROPUESTA DE METODOLOGÍA DE ANÁLISIS
CAPÍTULO III.- PROPUESTA DE METODOLOGÍA DE ANÁLISIS
III.1 INTRODUCCIÓN
Se supone que la estructura del dique descansa sobre una banqueta gra-
nular, libremente drenante, que, a su vez, apoya sobre el terreno natural for-
mado por suelos susceptibles de reducir su resistencia al aumentar el número
de ciclos de carga.
- Cálculos dinámicos.
- 87 -
Estudio dinámico de cimentaciones en obras portuarias
H2
Tiempo patrón =
Cv
K Em
Cv =
w
Los suelos que podrían ser conflictivos en los cimientos de los diques de
abrigo verticales, pueden tener módulos de deformación confinada, para el
rango de presiones de interés en el entorno 1 MPa < Em < 10 MPa.
Por otro lado, una indicación de la “longitud de drenaje” puede ser del or-
den de la mitad del ancho de la base del cajón. Un valor típico es del orden de
10 m.
0,1 m 1m
< Tiempo patrón <
K K
- 88 -
CAPÍTULO III.- PROPUESTA DE METODOLOGÍA DE ANÁLISIS
Para estar seguro de que no existe drenaje, en el otro extremo del pro-
blema con T = 105 s se tendría un coeficiente de permeabilidad del orden de
10-6 m/s.
- 89 -
Estudio dinámico de cimentaciones en obras portuarias
2
El valor (N1)60 es el índice N del ensayo de penetración estándar SPT, normalizado para una
presión efectiva de 1 atm y una energía efectiva del 60% de la nominal del equipo de hinca.
- 90 -
CAPÍTULO III.- PROPUESTA DE METODOLOGÍA DE ANÁLISIS
Como puede verse para un valor de (N1)60 (que en dicha figura aparece
como N 60) de 15 golpes se tiene que el valor de CRR crece con el contenido
- 91 -
Estudio dinámico de cimentaciones en obras portuarias
Por otro lado la resistencia dinámica, para otro número de ciclos de carga
distintos del nominal (N1)60 = 15 varía según se indica en Figura nº III.3.
- 92 -
CAPÍTULO III.- PROPUESTA DE METODOLOGÍA DE ANÁLISIS
Como se ve, el valor de CRR (definido para quince ciclos) es, en arenas,
claramente menor que la resistencia estática. Sin embargo, para suelos limo-
sos o limos, esa reducción de resistencia es menos notable; el número de ci-
clos no afecta de manera tan importante a la resistencia dinámica.
III.2 GEOMETRÍA
- 93 -
Estudio dinámico de cimentaciones en obras portuarias
Con mar en calma la resultante del peso del cajón y su excentricidad son:
- 94 -
CAPÍTULO III.- PROPUESTA DE METODOLOGÍA DE ANÁLISIS
Una vez definida la geometría, para realizar los cálculos es necesario dis-
cretizarla. En este caso esta discretización se ha efectuado mediante elemen-
tos rectangulares típicamente de 2,20 metros de anchura y 2,17 metros de
altura excepto en los elementos que representan el cajón en el que se han uti-
lizado elementos de, aproximadamente, el doble en cada una de las dimensio-
nes.
- 95 -
Estudio dinámico de cimentaciones en obras portuarias
(W' S)
E=
F
- 96 -
CAPÍTULO III.- PROPUESTA DE METODOLOGÍA DE ANÁLISIS
Parámetros resistentes.
Parámetros de deformabilidad.
Edómetros o presiómetros.
Parámetros de permeabilidad.
- 97 -
Estudio dinámico de cimentaciones en obras portuarias
Existen, dentro de los modelos del comportamiento dinámico del suelo, dos
que han sido desarrollados recientemente en la Escuela de Ingenieros de Ca-
minos, Canales y Puertos de la UPM, previamente citados en el Capítulo II de
esta tesis.
- 98 -
CAPÍTULO III.- PROPUESTA DE METODOLOGÍA DE ANÁLISIS
Se puede decir que el modelo que aquí se usará, que es el UBCS participa
de ambas aproximaciones. Siguiendo a M. Martín (2010) se utiliza la teoría de
la plasticidad generalizada y siguiendo a R. Galindo, se usa el programa FLAC,
que permite adoptar el UBCS a los resultados experimentales. El modelo que
se propone no reproducirá bien las deformaciones previas a la rotura, detalle
que resuelve excelentemente el modelo de R. Galindo, pero reproduce bien la
condición de rotura pues la resistencia dinámica del suelo, el valor de CRR,
puede ser un parámetro de entrada al modelo.
En efecto, uno de los probables usos del modelo UBCS, es el caso particu-
lar de los suelos arenosos limpios (menos del 5% del contenido en finos) que
podrían ser caracterizados únicamente en el valor del (N)60 del SPT. En ese
caso particular la relación entre el valor del CRR y dicho índice sería el que se
indica en Figura nº III.7.
- 99 -
Estudio dinámico de cimentaciones en obras portuarias
0.5
s'v =100 kPa
0.4
Resistencia cíclica (CRR)
0.3
0.2
0.1
0
0 5 10 15 20 25 30
Resistencia a la penetración dinámica (N)60
4000
3500
Kg p
3000
Número de módulo kg
2500
(adimensional)
2000
Kge
1500
1000
500
0
0 5 10 15 20 25 30
Resistencia a la penetración dinámica (N)60
- 100 -
CAPÍTULO III.- PROPUESTA DE METODOLOGÍA DE ANÁLISIS
Por otro lado el modelo UBCS estándar para arenas incluye un dato empíri-
co de enorme importancia que se refiere a la dependencia del valor de CRR
con la presión efectiva vertical inicial. Normalmente en suelos arenosos se su-
pone que ese factor adimensional es independiente de la presión de ensayo,
- 101 -
Estudio dinámico de cimentaciones en obras portuarias
pero la experiencia demuestra que dicho valor decrece con dicha presión. (Ver
Figura nº III.10).
(N)60=5
(N)60=10
(N)60=15
0.2
(N)60=20
CRR
0.1
0
0 100 200 300 400 500
Tension vertical efectiva inicial (kPa)
Como se ve, este efecto es muy notable en arenas limpias; el valor de CRR
puede disminuir del orden del 20 al 40% (dependiendo de la compacidad de
las arenas) al pasar de presiones efectivas iniciales de 1 atmósfera a presiones
de 4 atmósferas.
Este efecto es, según las experiencias publicadas por A. Soriano et al.
(2011) menos notable en suelos con mayores porcentajes de finos.
- 102 -
CAPÍTULO III.- PROPUESTA DE METODOLOGÍA DE ANÁLISIS
o
cuyo valor relativo es = y el valor correspondiente del ensayo normal,
' vo
CRR()
K () =
CRR (o)
En suelos arenosos limpios el valor de K que resulta del modelo UBCS se re-
coge en Figura nº III.11.
Como puede verse el valor de K puede ser ligeramente mayor que la uni-
dad para arenas densas, (N)60 > 20, pero es claramente menor que la unidad
para suelos menos densos. En este sentido la publicación de H. Patiño et al.
(2013) supone un gran avance al medir para suelos limosos normalmente
consolidados un valor medio del factor K que se indica en Figura nº III.12.
- 103 -
Estudio dinámico de cimentaciones en obras portuarias
K = CRR (a) /CRR (0)
0
0 0.2 0.4
Los parámetros adoptados para estos modelos son los que se presentan en
las siguientes tablas:
- 104 -
CAPÍTULO III.- PROPUESTA DE METODOLOGÍA DE ANÁLISIS
Cimiento
0 30 12500 15600 1,6
(efectivas)
Cimiento
0,32 · 'I 0 12500 15600 1,6
(corto plazo totales)
1,4 0,5 156 125 240 0,5 0,5 0,4 30 30 0,9 14 1,0
- 105 -
Estudio dinámico de cimentaciones en obras portuarias
60
s'vo=100 kPa
50
40
Tensión tangencial (kPa)
30
Mohr‐Coulomb (Su)
20 Mohr‐Coulomb sin drenaje
UBC sin drenaje
Mohr‐Coulomb drenado
10 UBC drenado
0
0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5
Deformacion (%)
- 106 -
CAPÍTULO III.- PROPUESTA DE METODOLOGÍA DE ANÁLISIS
60
s'vo=100 kPa
Mohr‐Coulomb (Su)
50
Mohr‐Coulomb sin drenaje
UBC sin drenaje
Mohr‐Coulomb drenado
40 UBC drenado
Tensión tangencial (kPa)
30
20
10
0
0 ‐20 ‐40 ‐60 ‐80 ‐100 ‐120
Tensión efectiva vertical (kPa)
- 107 -
Estudio dinámico de cimentaciones en obras portuarias
- 108 -
CAPÍTULO III.- PROPUESTA DE METODOLOGÍA DE ANÁLISIS
Fs Fs max
o
tan k n
L ks
- 109 -
Estudio dinámico de cimentaciones en obras portuarias
- 110 -
CAPÍTULO III.- PROPUESTA DE METODOLOGÍA DE ANÁLISIS
- 111 -
Estudio dinámico de cimentaciones en obras portuarias
- 112 -
CAPÍTULO III.- PROPUESTA DE METODOLOGÍA DE ANÁLISIS
- 113 -
Estudio dinámico de cimentaciones en obras portuarias
- 114 -
CAPÍTULO III.- PROPUESTA DE METODOLOGÍA DE ANÁLISIS
- 115 -
Estudio dinámico de cimentaciones en obras portuarias
W'
dw
E
H
O
S
- 116 -
CAPÍTULO III.- PROPUESTA DE METODOLOGÍA DE ANÁLISIS
1800
Mohr‐Coulomb
1600
UBC
1400
1200
1000
Empuje (kN)
800
600
400
200
0
0 2 4 6 8 10 12 14 16 18 20 22 24 26 28 30 32 34
Desplazamiento horizontal del cajón (cm)
- 117 -
Estudio dinámico de cimentaciones en obras portuarias
2000
Mohr‐Coulomb
1800 UBC
1600
1400
1200
Empuje (kN)
1000
800
600
400
200
0
0 2 4 6 8 10 12 14 16 18 20 22 24 26 28 30 32 34 36 38 40 42 44 46 48 50 52 54 56 58
Desplazamiento horizontal del cajón (cm)
2200
2000 Empuje Mohr
Empuje UBC
1800
1600
1400
Empuje (kN)
1200
1000
800
600
400
200
0
0 2 4 6 8 10 12 14 16 18 20 22 24 26 28 30 32 34 36 38 40 42 44 46 48
Desplazamiento horizontal del cajón (cm)
- 118 -
CAPÍTULO III.- PROPUESTA DE METODOLOGÍA DE ANÁLISIS
- 119 -
Estudio dinámico de cimentaciones en obras portuarias
σ σ σ σ
|σ | τ
2 2
- 120 -
CAPÍTULO III.- PROPUESTA DE METODOLOGÍA DE ANÁLISIS
Los resultados obtenidos para el caso de cimiento no drenado son los que
se presentan en la siguiente tabla:
- 121 -
Estudio dinámico de cimentaciones en obras portuarias
1500
1400
1300
1200
1100
1000
900
Empuje (kN)
800
700
Mohr‐Coulomb (efectivas)
600 Mohr‐Coulomb (Su)
UBC
500
400
300
200
100
0
0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15
Desplazamiento horizontal del cajón (cm)
- 122 -
CAPÍTULO III.- PROPUESTA DE METODOLOGÍA DE ANÁLISIS
1600
1400
1200
1000
Empuje (kN)
Mohr‐Coulomb (Efectivas)
800 Mohr‐Coulomb (Su)
UBC
600
400
200
0
0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28
Desplazamiento horizontal del cajón (cm)
2200
2000
1800
1600
1400
Empuje (kN)
1200
1000
800 Mohr‐Coulomb (efectivas)
Mohr‐Coulomb (Su)
UBC
600
400
200
0
0 2 4 6 8 10 12 14 16 18 20 22 24 26 28 30 32 34 36 38 40 42 44 46 48 50 52 54 56 58 60 62 64 66 68 70 72 74 76 78 80 82 84
Desplazamiento horizontal del cajón (cm)
- 123 -
Estudio dinámico de cimentaciones en obras portuarias
- 124 -
CAPÍTULO III.- PROPUESTA DE METODOLOGÍA DE ANÁLISIS
- 125 -
Estudio dinámico de cimentaciones en obras portuarias
- 126 -
CAPÍTULO III.- PROPUESTA DE METODOLOGÍA DE ANÁLISIS
0.9
0.8
k=10‐3 cm/s
k=10‐4 cm/s
0.7 k=10‐5 cm/s
Exceso de presión (Pmax/Pactual)
0.6
0.5
0.4
0.3
0.2
0.1
0
0 50 100 150 200 250 300
Tiempo (segundos)
Permeabilidad
-3
Tiempo 10 10-4 10-5
(s) cm/s cm/s cm/s
10 28,9 3,6 0,4
25 52,8 8,6 0,9
50 71,5 16,1 1,8
100 85,5 28,7 3,5
200 94,7 46,6 6,9
- 127 -
Estudio dinámico de cimentaciones en obras portuarias
Una vez establecidas las solicitaciones máximas para el caso estático, tan-
to en condiciones drenadas como sin drenaje, se procede a efectuar un cálculo
dinámico (cargas variables con el tiempo) con el modelo de comportamiento
del cimiento propuesto. La geometría utilizada para los cálculos dinámicos es
la misma que se empleó para realizar los cálculos estáticos. En el presente
apartado se resumen los resultados obtenidos en estos cálculos.
1
Carga unitaria
‐1
0 2 4 6 8 10 12
Tiempo (s)
- 128 -
CAPÍTULO III.- PROPUESTA DE METODOLOGÍA DE ANÁLISIS
Ese esquema se ha utilizado tanto para el empuje (E) como para la sub-
presión (S). Al igual que se hizo para los cálculos estáticos, se ha considerado
que la relación entre el empuje y la subpresión es E = S · 2,6.
En este caso se conoce que el punto de aplicación del empuje (E) sobre el
paramento varía en función de que el empuje se produzca con la cresta de la
ola o con el seno de la misma. Para reproducir esquemáticamente este efecto
se ha supuesto que los empujes negativos (empujes hacia el lado mar) tienen
un punto de aplicación situado a una cota menor que los empujes positivos.
Para las subpresiones se ha considerado siempre el mismo punto de aplica-
ción. Lo comentado en este párrafo se describe de manera gráfica en el es-
quema presentado en la Figura nº III.37.
- 129 -
Estudio dinámico de cimentaciones en obras portuarias
parían, durante el cálculo, menos del 10% de las presiones intersticiales, por
lo que podría realizarse el cálculo suponiendo que no existe drenaje.
Para cada una de estas alturas de banqueta se han realizado cuatro cálcu-
los cíclicos variando los empujes máximos (E) en cada cálculo y obteniendo,
para cada uno de esos empujes máximos el número de ciclos de carga en el
que se produce la rotura. En la Tabla XI se resumen los empujes máximos
considerados en cada uno de los cálculos efectuados.
- 130 -
CAPÍTULO III.- PROPUESTA DE METODOLOGÍA DE ANÁLISIS
- 131 -
Estudio dinámico de cimentaciones en obras portuarias
- 132 -
CAPÍTULO III.- PROPUESTA DE METODOLOGÍA DE ANÁLISIS
Número de ciclos
0 1 2 3 4 5 6 7 8
0.1
0.0
‐0.1
‐0.2
‐0.3
‐0.4
‐0.5
‐0.6
‐0.7 Movimiento vertical
Movmimiento horizontal
‐0.8
Movimiento (m)
‐0.9
‐1.0
‐1.1
‐1.2
‐1.3
‐1.4
‐1.5
‐1.6
‐1.7
‐1.8
‐1.9
‐2.0
- 133 -
Estudio dinámico de cimentaciones en obras portuarias
- 134 -
CAPÍTULO III.- PROPUESTA DE METODOLOGÍA DE ANÁLISIS
Estableciendo los ciclos de carga en los que se alcanza la rotura para cada
uno de los cálculos realizados se llega a los resultados que se muestran en la
Tabla XII (número de ciclos para alcanzar la rotura en función de la carga apli-
cada) y, de manera gráfica en la Figura nº III.45.
- 135 -
Estudio dinámico de cimentaciones en obras portuarias
- 136 -
CAPÍTULO III.- PROPUESTA DE METODOLOGÍA DE ANÁLISIS
1700
Banqueta 6 m
Banqueta 8 m
1500 Banqueta 10 m
1300
Empuje (kN/ml)
1100
900
700
500
1 10 100
Número de ciclos hasta rotura
- 137 -
Estudio dinámico de cimentaciones en obras portuarias
En ese cálculo se obtiene la rotura con 9 ciclos más de carga que en el ca-
so no drenado (o con poco drenaje). En la Figura nº III.46 se presenta una
comparación entre el resultado obtenido para el caso con mayor drenaje y el
efectuado con el cimiento poco drenante.
1100
Banqueta 6 m. 900 kN. Sin drenaje
Banqueta 6 m. 900 kN. Drenaje parcial
1000
900
800
700
600
500
1 10 100
- 138 -
CAPÍTULO III.- PROPUESTA DE METODOLOGÍA DE ANÁLISIS
- 139 -
Estudio dinámico de cimentaciones en obras portuarias
2500
2000
1500
Empuje (kN)
1000
M‐C Sin drenaje (efectivas)
M‐C Sin drenaje (Su)
500
UBC sin drenaje
E15
0
6 7 8 9 10
Altura de banqueta (m)
De los tres cálculos estáticos sin drenaje (equivalente a 1/2 ciclo de carga)
el que sistemáticamente conduce a mayores valores del empuje en rotura es
el correspondiente al modelo de Mohr-Coulomb. Esto se interpreta como una
consecuencia lógica del menor incremento de presiones intersticiales generado
cuando se compara con otro modelo cualquiera que incluya la degradación cí-
clica del terreno creada por los esfuerzos de corte de manera más explícita.
- 140 -
CAPÍTULO III.- PROPUESTA DE METODOLOGÍA DE ANÁLISIS
0.9
0.8
VALOR CRÍTICO E15/E0=0,65
E15 /E0
0.7
0.6
0.5
0.4
5 6 7 8 9 10
Espesor de banqueta (m)
- 141 -
Estudio dinámico de cimentaciones en obras portuarias
y elegir como situación límite la que conduzca a una situación más crítica. Esto
es, calcular Eo y E15 como se ha hecho en lo que precede y elegir como resis-
tencia dinámica del dique el valor:
y ese valor compararlo con la acción máxima del oleaje, Emax, para obtener el
coeficiente de seguridad correspondiente.
Er
F(temporal) =
E max
Espesor de banqueta
Nombre Descripción Resultado
6m 8m 10 m
- 142 -
CAPÍTULO III.- PROPUESTA DE METODOLOGÍA DE ANÁLISIS
- 143 -
CAPÍTULO IV
CONTRASTE CON UN CASO REAL
CAPÍTULO IV.- CONTRASTE CON UN CASO REAL
El proyecto original del dique era “en talud” pero al inicio de la construc-
ción se vio más conveniente modificar el proyecto para construir un dique ver-
tical reflejante.
- 147 -
Estudio dinámico de cimentaciones en obras portuarias
- 148 -
CAPÍTULO IV.- CONTRASTE CON UN CASO REAL
En la Figura nº IV.3 se indica la sección tipo del nuevo dique, hoy operati-
vo.
El terreno donde apoya el dique de abrigo está formado por los suelos del
delta del Llobregat que, en general, están normalmente consolidados o sólo
ligeramente sobreconsolidados.
- 149 -
Estudio dinámico de cimentaciones en obras portuarias
- 150 -
CAPÍTULO IV.- CONTRASTE CON UN CASO REAL
Figura nº IV.4. Sección tipo del dique de la nueva bocana del puerto
de Barcelona, el día de la rotura
El proyecto de construcción del dique indica que, de acuerdo con los datos
existentes y en base a la experiencia de los proyectistas, en esta zona, las
“arenas superiores” tienen cierto porcentaje de finos (15% a 40%). Su ángulo
de rozamiento, supuesto para los cálculos de estabilidad, es = 34º. La hu-
medad saturada de estos suelos queda próxima al 25% y el peso específico
relativo de sus partículas es G = 2,7. Esto es, tienen un índice de poros pró-
ximo a eo 0,68 y una densidad saturada sat 2 gr/cm3.
- 151 -
Estudio dinámico de cimentaciones en obras portuarias
Con estos nuevos datos se pudo precisar que las arenas mencionadas tie-
nen un contenido en finos cuyo valor medio es 15% y cuyo valor máximo es el
50% y que su ángulo de rozamiento es = 35º.
- 152 -
CAPÍTULO IV.- CONTRASTE CON UN CASO REAL
- 153 -
Estudio dinámico de cimentaciones en obras portuarias
- 154 -
CAPÍTULO IV.- CONTRASTE CON UN CASO REAL
- 155 -
Estudio dinámico de cimentaciones en obras portuarias
Los ensayos realizados indican que los finos son arcillas de baja plastici-
dad. Los índices de plasticidad parecen quedar en el rango 5% < IP < 20%.
Ver Figura nº IV.8.
Como puede verse el valor del índice SPT de las zonas menos arenosas es,
como media del orden de N (SPT) = 15 mientras que para las zonas más are-
nosas ese índice es del orden del doble del indicado.
- 156 -
CAPÍTULO IV.- CONTRASTE CON UN CASO REAL
Coeficientes de permeabilidad
Las muestras más arenosas parecen tener una permeabilidad mayor pero
no mayor que 4 10-7 cm/s para cargas moderadas.
Ángulo de rozamiento
- 157 -
Estudio dinámico de cimentaciones en obras portuarias
Los ensayos de corte simple indican los resultados que se recogen en Figu-
ra nº IV.11 (para una presión de consolidación aproximada de 350 kPa). Para
otras presiones los resultados se representan gráficamente en Figura nº IV.12.
Deformabilidad
- 158 -
CAPÍTULO IV.- CONTRASTE CON UN CASO REAL
250
150
100
50
0
20 25 30 35 40 45
1 0.5
ARENAS LIMOS
0.9 0.45
Resistencia al corte sin drenaje relativo, rot /' v
Resistencia al corte sin drenaje relativo, rot /' v
0.8 0.4
0.7 0.35
0.6 0.3
0.5 0.25
0.4 0.2
0.3 0.15
0.2 0.1
0.1 0.05
0 0
10 15 20 25 30 35 30 35 40 45 50 55 60 65
Humedad natural, w (%) Humedad natural, w (%)
- 159 -
Estudio dinámico de cimentaciones en obras portuarias
Figura nº IV.13.- Gráficos de ajuste de la deformabilidad en la interpretación de los ensayos triaxiales CU según la
teoría del modelo hiperbólico
- 160 -
CAPÍTULO IV.- CONTRASTE CON UN CASO REAL
Resistencia dinámica
0.3
SC02‐SH2
SC02‐TP2
SC03‐SH2
0.25
SC03‐SH4
0.2
∆t/s'v0
0.15
0.1
0.05
Nota: La muestra SC02‐SH2 es la más arenosa
0
1 10 100 1000 10000
Numero de ciclos para alcanzar la rotura
- 161 -
Estudio dinámico de cimentaciones en obras portuarias
En el manual de uso del modelo (M.H. Beaty y P.M. Byrne, 2011) se esta-
blecen unos parámetros generales para estos datos (en arenas) basados en
correlaciones con el resultado obtenido en ensayos SPT ((N1)60) según las ex-
presiones siguientes:
Kb 0,7 K Ge
0,15
R f min (0,99; 1,1 (N1 )60 )
Para los exponentes de los módulos se propone el valor 0,5 excepto para
el exponente del módulo plástico para el que se propone el valor de 0,4.
- 162 -
CAPÍTULO IV.- CONTRASTE CON UN CASO REAL
4000
3500
3000
Número de módulo kg
2500
(adimensional)
2000
Kg e
1500
1000
500
0
0 5 10 15 20 25 30
Resistencia a la penetración dinámica (N)60
1.2
0.8
0.6
Rf
0.4
0.2
0
0 5 10 15 20 25 30 35
Nspt
- 163 -
Estudio dinámico de cimentaciones en obras portuarias
Tal como se observa en las figuras anteriores, los autores, para arenas con
un golpeo corregido en el ensayo SPT comprendido entre 15 y 20 presentarían
unos números de módulo de deformación tangencial (tanto elástico como plás-
tico) del orden de 1100 (adimensional) y una relación de rotura en torno a
0,7 (adimensional).
Vx Vx
x0
y0
U0
y
- 164 -
CAPÍTULO IV.- CONTRASTE CON UN CASO REAL
Para efectuar los cálculos cíclicos con el modelo presentado en la figura an-
terior el procedimiento es establecer una velocidad en los nodos superiores
suficientemente reducida e ir controlando (en cada paso de cálculo) el valor de
la tensión tangencial en el elemento. Al alcanzar la tensión tangencial deseada
se invierte el signo de las velocidades y, de este modo se completan los ciclos.
De este modo, con los parámetros para las arenas y los limos resumidos
en la Tabla XVI se obtienen los resultados para la resistencia que muestran en
la Figura nº IV.18.
Arenas 1,6 0,4 1368 957,6 2098,8 0,5 0,5 0,4 35 35 0,8 20 1,0
Limos 1,6 0,4 624 500,0 960,0 0,5 0,5 0,4 30 30 0,9 14 1,0
- 165 -
Estudio dinámico de cimentaciones en obras portuarias
0.3
Arenas
Limos
0.25
0.2
t/s 'v0
0.15
0.1
0.05
0
1 10 100 1000 10000
Numero de ciclos para alcanzar la rotura
SC3‐SH4
0.9 SC02‐SH2
SC02‐TP2
SC03‐SH2
0.8 Limos (ensayo numérico)
Arenas (ensayo numérico)
0.7
0.6
s'v/s' vo
0.5
0.4
0.3
0.2
0.1
0
0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20
Número de ciclos
- 166 -
CAPÍTULO IV.- CONTRASTE CON UN CASO REAL
G
Gmax Gmax
G1
c G1
- 167 -
Estudio dinámico de cimentaciones en obras portuarias
Gmax v2s
- 168 -
CAPÍTULO IV.- CONTRASTE CON UN CASO REAL
- 169 -
Estudio dinámico de cimentaciones en obras portuarias
350
s'm=100 kPa
s'm=200 kPa
s'm=300 kPa
300
s'm=400 kPa
s'm=500 kPa
250
200
Gmax (MPa)
150
100
50
0
0 5 10 15 20 25 30
(N 1)60
- 170 -
CAPÍTULO IV.- CONTRASTE CON UN CASO REAL
180
Ohta y Goto(1976). Seed et al.(1986)
160 s'm=100 kPa
Imai y Tonouchi(1982)
140
120
Gmax (MPa)
100
80
60
40
20
0
0 5 10 15 20 25 30
N1 60
- 171 -
Estudio dinámico de cimentaciones en obras portuarias
G
K , PI ( 'm )m ( ,PI) mo
Gmax
Donde:
0,000102 n (PI)0,492
K (, PI) = 0,5 1 tanh ln
0,000556 0,4
m (, PI) - mO = 0,272 1 tanh ln exp (0,0145 PI1,3 )
- 172 -
CAPÍTULO IV.- CONTRASTE CON UN CASO REAL
0,0 for PI 0
3,37 10 6 PI1,404 for 0 PI 15
n (PI) =
7,0 10 7 PI1,976 for 15 PI 70
2,7 10 5 PI1,115 for PI 70
0.9
Confinamiento 1 kPa
0.8
Confinamiento 50 kPa
Confinamiento 100 kPa
0.7 Confinamiento 200 kPa
Confinamiento 300 kPa
Confinamiento 400 kPa
0.6
G/Gmax
0.5
0.4
0.3
0.2
0.1
0
1.E‐06 1.E‐05 1.E‐04 1.E‐03 1.E‐02 1.E‐01 1.E+00
Amplitud de deformación tangencial cíclica( g)
- 173 -
Estudio dinámico de cimentaciones en obras portuarias
0.9
Confinamiento 1 kPa
0.8
Confinamiento 100 kPa
Confinamiento 200 kPa
0.7 Confinamiento 400 kPa
0.6
G/Gmax
0.5
0.4
0.3
0.2
0.1
0
1.E‐06 1.E‐05 1.E‐04 1.E‐03 1.E‐02 1.E‐01 1.E+00
Amplitud de deformación tangencial cíclica( g)
IP=0
0.9 IP=10
IP=20
IP=30
0.8 IP=40
IP=50
0.7
0.6
G/Gmax
0.5
0.4
0.3
0.2
0.1
0
1.E‐06 1.E‐05 1.E‐04 1.E‐03 1.E‐02 1.E‐01 1.E+00
Amplitud de deformación tangencial cíclica( g)
- 174 -
CAPÍTULO IV.- CONTRASTE CON UN CASO REAL
0.9
Limos
Arenas
0.8
IP=0 Conf. 100 kPa (Ishibashi y Zang (1993))
IP=20 Conf 100 kPa. (Ishibashi y Zang (1993))
0.7
0.6
G/Gmax
0.5
0.4
0.3
0.2
0.1
0
0.000001 0.00001 0.0001 0.001 0.01 0.1 1
Amplitud de deformación tangencial (g)
- 175 -
Estudio dinámico de cimentaciones en obras portuarias
- 176 -
CAPÍTULO IV.- CONTRASTE CON UN CASO REAL
35
IP=0
IP=10
30 IP=20
IP=30
IP=40
IP=50
25
Amortiguamiento (%)
20
15
10
0
0,000001 0,00001 0,0001 0,001 0,01 0,1 1
Amplitud de deformación tangencial cíclica( g)
35
Confinamiento 1 kPa
30 Confinamiento 100 kPa
Confinamiento 300 kPa
Confinamiento 400 kPa
Confinamiento 200 kPa
25
Confinamiento 50 kPa
Amortiguamiento (%)
20
15
10
0
1,E‐06 1,E‐05 1,E‐04 1,E‐03 1,E‐02 1,E‐01 1,E+00
Amplitud de deformación tangencial cíclica( g)
- 177 -
Estudio dinámico de cimentaciones en obras portuarias
50
Limos
Arenas
45 IP=0 Confinamiento 100 kPa (Ishibashi y Zang (1993))
IP=20 Confinamiento 100 kPa (Ishibashi y Zang (1993))
40
35
30
Amortiguamiento (%)
25
20
15
10
0
0.00001 0.0001 0.001 0.01 0.1
Deformación tangencial (g)
Por último, como ya se indica en el manual del modelo escrito por sus au-
tores (Beaty y Byrne 2011) las deformaciones que se registran en los ensayos
cíclicos son claramente mayores a las que se registran en los ensayos numéri-
cos efectuados con el modelo UBC. Para ilustrar este extremo se incluye la
Figura nº IV.30 en la que se comparan la curva tensión tangencial versus de-
- 178 -
CAPÍTULO IV.- CONTRASTE CON UN CASO REAL
- 179 -
Estudio dinámico de cimentaciones en obras portuarias
- 180 -
CAPÍTULO IV.- CONTRASTE CON UN CASO REAL
Blanes
Tortosa
5
Llobregat
Roses
4
Altura de ola significante Hs (m)
0
0 3 6 9 12 15 18 21 24 27 30 33 36 39 42 45 48 51 54 57 60 63 66 69 72
Tiempo (horas). Inicio 0 h del día 10 de noviembre de 2001
- 181 -
Estudio dinámico de cimentaciones en obras portuarias
El peso sumergido del cajón, estimado por los proyectistas para la situa-
ción del 10-11 noviembre, poco antes de hundirse es:
La acción del temporal, según esos mismos técnicos alcanza, para la altura
de ola significante de Hs = 4,68 m y siguiendo las recomendaciones del méto-
do de Goda los siguientes esfuerzos:
E = 128 t/m
S = 56 t/m
- 182 -
CAPÍTULO IV.- CONTRASTE CON UN CASO REAL
Rv 376
comp = = 32,1 t/m2
B *
2 5,86
B* = 11,7 m
Inclinación de la resultante:
128
tg = = 0,34 ( = 18,8º)
376
a) Deslizamiento
Fd =
tg
Fd = 1,8
b) Vuelco
- 183 -
Estudio dinámico de cimentaciones en obras portuarias
1
MEST = 432 19,6 = 4233 mt/ml
2
2
MDES = 128 101 + 56 19,6 = 2.025 mt/ml
3
y resultaría:
4.233
Fv = = 2,1
2.025
1 * 1 432
B (min) = = 4,32 m
2 2 50
Con este cambio, los momentos respecto al eje de giro retranqueado se-
rían:
1
MEST = 432 19,6 4,32 = 2.367 mt/ml
2
2
MDES = 128 10,1 + 56 19,6 4,32 = 1782 mt/ml
3
2.367
Fv = = 1,33
1.782
c) Hundimiento
- 184 -
CAPÍTULO IV.- CONTRASTE CON UN CASO REAL
432
’v = = 22 t/m2 (banqueta)
19,6
432
’v = + 3 1,1 = 22,4 t/m2
19,6 3
su = 0,40 · ’v
y en consecuencia:
20,6
Fh (sin drenaje) = = 0,9
22,4
- 185 -
Estudio dinámico de cimentaciones en obras portuarias
1
v (hundimiento, con drenaje) = ’ B* N (N según ROM 0.5-94)
2
v = 1 5,86 22 = 129
y en consecuencia:
129
Fh (con drenaje) = =4 (sin temporal)
32,1
y en consecuencia:
- 186 -
CAPÍTULO IV.- CONTRASTE CON UN CASO REAL
IV.8.1 Geometría
- 187 -
Estudio dinámico de cimentaciones en obras portuarias
- 188 -
CAPÍTULO IV.- CONTRASTE CON UN CASO REAL
- 189 -
Estudio dinámico de cimentaciones en obras portuarias
Escollera
0 40 20830 9615 1,6 0,4
(banqueta)
Arenas
0 35 83333,3 38461,5 1,6 0,4
inferiores
- 190 -
CAPÍTULO IV.- CONTRASTE CON UN CASO REAL
Para los otros dos materiales (arenas superiores y limos) se han realizado
tres (3) hipótesis (casos) distintas sobre el modelo de comportamiento y las
características de los materiales para analizar las dos situaciones que se calcu-
lan con el modelo estático: fondeo del cajón y cálculo pseudoestático. A conti-
nuación se presentan las distintas hipótesis consideradas.
- 191 -
Estudio dinámico de cimentaciones en obras portuarias
En este caso se considera que tanto los limos como las arenas tienen un
modelo de comportamiento elasto-plástico con criterio de rotura de Mohr-
Coulomb y una resistencia al corte basada en los parámetros obtenidos de los
ensayos de laboratorio disponibles. Las características de este caso son las que
se resumen en la Tabla XX.
En este tercer (y último caso) se ha considerado que tanto las arenas como
los limos responden al modelo de comportamiento UBC (descrito en capítulos
previos). Las características de los materiales para este caso son las que se re-
sumen en la Tabla XXI y han sido deducidas de los ensayos de laboratorio dis-
ponibles.
- 192 -
CAPÍTULO IV.- CONTRASTE CON UN CASO REAL
Arenas 1,6 0,4 342 239,4 524,7 0,75 0,75 0,4 35 35 0,8 20 1,0
Limos 1,6 0,4 156 125,0 240,0 1,0 1,0 0,4 30 30 0,9 14 1,0
- 193 -
Estudio dinámico de cimentaciones en obras portuarias
1) Fondeo del cajón sin drenaje para cada una de las hipótesis considera-
das. Para efectuar este cálculo, partiendo de la situación inicial, se in-
troduce en el modelo el cajón y se impide el drenaje en las arenas su-
periores y en los limos, considerando tanto la escollera como las are-
nas inferiores libremente drenantes. Para cada uno de los casos se cal-
cula el coeficiente de seguridad obtenido (una vez fondeado el cajón)
mediante la técnica de reducción de los parámetros de resistencia al
corte de los materiales.
- 194 -
CAPÍTULO IV.- CONTRASTE CON UN CASO REAL
Situación inicial
- 195 -
Estudio dinámico de cimentaciones en obras portuarias
- 196 -
CAPÍTULO IV.- CONTRASTE CON UN CASO REAL
- 197 -
Estudio dinámico de cimentaciones en obras portuarias
- 198 -
CAPÍTULO IV.- CONTRASTE CON UN CASO REAL
- 199 -
Estudio dinámico de cimentaciones en obras portuarias
- 200 -
CAPÍTULO IV.- CONTRASTE CON UN CASO REAL
- 201 -
Estudio dinámico de cimentaciones en obras portuarias
- 202 -
CAPÍTULO IV.- CONTRASTE CON UN CASO REAL
- 203 -
Estudio dinámico de cimentaciones en obras portuarias
Aunque en principio pudiera pensarse que las arenas drenan algo durante
el temporal de cálculo, las estimaciones del coeficiente de permeabilidad con-
duce, en todo caso a K < 10-3 cm/s y en consecuencia se cree que el drenaje
tendrá un efecto imperceptible. Únicamente la banqueta de escollera se su-
pondrá drenante en estos cálculos.
Se supone que, como quiera que desde el fondeo y relleno de celdas trans-
currieron unos dos meses y medio, que en el momento de ocurrencia del tem-
poral no existían en el cimiento presiones intersticiales relevantes. En este
cálculo seudoestático sólo se consideran las presiones intersticiales generadas
por la acción del oleaje sobre los cajones.
Para conocer las cargas que producen la rotura, para cada modelo (casos
a, b y c comentados anteriormente) se introduce en el modelo una carga hori-
zontal a 10 metros de la base del cajón (cota -7,5) que simula el empuje y
una carga vertical positiva en la base del cajón que simula la subpresión (ver
Figura nº IV.33). Para cada carga se calcula el desplazamiento del cajón obte-
niéndose las curvas carga-desplazamiento que se muestran en la Figura nº
IV.49.
- 204 -
CAPÍTULO IV.- CONTRASTE CON UN CASO REAL
700
600
500
400
Empuje (kN)
Caso a) M‐C. Su
300
Caso b) M‐C Sin drenaje
Caso c) UBC sin drenaje
200
100
0
0 2 4 6 8 10 12 14 16 18 20
Desplazamiento (cm)
Del gráfico anterior se obtiene que el empuje máximo para el caso a) re-
sultaría ser de 517 kN, 591 para el caso b) y, por último de 595 kN para el
caso c). Teniendo en cuenta que el empuje que se considera para el temporal
es de 1280 kN, se obtienen los siguientes coeficientes de seguridad:
Las formas de rotura que se obtienen son las indicadas en Figuras nos
IV.50, IV.51 y IV.52. Tal como se observa en esas figuras, en el caso a) la lí-
nea de rotura profundiza más en el cimiento mientras que la forma de rotura
en los otros dos casos (b y c) son muy similares.
- 205 -
Estudio dinámico de cimentaciones en obras portuarias
- 206 -
CAPÍTULO IV.- CONTRASTE CON UN CASO REAL
- 207 -
Estudio dinámico de cimentaciones en obras portuarias
60
40
Sobrepresión intersticial (kPa)
20
‐20
‐40
0 2 4 6 8 10 12 14 16 18
Desplazamiento (cm)
- 208 -
CAPÍTULO IV.- CONTRASTE CON UN CASO REAL
- 209 -
Estudio dinámico de cimentaciones en obras portuarias
E1
10 m E2
6,5 m
Las cargas introducidas varían en el tiempo según una onda sinusoidal con
un periodo de 12 segundos y una amplitud variable (fija en cada caso de car-
ga). En la Figura nº IV.56 se muestra el esquema de la carga cíclica unitaria
introducida en el modelo.
- 210 -
CAPÍTULO IV.- CONTRASTE CON UN CASO REAL
E1
Carga unitaria
E2
‐1
0 1 2 3 4 5 6 7 8 9 10 11 12
Tiempo (segundos)
Ese esquema se ha utilizado tanto para el empuje (E) como para la sub-
presión (S). Al igual que se hizo para los cálculos estáticos, se ha considerado
que la relación entre el empuje y la subpresión es E = S · 2,6.
- 211 -
Estudio dinámico de cimentaciones en obras portuarias
Número de ciclos
0 1 2 3 4 5 6 7 8 9 10 11 12
0.0
‐0.1
‐0.2
‐0.3
‐0.4
‐0.5
‐0.6
‐0.7
Asiento del cajón (m)
‐0.8
‐0.9
‐1.0
‐1.1
‐1.2
‐1.3
‐1.4
‐1.5
‐1.6
‐1.7
- 212 -
CAPÍTULO IV.- CONTRASTE CON UN CASO REAL
- 213 -
Estudio dinámico de cimentaciones en obras portuarias
- 214 -
CAPÍTULO IV.- CONTRASTE CON UN CASO REAL
Número de ciclos
0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18
0.0
‐0.1
‐0.2
‐0.3
‐0.4
Asiento del cajón (m)
‐0.5
‐0.6
‐0.7
‐0.8
‐0.9
‐1.0
‐1.1
‐1.2
- 215 -
Estudio dinámico de cimentaciones en obras portuarias
- 216 -
CAPÍTULO IV.- CONTRASTE CON UN CASO REAL
Con la menor de las cargas cíclicas con las que se ha efectuado el cálculo
se obtienen los asientos en función del número de ciclos de carga que se pre-
senta en la Figura nº IV.65.
En este caso los incrementos de asiento del cajón entre ciclos se van in-
crementando paulatinamente, siendo este incremento de orden milimétrico
entre los ciclos 1 y 19, de orden centimétrico entre los ciclos 20 y 29 y de or-
den decimétrico a partir del ciclo nº 30. Considerando como criterio de rotura
un asiento del cajón de un (1) metro, se obtiene que la rotura se produce a
partir del ciclo 31.
- 217 -
Estudio dinámico de cimentaciones en obras portuarias
Número de ciclos
0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33
0.0
‐0.1
‐0.2
‐0.3
‐0.4
Asiento del cajón (m)
‐0.5
‐0.6
‐0.7
‐0.8
‐0.9
‐1.0
‐1.1
‐1.2
- 218 -
CAPÍTULO IV.- CONTRASTE CON UN CASO REAL
Al igual que sucedía para los casos anteriores de carga, en la zona del ci-
miento situada bajo el cajón se producen sobrepresiones intersticiales debido
a la actuación de la carga cíclica. En la Figura nº IV.67 se presenta la distribu-
ción de estas presiones intersticiales cuando el modelo se encuentra ya en ro-
tura.
Al igual que para los casos anteriores, se ha preparado una figura (Figura
nº IV.68) en la que se muestra la evolución de las sobrepresiones intersticiales
en la zona central del modelo.
- 219 -
Estudio dinámico de cimentaciones en obras portuarias
- 220 -
CAPÍTULO IV.- CONTRASTE CON UN CASO REAL
IV.9.5 Resumen
Número de ciclos
0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32
0.0
‐0.1
‐0.2
‐0.3
‐0.4
Asiento del cajón (m)
‐0.5
Carga 500 kN/ml
Carga 400 kN/ml
‐0.6 Carga 300 kN/ml
‐0.7
‐0.8
‐0.9
‐1.0
Dado que se ha adoptado como criterio de rotura que la base del cajón
presente un asiento de un (1) metro, se tendría que el cajón soportaría 11
ciclos para la carga de 500 kN/ml, 17 ciclos para la carga de 400 kN/ml y 31
ciclos para la carga e 300 kN/ml. Estos datos se presentan de manera gráfica
en la Figura nº IV.70. Interpolando en ese gráfico se obtiene que la carga que
produciría la rotura en 15 ciclos sería de 428 kN/ml.
- 221 -
Estudio dinámico de cimentaciones en obras portuarias
600
500
400
Empuje (kN/ml)
300
200
100
0
10 100
Número de ciclos
F = Er/Emax = 0,46
- 222 -
CAPÍTULO V
CONCLUSIONES Y FUTURAS LÍNEAS DE
INVESTIGACIÓN
CAPÍTULO V.- CONCLUSIONES Y FUTURAS LÍNEAS DE INVESTIGACIÓN
- 225 -
Estudio dinámico de cimentaciones en obras portuarias
Tal como se ha comentado, para realizar la tesis se han tenido que investi-
gar distintos aspectos que influyen en los resultados obtenidos. En general en
cada aspecto estudiado se ha concluido que es necesario realizar investigacio-
nes adicionales:
- 226 -
CAPÍTULO V.- CONCLUSIONES Y FUTURAS LÍNEAS DE INVESTIGACIÓN
- 227 -
Estudio dinámico de cimentaciones en obras portuarias
- 228 -
CAPÍTULO VI
BIBLIOGRAFÍA
CAPÍTULO VI- BIBLIOGRAFÍA
Allsop, W, Vicinanza, D and McKenna, J.E. (1996) " Wave forces on vertical
and composite breakwaters" Strategic Research Report SR 443, HR Wall-
ingford, Wallingford, pp 1-94.
- 231 -
Estudio dinámico de cimentaciones en obras portuarias
Baldi, G., Bellotti, R., Ghionna, V., Jamiolkowsky, M., Marchetti, S., and
Pasquelini, E. (1986). “Flat dilatometer tests in calibration chambers”.
Proceedings, In Situ ’86, Geotechnical Special Publication 6, ASCE, New
York, pp. 431-446.
Beaty, M.H. and Byrne P.M. (2011). “UBCSAND constitutive model: Version
904aR”. Documentation Report: UBCSAND Constitutive Model on Itasca
UDM Web Site.
Been, K. and Jefferies, M. G. (1985). "A state parameter for sands". Géotech-
nique, 35(2), 99-112.
Blackmore, P.A and Hewson, P.J. (1984). "Experiments on full-scale wave im-
pact pressures". Coastal Engineering Volume:08 Issues:4, pp 331-346.
British Standar BS6349 (1990) Part 7: Guidelines for design and construction
of breakwaters. British Standar Institution. London
- 232 -
CAPÍTULO VI- BIBLIOGRAFÍA
Chiu, Y, Lin, J., Chang, S., Lin, Y. and Chen, C. (2007). "An experimental
study of wave forces on vertical breakwater". Journal of Marine Science
and Technology, Vol 15, No 3, pp 158-170.
CIRIA, CUR, CETMEF (2007). The rock manual. The use of rock in hydraulic
engineering (2nd edition). C683, CIRIA, London.
Cuéllar V., Bazant, Z.P., Krizek, R.J. y Silver, M.L. (1977). “Densificación e
histéresis de la arena en estados cíclicos de corte simple”. Journal of Ge-
otechnical and Geoenvironmental Engineering. May. 1997.
Dafalias, Y.F, Herrmann, L.R. and DeNatale, J.S. (1982). “The bounding sur-
face plasticity model for isotropic cohesive soils and its application at the
Grenoble Workshop”. Proc. Int. Workshop on Constitutive Relations for
Soils. Grenoble, 273-287.
De Groot, M.B., Andersen, K.H., Burcharth, H.F., Ibsen, L.B., Kortenhaus, A.,
Lundgren, H., Magda, W., Oumeraci, H. and Richwien, W. (1996). “Foun-
dation design of caisson breakwaters”. Publ. 198. Norwegian Geotech-
nical Institute, Oslo, Norway.
- 233 -
Estudio dinámico de cimentaciones en obras portuarias
De Groot, M.B., Bolton, M.D. Foray P., et al. (2006). “Physics of liquefaction
phenomena around marine structures”. ASCE Journal of Waterway, Port,
Coastal and Ocean Engineering, 132 (4), pp. 227-243.
Dinesh, S.V., Sitharam, T.G., and Vinod, J.S. (2004). “Dynamic properties and
liquefaction behavior of granular materials using discrete element me-
thod”. Current Science, Special Section: Geotechnics and earthquake
hazards, Vol. 87, No. 10.
Finn, W.D.L., Lee, K.W. and Martin, G.R. (1977). “An effective stress model for
liquefaction”. Journal of the Geotechnical Engineering Division. ASCE, Vol.
103, No. GT6, pp. 517-533.
FLAC User’s Manual, version 7.0 (2011). Itasca Consulting Group, Inc, Minea-
polis, Minessota.
Galmés Giralt, J.L. (2005). "Diseño innovador de diques verticales con geome-
trías colíndricas. Anejo 2.- Cálculo de estabilidad estructural de diques
verticales". Revisión histórica y estado del arte.
- 234 -
CAPÍTULO VI- BIBLIOGRAFÍA
Goda, Y. (1974). "A new method of wave pressure calculation for the design of
composite breakwater". International Conference on Coastal Engineering.
Copenhagen, , pp 1702-1720.
Green, T (1989). "Wave forces on vertical walls". Department of Civil and En-
vironmental Engineering. University of Wisconsin.
Hardin, B.O. y Black, A.M. (1966). “Sand stiffness under various triaxial
stresses”. Journal of the Soil Mechanics and Foundations Division, ASCE,
Vol. 92 (SM2), pp. 27-42.
Hardin, B.O. y Dmevich, V.P. (1972a). “Shear modulus and damping in soils:
measurement and parameter effects”. Journal of the Soil Mechanics and
Foundation, ASCE, 98 (SM6), pp. 603-624.
- 235 -
Estudio dinámico de cimentaciones en obras portuarias
Hardin, B.O. y Dmevich, V.P. (1972b). “Shear modulus and damping in soils:
design equation and curves”. ASCE, Journal of the Soil Mechanics and
Foundation, 98 (SM7), pp. 667-692.
Hiroi, I (1919). "On a method of estimating the force of walls". Journal of Col-
lege of Engineering, University of Tokio.
Ishibashi, I. and Zhang, X. (1993). “Unified dynamic shear moduli and damp-
ing ratios of wide strain range”. Soils and Foundations, Vol. 33, No. 1, pp.
182-191.
Ito, Y. and Tanimoto, K. (1971). "A new method of numerical analysis of wave
propagation and its application to wave height distribution along struc-
tures". Proceedings 18th Japanese Conferencia on Coastal Engineering, pp
87-90 (in Japanese).
Iwan, W.D. (1967). “On a class of models for the yielding behavior of continu-
ous and composite systems”. J. Applied Mechanics, ASME 34, 612-617.
Iwasaki, T., Tatsouka, F., and Takagi Y. (1978). “Shear modulus of sands un-
der torsional shear loading”. Soils and Foundations, Vol. 18, no. 1, pp.
39-56.
- 236 -
CAPÍTULO VI- BIBLIOGRAFÍA
Jafarian, Y., Alielahi, H., Abdollahi, A. and Vakili, R. (2010). “Seismic numeri-
cal simulation of breakwater on a liquefiable layer: Iran LNG port”. Elec-
tron Journal Geotech. Eng. 15, pp. 1-11.
Jefferies, M.G. (1993). “Nor-Sand: a simple critical state model for sand”.
Geotechnique, 43(1), 91-103.
Jeng, D.S., Cha, D.H., Lin, Y.S. y Hu, P.S. (2001). “Wave-induced pore pres-
sure around a composite breakwater”. Ocean Eng. 28(10). 1413-1432.
Kaliakin, V.N. and Dafalias, Y.F. (1989). “Simplifications to the bounding sur-
face model for cohesive soils”. Int. J. Numer. Analyt. Meth. Geomech. 13,
91-100.
Kim, M.K. and Lade, P.V. (1988). “Single Hardening Constitutive Model for
Frictional Materials, I., Plastic Potential Function”. Computers and Ge-
otechnics, 5(4), 307-324.
Klisinski, M., Alawi, M.M., Sture, S., Ko, H.-Y. and Muir Wood, D. (1987).
“Elastoplastic model for sand based on fuzzy sets”. Proc. Int. Workshop
on Constitutive Equations for Granular Non-Cohesive Soil. A. Sada and G.
Bianchini (eds.). Case Western Reserve University, Cleveland, 325-347.
Kokusho, T. (1980). “Cyclic triaxial test of dynamic soil properties for wide
strain range”. Soils and Foundations, Vol. 18, No. 1, pp. 39-56.
Kokusho, T., Yoshida, Y., and Esashi, Y. (1982). “Dynamic properties of soft
clay for wide strain range”. Soils and Foundations, Vol. 22, No. 4, pp. 1-
18.
- 237 -
Estudio dinámico de cimentaciones en obras portuarias
Lade, P.V. (2005). “Overview of constitutive models for soils”. ASCE Geotech-
nical Special Publication No. 128. Soil Constitutive Models: Evaluation,
Selection and Calibration, Edited by J.A. Yamamuro and V.N. Kaliakin, pp.
1-34.
Lade, P.V. and Kim, M.K. (1988a). “Single Hardening Constitutive Model for
Frictional Materials, II. Yield Criterion and Plastic Work Contours”. Com-
puters and Geotechnics, 6(1), 13-29.
Lade, P.V. and Kim, M.K. (1988b). “Single Hardening Constitutive Model for
Frictional Materials, III. Comparisons with Experimental Data”. Comput-
ers and Geotechnics, 6(1), 31-47.
- 238 -
CAPÍTULO VI- BIBLIOGRAFÍA
Martín, G.R., Finn, W.D.L., and Deed, H.B. (1975). “Fundamentals of Liguefac-
tion under Cyclic Loading”. J. Geotech., Div. ASCE, 101 (GT5), pp. 423-
438.
Martín, M. (2010). “Sobre la respuesta dinámica del terreno bajo la acción del
oleaje en cajones fondeados en suelos arcillosos”. Tesis doctoral.
Martinelli, L., Voortman, H.G., Lamberti, A and Vrijling, J.k (2001) "Hazard
Analysis of dynamically loaded caisson breakwaters". Coastal Engineering
2000: pp 1636-1649.
Mayne, P.W. and Rix, G.J. (1993). “Gmax·qc relationship for clays”. Geotechni-
cal Testing Journal, ASTM. Vol. 16, No. 1, pp. 54-60.
Mogami, T., and Zubu, K. (1953). “The behaviour of soil during vibration”.
Proceedings, 3rd International Conference on Soil Mechanics and Founda-
tion Engineering, Zurich, Vol. 1, pp. 152-155.
Moradi, G., Sotoubadi, M.H. and Khatibi, B.R. (2015). “The influence of over-
burden pressure on liquefaction potential”. Turkish Journal of Engineering
and Environmental Sciences.
- 239 -
Estudio dinámico de cimentaciones en obras portuarias
Nabili, S. Jafarian, Y. and Baziar, M.H. (2008). “Seismic pore water pressure
generation models: Numerical evaluation and comparison”. 14th World
Conference on Earthquake Engineering. October 12-17, Beijing, China.
Negro V., Varela, O., García, J. y Santos, J. (2001). "Diseño de diques vertica-
les". Colegio de Ingenieros de Caminos, Canales y Puertos. Colección
Seinor.
Oumeraci, H., Allsop, N.W.H., de Groot, M.B., Crouch, R.S. and Vrijling, J.k.
(2000). "Probabilistic design methods for vertical breakwaters (PROV-
ERBS)". Coastal Structures, 99. pp 585,594.
Papadimitriou, A.G., Bouckovalas, G.D. (2002). “Plasticity model for sand un-
der small and large cyclic strains: a multiaxial formulation”. Soil Dynam-
ics and Earthquake Engineering, 22. pp. 191-204.
- 240 -
CAPÍTULO VI- BIBLIOGRAFÍA
Pestana, J.M., Whittle, A.J. and Salvati, L.A. (2002a). “Evaluation of a consti-
tutive model for clays and sands: Part I - sand behaviour”. Int. J. Numer.
Analyt. Meth. Geomech., 26(11), 1097-1121.
Puertos del Estado (2009). “ROM 1.0-09: Recomendaciones del diseño y eje-
cución de las Obras de Abrigo (Parte 1º. Bases y factores para el proyec-
to. Agentes climáticos)”.
Pyke, R. (1979). “Nonlinear soil models for irregular cyclic loadings”. Journal
of Geotechnical and Geoenvironmental Engineering, June 1979.
Rix, G.J. and Stokoe, K.H. (1991). “Correlation of initial tangent modulus and
cone penetration resistance”. Calibration Chamber Testing, International
Symposium on Calibration Chamber Testing, A.B. Huand, ed. Elsevier
Publishing. New York, pp. 351-362.
Roscoe, K.H., and Burland, J.B. (1968). “On the generalized stress-strain be-
haviour of ‘wet’ clay”. In Engineering plasticity (eds. J. Heyman and F.A.
Leckie). Cambridge University Press. pp. 535-609.
- 241 -
Estudio dinámico de cimentaciones en obras portuarias
Sainflou (1928). "Essai sur les digues maritimes verticales". Ann. Ponts et
Chausses.
Schofield, A. and Wroth, C.P. (1968). “Critical State Soil Mechanics”. McGraw-
Hill.
Seed, H.B. and Idriss, I.M. (1970). “Soil moduli and damping factors for dy-
namic response analyses”. Report EERC 70-10, Earthquake Engineering
Research Center, University of California, Berkeley.
Seed, H.B., Wong, R.T., Idriss, I.M., and Tokimatsu, K. (1986). “Moduli and
damping factors for dynamic analyses of cohesionless soils”. Journal of
Geotechnical Engineering, ASCE, Vol. 112, nº. GT11, pp. 1016-1032.
- 242 -
CAPÍTULO VI- BIBLIOGRAFÍA
Sun, J.I., Golesorkhi, R., and Seed, H.B. (1988). “Dynamic moduli and damp-
ing ratios for cohesive soils”. Report no. EERC-88/15. Earthquake Engi-
neering Research Center, University of California, Berkeley.
Wang, Z.L., Dafalias, Y.F., and Shen, C.K. (1990). “Bounding surface hypo-
plasticity model for sand”. Engrg. Mech., ASCE, 116 (5), pp. 983-1001.
Wang, Z.L., Makdisi F.I., and Egan, J. (2006). “Practical Applications of a Non-
linear Approach to Analysis of Earthquake-Induced Liquefaction and De-
formation of Earth Structures”. Journal of Soil Dynamics & Earthquake
Engineering. v. 26, no. 2-4, pp. 231-252.
Zen, K., Umehara, Y. and Hamada, K. (1978).” Laboratory tests and in-situ
seismic survey on vibratory shear modulus of clayey soils with different
plasticities”. Proceedings, Fifth Japan Earthquake Engineering Symposi-
um, Tokyo, pp. 721-728.
- 243 -
APÉNDICE I
ENSAYOS DE LABORATORIO
APÉNDICE I.- ENSAYOS DE LABORATORIO
ÍNDICE
Págs.
1.- IDENTIFICACIÓN
L. Atterberg Granulometría
G
Nombre Prof. % Finos
LL LP IP % Finos (t/m³)
< 2 mm
S1-MI3 42 33,7 21,1 12,6 70,9 34 2,79
-1-
Estudio dinámico de cimentaciones en obras portuarias
-2-
APÉNDICE I.- ENSAYOS DE LABORATORIO
Humedad Humedad
Nombre
Inicial (%) Final (%)
S1-MI3 23,9 18,35
S2-MI1 36,8 23,1
S2-TP2 41,2 33,8
S3-SH2 31,4 23,4
S3-SH5 21,1 19,9
1.2
0.8
Índice de poros (e)
0.6
0.4
S1-MI3
S2-MI1
S2-TP2
S3-SH2
0.2 S3-SH5
0
1 10 100 1000
Tensión vertical (kPa)
-3-
Estudio dinámico de cimentaciones en obras portuarias
30
S1-MI3
S2-MI1
S3_SH2
S3-SH4
28
S3-SH5
26
Humedad final (%)
24
22
20
0 50 100 150 200 250 300 350
Tensión vertical (kPa)
-4-
APÉNDICE I.- ENSAYOS DE LABORATORIO
-5-
Estudio dinámico de cimentaciones en obras portuarias
S1-MI3
60
50
Tensión tangencial (kPa)
40
30
20
10
0
200 300 400 500
Tensión normal (kPa)
S1-SP3
160
140
120
Tensión tangencial (kPa)
100
80
60
40
20
0
200 300 400 500
Tensión normal (kPa)
-6-
APÉNDICE I.- ENSAYOS DE LABORATORIO
S2-MI1
60
50
Tensión tangencial (kPa)
40
30
20
10
0
200 300 400 500
Tensión normal (kPa)
S2-SH2
300
250
Tensión tangencial (kPa)
200
150
100
50
0
200 300 400 500
Tensión normal (kPa)
-7-
Estudio dinámico de cimentaciones en obras portuarias
S2-TP2
160
140
120
Tensión tangencial (kPa)
100
80
60
40
20
0
200 300 400 500
Tensión normal (kPa)
S3-SH2
100
90
80
Tensión tangencial (kPa)
70
60
50
40
30
20
10
0
200 300 400 500
Tensión normal (kPa)
-8-
APÉNDICE I.- ENSAYOS DE LABORATORIO
S3-SH4
140
120
100
Tensión tangencial (kPa)
80
60
40
20
0
200 300 400 500
Tensión normal (kPa)
S3-SH5
250
200
Tensión tangencial (kPa)
150
100
50
0
200 300 400 500
Tensión normal (kPa)
-9-
Estudio dinámico de cimentaciones en obras portuarias
S3-SPT4
160
140
120
Tensión tangencial (kPa)
100
80
60
40
20
0
200 300 400 500
Tensión normal (kPa)
S3-SPT6
250
200
Tensión tangencial (kPa)
150
100
50
0
200 300 400 500
Tensión normal (kPa)
- 10 -
APÉNDICE I.- ENSAYOS DE LABORATORIO
0.80
0.70
0.60
S1-MI3
S1-SP3
Arenas t/s'vo=0,50
0.50 S2-MI1
S2-SH2
t/s'v0
S2-TP2
0.40
S3-SH2
S3-SH4
0.30 S3-SH5
Limos t/s'vo=0,25
S3-SPT4
S3-SPT6
0.20
Series11
Series12
0.10
0.00
200 300 400 500
s'v (kPa)
Asiento N
/’ v0 ru max (%)
(mm) (rotura)
0,10 0,95 8 0,6 350
S2-SH2 0,15 0,91 10 0,6 68
0,20 0,88 10 0,5 19
0,25 0,65 15 2,4 7
S2-TP2 0,15 0,88 15 1,9 450
0,20 0,82 15 1,6 45
0,10 0,81 3 4,8 1300
S3-SH2 0,15 0,81 15 3,4 37
0,25 0,70 15 3,8 8
0,10 0,66 13 2,4 1300
S3-SH4 0,15 0,86 15 2,8 76
0,20 0,71 15 2,6 9
- 11 -
Estudio dinámico de cimentaciones en obras portuarias
S2-SH2
0.3
0.2
t/s'v
0.1
0
1 10 100 1000
Número de ciclos para rotura
S2-TP2
0.3
CRR=0,225
0.2
t/s'v
0.1
0
1 10 100 1000
Número de ciclos para rotura
- 12 -
APÉNDICE I.- ENSAYOS DE LABORATORIO
S3-SH2
0.3
CRR=0,175
0.2
t/s 'v
0.1
0
1 10 100 1000 10000
Número de ciclos para rotura
S3-SH4
0.3
CRR=0,185
0.2
t/s'v
0.1
0
1 10 100 1000 10000
Número de ciclos para rotura
- 13 -
Estudio dinámico de cimentaciones en obras portuarias
0.3
S2-SH2
S2-TP2
S3-SH2
S3-SH4
0.2
t/s'v
0.1
0
1 10 100 1000 10000
Número de ciclos para rotura
- 14 -
0,0 350 350
CONSOLIDACION - CORTE SIMPLE ESTATICO
0,1 300 CORTE SIMPLE ESTATICO 300
τ
250 250
0,2
200 u 200
0,3
150 150
0,4 100 100
CORTE SIMPLE ESTATICO
Deformación, mm .
0,5 50 50
Tensión Tangencial y
Presión de Poros, kPa.
Tensión Tangencial, (kPa) .
0,6 0 0
1 10 100 1000 10000 0 5 10 15 20 25 30 0 100 200 300 400 500
0,0 50 0,250
CONSOLIDACION CORTE SIMPLE CICLICO 40 0,200
0,1 Máx.
30
0,2 0,150
20 Cíclica
10 0,100
0,3
0 0,050
0,4 -10
Unitaria
0,000 Permanente
0,5 -20
-0,050 Mín
-30
Defornación, mm .
0,6 -0,100
-40
Deformación Tangencial .
0,7 -50 -0,150
0,10 30
2,0
20
Amortiguamiento .
0,05 1,0 10
Amortiguamiento .
0,05 1,0 Rampa Kpa 0,00
Deformación, mm .
0,5 50 50
Tensión Tangencial y
Presión de Poros, kPa.
Tensión Tangencial, (kPa) .
0,6 0 0
1 10 100 1000 10000 0 5 10 15 20 25 30 0 100 200 300 400 500
0,0 80 0,300
CONSOLIDACION CORTE SIMPLE CICLICO 60
0,1 0,200 Cíclica
0,2 40 0,100
20 Máx.
0,3 0,000
0 Permanente
0,4 -0,100
Unitaria
-20
0,5 -40 -0,200
Defornación, mm .
0,6 -60 -0,300 Mín
Deformación Tangencial .
0,7 -80 -0,400
0,10 30
2,0
20
Amortiguamiento .
0,05 1,0 10
Amortiguamiento .
0,05 1,0 Rampa Kpa 0,00
Deformación, mm .
0,5 50 50
Tensión Tangencial y
Presión de Poros, kPa.
Tensión Tangencial, (kPa) .
0,6 0 0
1 10 100 1000 10000 0 5 10 15 20 25 30 0 100 200 300 400 500
Unitaria
-40 -0,050
0,5
-60
Defornación, mm .
0,6 -0,100
-80 Mín
Deformación Tangencial .
0,7 -100 -0,150
0,30 7,0 90
6,0 80
0,25
70
5,0
0,20 60
4,0 50
0,15
3,0 40
σ 'v/σ 'ov, %
0,10 30
2,0
20
Amortiguamiento .
0,05 1,0 10
Amortiguamiento .
0,05 1,0 Rampa Kpa 0,00
Deformación, mm .
50 20
Tensión Tangencial y
Presión de Poros, kPa.
1,6
Unitaria
-0,050 Permanente
-50
-0,100
2,0
Defornación, mm .
-100 -0,150
Mín
Deformación Tangencial .
2,5 -150 -0,200
0,286 3,5 90
0,284 80
0,282 3,0
70
0,280 2,5
0,278 60
0,276 2,0 50
0,274 1,5 40
0,272
σ 'v/σ 'ov, %
30
0,270 1,0
20
0,268
Amortiguamiento .
0,5 10
0,266
Amortiguamiento .
0,268 0,5
0,266 Rampa Kpa 0,00
Deformación, mm .
50 20
Tensión Tangencial y
Presión de Poros, kPa.
1,6
0,0 80 0,200
CONSOLIDACION CORTE SIMPLE CICLICO 60 0,150 Cíclica
0,5 0,100
40 Máx.
1,0 20 0,050
0 0,000 Permanente
1,5
Unitaria
-20 -0,050
-40 -0,100
2,0
Defornación, mm .
-60 -0,150
Mín
Deformación Tangencial .
2,5 -80 -0,200
0,10 2,0 30
20
Amortiguamiento .
0,05 1,0
10
Amortiguamiento .
0,05 1,0
Rampa Kpa 0,00
Deformación, mm .
50 20
Tensión Tangencial y
Presión de Poros, kPa.
1,6
Unitaria
1,2 -0,050 Permanente
-40
1,4 -0,100
-60
Defornación, mm .
1,6 -80 -0,150
Mín
Deformación Tangencial .
1,8 -100 -0,200
0,10 2,0 30
20
Amortiguamiento .
0,05 1,0
10
Amortiguamiento .
0,05 1,0
Rampa Kpa 0,00
Deformación, mm .
Tensión Tangencial y
3,5 10
0,0 40 0,040
CONSOLIDACION CORTE SIMPLE CICLICO Máx.
30 0,030
1,0
20 0,020
2,0 Cíclica
10
0,010
3,0 0 Permanente
0,000
Unitaria
4,0 -10
-20 -0,010
Defornación, mm .
5,0 -0,020
-30 Mín
Deformación Tangencial .
6,0 -40 -0,030
1,5 30
0,05 1,0 20
Amortiguamiento .
0,5 10
Amortiguamiento .
0,05
0,5 Rampa Kpa 0,00
Deformación, mm .
Tensión Tangencial y
3,5 10
0,0 60 0,200
0,5 CONSOLIDACION CORTE SIMPLE CICLICO 0,150
40 Máx.
1,0 20 0,100
Cíclica
1,5 0,050
0
2,0 0,000
-20
Unitaria
2,5 -0,050 Permanente
3,0 -40 -0,100
Defornación, mm .
3,5 -60 -0,150 Mín
Deformación Tangencial .
4,0 -80 -0,200
0,10 1,5 30
1,0 20
Amortiguamiento .
0,05
0,5 10
Amortiguamiento .
0,05
0,5 Rampa Kpa 0,00
Deformación, mm .
Tensión Tangencial y
3,5 10
0,0 80 0,200
0,5 CONSOLIDACION CORTE SIMPLE CICLICO 60 0,150 Máx.
40 0,100 Cíclica
1,0
1,5 20 0,050
0
2,0 0,000
-20
Unitaria
2,5 -0,050
-40 Permanente
3,0 -60 -0,100
Defornación, mm .
3,5 -80 -0,150 Mín
Deformación Tangencial .
4,0 -100 -0,200
0,31 3,5 90
0,30 80
3,0
0,30 70
0,29 2,5
60
0,29 2,0 50
0,28
1,5 40
0,28
σ 'v/σ 'ov, %
0,27 30
1,0
0,27 20
Amortiguamiento .
0,5 10
0,26
Amortiguamiento .
0,5 Rampa Kpa 0,00
0,26
Módulo de Rigidez, MPa .
0,26 0,0 Tensión tangencia cíclica +/- kPa 68,00
-0,080 -0,070 -0,060 -0,050 -0,040 -0,030 -0,020 -0,010 0,000 -0,080 -0,060 -0,040 -0,020 0,000 Archivos Originales Descripción Arcilla limosa marrón, consistencia blanda.
SSE_SC03_SH2_350 Observaciones
Deformación Tangencial Unitaria Acumulada Deformación Tangencial Unitaria Acumulada SSC_SC03_SH2_340_0_68_EQU
SSC_SC03_SH2_340_0_68_CYC
0,0 250 140
CONSOLIDACION - CORTE SIMPLE ESTATICO u
0,2 CORTE SIMPLE ESTATICO 120
200
0,4 100
150
0,6 τ 80
100
0,8 60
50
1,0 40
0 CORTE SIMPLE ESTATICO
Deformación, mm .
1,2 20
Tensión Tangencial y
Presión de Poros, kPa.
Tensión Tangencial, (kPa) .
1,4 -50 0
1 10 100 1000 10000 0 5 10 15 20 25 30 0 50 100 150 200 250 300 350 400
0,0 50 0,015
CONSOLIDACION CORTE SIMPLE CICLICO 40 Máx.
Cíclica
0,5 30 0,010
20 Permanente
1,0 10 0,005
0
1,5 -10 0,000
Unitaria
-20
2,0 -30 -0,005
Defornación, mm .
-40 Mín
Deformación Tangencial .
2,5 -50 -0,010
0,06 1,5 30
0,04 1,0 20
Amortiguamiento .
0,02 0,5 10
Amortiguamiento .
0,02 Rampa Kpa 0,00
Deformación, mm .
1,2 20
Tensión Tangencial y
Presión de Poros, kPa.
Tensión Tangencial, (kPa) .
1,4 -50 0
1 10 100 1000 10000 0 5 10 15 20 25 30 0 50 100 150 200 250 300 350 400
0,0 80 0,150
CONSOLIDACION CORTE SIMPLE CICLICO 60 Cíclica
0,5 0,100
40 0,050 Máx.
1,0
20
1,5 0,000
0 Permanente
2,0 -0,050
Unitaria
-20
2,5 -40 -0,100
Defornación, mm .
3,0 -60 -0,150
Mín
Deformación Tangencial .
3,5 -80 -0,200
0,10 2,0 30
20
Amortiguamiento .
0,05 1,0
10
Amortiguamiento .
0,05 1,0
Rampa Kpa 0,00
Deformación, mm .
1,2 20
Tensión Tangencial y
Presión de Poros, kPa.
Tensión Tangencial, (kPa) .
1,4 -50 0
1 10 100 1000 10000 0 5 10 15 20 25 30 0 50 100 150 200 250 300 350 400
Unitaria
2,0 -40 -0,100
-60
Defornación, mm .
2,5 -0,150
-80 Mín
Deformación Tangencial .
3,0 -100 -0,200
0,30 4,0 90
3,5 80
0,25
3,0 70
0,20 60
2,5
50
0,15 2,0
40
1,5
σ 'v/σ 'ov, %
0,10 30
1,0 20
Amortiguamiento .
0,05 0,5 10
Amortiguamiento .
0,05 0,5 Rampa Kpa 0,00
ÍNDICE
Págs.
1.- INTRODUCCIÓN
-1-
Estudio dinámico de cimentaciones en obras portuarias
Dado que los resultados presentan aspectos similares en los 6 cálculos rea-
lizados, se comentan en lo que sigue los correspondientes a la relación
/σ'v0 = 0,20 para cada uno de los materiales.
-2-
APÉNDICE II.- CÁLCULOS NUMÉRICOS
80
60
40
Tensión tangencial (kPa)
20
-20
-40
-60
-80
-3 -2.5 -2 -1.5 -1 -0.5 0 0.5 1 1.5 2
Deformación tangencial (%)
80
60
40
Tensión tangencial (kPa)
20
-20
-40
-60
-80
0 50 100 150 200 250 300 350 400
Tensión efectiva vertical (kPa)
-3-
Estudio dinámico de cimentaciones en obras portuarias
1.5
1
Deformación tangencial (%)
0.5
-0.5
-1
0 20000 40000 60000 80000 100000 120000 140000
Pasos de cálculo
-4-
APÉNDICE II.- CÁLCULOS NUMÉRICOS
' v
ru 1
' v0
0.9
0.8
0.7
0.6
0.5
ru
0.4
0.3
0.2
0.1
0
0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20
Número de ciclos
En la figura anterior se han puesto, a partir del 15º ciclo de carga los valo-
res máximo (azul) y mínimo (rojo) entre los que varía el valor de ru en los ci-
clos cuando el modelo se encuentra en rotura.
-5-
Estudio dinámico de cimentaciones en obras portuarias
0.9
0.8
0.7
0.6
G/Gmax
0.5
0.4
0.3
0.2
0.1
0
0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20
Número de ciclos
Al igual que en el caso anterior se ha puesto, a partir del 15º ciclo de car-
ga, los valores máximo (azul) y mínimo (rojo) entre los que varía el valor de
G/Gmax en los ciclos cuando el modelo se encuentra en rotura.
Al igual que en el caso del material arenoso, se han realizado ensayos con
modelo numérico con las características de los materiales adoptadas para los
materiales limosos para las mismas condiciones de tensión inicial y variación
de tensión tangencial.
-6-
APÉNDICE II.- CÁLCULOS NUMÉRICOS
80
60
40
Tensión tangencial (kPa)
20
-20
-40
-60
-80
-10 -8 -6 -4 -2 0 2 4 6 8 10
Deformación tangencial (%)
80
60
40
Tensión tangencial (kPa)
20
-20
-40
-60
-80
0 50 100 150 200 250 300 350 400
Deformación tangencial (%)
-7-
Estudio dinámico de cimentaciones en obras portuarias
2.5
2
Deformación tangencial (%)
1.5
0.5
-0.5
0 20000 40000 60000 80000 100000 120000 140000 160000 180000
Pasos de cálculo
-8-
APÉNDICE II.- CÁLCULOS NUMÉRICOS
0.9
0.8
0.7
0.6
0.5
ru
0.4
0.3
0.2
0.1
0
0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15
Número de ciclos
En la figura anterior se han puesto, a partir del ciclo de carga nº 11, los
valores máximo (azul) y mínimo (rojo) entre los que varía el valor de ru en los
ciclos cuando el modelo se encuentra en rotura.
-9-
Estudio dinámico de cimentaciones en obras portuarias
0.9
0.8
0.7
0.6
G/Gmax
0.5
0.4
0.3
0.2
0.1
0
0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20
Número de ciclos
Al igual que en los casos anteriores se ha puesto, a partir del ciclo de carga
número 11, los valores máximo (azul) y mínimo (rojo) entre los que varía el
valor de G/Gmax en los ciclos cuando el modelo se encuentra en rotura.
- 10 -
APÉNDICE II.- CÁLCULOS NUMÉRICOS
ya que, para ilustrar los resultados se ha considerado que ese caso podría re-
sultar representativo.
Para esta carga cíclica se obtiene que la rotura se produce a partir del ciclo
de carga nº 20. En la Figura nº 11 se presenta el desplazamiento vertical del
cajón en función del número de ciclos obtenidos en el cálculo.
- 11 -
Estudio dinámico de cimentaciones en obras portuarias
Número de ciclos
0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25
0
-0.2
-0.4
-0.6
Asiento del cajón (m)
-0.8
-1
-1.2
-1.4
-1.6
-1.8
-2
Figura nº 11.- Desplazamientos verticales del cajón en función del
número de ciclos. Banqueta 6 m de altura. Empuje ±600 kN/ml
- 12 -
APÉNDICE II.- CÁLCULOS NUMÉRICOS
- 13 -
Estudio dinámico de cimentaciones en obras portuarias
180
160
140
Sobrepresiones intersticiales (kPa)
120
100
80
60
40
20
0
0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25
Numero de ciclos
10
0
Tensión tangencial (kPa)
-10
Promedio
-20
-30
-40
0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25
Numero de ciclos
- 14 -
APÉNDICE II.- CÁLCULOS NUMÉRICOS
Para esta carga cíclica se obtiene que la rotura, teniendo en cuenta el cri-
terio adoptado de que el movimiento vertical en la base del cajón sea de 1
metro, se produce a partir del ciclo de carga nº10. En la Figura nº 16 se pre-
senta el desplazamiento vertical del cajón en función del número de ciclos ob-
tenidos en el cálculo.
Número de ciclos
0 1 2 3 4 5 6 7 8 9 10 11
0
-0.2
-0.4
-0.6
Asiento del cajón (m)
-0.8
-1
-1.2
-1.4
-1.6
-1.8
-2
Figura nº 16.- Desplazamientos verticales del cajón en función del
número de ciclos. Banqueta 6 m de altura. Empuje ±800 kN/ml
- 15 -
Estudio dinámico de cimentaciones en obras portuarias
- 16 -
APÉNDICE II.- CÁLCULOS NUMÉRICOS
- 17 -
Estudio dinámico de cimentaciones en obras portuarias
140
120
Sobrepresiones intersticiales (kPa)
100
80
60
40
20
0
0 1 2 3 4 5 6 7 8 9 10 11
Numero de ciclos
10
0
Tensión tangencial (kPa)
-10
-20
-30
-40
0 1 2 3 4 5 6 7 8 9 10 11
Numero de ciclos
- 18 -
APÉNDICE II.- CÁLCULOS NUMÉRICOS
Para esta carga cíclica se obtiene que la rotura, teniendo en cuenta el cri-
terio adoptado de que el movimiento vertical en la base del cajón sea de 1
metro, se produce a partir del ciclo de carga nº5. En la Figura nº 21 se pre-
senta el desplazamiento vertical del cajón en función del número de ciclos ob-
tenidos en el cálculo.
Número de ciclos
0 1 2 3 4 5 6
0
-0.2
-0.4
-0.6
Asiento del cajón (m)
-0.8
-1
-1.2
-1.4
-1.6
-1.8
-2
Figura nº 21.- Desplazamientos verticales del cajón en función del
número de ciclos. Banqueta 6 m de altura. Empuje ±1000 kN/ml
- 19 -
Estudio dinámico de cimentaciones en obras portuarias
- 20 -
APÉNDICE II.- CÁLCULOS NUMÉRICOS
- 21 -
Estudio dinámico de cimentaciones en obras portuarias
140
120
Sobrepresiones intersticiales (kPa)
100
80
60
40
20
0
0 1 2 3 4 5 6
Numero de ciclos
20
10
Tensión tangencial (kPa)
-10
-20
-30
-40
-50
0 1 2 3 4 5 6
Numero de ciclos
- 22 -
APÉNDICE II.- CÁLCULOS NUMÉRICOS
Para esta carga cíclica se obtiene que la rotura, con el criterio adoptado
definido anteriormente, se produce a partir del ciclo de carga nº45. En la Figu-
ra nº 26 se presenta el desplazamiento vertical del cajón en función del núme-
ro de ciclos obtenidos en el cálculo.
Número de ciclos
0 5 10 15 20 25 30 35 40 45
0
-0.1
-0.2
-0.3
Asiento del cajón (m)
-0.4
-0.5
-0.6
-0.7
-0.8
-0.9
-1
Figura nº 26.- Desplazamientos verticales del cajón en función del
número de ciclos. Banqueta 8 m de altura. Empuje ±700 kN/ml
- 23 -
Estudio dinámico de cimentaciones en obras portuarias
- 24 -
APÉNDICE II.- CÁLCULOS NUMÉRICOS
- 25 -
Estudio dinámico de cimentaciones en obras portuarias
250
200
Sobrepresiones intersticiales (kPa)
150
100
50
0
0 5 10 15 20 25 30 35 40 45
Numero de ciclos
30
20
10
Tensión tangencial (kPa)
-10
Promedio
-20
-30
-40
0 5 10 15 20 25 30 35 40 45
Numero de ciclos
- 26 -
APÉNDICE II.- CÁLCULOS NUMÉRICOS
Para esta carga cíclica se obtiene que la rotura, con el criterio adoptado
definido anteriormente, se produce a partir del ciclo de carga nº 29.
Número de ciclos
0 5 10 15 20 25 30
0
-0.1
-0.2
-0.3
Asiento del cajón (m)
-0.4
-0.5
-0.6
-0.7
-0.8
-0.9
-1
Figura nº 31.- Desplazamientos verticales del cajón en función del
número de ciclos. Banqueta 8 m de altura. Empuje ±800 kN/ml
- 27 -
Estudio dinámico de cimentaciones en obras portuarias
- 28 -
APÉNDICE II.- CÁLCULOS NUMÉRICOS
250
200
Sobrepresiones intersticiales (kPa)
150
100
50
0
0 5 10 15 20 25 30
Numero de ciclos
30
20
10
Tensión tangencial (kPa)
-10
Prome…
-20
-30
-40
0 5 10 15 20 25 30
Numero de ciclos
- 29 -
Estudio dinámico de cimentaciones en obras portuarias
Para esta carga cíclica se obtiene que la rotura, con el criterio adoptado, se
produce a partir del ciclo de carga nº 14.
Número de ciclos
0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15
0
-0.1
-0.2
-0.3
Asiento del cajón (m)
-0.4
-0.5
-0.6
-0.7
-0.8
-0.9
-1
Figura nº 36.- Desplazamientos verticales del cajón en función del
número de ciclos. Banqueta 8 m de altura. Empuje ±1000 kN/ml
- 30 -
APÉNDICE II.- CÁLCULOS NUMÉRICOS
- 31 -
Estudio dinámico de cimentaciones en obras portuarias
180
160
140
Sobrepresiones intersticiales (kPa)
120
100
80
60
40
20
0
0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15
Numero de ciclos
30
20
10
Tensión tangencial (kPa)
-10
Promedio
-20
-30
-40
0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21
Numero de ciclos
- 32 -
APÉNDICE II.- CÁLCULOS NUMÉRICOS
Para esta carga cíclica se obtiene que la rotura, con el criterio adoptado, se
produce a partir del ciclo de carga nº9.
Número de ciclos
0 1 2 3 4 5 6 7 8 9 10
0
-0.1
-0.2
-0.3
Asiento del cajón (m)
-0.4
-0.5
-0.6
-0.7
-0.8
-0.9
-1
Figura nº 41.- Desplazamientos verticales del cajón en función del
número de ciclos. Banqueta 8 m de altura. Empuje ±1200 kN/ml
- 33 -
Estudio dinámico de cimentaciones en obras portuarias
- 34 -
APÉNDICE II.- CÁLCULOS NUMÉRICOS
180
160
140
Sobrepresiones intersticiales (kPa)
120
100
80
60
40
20
0
0 1 2 3 4 5 6 7 8 9 10
Numero de ciclos
30
20
10
Tensión tangencial (kPa)
-10
-20 Promedio
-30
-40
-50
0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16
Numero de ciclos
- 35 -
Estudio dinámico de cimentaciones en obras portuarias
Para este caso de carga se obtiene la rotura a partir del ciclo de carga
nº40. A continuación se presentan las mismas figuras utilizadas en los aparta-
dos anteriores para resumir los resultados obtenidos.
Número de ciclos
0 5 10 15 20 25 30 35 40 45
0
-0.1
-0.2
-0.3
Asiento del cajón (m)
-0.4
-0.5
-0.6
-0.7
-0.8
-0.9
-1
Figura nº 46.- Desplazamientos verticales del cajón en función del
número de ciclos. Banqueta 10 m de altura. Empuje ±900kN/ml
- 36 -
APÉNDICE II.- CÁLCULOS NUMÉRICOS
- 37 -
Estudio dinámico de cimentaciones en obras portuarias
250
200
Sobrepresiones intersticiales (kPa)
150
100
50
0
0 5 10 15 20 25 30 35 40 45
Numero de ciclos
10
0
Tensión tangencial (kPa)
-10
Promedio
-20
-30
-40
0 5 10 15 20 25 30 35 40 45
Numero de ciclos
- 38 -
APÉNDICE II.- CÁLCULOS NUMÉRICOS
Para este caso de carga se obtiene la rotura a partir del ciclo de carga
nº30. A continuación se presentan las mismas figuras utilizadas en los aparta-
dos anteriores para resumir los resultados obtenidos. El punto en el que se
representa la evolución de las presiones intersticiales (punto 1 en los gráficos)
se ha tomado, para este y los siguientes casos de empuje, (siempre en el eje
del cajón) a 15 metros de profundidad con respecto al contacto banqueta-
cimiento. Si bien en ese punto no se alcanzan las máximas sobrepresiones, si
se observa la evolución de las sobrepresiones en función del número de ciclos,
que es similar a la que se produce en la zona de contacto banqueta-cimiento
(para el eje del cajón).
Número de ciclos
0 5 10 15 20 25 30 35
0
-0.1
-0.2
-0.3
Asiento del cajón (m)
-0.4
-0.5
-0.6
-0.7
-0.8
-0.9
-1
Figura nº 51.- Desplazamientos verticales del cajón en función del
número de ciclos. Banqueta 10 m de altura. Empuje ±1000 kN/ml
- 39 -
Estudio dinámico de cimentaciones en obras portuarias
- 40 -
APÉNDICE II.- CÁLCULOS NUMÉRICOS
90
80
70
Sobrepresiones intersticiales (kPa)
60
50
40
30
20
10
0
0 5 10 15 20 25 30
Numero de ciclos
20
10
Tensión tangencial (kPa)
-10
-20 Promedio
-30
-40
-50
0 5 10 15 20 25 30
Numero de ciclos
- 41 -
Estudio dinámico de cimentaciones en obras portuarias
Para este caso de carga se obtiene la rotura a partir del ciclo de carga
nº20. A continuación se presentan las mismas figuras utilizadas en los aparta-
dos anteriores para resumir los resultados obtenidos.
Número de ciclos
0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21
0
-0.1
-0.2
-0.3
Asiento del cajón (m)
-0.4
-0.5
-0.6
-0.7
-0.8
-0.9
-1
Figura nº 56.- Desplazamientos verticales del cajón en función del
número de ciclos. Banqueta 10 m de altura. Empuje ±1250 kN/ml
- 42 -
APÉNDICE II.- CÁLCULOS NUMÉRICOS
- 43 -
Estudio dinámico de cimentaciones en obras portuarias
80
70
60
Sobrepresiones intersticiales (kPa)
50
40
30
20
10
0
0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20
Numero de ciclos
30
20
10
Tensión tangencial (kPa)
-10
Promedio
-20
-30
-40
-50
0 5 10 15 20
Numero de ciclos
- 44 -
APÉNDICE II.- CÁLCULOS NUMÉRICOS
Para este caso de carga se obtiene la rotura a partir del ciclo de carga nº
14. A continuación se presentan las mismas figuras utilizadas en los apartados
anteriores para resumir los resultados obtenidos.
Número de ciclos
0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15
0
-0.1
-0.2
-0.3
Asiento del cajón (m)
-0.4
-0.5
-0.6
-0.7
-0.8
-0.9
-1
Figura nº 61.- Desplazamientos verticales del cajón en función del
número de ciclos. Banqueta 10 m de altura. Empuje ±1500 kN/ml
- 45 -
Estudio dinámico de cimentaciones en obras portuarias
- 46 -
APÉNDICE II.- CÁLCULOS NUMÉRICOS
100
90
80
Sobrepresiones intersticiales (kPa)
70
60
50
40
30
20
10
0
0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15
Numero de ciclos
30
20
10
Tensión tangencial (kPa)
-10
-20 Promedio
-30
-40
-50
-60
0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15
Numero de ciclos
- 47 -
APÉNDICE III
NOTAS SOBRE EL DAÑO ACUMULADO
APÉNDICE III.- NOTAS SOBRE EL DAÑO ACUMULADO
ÍNDICE
Págs.
El “daño” producido por una ola es tanto menor cuanto menor es su altura.
Ese daño se puede calificar con el número de olas N(H) que, en caso de ocurrir
con amplitud constante, produciría la ruina del dique.
Ho
H15
Altura de ola estacionaria
que rompe al dique, H
N = 15
·H0
·H0
-1-
Estudio dinámico de cimentaciones en obras portuarias
Ho H
lg10 N log10 15 Ho ≥ H ≥ H o
Ho H15
Ho H
lg10 D(1) log10 15
Ho H15
H H15
= ; = <1 , <1
Ho Ho
1
lg10 D 1,18
1
1
lg10 D(15) 1,18 1 1,18
1 1
-2-
APÉNDICE III.- NOTAS SOBRE EL DAÑO ACUMULADO
2 H
f(x) = 2 e x ; x=
Hrms
-3-
Estudio dinámico de cimentaciones en obras portuarias
1000
100
Número de olas >H
10
1
0 0.2 0.4 0.6 0.8 1 1.2
H/Hmax
f(x)
Daño (temporal) = N
dx
-4-
APÉNDICE III.- NOTAS SOBRE EL DAÑO ACUMULADO
100
10
b=0,6
Daño acumulado durante el temporal
b=0,8
0.1
0.01
Temporalde 1000 olas
0.001
0.8 1 1.2 1.4 1.6 1.8 2
Valor de H0 /Hmax
lg D(15) = 1,18
1
-5-
Estudio dinámico de cimentaciones en obras portuarias
b=0,6
Valor de r que produce el daño del temporal
Valor propuesto
b=0,8 r=0,65
0.5
0
1 1.5 2
Coeficiente de seguridad, F
Como puede verse el valor = 0,65, propuesto en esta tesis puede ser
adecuado.
Se verifica en este escrito que esto es cierto bajo las siguientes hipótesis
simplificativas:
-6-
APÉNDICE III.- NOTAS SOBRE EL DAÑO ACUMULADO
d) Se supone que existe una relación lineal entre la altura de ola y el valor
del empuje sobre el dique.
-7-