1

III.- MARCO TEÓRICO

socavación

Rocha (1998) menciona que. “En general se denomina erosión, socavación o

degradación a la remoción que efectúa el agua de las partículas solidas constituyentes de un

lecho fluvial y cuyo resultado es una profundización del cauce la erosión puede ser local o

generalizada” (p.201).

Aguirre (2006) afirma que. “El fenómeno de socavación se relaciona con dos de

los problemas más complejos de la mecánica de fluidos, como son la mecánica de

transportes de sedimentos y· la capa límite de tridimensional” (p. 45).

Martinez (2007) menciona que. “en el estudio de la socavacion de obras civiles, el

problema principal es el colapso de obras hidraulicas por las grandes recargas hidraulicas

de los rios, dando lugar perdida de vidas y economicas” (p.93).

Rodriguez (2010) menciona que. “socavación en cauces, resulta fundamental

evaluar en forma pormenorizada la profundización del lecho producida por las estructuras

que se implantan y alteran las condiciones del flujo en un cauce natural” (p.26).

Trujillo (2015) menciona que. “La socavación es el resultado de la acción erosiva

del flujo de agua que arranca y acarrea material de lecho y de las bancas de un cauce,

convirtiéndose, por ejemplo, en una de las causas más comunes de falla en puentes” (p.38).

Enciso (2015) menciona que. “En general se denomina erosión, socavación o

degradación a la remoción que efectúa el agua de las partículas sólidas constituyentes de un

lecho fluvial y cuyo resultado es una profundización del cauce la erosión puede ser local o

generalizada” (p.8).
 2

Venegas (2011) menciona que. “La socavación, como fenómeno hidráulico, se

origina en movimientos vorticosos que ocurren al pie de determinadas estructuras

hidráulicas, como un pilar y estribo de puente, o en la descarga de un vertedero.” (p.14).

Chirino (2017) menciona que. “La socavación es un fenómeno natural, que se debe

al aumento de la velocidad del flujo, a causa del incremento del caudal, las partículas que se

encuentran en el fondo del rio entran en inestabilidad. Esto provoca que el suelo que forma

el cauce del rio no pueda soportar ninguna carga” (p.7).

3.1 Tipos de socavación:

Los tipos de socavación mas importante que se consideran en los estudios son los

siguientes:

 Socavación general.

 Socavación transversal.

 Socavación en curvas.

 Socavación local.

 Socavación aguas debajo de grandes presas.

3.1.1. Socavación general.

Rocha (1998) menciona que. “La erosión generalizada, o degradación, se produce a

lo largo y ancho de un cauce fluvial. La degradación del cauce es extensa y progresiva”

(p.202).

Apaclla (2014) afirma que.

Consiste en una disminución generalizada del nivel del fondo como consecuencia

del incremento de la capacidad de flujo para transportar material durante las

 3

avenidas. También se debe al desequilibrio ocasionado por la diferencia de

capacidad de transporte entre dos secciones consecutivas que definen a un tramo;

aquí las erosiones tienen lugar cuando el volumen de sedimento que sale por la

sección de aguas abajo en un cierto periodo es mayor que el que entra al tramo por

la sección aguas arriba en el mismo periodo. Este fenómeno es un proceso natural

que puede ocurrir a todo lo largo del rio y no es provocado por factores humanos.

APACLLA R. (2014). Hidráulica fluvial. Fondo editorial Universidad Nacional

Agraria la Molina, lima – Perú.

3.1.2 Socavación transversal.

Apaclla (2014) menciona que.

La reducción del ancho de la sección de un cauce, sea ocasionada en forma natural

o como consecuencia de una obra, es compensada con un incremento en la

profundidad hasta el punto en el cual se alcanza la capacidad necesaria en la

sección; este incremento de la profundidad del cauce se conoce como socavación

transversal. (p.170)

3.1.3 Socavación en curvas.

Apaclla (2014) afirma que.

Se genera en el fondo cercano al lado cóncavo de la curva o talud exterior debido al

flujo helicoidal que se forma cuando el rio cambia de dirección. Esto se debe a la

sobreelevación del nivel del agua en esta zona producida por la fuerza centrífuga. Aun

cuando no sea inducida por el hombre debe tomarse en cuenta que puede incrementarse al
 4

proteger, con enrocados o protecciones marginales, las curvas para estabilizar sus

márgenes. (p.170)

Apaclla (2014) menciona que

Generalmente la profundidad máxima se observa en el tramo de la mitad de la curva

hacia aguas abajo; sin embargo (y por su posición), depende de la dirección del flujo a la

entrada y su trayectoria a lo largo de la curva. Para problemas prácticos se considera que

puede existir en cualquier sección de la curva. (p.170)

3.1.4 Socavación local.

Apaclla (2014) afirma que.

Se presenta al pie de las estructuras interpuestas a las corrientes, sumergidas o que

emergen de la superficie del agua, como resultado de la deflexión de las líneas de flujo, la

turbulencia y la vorticidad provocada por la presencia del obstáculo. Como ejemplo, se

tiene los pilares, estribos, espigones, etc. Esta clase de erosión puede manifestarse en dos

tipos de estructuras: en las rodeadas por el flujo, como es el caso de los pilares de un

puente, y en estructuras pegadas a la margen, como en caso de estribos y espigones. (p.170)

3.1.5 Socavación de aguas debajo de grandes presas.

Apaclla (2014) afirma que.

Consiste en descenso gradual del fondo inducido por la interrupción del transporte

de sedimentos de aguas arriba, ocasionada, a su vez, por la presencia de la cortina de la

presa, en donde se retiene la mayor parte del sedimento que entra al vaso, lo anterior

permite el incremento de la capacidad de transporte de las partículas del fondo en el tramo

aguas abajo, las cuales, además, no pueden ser reemplazadas. (p.171)

 5

Pareja (2010) afirma que, “Esta sección presenta los parámetros que tienen

influencia en la socavación de materiales o mezclas no uniformes y se deriva de relaciones

para cuantificar las profundidades de socavación.”(p.64)

3.2 Estimación de la socavación general

Apaclla (2014) menciona que.

La estimación de la magnitud de la socavación general es muy importante cuando

se pretende construir o colocar obras y equipos cercanos, o bien cruza una obra por el

cauce. La intención es estimar la sección de máxima erosión correspondiente a un gasto de

diseño, de tal forma que, al construir la obra, tal sección ojo tal obra no afecte ni sea

afectado seriamente. (p.171)

Apaclla (2014) afirma que.

Para el cálculo de socavación general se recomienda utilizar el método de Lischtvan

y Lebediev, el cual se basa en determinar la condición de equilibrio entre la velocidad

media de la corriente y la velocidad media del flujo necesaria para erosionar el material que

forma el cauce. Es aplicable tanto para materiales del subsuelo con distribución homogénea

o heterogénea. (p.171)

La condición de equilibrio está dada por:

𝑈𝑒 = 𝑈𝑟

Donde:

𝑈𝑒 = velocidad media de la corriente necesaria para erosionar el material del fondo (inicio

de arrastre), en m/s.

𝑈𝑟 = velocidad media “real” de la corriente para el gasto de diseño, en m/s.

Para poder aplicar el método se requiere los datos siguientes:

 6

 El gasto de diseño, Qd y su periodo de retorno.

 Elevación del agua en la sección de estudio para el gasto de diseño

 El perfil de la sección transversal o del tramo en estudio obtenido durante el estiaje

anterior.

 Estratigrafía del material bajo el fondo.

 Si el suelo es granular, la granulometría del material del fondo y su diámetro medio,

Dm.

 Si el suelo es cohesivo, el peso volumetrico seco 𝛾𝑠 , que es cociente del peso de la

muestra seca entre el volumen de la muestra inalterada al recolectarse.

Cuando la sección esta bajo un puente se requiere además la siguiente información:

 Geometría del puente y su cimentación.

 Dirección de las líneas de flujo.

3.2.1 Calculo de la velocidad media real (ur)

Apaclla (2014) menciona que. “La hipótesis fundamental consiste en suponer que el caudal

que pasa por cualquier franja de ancho fijo de la sección permanece constante durante el

proceso de la erosión” (p.172).

La variación de la velocidad media del flujo 𝑈𝑟 , como función de la profundidad en un

punto de la sección puede ser obtenida mediante el análisis de una franja vertical de ancho

∆𝐵 (ver fig. 3.1)

 7

Figura 1. Sección transversal. - Variables utilizados para determinar Ur y Ue

La ecuación de Manning aplicada a una franja vertical, de ancho ∆𝐵, para el caudal

correspondiente ∆𝑄 es:

1 1/2 5/3
∆𝑄 = 𝑈∆𝐴 = 𝑆 𝐷0 ∆𝐵
𝑛

Donde:

S = gradiente hidraulico

𝐷0 = profundidad inicial de la franja analizada, medida sobre su eje vertical entre el nivel

del agua correspondiente al caudal de diseño y el nivel del fondo correspondiente al estiaje,

en m.

B = ancho local de la sección. Para encontrar B se traza una línea perpendicular a las líneas

de corriente proyectada sobre la superficie libre del agua. Asi se toma en cuenta además el

esviajamiento de la corriente.
 8

Ya que la rugosidad n y el gradiente hidraulico se consideran constantes en la sección, el

1
producto 𝑛 𝑆 1/2 es constante para cualquier punto y se designa con la letra, luego:

∆𝑄 = 𝛼𝐷0 5/3 ∆𝐵

El valor de se puede ser expresado en forma general como una función de la profundidad

media dm del flujo en la sección antes de la socavación, la velocidad media en la sección

𝑈𝑟 , y el caudal de diseño Qd:

1 1/2 5/3
𝑄𝑑 = 𝑆 𝐷𝑚 𝐵𝑒
𝑛

Donde:

Be = ancho efectivo en la sección, descontado todos los obstáculos, en metro. Para

encontrar Be se traza una línea perpendicular a las líneas de corriente, sobre esta línea se

proyecta todos los obstáculos y Be es la suma de todos los espacios libres. Así se toma en

cuenta el esviajamiento de la corriente (ver figura 3.2). es muy importante conocer Be

cuando existe una obra en la sección de estudio, cuando no existe obstáculos, Be = B.

Dm = profundidad media de la sección que resulta de dividir su área hidráulica por el ancho

efectivo Be en metros: A/Be.

Figura 2. Determinación del ancho efectivo Be

 9

Apaclla (2014) menciona que “considerando la turbulencia que el flujo produce cerca de

los pilares y estribos, es necesario aplicar en la fórmula el coeficiente de contracción 𝜇:”

(p.174)

1 1/2 5/3
𝑄𝑑 = 𝑆 𝐷𝑚 𝐵𝑒
𝑛

Apaclla (2014) menciona que “𝜇 es un coeficiente que considera el efecto de contracción

producido por los pilares y los estribos. Se puede calcular como función de la velocidad

media del flujo y el claro entre los pilares y estribos L utilizando la siguiente expresión:”

(p.174)

0.387𝑈𝑜
𝜇 =1−
𝐿

𝑄𝑑
𝑈𝑜 =
𝐴

Al sustituir en la ecuación:

𝑄𝑑 = 𝜇𝛼𝐷𝑚 5/3 𝐵𝑒

𝑄𝑑
𝛼=
𝜇𝐷𝑚 5/3 𝐵𝑒

Apaclla (2014) afirma que “para una franja de ancho 𝛥𝐵, cuando 𝑑𝑜 aumenta a cualquier

valor de profundidad _ de socavación, 𝑑𝑠 , la velocidad decrece a un nuevo valor 𝑈𝑟 y 𝛥𝑄

se puede expresar como función de la velocidad y la profundidad:” (P.174)

𝛥𝑄 = 𝑈𝑟 𝐷𝑠 𝛥𝐵

Apaclla (2014) menciona que “𝐷𝑠 es la profundidad después de producirse la socavación

del fondo, medida desde la superficie libre del agua al pasar la avenida hasta el fondo

de la sección erosionada.” (p.174)

 10

𝛼𝐷0 5/3
𝑈𝑟 =
𝐷𝑠

Apaclla (2014) menciona que “el valor de 𝑈𝑟 propuesto en la ecuación () es válido

mientras el ancho de · la sección permanezca constante durante el tiempo de la avenida, y

el fondo descienda uniformemente en todo el ancho de la sección; la primera condición

generalmente se cumple, mientras que la segunda ocurre sólo en los cauces de ríos donde

todo el ancho del cauce principal se encuentra cubierto de agua durante la época de sequía.”

(p.174)

Apaclla (2014) dice que “esta es una de las deficiencias del método, ya que, para

cumplir completamente la condición establecida, todos los tramos de la sección deben

presentar una resistencia similar a la erosión; incluso para todos los estratos descubiertos.

Se pueden mencionar dos ejemplos sencillos donde esta condición no se cumple totalmente.”

(p.174)

 Apaclla (2014) menciona que “el primero es cuando el canal principal está cubierto de

vegetación, excepto en la zona donde fluye el agua con gastos bajos. Al presentarse la

avenida se alcanzan mayores velocidades en la zona desprotegida que en la zona con

vegetación, por lo tanto, será socavado con mayor rapidez incrementándose también el

gasto en esta zona.” (p.175)

 Apaclla (2014) dice que “el segundo ejemplo es cuando una sección está compuesta por

una parte de arcilla y otra de arena. Durante la avenida, la socavación será mayor en la

parte arenosa, ya que se requiere de 1nayor tie1npo para erosionar la parte arcillosa. Sin

embargo, cuando sucede algo similar, al1nenos es posible conocer que se tendrá una

socavación mayor en las zonas más débiles, aun cuando no sea posible evaluarla

exactamente.” (p.175)
 11

Apaclla (2014) dice que “una particularidad del primer ejemplo se presenta en los ríos donde

el ancho del cauce durante la época de sequía es menor que en la época de lluvias; aparecen

islas y zonas descubiertas sobre las que crece vegetación, lo cual incrementa la rugosidad y

disminuye la velocidad. Como resultado se tiene que en la época de avenidas la corriente

principal tiende a fluir sobre el cauce de la época de sequía; cuando esto sucede, la velocidad

del flujo se incrementa en esa parte de la sección. En estas zonas, con velocidades más altas,

se produce también mayor socavación que la esperada.” (p.175)

3.2.2 Cálculo de la velocidad media 𝑼𝒆

Apaclla (2014) menciona que “la velocidad mínima necesaria para arrastrar los materiales

erosionados, 𝑈𝑒 , depende de la naturaleza de los mismos.” (p.175)

 Suelos granulares o no cohesivos

𝑆𝑖: 0.05 𝑚𝑚 < 𝐷𝑚 < 2.6 𝑚𝑚

0.03
𝑈𝑒 = 4.7𝛽𝑑84 0.28 𝐷𝑠 0.322/𝑑𝑚

𝑆𝑖: 2.6 𝑚𝑚 < 𝐷𝑚 < 182 𝑚𝑚

0.092
𝑈𝑒 = 4.7𝛽𝑑84 0.28 𝐷𝑠 0.223/𝑑𝑚

𝑆𝑖: 182 𝑚𝑚 < 𝐷𝑚 < 1000 𝑚𝑚

0.187
𝑈𝑒 = 4.7𝛽𝑑84 0.28 𝐷𝑠 0.191/𝑑𝑚

 Suelos cohesivos
0.725
𝑈𝑒 = 0.000173𝛽𝛾𝑆1.18 𝐷𝑠 66.28/𝛾𝑠

𝛽 = 0.8416 + 0.03342𝐿𝑛(𝑇)

Donde:
 12

𝛽 ∶ Coeficiente que toma en cuenta el periodo de retorno T, correspondiente al

gasto de diseño.

𝛾𝑠 : Peso volumétrico seco en ton/m³.

𝐷𝑠 : Profundidad de socavación medida sobre el eje vertical de la franja analizada.

𝑑𝑚 : Diámetro medio del material de fondo

Apaclla (2014) menciona que “la velocidad 𝑈𝑒 , calculada con la ecuación () es muy similar

a la que se requiere para el inicio del arrastre de las partículas.” (p.175)

3.2.3 Cálculo de la socavación, 𝑫𝒔 , para suelos homogéneos

Apaclla (2014) menciona que “una vez conocido el tipo de suelo que existe en el sitio y

suponiendo que la rugosidad es constante en toda la sección, la profundidad hasta la que

llegará la socavación general se obtiene al igualar el valor de la velocidad media 𝑈𝑟 con el

valor de la velocidad 1nínima necesaria para arrastrar los materiales erosionados 𝑈𝑒 .” (p.176)

 Para suelos granulares

𝛼𝐷0 5/3 0.03

= 4.7𝛽𝑑84 0.28 𝐷𝑠 0.322/𝑑𝑚
𝐷𝑠

𝛼𝐷0 5/3 0.092

= 4.7𝛽𝑑84 0.28 𝐷𝑠 0.223/𝑑𝑚
𝐷𝑠

𝛼𝐷0 5/3 0.187

= 4.7𝛽𝑑84 0.28 𝐷𝑠 0.191/𝑑𝑚
𝐷𝑠

 Para suelos cohesivos

𝛼𝐷0 5/3 0.725

= 0.000173𝛽𝛾𝑆1.18 𝐷𝑠 66.28/𝛾𝑠
𝐷𝑠

Apaclla (2014) dice que “al despejar la profundidad de socavación 𝐷𝑠 , se obtendrán las

relaciones correspondientes para su cálculo.” (p.176)

 13

𝑆𝑖: 0.05 𝑚𝑚 < 𝐷𝑚 < 2.6 𝑚𝑚

𝑑𝑚 0.03 /(0.322+𝑑𝑚 0.03 )

𝛼𝐷0 5/3
𝐷𝑠 = [ ]
4.7𝛽𝑑84 0.28

𝑆𝑖: 2.6 𝑚𝑚 < 𝐷𝑚 < 182 𝑚𝑚

𝑑𝑚 0.092 /(0.223+𝑑𝑚 0.092 )

𝛼𝐷0 5/3
𝐷𝑠 = [ ]
4.7𝛽𝑑84 0.28

𝑆𝑖: 182 𝑚𝑚 < 𝐷𝑚 < 1000 𝑚𝑚

𝑑𝑚 0.187 /(0.191+𝑑𝑚 0.187 )

𝛼𝐷0 5/3
𝐷𝑠 = [ ]
4.7𝛽𝑑84 0.28

Para suelos cohesivo

𝛾𝑚 0.725 /(66.28+𝛾𝑚 0.725 )

5780𝛼𝐷0 5/3
𝐷𝑠 = [ ]
𝛽𝛾𝑚 1.18

Apaclla (2014) afirma que “se considera que existe rugosidad uniforme cuando la misma

sección es utilizada por el flujo en las épocas de estiaje y d avenidas por no haber cauces de

inundación o bien cuando estos últimos son pequeños y no crecen en ellos plantas o

arbustos.” (p.177)

Apaclla (2014) menciona que “las ecuaciones (- y -) se aplican en las franjas verticales en

que se divide la sección (ver Fig. 7.3); en cada una se obtiene una profundidad de socavación

𝑑𝑠 función de la profundidad inicial 𝑑0 . Al unir todos los puntos obtenidos se obtiene el perfil

teórico de la sección erosionada.” (p.177)

Apaclla (2014) menciona que “para el cálculo de la socavación general se presupone una

condición inicial de equilibrio, por lo que la socavación calculada ocurre cuando se presenta

el gasto de diseño.” (p.177)

 14

Apaclla (2014) dice que “si el tramo de río está sujeto a procesos de sedimentación y erosión,

el cálculo de la erosión será válido sólo para el periodo de avenidas que sigue al periodo seco

durante el que se realizó el levantamiento topográfico de la sección transversal utilizada en

el cálculo.” (p.177)

3.2.4 Cálculo de la socavación, 𝑫𝒔 , para suelos heterogéneos

Apaclla (2014) menciona que “un suelo heterogéneo está integrado por dos o más tipos de

materiales muy diferentes entre sí, conformados en capas. Puede consistir de un perfil

estratigráfico de capas intercaladas de material cohesivo y granular con diferentes

densidades, distribuciones de diámetro, etcétera. No importa el tipo de estratificación, la-

profundidad de equilibrio puede ser calculada aplicando las ecuaciones (- y -) en forma

analítica por el método de prueba y error o utilizando el método semigráfico.” (p.177)

Figura 3. Sección para el cálculo de ds, para suelos heterogéneos.

[Link]. Método de prueba y error

De la distribución estratigrafía de los materiales bajo una vertical, se escoge el

estrato superior y de acuerdo con la naturaleza del material se aplica la formula

correspondiente para determinar dS. Si la profundidad obtenida resulta bajo el nivel del
 15

límite inferior del estrato, se escoge el segundo estrato y se repite en forma similar el

cálculo anterior, considerando la formula correspondiente al tipo de suelo del segundo

estrato. En tanteo en el que la profundidad dS calculaba este dentro del espesor del estrato

en estudio, se habrá obtenido la dS buscada. El cálculo se debe hacer en forma ordenada,

iniciándolo por el estrato superior y continuando hacia los más profundos cuando la dS

calculada se encuentra a un nivel superior al del estrato analizado, significa que ese

material es muy resistente a la erosión, en ese caso, la profundidad de erosión corresponde

al nivel de la frontera superior del estrado considerado, en todos los tanteos se utiliza el

mismo valor de do. (Apaclla, 2014, p.178)

Figura 4. Calculo de la socavación para suelos heterogéneos. Método de prueba y error

[Link]. Método semigráfico

Apaclla (2014) menciona que “Para cada franja vertical se calculan las velocidades

Ur y Ue a diferentes profundidades. Los resultados obtenidos se grafican en un sistema

coordenado, con la velocidad sobre el eje horizontal y las profundidades sobre el eje
 16

vertical. El punto de intersección de ambas curvas determina la profundidad de equilibrio

de la socavación y la velocidad media correspondiente.” (Apaclla, 2014, p.178)

3.2.5. Calculo de la socavación general cuando la rugosidad no es uniforme en la sección

Apaclla (2014) menciona que “Cuando existen dos o más zonas con diferente

rugosidad a lo ancho de la sección en estudio, el procedimiento de cálculo es semejante,

con la diferencia de trabajar en forma aislada con cada franja y que para cada una debe

calcularse el coeficiente a, correspondiente.”

𝑄𝑑𝑖
𝛼𝑖 = (7.17)
𝜇𝑖 𝐷𝑖 5/3 𝐵𝑒𝑖

Figura 5. Calculo de la socavación para suelos heterogéneos. Método semigráfico.

Apaclla (2014) menciona que “Generalmente este problema se presenta cuando se

tiene un cauce de avenidas extenso y cubierto de vegetación, aunque también se presenta en

caso de suelos heterogéneos. El subíndice i de las variables que intervienen en ella se

refiere al valor para cada franja en que se divide la sección en estudio.”

Para calcular el gasto de cada franja se aplica la expresión siguiente:

 17

𝑄 𝐴 𝐶𝑖 √𝐷𝑖
𝑄𝑑𝑖 = ∑𝑛 𝑑(𝐴𝑒𝑖 (7.18)
𝑖=1 𝑒𝑖 𝐶𝑖 √𝐷𝑖 )

𝐷𝑖 1/6
𝐶𝑖 = (7.19)
𝑛𝑖

Donde:

Qdi = gasto que pasa por cada franja de la sección en m3/s

Aei = área hidráulica efectiva antes de la erosión de la sección en estudio, es decir, el área

hidráulica total menos el área proyectada por los obstáculos en el plano perpendicular al

flujo.

Qd = gasto total de diseño en m3/s

Ci = coeficiente de rugosidad de Chezy de cada franja.

Di = profundidad media en cada franja en m.

ni = coeficiente de rugosidad de Manning de cada franja.

3.2.6. Efecto en la socavación cuando el flujo tiene una alta concentración d sedimentos en

suspensión

Cuando las condiciones aguas arriba favorecen el transporte de materiales finos en

suspensión, existe una reducción en la profundidad de socavación calculada para una

determinada velocidad media. Esto se puede atribuir a que se requiere cierta turbulencia

para levantar las partículas del fondo. La turbulencia, a su vez, es función de la velocidad

del flujo dividida por la viscosidad cinemática del fluido. Cuando se tiene una gran

cantidad de material en suspensión, su peso específico y viscosidad se incrementan, lo cual

reduce el grado de turbulencia del flujo. (Apaclla, 2014, p. 180)

Apaclla (2014) menciona que “Por lo tanto, para un cierto flujo con sólidos en suspensión

es necesario incrementar la velocidad media con el fin de obtener el mismo grado de

socavación que con agua limpia. Lo anterior se logra introduciendo el coeficiente ѱ como
 18

producto este coeficiente depende del valor del peso específico del agua con sedimentos

suspendidos γ.”

Se puede calcular con la siguiente expresión.

𝛾 2
𝑚
ѱ = 0.38 + (1272) (7.20)

La profundidad de socavación general para material granular se obtiene utilizando las

ecuaciones siguientes:

Si 0.05 𝑚𝑚 < 𝑑𝑚 < 2.6 𝑚𝑚

0.03 )
𝑑𝑚 0.03/(0.322+𝑑𝑚
𝛼𝐷𝑜 5/3
𝐷𝑠 = [4.7𝛽ѱ𝑑 0.28 ] (7.21a)
𝑚

Si 2.6 𝑚𝑚 < 𝑑𝑚 < 182 𝑚𝑚

0.092 )
𝑑𝑚 0.092/(0.223+𝑑𝑚
𝛼𝐷𝑜 5/3
𝐷𝑠 = [4.7𝛽ѱ𝑑 0.28 ] (7.21b)
𝑚

Si 182 𝑚𝑚 < 𝑑𝑚 < 1000 𝑚𝑚

0.187 )
𝑑𝑚 0.187/(0.191+𝑑𝑚
𝛼𝐷𝑜 5/3
𝐷𝑠 = [4.7𝛽ѱ𝑑 0.28 ] (7.21c)
𝑚

Para suelos cohesivos

0.725 )
𝑑𝑚 0.725/(66.28+𝛾𝑚
5780𝛼𝐷𝑜 5/3
𝐷𝑠 = [ ] (7.22)
𝛽ѱ𝛾𝑚 1.18

3.3. Estimación de la socavación transversal

Apaclla (2014) menciona que “Se puede calcular con el método para la socavación general,

ya que se toman en cuenta las reducciones producidas dentro del valor del área efectiva, Ae
 19

y el ancho efectivo Be. Sin embargo, para tener una idea aproximada del valor de la

socavación transversal se puede utilizar la fórmula de Straub.”

Partiendo de que el gasto solido total por cualquier sección del tramo de rio en estudio debe

ser el mismo, es decir:

𝐺𝐵𝑂 = 𝐺𝐵1 (7.23)

Donde GB es el gasto solido transportado por el fondo y los subíndices 0 y 1 corresponden

a las secciones normal y reducida, respectivamente.

Apaclla (2014) menciona que “Al aplicar en la ecuación (7.11) para material granular, la

fórmula de Meyer, Peter y Müller para calcular el transporte de sedimentos, y se consideran

las condiciones hidráulicas durante la ocurrencia de la avenida.”

Se obtiene la siguiente ecuación:

𝐵 2/3
𝐷𝑚1 = 𝐷𝑚0 (𝐵𝑜 ) (7.24)
1

Donde:

Dmi = profundidad media de la sección analizada.

Bi = ancho efectivo de la sección analizada.

El subíndice 1 corresponde a la sección reducida.

El subíndice 0 corresponde a los valores en una sección inalterada localizada aguas arriba.

𝑆 𝐵 2/3
𝐷𝑚1 = 𝐷𝑚0 𝑆𝑜 ( 𝐵𝑂 ) (7.25)
1 1

La fórmula de Straub debe ser aplicada en fondos arenosos con una distribución

homogénea del material. Dada su simplicidad, es conveniente utilizarla para estimaciones

preliminares. Es útil para calcular la erosión en una sección reducida de ancho fijo, ya que

se supone que el ancho de la superficie libre permanece constante; es decir, que no se tiene

erosión lateral. (Apaclla, 2014, p.181)

 20

El fundamento del método de Straub es suponer que el transporte de sedimentos entre una

sección no alterada aguas arriba y la sección de ancho reducido es constante. Como

acciones preventivas de una erosión intensa localizada en las secciones de los puentes, es

conveniente mantener el cauce limpio de vegetación en una longitud al menos igual a una o

dos veces el ancho del cauce (2Be) hacia aguas arriba y a la mitad (Be/2) hacia aguas abajo.

(Apaclla, 2014, p.181)

3.4 Estimación de la socavación en curvas

Existen varias formas de resolver este problema. Si se dispone a la sección transversal en

estiaje, se puede aplicar el método de la socavación general ya que este indicara mayores

profundidades en el lado exterior de las curvas. Además, su aplicación presenta la ventaja de

que permite evaluar la profundidad máxima y la forma aproximada de la sección transversal

ya socavada.

Si no se tiene la sección de estiaje en las profundidades medias, se calculan las características

de la curva en planta, tales como el radio de curvatura r, medido al centro del cause, y el

ancho de la superficie libre B. Maza Albarez, propone la siguiente formula proporcionado

por Altunin:

𝐵 2.24
𝐷𝐶 = 1.9 ( ) 𝐷𝑚
𝑟

Donde:

Dm= Profundidad media en el tramo recto, situado aguas arriba de la curva, en metros.

Dc= Profundidad media a lo largo de la zona erosionada de la curva, en metros.

Además, en un momento dado y en una de las secciones de la curva extradós, la profundidad

máxima, Dc max., puede ser determinada con:

𝑟
Si 2.0 ≤ 𝐵 ≤ 5.91
 21

0.38𝑟
𝐷𝑐 = (3.73 − ) 𝐷𝑚
𝐵
𝑟
Si 𝐵 ≥ 5.91

1.26𝐵
𝐷𝑐 = (1.27 − ) 𝐷𝑚
𝑟

Existen otros métodos para evaluar las profundidades medias en curvas y meandros. La

máxima profundidad puede ocurrir en cualquier zona a lo largo de la curva, pero no en toda

al mismo tiempo.

3.5 Estimación de la socavación local en pilares de puente

Cuando un puente cruza un rio en una zona donde no es factible alcanzar un manto rocoso

en el que apoyar los pilares y estribos, el principal problema que se presenta (tanto en

proyecto como en mantenimiento) es el conocimiento de las erosiones locales que sufre el

fondo del cauce, que pueden ser de tal magnitud que lleguen a alcanzar la base de los pilares

y provocar la falla total de la estructura.

En este caso a sido estudiado por muchos investigadores, entre los que se encuentran: Jarocki

(1961), Maza y Sanchez (1963), Breusers (1967), Shen-Shneider y Karaki (1969) y Melville

(1976) entre otros.

El primero en de esos criterios es el propuesto por Lauren y Toch de acuerdo a los estudios

que realizaron en la universidad de Lowa y que después ha sido ampliado por Maggiolo,

Romitta y Souza Pinto entre otros. el segundo criterio que se debe a Varoslav es el resultado

de mediciones hechas en varios puentes de la unión soviética y ha sido corroborados por las

investigaciones de bata Andreu y otros investigadores para Lawrense Touch y su escuela, la

socavación depende fundamentalmente del tirante y ancho del Pilar Y en segundo término
 22

de la forma de este; si bien este parámetro no es significativo cuando la corriente hasta K el

pilar en forma oblicua en cambio la velocidad media de la corriente y el tamaño de los granos

del fondo no son factores que sean tomados en cuenta para determinar la profundidad de la

socavación con este criterio.

Por el contrario, para dialogar la socavación depende del cuadro de la velocidad media de la

corriente En primer término y tanto el tamaño de los granos del fondo como la geometría del

Pilar, aún en el caso de ataque oblicuo de la corriente son parámetros significativos de la

valuación de la profundidad de la socavación.

Resultados experimentales han mostrado que para un tirante y Pilar dados la profundidad de

la socavación concuerda bastante bien con la variación que indica hallazgos al variar la

velocidad con tal de la relación H sobre b/h tirante aguas arriba Pilar y de ancho Pilar sea

mayor que 1.5 y en tanto no se ha hallado una curva de la orense Touch; Esto es existe un

valor límite para la velocidad Más allá del cual la socavación no progresa Mientras no varía

en tirantes todo hace suponer entonces que con el criterio de lawrense Touch se determina el

intervalo de aplicabilidad de criterio dejar hablar de suerte que Los criterios antagónicos en

apariencia resultan ser complementarios.

El desconocimiento de este hecho puede conducir Sin embargo a resultados absurdos en

efecto en un río de montaña con velocidad muy grande el criterio de Yardas puede Conducir

a la predicción de socavaciones muy fuertes que no ocurrirán en la naturaleza en tanto que

un río de planicie de criterio lawrense es que es útil sólo en condiciones extremas conducirá

también a prever socavaciones muy grandes que nunca llegarán a presentarse.

Aplicar los métodos de cálculo de la socavación local es conveniente conocer el proceso con

la finalidad de tener una idea clara de la formación y el desarrollo de la socavación al pie de

los pilares de un puente.

 23

3.5.1 Proceso de la socavación local en pilares

Las siguientes observaciones fueron realizadas en un canal de laboratorio de hidráulica del

instituto de ingeniería de la Universidad Nacional Autónoma de México (UNAM).

Parámetros que influyen en la profundidad de socavación Al pie de un Pilar de fuentes se

pueden agrupar de acuerdo con su naturaleza en 5 diferentes clases:

 Parámetros hidráulicos

Velocidad media de flujo

Profundidad de flujo de la cara del Pinar

Distribución de velocidades

Dirección del flujo con respecto a los trajes de Los pilares

 Parámetros del material del fondo

Diámetro del material.

Curva granulométrica

Forma del grano

Grado de cohesión

Peso específico sumergido

Espesor de las capas del subsuelo

 Parámetros geométricos de Los pilares

Ancho del pilar.

Relación entre la longitud y el ancho del Pilar

Sección transversal del Pilar

 Características que dependen de la localización del puente

Contracción de la sesión del puente

 24

Radiación de la curva del tramo del río

Obras de control de fluido ya sea aguas arriba o aguas abajo

 Parámetro de tiempo

Duración del pico de la avenida

Tiempo requerido para remover el material y alcanzar una condición estable.

Debido a la forma en que Se realizaron los estudios se despreciaron los factores 10, 14, 16,

17 y 18 así también se desarrollaron experimentos con tres tipos de material granular y 4

tipos de Pilares y diferentes ángulos de incidencia en cada caso para un determinado material

tipo de Pilar y su localización sólo se varió la profundidad y la velocidad del flujo

[Link] Pilares rectangulares alineadas con el flujo

Comprender el mecanismo de la socavación se presenta un resumen de lo que sucede bajo

condiciones ideales es decir manteniendo una profundidad constante en el modelo mientras

Se incrementa en forma progresiva la velocidad.

Flujo alcanza una cierta velocidad media la erosión pinilla en las dos esquinas de la cara de

aguas arriba y el material erosionado se deposita a los lados conforme Se incrementa la

velocidad la erosión en las esquinas continuas alcanzar la misma profundidad a lo ancho de

la cara de aguas arriba al aumentar aún más la velocidad se inicia el transporte del fondo de

aguas arriba y se deposita en la zona de erosión y el tirante del flujo se mantiene constante y

Se incrementa paulatinamente la velocidad la erosión aumenta hasta que se alcanza la

estabilidad de socavación en material erosionado que al principio se deposita los lados es

desplazado aguas abajo debido al incremento de la velocidad es desplazamiento persiste hasta

que alcanza al otro extremo del Pilar desde donde el material es transportado por el flujo con
 25

las velocidades altas siempre habrá suspensión de material en este extremo, y en algunos

casos deposito a poca distancia aguas abajo del pilar.

Fórmula profundidad de erosión varía la estabilización del proceso es interferida por la

formación de rizos y donas cuando la cresta de una dona se aproxima en la zona de erosión

la profundidad Se incrementa y alcanzando un máximo cuando la cresta se encuentra en la

frontera de la zona de saturación este valor se mantiene constante hasta que media luna ha

sido desplazada aguas abajo decreciendo después hasta que pasa a la siguiente ronda los

experimentos realizados en el instituto de ingeniería mostraron que se requerían de 2 a 5

horas para obtener el equilibrio de socavación.

3.5.1.2Pilar rectangular esviajado con respecto a la dirección del flujo

Se inicia en la esquina exterior de la cara de aguas arriba y a medida que aumenta la velocidad

del flujo se hace uniforme frente a esa cara para velocidades altas con mucho transporte de

sedimentos la máxima socavación ocurre en la esquina de aguas

Abajo opuesta a aquella en la que se inició, es decir, en la esquina de aguas abajo expuesta

al flujo y no protegida por el cuerpo del pilar. En un pilar esviajado, las magnitudes de las

erosiones pueden ser dos o más veces mayores que para pilares alineados, por lo que es

conveniente evitar esta condición. Por otro lado, como la máxima erosión se presente aguas

abajo del pilar, se observa que cuando hay grandes avenidas, los pilares esviajados fallan,

se inclinan hacia aguas abajo; esto se interpreta erróneamente, al suponer que el empuje de

la corriente es lo que las hace fallar. (pág. 186)

 26

Figura 6. Proceso de socavación alrededor de pilares rectangulares

[Link] Pilares circulares y rectangulares con frentes redondeados

En estos casos la erosión se desarrolla en forma similar al de pilares rectangulares, pero con

ciertas peculiaridades: la erosión se inicia en dos regiones aproximadamente a 65° a cada

lado del eje del pilar (ver figuras). El resto del proceso es el mismo, excepto por las

diferencias registradas de las velocidades a las cuales los dos conos de erosión llegan a

juntarse, y a las que el mismo grado de erosión ocurre en el pilar y lleva a la máxima

erosión. Aguas debajo de pilares, las velocidades se incrementan en un grado mayor que en

las rectangulares, formando un vórtice con gran cantidad de material en suspensión.

Cuando no se puede garantizar que la dirección del flujo se mantendrá fija, conviene

utilizar pilares circulares. (pág. 186)

 27

Figura 7. Proceso de socavación alrededor de una pilar circular

 28

Figura 8. Diferentes etapas de los procesos de socavación en pilares.

3.5.2 Formulas de socavación local al pie de pilares de puente

[Link] Método de Laursen y Toch

Está basado principalmente en las experiencias realizadas por los investigadores en el

Instituto de Investigación Hidráulica de lowa, fue confinado con mediciones realizadas en

un puente sobre el rio Shunk por P. G. Hubbard del mismo laboratorio.

 29

Figura 9. Relación entre la erosión relativa y la profundidad relativa

Los autores de estos trabajos distinguen dos cosas generales; uno cuando la corriente incide

paralelamente al eje de los pilares y otro cuando forma un cierto ángulo con el mismo.

Cuando la mayor dimensión transversal del pilar está alineado con el flujo, la socavación

puede expresarse por:

𝑆0 = 𝐾1 𝐾2 𝑏

Donde:

𝑆0 = profundidad de la socavación, a partir del fondo

𝐾1 = coeficiente que depende de la relación tirante entre ancho del pilar y que se

encuentra en la gráfica de la figura.

𝐾2 = coeficiente que depende de la forma de la nariz del pilar y que se encuentra en la

tabla.

𝑏= ancho del pilar.

Como puede observarse, para Laursen y Toch la socavación depende únicamente del

tirante, ancho del pilar y de la forma de este, sin tomar en cuando la velación, ni el diámetro

del material del fondo. Este se considera únicamente arenoso, por lo que el método no es

aplicable si existen boleos en el cauce.

 30

En el caso de incidir oblicuamente la corriente y formar un ángulo con el eje del pilar, la

socavación puede determinarse con la expresión:

𝑆0 = 𝐾1 𝐾3 𝑏

Donde:

𝐾3 = coeficiente que depende del ángulo y de la relación a/b, el cual se determina con la

ayuda de la Fig.

En este caso la socavación no depende de la forma de la nariz del pilar. Laursen y Toch

realizaron sus observaciones fijándose en la socavación máxima que se puede presentar

para un tirante dado de la corriente. Observaron que sin variar el tirante y a pesar de

aumentar considerablemente la velocidad de la corriente, la socavación no progresaba.

Figura 10. Coeficiente de correlación cuando existe un ángulo de incidencia entre el eje

del pilar y la corriente.

 31

Figura 11. Coeficiente de correlación que depende de la forma del pilar. Aplicable a

pilares orientados según la corriente.

Según parece, la mayor socavación es lo que les preocupaba y no dan ningún criterio para

el caso en que no exista arrastre en el fondo o en que el arrastre sea menor con un fondo

con rizos o dunas pequeñas. El valor de esa máxima socavación obtenida no se ve afectado

por el diámetro del material del fondo, mientras se trate de arenas. Para gravas, no aclaran

si su gráfica puede usarse o no, pero en baleos definitivamente no es válida. (Apaclla, 2014,

p.190)

[Link] Método de Yaroslavtziev

Apaclla (2014) dice que “este investigador distingue dos casos, uno cuando el fondo del

cauce está formado por materiales no cohesivos y otro cuando está formado por materiales

cohesivos”. (p.190)

a. Método para suelos granulares

 32

Apaclla (2014) indica que “la expresión propuesta por Yaroslavtziev fue obtenida al

parecer de la observación directa en varios puentes de la Unión Soviética y es la siguiente”:

(p.190)

2
𝑆0 = 𝐾𝑓 𝐾𝑏 (𝑒 + 𝐾𝐻 ) − 30𝑑
𝑔

Donde:

𝑆0 = Profundidad de socavación, en m

𝐾𝑓 = Coeficiente que depende, en general, de la forma de la nariz del pilar y del ángulo de

incidencia entre la corriente y el eje de la misma. Se encuentra en la Fig.11.

Kv = Coeficiente definido por la expresión:

𝑣2
𝐾𝑣 = −0.28𝑠 √ 𝑔 − 30𝑑 = El cual se encuentra en la Fig. 12.

V= velocidad media de la corriente aguas arriba del pilar, después de producirse la erosión

general, en m/s.

g= 9.81m/s2

b1= Proyección de un plano perpendicular a la corriente, de la sección del pilar. Cuando el

ángulo de incidencia vale 0°, b1 es igual al ancho b del pilar.

e= Coeficiente de corrección, cuyo valor depende del sitio en donde están

colocados los pilares; vale 0,6 si se encuentran en el cauce principal y 1,0 para las

construidas en el cauce de avenidas.

KH= Coeficiente que toma en cuenta la profundidad de la corriente, definido por la

expresión:
𝐻
𝑙𝑜𝑔𝐾𝐻 = 0.17 − 0.35 𝑏 y que puede encontrarse además con ayuda de la curva Fig. 13.
1
 33

H= tirante de la corriente frente al pilar. Este valor es el obtenido al presentarse una

avenida.

d= diámetro en m de las partículas más gruesas que forman en el fondo y está representado

aproximadamente por el d85de la curva granulométrica.

Figura 11: Valores de Kf y b1

Fuente: Apaclla Nalvarte (2014). Hidráulica fluvial. Lima – Perú. Pág. 191.

El valor del ángulo de incidencia entre la corriente y los pilares es tomado en cuenta en el

valor de b1, tal y como puede observarse en la Fig. 11, el cual, a su vez, interviene rn la

valuación de KH y KV. Además, afecta directamente al valor del coeficiente K1, como se

observa en esas mismas figuras, excepto en los pilares rectangulares y circulares. (Apaclla,

2014, p.192)
 34

Yaroslavtziev hace hincapié en que, en vista de que el enviajamiento de la corriente influye

considerablemente en la erosión, puede resultar que para un caudal de agua menor, pero

que incida con el ángulo máximo, la erosión local llegue a ser mayor que para las

condiciones de gasto máximo con el ángulo menor. (Apaclla, 2014, p.192)

Yaroslavtziev advierte, además, que su fórmula puede conducir a errores en los casos en

que la relación H/b1 sea menor de 2 y el pilar este inclinado respecto a la corriente y añade

también que los valores con ella obtenidos en esas condiciones son menores que los que

realmente se presentan. (Apaclla, 2014, p.192)

Figura 11: Valores de Kf y b1

Fuente: Apaclla Nalvarte (2014). Hidráulica fluvial. Lima – Perú. Pág. 192.
 35

Figura 12. Cálculo de Kv (Método de Yaroslavtziev)

Fuente: Apaclla Nalvarte (2014). Hidráulica fluvial. Lima – Perú. Pág. 192.

Figura 13. Cálculo de KH (Método de Yaroslavtziev)

Fuente: Apaclla Nalvarte (2014). Hidráulica fluvial. Lima – Perú. Pág. 193.

b) Método de Yaroslavtziev para suelos cohesivos.

La expresión utilizada es la misma que para suelos granulares y permite dar un resultado

aproximado mediante la apreciación de la resistencia a la erosión del suelo cohesivo en

comparación con la resistencia a la erosión del suelo granular. Este es tomado en cuenta en

el segundo término (30d) en donde se considera un diámetro d, equivalente para los suelos

cohesivos tal como te muestra en la tabla 2. (Apaclla, 2014, p.192)

 36

Tabla 2.

Diámetro equivalente en suelos granulares

Peso volumétrico del Dimensiones del diámetro equivalente en suelos granulares

material seco en ton/m3 (cm)

Arcillas y suelos Suelos mediante Suelos de aluvión y

altamente plásticos plásticos arcillas margosas

<1.2 1 0.5 0.5

1.2-1.6 4 2 2

1.6-2.0 8 8 3

2.0-2.5 10 10 6

Fuente: Apaclla Nalvarte (2014). Hidráulica fluvial. Lima – Perú. Pág. 193.

[Link] Experimentos en el laboratorio de Poona, India

En la estación central de investigaciones de agua y fuerza de Poona, India, se realizaron una

serie de ensayos con el fin de encontrar la profundidad de socavación para un solo pilar

colocado en el centro de la corriente, en donde el material de fondo estaba formado por

arena uniforme de un diámetro medio de 0.29 mm al concluir el estudio se encontró la

siguiente expresión.
0.78
𝑆𝑇 𝑞 2/3
= 2.30 ( ) … … … … . (7.32)
𝐻 𝑏

Donde:

𝑆𝑇 = Nuevo valor del tirante después de la socavación.

𝑏 = Ancho del pilar, en metros.

𝑞 = Gasto unitario en m3/s/m.

 37

“El ángulo de incidencia es tomado en cuenta al considerar en lugar de b, la proyección del

pilar según un plano normal a la dirección del flujo, es decir, bl de Yaroslavtziev.”

(Apaclla, 2014, pág.194)

[Link] Método de la división de investigación de la Facultad de Ingeniería de la UNAM

Los estudios experimentales se llevaron a cabo en tres canales y en una instalación para

estudios aerodinámicos. En el primer canal se estudió la forma como se inicia y prosigue la

socavación, líneas de corriente, etc., y principalmente se realizó la comparación entre los

métodos de Yaroslavtziev y la Laursen – Toch. En el segundo canal se trató de verificar

para otras condiciones, las modificaciones propuestas para la determinación de la

socavación obtenida del primer modelo. En el tercero se estudió someramente la

socavación en gravas y protecciones formadas con pedraplenes y, por último, en las

instalaciones para estudios aerodinámicos se probaron distintos tipos de protección. El

trabajo fue desarrollado por José A. Maza Alvares como tesis de maestría, bajo la dirección

del profesor José L. Sánchez Bribiesca.

En las figuras 7.14, 7.15 y 7.16 se muestran las observaciones realizadas respecto a las

trayectorias de las partículas tanto liquidas como solidas alredor del obstáculo, forma de

iniciarse el proceso erosivo y los avances sucesivos de la socavación. (pág. 194)

 38

Figura 14. Diagrama para determinar la socavación alrededor de pilares rectangulares

Es importante aclarar que todos los valores de la socavación registrados con los que se

trabajó corresponden a los máximos observados, los cuales no siempre se presentan en el

mismo sitio del pilar. Esa ubicación depende de la velocidad de la corriente y de si el pilar

está o no esviajado. (pág. 195)

 39

Figura 15. Diagrama para determinar la socavación alrededor de pilares rectangulares con

frentes redondeados.

Así la magnitud de la socavación cuando la velocidad es pequeña y apenas se inicia la

erosión, es medida en la esquina de los pilares rectangulares y en las zonas laterales a los

65° respecto al eje en las circulares. Para velocidades mayores y siempre que el ángulo de

incidencia sea de cero grados, la magnitud máxima de la erosión es medida en el frente del

pilar cualquiera que sea su forma.(pág. 196)

 40

Figura 16. Diagrama para determinar la socavación alrededor de pilares circulares.

Los parámetros mostrados en las figuras tienen el siguiente significado:

d = Tirante de la corriente frente al pilar es una zona no afectada por la erosión local,

obtenida después de calcular la socavación general, transversal o en curvas. Es

decir, las erosiones locales se calculan a parir de las condiciones de la corriente, una

vez que se han producido las erosiones que afectan al fondo del cauce. (pág. 197)

𝑏1 = Proyección del pilar en un plano perpendicular a la corriente. Cuando el flujo y el

eje longitudinal del pilar están alineados b1 es igual al ancho del pilar b.

∅ = Ángulo entre la dirección del flujo y el eje longitudinal del pilar.

𝐹 = Numero de Froude del escurrimiento aguas arriba del pilar.

 41

𝑈 = Velocidad media del flujo inmediatamente aguas arriba de los pilares (no se

considera la presencia de los pilares), después de que el fondo ha sido erosionado,

sin tomar en cuenta la erosión local.

𝑆 = Erosión local medida desde el fondo no afectado por esta erosión.

𝑆𝑇 = Erosión local medida desde la superficie libre del agua.

𝑓𝑐 = Coeficiente por el que hay que afectar al número de Froude, 𝐹,cuando≠ 0° y 𝐹 >

0.25

La erosión local que se obtiene con ayuda de estas figuras es la máxima que puede

presentarse. Este método es útil solamente para fondos compuestos de arena y grava. Los

diagramas han sido desarrollados para pilares con secciones de forma rectangulares,

rectangulares con frentes redondeados, y circulares. Son también útiles cuando los ejes de

los pilares forman un cierto ángulo con la dirección del flujo. Para pilares esviajados el

número de Froude debe corregirse por el factor fc si este es mayor 0.06; el valor de fc se

indica en las figuras 7.14y 7.15 como una función de (ángulo que forma los ejes de los

pilares con la dirección del flujo)

[Link] Formula de Laursen

La fórmula de laursen contiene la profundidad de socavación como una variable implícita:

1.69
𝑏 𝐷𝑠 1𝐷
= 5.5 [( + 1) − 1] … … (7.33)
𝑑0 𝐷𝑜 11.5𝐷

Donde:

Ds = profundidad de socavación por debajo de la elevación del fondo del cauce y está

en forma implícita.

b = ancho del pilar.

Do = tirante medio del flujo aguas arriba del pilar.

 42

La profundidad de socavación DS está implícita.

[Link] Formula de Neil

Neil usando la ecuación de Laursen, diseño curvas para obtener la siguiente formula

explícita.

𝐷𝑠 𝐷𝑠
= 1.5( )0.3
𝑏 𝑏

[Link] Formula de Shen

Shen utilizo el número de Froude en la predicción de la profundidad de socavación.

𝐷𝑠 𝐷0
= 3.4(𝐹0 )2/3 ( )1/3
𝑏 𝑏
𝑈0
Donde 𝐹0 = 1/2 es el número de Froude basado en la velocidad medias aguas arriba y el
𝑔𝐷0

tirante 𝐷0 .

[Link] Universidad de colorado

La universidad de colorado desarrollo la siguiente formula reportada en el manual del

federal Highwas Administration, mostrada en la fig. N° 17

𝐷𝑠 𝑏
= 2.2( )0.65 (𝐹0 )0.43 …… (7.36)
𝐷0 𝐷0
 43

Figura. 17: formula del SCU para pilares circulares. La constante 2.0 cambia a 2.2 para

pilares rectangulares.

[Link] Formula de Jai y Fischer

Jain y Fischer (1979) estudiaron la socavación alreedor d pilares circulares para números de

Froude altos y propusieron las siguientes formulas.

𝐷𝑠 𝐷00.5
= 2.0(𝐹0 − 𝐹𝐶 )0.25 𝑝𝑎𝑟𝑎 𝐹0 − 𝐹𝑐 > 0.2 … … … (7.37)
𝑏 𝑏
𝐷𝑠 𝐷0
= 1.84( )0.3 𝐹 0.25 𝑝𝑎𝑟𝑎 𝐹0 > 𝐹𝐶 … … (7.38)
𝑏 𝑏

Donde Fc es el número de Froude crítico para el movimiento inminente del sedimento. En

la región 0<F0<Fc<0.2 se usa el mayor valor de ambas fórmulas.

El procedimiento para calcular Fc es el siguiente:

1. Determinar τc del diagrama de Shield basados en la estimación del diámetro del

material de fondo.

2. Calcular el espesor de la supcapa laminar ẟ y obtener la relación ẟ/d50.

3. Seleccionar el factor de corrección X en la distribución de velocidad logarítmica.

4. Calcule la velocidad media critica de la ecuación:

 44

𝜏𝑐 11.02𝐷0 𝑋
𝑈𝐶 = 2.5( )1/2 ln ( )
𝜌 𝑑50
𝑈
Calcular FC de: (𝑔𝐷 𝐶)1/2
0

3.5.3 Socavación en estribos y espigones

Aunque este tipo de erosión ha sido también estudiada por varios autores, entre quienes

pueden citarse a Laursen, Shen, Veiga Da Cunha, es la más difícil de determinar y estudiar.

El problema se da por la cantidad de parámetros que es necesario tomar en cuenta y porque

varían notablemente de un punto a otro, como es el caso del ancho total del estribo con

respecto al primer claro del puente, el ángulo de ataque de la corriente, el gasto teórico

interceptado por el estribo, sus taludes y la forma e inclinación de las esquinas cuando la

cara del estribo es vertical, además de considerar la velocidad y el tirante de la corriente y

las propiedades del material de fondo. El figura 7.18 se presenta un esquema de la

disposición de estribos. Para valuar esta erosión, los autores proponen las relaciones que a

continuación se indica y que fueron obtenidas a partir de datos de Artamonof y Veiga Da

Cunha, principalmente” (Apaclla, 2014, pág.200).

. La socavación al pie de un estribo, medida desde la superficie libre del agua, está dada

por:

𝑑𝑒 = 𝐾𝛼 𝐾𝑞 𝑑𝑎 … . . (7.39)

Donde:

de = tirante del flujo en el extremo del estribo considerando la erosión local.

da = tirante del flujo aguas arriba del estribo en una zona no afectada por la erosión del

estribo, pero que incluye la erosión general, transversal, en curvas y otro tipo erosión que

afecta el fondo.

Kα = coeficiente que depende del ángulo α; su valor se obtiene con la relación.

 45

α = ángulo medido aguas abajo del eje del puente y formado entre ese eje y la dirección del

flujo.

Kk = coeficiente que depende del talud, k, del extremo del estribo; se obtiene con ayuda de

la expresión:

𝐾𝛼 = 0.782𝑒 0.027𝛼

Si 0.001 ≤ k ≤ 1.3

0.85
Si 1.5 ≤ k ≤ 3 𝐾𝑘 = 𝐾0.0237

Kq = coeficiente que depende del cociente de Q1 , Q, donde Q1 es el gasto teorico que

podría pasar atravez del área ocupada por el estribo si este no existiera y Q es el gasto total

en el rio.

𝑘
𝐾𝑘 =
2.8𝑘 − 2.34

Su valor se obtiene de:

𝑄1
𝐾𝐾 = 4.429 + 1.063𝑙𝑛
𝑄

De la misma forma puede ser aplicada a espigones, siendo necesario saber si están

construidos en ambas márgenes y opuesto uno al otro, en este caspo, la socavación puede

ser reducida en un 755.

𝑑𝑒 = 0.75𝐾𝛼 𝐾𝐾 𝐾𝑞 𝑑𝑎

Aun cuando no existe forma de determinar la socavación cuando el espigón está cubierto de

agua, y ya que el método propuesto da los valores máximos posibles, es aconsejable

considerar Q1 como el flujo teórico máximo interceptado por el espigón hasta su cresta. El

error máximo entre los resultados de las ecuaciones propuestas y los datos de artamovov en

menor que 1.8%.” (Apaclla, 2014, pág.201)

 46

3.6.

Apaclla Nalvarte (2014) Quien indica

Este tipo de socavación es generada por las descargas de la presa casi libres de sedimentos;

esto ocasiona que las partículas que se erosionadas y transportadas de las primeras

secciones cercanas a la presa no sean remplazadas por otras que proceden de aguas arriba,

ello ocurre en el primer tramo. Este tramo inicia en la primera seccione inmediatamente

aguas debajo de la obra de toma o vertedor que no es afectada por la socavación local. La

erosión es mayor en las primeras secciones y disminuye paulatinamente a medida que se

alejan hacia aguas abajo, hasta una sección que se puede considerar no alterada o estable,

con respecto a este tipo de erosión.

Conforme las descargas continúan, el tramo erosionado se alarga debido a que las descargas

del rio remueven y transportan el sedimento del fondo. Cuando el material del fondo

contiene boleos tiende a acorazarse, es decir, quedan en la superficie las partículas de

mayor diámetro, las cuales protegen a las que se encuentran debajo.

Se observa también que aguas debajo de la sección no alterada, el material del fondo y la

pendiente varían con respecto a su condición original, como resultado de los cambios

morfológicos producidos por las modificaciones del gasto formativo: primero cambia el

tirante y después el ancho. La pendiente cambia lentamente durante un periodo de tiempo

muy grande. Al tramo después de la sección estable se le llama segundo tramo del río aguas

debajo de la presa (ver Fig. 19).

La erosión en el primer tramo se lleva a cabo lenta, pero continuamente durante la vida de

la presa, mientras se presenten descargas que puedan transportar material del fondo, de tal

forma que incrementa su longitud y desciende su fondo en forma gradual y continua (ver

Fig. 19).
 47

Figura 18. Esquema de la disposición de estribos.

Fuente 1
 48

Figura 19. Esquema del proceso de socavación aguas debajo de grandes presas.

Fuente 2

La erosión aguas debajo de grandes embalses depende principalmente de las descargas, su

variación y duración, de las características del material del cauce, y por consecuencia, del

transporte de sedimentos producido en el tramo no alterado o segundo tramo.

A pesar de su complejidad, el cálculo del proceso de socavación puede hacerse fácilmente

siempre y cuando el material en el cauce y en la profundidad afectada por la erosión sea el

mismo en todo el primer tramo.

Este proceso se vuelve más complicado, no solo cuando la erosión descubre materiales con

diferentes granulometrías, sino también cuando se produce acorazamiento del cauce. Si en

toda la profundidad del primer tramo el material es el mismo, pero con granulometría

extendida o bien graduada, 𝜎g>3, este puede llegar a acorazarse y debe considerarse una

condición especial, ya que solo descargas que la exceden pueden destruir el acorazamiento

y producir transporte de sedimentos. (p.203)

3.6.1. Método de solución

El método de cálculo propuesto por Maza ofrece varias ventajas:

 Tanto el proceso erosivo y los cambios morfológicos del canal están considerados

en el mismo tramo; es decir, el ancho de la sección varia simultáneamente con la

pendiente, la profundidad y el transporte de sedimentos.

 No se presentan inestabilidades de los métodos numéricos de solución.

 49

 Es más fácil de utilizar por los ingenieros proyectistas que la solución simultánea de

las ecuaciones diferenciales de continuidad, energía y continuidad del transporte de

sedimentos.

 Para la cuantificación de la erosión aguas abajo de la presa se utiliza la ecuación de

Maza y Cruickshank propuesta en 1973 para el estudio de estabilidad de canales.

(p.204)

Además, considera la estabilidad dinámica del tramo en estudio, para lo cual se utilizan las

fórmulas de iseño para cauces estables deducidas de las teorías de Cruickshank y Maza,

Meyer, Peter y Muller y Gluschkov:

2
𝑏 (2.583𝑚)+1.526 (𝑄𝐵 3 + 0.047(𝜀𝑏)2/3 = 𝐸 2/3 𝐾 2.588 (𝑄/𝛼)2.193 …………… (7.41)

𝑑 = 𝑏 𝑚 /𝐾 …………………………………………… (7.42)

𝑄 (2+3𝑚/3𝑤)
𝑄𝐵 2/3 = 𝐸 2/3 𝐾 0.0893/𝑊 (𝛼 ) 𝑆 (1.178𝑚+0..696)/𝑤 −

𝑄 (2+3𝑚/3𝑤)
0.047𝜀 2/3 (𝛼 ) 𝐾 1.089/𝑊 𝑆 −0.304𝑤 …………………………………. (7.42)

Donde:

b = ancho de la superficie libre del agua en m.

d = tirante medio. Es la relación entre el área de la sección transversal y el ancho de la

superficie libre del agua.

S = pendiente hidráulica.

K= coeficiente de forma cuyo valor depende de la facilidad de erosión del material del

cauce, para canales aluviales su valor medio es de 10.

 50

m = exponente función de la pendiente del tramo y las características del material. Varia de

0,5 a 1, su valor medio para cauces aluviales es de 0.7.

Q = gasto formativo o dominante.

QB = gasto solido de fondo.

W= 1+1,634 m

1/2
𝜀 = 8(gΔ𝐷𝑚 3 ) ………………………………………………… (7.43)

n᾿ 1.5 1
𝑁𝑛 = ( 𝑛 ) …………………………………………………… (7.46)
Δ𝐷𝑚

𝐸 = 𝜀𝑁𝑛 1.5 …………………………………………………………. (7.4)

0.04
𝛽 = Δ2 𝑔0.5𝐷 ………………………………………………………. (7.4)
𝑚

7.85𝜔50
⍺=𝐷 0.634 Δ0.456 ……………………………………………………. (7.4)
𝑚

n = rugosidad total de la sección. Se obtiene de la fórmula de Manning.

n᾿ = rugosidad debida a las partículas. Se determina con las expresiones de Meyer y Peter.

ω = velocidad de caída de las partículas de diámetro D50

Δ = densidad relativa de las partículas sumergidas.

Dm= diámetro medio del material del fondo en metros.

De acuerdo con Stebbings y Friedkin, quienes realizaron estudios en canales de laboratorio

observando los cambios de ancho y profundidad en secciones con diferente transporte de

sedimentos, encontraron que el exponente m es en todos los casos es igual a 0.7, mientras

que el coeficiente k varia de 8,22 a 25,82 para el mínimo y máximo transporte. Así, cuando

se estudia este tipo de socavación y no se tienen datos para calibrar las ecuaciones, se

recomiendan los siguientes valores:

 51

K= 10 para la sección inmediatamente aguas debajo de la presa donde el transporte

sedimentos es casi nulo.

K= 18 para la última sección de agua abajo del tramo sujeto a erosión, debido a que no se

conoce con certeza que el transporte en él sea el máximo.

𝑥
𝐾𝑖 = 𝐾𝑐 + (𝐾𝑟 − 𝐾𝑐 )√1 − ( 𝐿𝑖 )2 ………………………………………… (7.4)

El sistema de ejes y numeración se muestra en el perfil de la Fig.7.20. El significado de las

variables es:

Kc= coeficiente de la primera sección aguas arriba del primer tramo, kc ≈10.

Kr= coeficiente de la última sección aguas abajo del primer tramo, kc ≈18.

L= longitud del primer tramo o distancia entre las secciones arriba mencionadas la

cual se incrementa con el tiempo.

Xi = distancia entre la sección final y aquella para la cual k, será obtenido.

Es importante mencionar los motivos que dificultan la solución exacta del problema, la cual

es descrita por Apaclla Nalvarte (2014)

1. No existen mediciones completas y sistemáticas tomadas durante varios años para

comparar los resultados analíticos con los reales con el fin calibrar el modelo.

2. La cantidad de variables influye en la cuantificación de la erosión y cambios

morfológicos, aun cuando algunas deben ser obtenidas de datos reales, si estos no

existen, su valor analítico es aceptado sin ajuste. Este es el caso del exponente m y el

coeficiente K.
 52

3. Para cuantificar la socavación aguas debajo de la presa se considera la existencia de un

mismo material a todo lo ancho y en toda la profundidad erosionada de la sección del

cauce, lo cual no es necesariamente cierto.

4. Cuando se presenta poco transparente de sedimentos se tiene a formar un acorazamiento

del cauce, tendencia que es diferente proporcional a la desviación estándar de la

granulometría. Este tipo de fenómenos no han sido estudiados con detalle en la

naturaleza, pero se sabe de su existencia, como el que se produce durante el proceso de

acorazamiento del centro hacia las márgenes hasta alcanzar la superficie del agua.

Este método permite estimar, para el primer tramo, las alteraciones producidas por la

presencia de la presa, tales como la erosión en el cauce, los cambios en sus secciones, la

modificación de la pendiente y el tiempo en el cual ocurren. (p.206)

3.6.2. Cálculos preliminares

Los cálculos se realizan tomando en cuenta las hipótesis mostradas anteriormente en las

figuras 7.19 y 7.20 el cual Apaclla Nalvarte menciona:

1. Es necesario conocer el gasto dominante o formativo, el perfil y las secciones de

interés del rio agua abajo de la futura presa hasta una longitud 50km

aproximadamente, de esta forma K, m y α para este tramo se podrán determinar sin

suponerlos. Si es posible calibrar estos parámetros ajustándolos a observaciones

reales, el cálculo de la socavación será más preciso.

2. Al final del primer tramo, el canal no es afectado por esta socavación. Aguas abajo

de este punto, la pendiente hidráulica permanece igual a la original Sr. Entonces,

S=Sr.
 53

3. Se debe determinar el nuevo gasto formativo de la presa. Cuando se trata de una

planta hidroeléctrica puede ser considerado el gasto medio máximo diario. Si el

vertedor presenta descargas durante largos periodos de tiempo, el gasto asociado

con esas descargas puede ser usado como el gasto formativo.

4. El ancho Br y la profundidad r, y el nuevo gasto formativo Q en el segundo tramo se

obtiene con Sr. Estos parámetros corresponden a la última sección del primer tramo.

A continuación, se describe el procedimiento a seguir:

[Link] calcula QB con la ecuación (7.43) haciendo –S=Sr.

[Link] QB se obtiene Br, con la ayuda de la ecuación 7.41

[Link] Br, se obtiene dr con la ayuda de la ecuación 7.42

Donde:

El subíndice r es el término asociado con la selección final del primer tramo (inicio del

tramo de rio sin erosión y bajo estabilidad dinámica)

El subíndice es el término asociado con la selección inicial del primer tramo (condición

crítica de erosión).

5. Para la primera sección del primer tramo, el transporte de sedimentos se considera

nulo, lo que significa que la pendiente hidráulica S=Sc. Esto se calcula con la

ecuación 7.43, con QB=0, la cual toma la siguiente forma:

𝑊
0.047 1/𝑊 𝛼 1.178𝑚+1
𝑆𝑐 = [ 𝐾 (𝑄 ) ]……………………………………(7.45)
𝑁

6. Para hacer el cálculo, el primer tramo se divide en N segmentos iguales. En cada

segmento se aplican las ecuaciones, 7.41 a 7.44. Para cada tramo existen cuatro

cantidades desconocidas: b, S, d y QB, y solo se tienen tres ecuaciones, por lo que

deberán ser aceptadas alguna de las hipótesis planteadas anteriormente.

 54

7. Si en el primer tramo del rio las márgenes están formadas de roca o son muy

resistentes, se pierde el grado de libertad, ya que b se considera constante durante

todo el proceso de erosión, entonces solo se requieren dos ecuaciones (7.41 y 7.43)

para resolver el problema:

2 1.387
𝛼𝛽 0.849 1 𝑄 3
𝑆 =(𝑄) [( 𝜀𝐵𝐵 ) + 0.047] ………………… (7.46)
𝑁 1.387

𝑄 0.849 𝑁 0.387
𝑑 = (𝛼𝐵) [(𝑄 2/3 +0.047
] …………………………………. (7.47)
𝐵 /𝜀𝐵)

Se tienen dos ecuaciones con tres cantidades desconocidas: QB, d y S, por lo que el

procedimiento a seguir es similar al explicado para tres grados de libertad. (p.206,207)

3.6.3. Hipótesis

Cuando existen más variables desconocidas que ecuaciones, se pueden establecer las

siguientes hipótesis:

a) El transporte de sedimentos varía de acuerdo con una ley preestablecida desde cero

en la primera sección hasta QBr en la última sección del primer tramo. Estas

secciones se encuentran a una distancia Lr medida a lo largo del cauce, la cual

crece con el tiempo. Como una primera estimación, se puede aplicar una variación

lineal, pero cualquier otra puede utilizarse siempre y cuando se cumpla el principio

de continuidad del sedimento para cada tramo. Si se conoce la variación real, esta

debe ser aplicada.

b) El coeficiente K también depende del transporte de sedimento. Para material de

fondo se han obtenido en forma experimental los valores de K=8.22 para 𝑄𝐵 = 0 y

𝐾𝑚𝑎𝑥 = 22.82 para 𝑄𝐵 máximo. Considerando lo anterior, K puede hacerse variar

desde 𝐾𝑐 =10 en la primera sección hasta K=Kr, en la última. Si no se tiene

 55

información para determinar Kc, es recomendable iniciar con Kr=18. Para iniciar el

cálculo se recomienda utilizar la ecuación 7.44

Figura 20. Consideraciones y variables en el cálculo

Fuente. Rocha Felices, Hidráulica Fluvial, p.208, 1998(Maza Álvarez,

México, 1987)

Figura 21. Volumen de control de la erosión entre las secciones i, i+1 y los instantes j, j+1
 56

Fuente. Rocha Felices, Hidráulica Fluvial, p.208, 1998 (Maza Álvarez, México, 1987)

Bibliografía

Aguirre, J. (2006). Hidraulica fluvial. Venezuela: Universidad los Andes.

Apaclla, R. (2014). Hidraulica Fluvial. lima: Universidad Nacional Agraria la Molina, lima

– Perú.

Chirino, F. (2017). Estudio de socavación de la subestructura de puente. Lima - Peru:

Universidad Cesar Vallejo.

Enciso, P. (2015). Estimacion de S ocavacion en el Puente. Ayacucho: Universidad

Nacional San Cristobal de Huamanga.

Martinez, G. (2007). Estudio de socavacion en Obras Civiles y su problematica. Mexico:

Instituto Politécnico Nacional.

Pareja, P. (2010). Diseño de Pequeños Puentes. lima : Universidad Nacional dde

Ingenieria.

Rocha, A. (1998). Hidraulica Fluvial. Lima: Universidad Nacional de Ingenieria.

Rodriguez, H. (2010). Socavacion en cauces. Colombia: Escuela Colombiana de Ingeniería.

Trujillo, H. (2015). Tesis para Optar el titulo de Ingeniero Civil. Trujillo: Universidad

Privada Antenor Orrego .

Venegas, Z. (2011). Socavacion al pie de Estribos de Puentes. Lima - Peru: Universidad

Nacional de Ingenieria Lima - Peru.