Saltar al contenido

Inicio de arrastre y transporte de fondo: Shields y Meyer-Peter-Müller aplicados

Pablo Rojas4 min de lectura

En este artículo
  1. Primera pregunta: la tensión de corte crítica
  2. Segunda pregunta: cuánto transporta
  3. Acorazamiento: por qué el lecho se defiende solo
Lecho de río de gravas y bolones bajo agua clara y poco profunda, con el flujo visible sobre las piedras

Dos preguntas gobiernan casi todo el diseño fluvial: ¿se mueve el lecho? y, si se mueve, ¿cuánto material transporta? De la primera dependen la estabilidad de un canal no revestido, el acorazamiento bajo un puente y la vida útil de una protección. De la segunda, la sedimentación de un embalse, la degradación aguas abajo de una presa y el equilibrio de un cauce intervenido.

Primera pregunta: la tensión de corte crítica

El flujo ejerce sobre el lecho una tensión de corte media:

τ0=γw R S\tau_0 = \gamma_w\, R\, S

con RR el radio hidráulico y SS la pendiente de energía. En cauces anchos, R≈yR \approx y y basta el tirante. El grano empieza a moverse cuando τ0\tau_0 supera la tensión crítica τc\tau_c, que se expresa adimensionalmente mediante el parámetro de Shields:

θ=τ0(γs−γw) D\theta = \frac{\tau_0}{(\gamma_s - \gamma_w)\, D}

El diagrama de Shields relaciona el valor crítico θc\theta_c con el número de Reynolds del grano, Re∗=u∗D/νRe_* = u_* D/\nu, donde u∗=τ0/ρu_* = \sqrt{\tau_0/\rho} es la velocidad de corte. Para el rango de interés práctico —arenas gruesas, gravas y bolones, con Re∗>400Re_* \gt 400— el diagrama se aplana y se puede usar directamente:

θc≈0,047 (valor de Meyer-Peter-Mu¨ller)oθc≈0,06 (Shields original)\theta_c \approx 0{,}047 \ \text{(valor de Meyer-Peter-Müller)} \qquad \text{o} \qquad \theta_c \approx 0{,}06 \ \text{(Shields original)}

La diferencia entre 0,047 y 0,06 no es menor: define un umbral de movimiento un 28 % distinto. En proyecto conviene usar 0,047 cuando se busca ser conservador respecto del movimiento (por ejemplo, verificar que un canal no se erosione) y 0,06 cuando se busca ser conservador respecto de la estabilidad de una protección.

Despejando, la tensión crítica y el tamaño estable resultan:

τc=θc (γs−γw) D50⟹D50,estable=γw R Sθc (γs−γw)\tau_c = \theta_c\,(\gamma_s - \gamma_w)\, D_{50} \qquad \Longrightarrow \qquad D_{50,\text{estable}} = \frac{\gamma_w\, R\, S}{\theta_c\,(\gamma_s - \gamma_w)}

Ejemplo 1 — ¿se mueve?

Cauce de grava con D50=40D_{50} = 40 mm, tirante y=1,20y = 1{,}20 m, pendiente S=0,004S = 0{,}004, γs=26\gamma_s = 26 kN/m³.

  • τ0=9,81×1,20×0,004=0,0471\tau_0 = 9{,}81 \times 1{,}20 \times 0{,}004 = 0{,}0471 kN/m² =47,1= 47{,}1 N/m²
  • τc=0,047×(26−9,81)×0,040=0,047×16,19×0,040=0,0304\tau_c = 0{,}047 \times (26 - 9{,}81) \times 0{,}040 = 0{,}047 \times 16{,}19 \times 0{,}040 = 0{,}0304 kN/m² =30,4= 30{,}4 N/m²

Como τ0>τc\tau_0 \gt \tau_c, el lecho está en movimiento. El parámetro de Shields real vale θ=47,1/(16190×0,040)=0,0727\theta = 47{,}1/(16190 \times 0{,}040) = 0{,}0727, es decir, un 55 % sobre el umbral.

Segunda pregunta: cuánto transporta

La ecuación de Meyer-Peter-Müller (1948) sigue siendo la referencia para transporte de fondo en lechos de grava. En forma adimensional:

qb∗=8 (θ−θc)3/2,qb∗=qb(γsγw−1)g D3q_b^* = 8\,(\theta - \theta_c)^{3/2}, \qquad q_b^* = \frac{q_b}{\sqrt{\left(\dfrac{\gamma_s}{\gamma_w} - 1\right) g\, D^3}}

donde qbq_b es el caudal sólido de fondo por unidad de ancho, en volumen de sólidos por segundo y por metro. Tres observaciones que cambian cómo se interpreta el resultado:

  • El exponente 3/2 sobre el exceso de tensión hace el transporte brutalmente no lineal. Duplicar el exceso de tensión casi triplica el caudal sólido. Por eso el balance de sedimentos anual de un río lo definen unos pocos días de crecida, no el régimen medio.
  • Cerca del umbral el resultado es muy incierto. Con θ\theta apenas sobre θc\theta_c, un error del 10 % en θc\theta_c produce errores de cientos por ciento en qbq_b. Ahí la ecuación entrega orden de magnitud, no cifras.
  • MPM da transporte de fondo, no total. En cauces arenosos la carga en suspensión suele ser 10 a 100 veces mayor. Para eso hay que ir a formulaciones de carga total (Engelund-Hansen, Ackers-White, van Rijn).

Ejemplo 2 — cuánto

Con los datos anteriores (θ=0,0727\theta = 0{,}0727, θc=0,047\theta_c = 0{,}047, D=0,040D = 0{,}040 m):

qb∗=8 (0,0727−0,047)3/2=8×(0,0257)3/2=8×0,00412=0,0330q_b^* = 8\,(0{,}0727 - 0{,}047)^{3/2} = 8 \times (0{,}0257)^{3/2} = 8 \times 0{,}00412 = 0{,}0330

El factor de escala: (2,65−1)×9,81×0,0403=1,65×9,81×6,4×10−5=1,036×10−3=0,0322\sqrt{(2{,}65 - 1) \times 9{,}81 \times 0{,}040^3} = \sqrt{1{,}65 \times 9{,}81 \times 6{,}4\times10^{-5}} = \sqrt{1{,}036\times10^{-3}} = 0{,}0322 m²/s.

qb=0,0330×0,0322=1,06×10−3 m3/s por metro de anchoq_b = 0{,}0330 \times 0{,}0322 = 1{,}06\times10^{-3}\ \text{m}^3/\text{s por metro de ancho}

En un cauce de 25 m de ancho: 0,02660{,}0266 m³/s de sólidos. Con porosidad de 0,35, eso equivale a unos 0,041 m³/s de volumen de lecho, o 3.500 m³ por día si la crecida se mantiene. Un embalse aguas abajo recibiría ese material completo.

Acorazamiento: por qué el lecho se defiende solo

En un lecho de granulometría extendida, el flujo se lleva primero las partículas finas y deja en superficie las gruesas. Se forma una coraza de una o dos capas de espesor cuyo D50D_{50} puede ser dos o tres veces el del material subyacente. Una vez formada, el lecho resiste caudales bastante mayores que los que lo erosionaron inicialmente.

Dos consecuencias prácticas:

  • Un cauce natural aparentemente estable puede estarlo solo por la coraza. Si una crecida excepcional la rompe, la erosión posterior es acelerada porque debajo hay material mucho más fino.
  • Aguas abajo de una presa llega agua limpia ("agua hambrienta") que erosiona sin reponer. El lecho se degrada hasta acorazarse, y en el camino baja la cota del cauce, descalza fundaciones de puentes y desconecta las tomas de riego. Es un efecto que puede propagarse decenas de kilómetros y durar décadas.

Relacionado: pérdida de capacidad de embalses por sedimentación, sedimento en suspensión y la calculadora del número de Froude.