CAPITULO III : FENOMENO DE ACORAZAMIENTO

3.1 ACORAZAMIENTO DE UN CAUCE Los sedimentos de granulometra extendida o suelos bien graduados, g > 2, son los que revelan en su curva granulomtrica un amplio intervalo de tamaos de partculas y cantidades apreciables de cada tamao. La caracterstica ms importante de los sedimentos no cohesivos de granulometra extendida, cuando estn sujetos a la accin de un flujo de agua, es su capacidad para desarrollar una coraza o armadura en la superficie del lecho, como un empedrado. Para que ocurra esto, es necesario que el flujo pueda arrastrar los granos finos del material del cauce, pero no los ms gruesos: las partculas pequeas, incapaces de resistir el paso del agua, van siendo lavadas o transportadas por la corriente, lo que hace que afloren partculas de mayor tamao que s pueden resistir a las fuerzas de arrastre y sustentacin debidas al flujo; de este modo, por la permanencia de las partculas gruesas, se va formando una especie de empedrado o coraza en la superficie del lecho, que protege del arrastre o acarreo al material fino subyacente. El proceso de acorazamiento es fundamentalmente un problema de inicio de arrastre o de transporte incipiente; dicho proceso debe tenerse en cuenta al tratar con sedimentos no cohesivos de granulometra no uniforme y extendida, ya que en este tipo de sedimentos no se tiene una condicin nica o situacin promedio de transporte incipiente, como acontece con material uniforme, sino que pueden darse varias condiciones de transporte incipiente: Si las partculas de dimetro mnimo estn a punto de ser movidas por el flujo, se est en la condicin llamada punto inferior de transporte incipiente; por el contrario, cuando el flujo est a punto de mover las partculas de dimetro mximo, se alcanza entonces el punto superior de transporte incipiente. Entre ambos puntos o situaciones extremas de transporte incipiente, dependiendo de las caractersticas del flujo y de las del material del cauce, tambin podrn suceder otras situaciones intermedias de transporte incipiente. Por tanto, entre los puntos inferior y superior de transporte Incipiente es donde ocurre el fenmeno de acorazamiento y, por consiguiente, podrn darse distintos grados de acorazamiento, segn los tamaos de las partculas que permanezcan formando la coraza o armadura del cauce. Sin embargo, dado el carcter aleatorio de la turbulencia en el flujo de agua, siempre existir la posibilidad de que cualquier grano sea arrastrado por el flujo. Bajo igualdad de condiciones, la probabilidad de que una partcula no sea arrastrada y forme parte de la coraza, q, ser mayor para las partculas gruesas que para las finas. Gessler obtuvo experimentalmente dicha probabilidad, analizando el material erosionado y el que formaba coraza; encontr que la variable aleatoria c/0 segua una distribucin gaussiana o normal de probabilidades, con media = 1 y desviacin estndar = 0.57, fig 3.1. Por tanto, si se estandariza la variable aleatoria, escribirse Z = c/0 - 1 0.57 (3.1)

Al igual que la fig 3.1, la ec 3.1 permite tambin valuar, con ayuda de la tabla 3.1, la probabilidad q que tienen las partculas de un tamao especfico de no ser arrastradas por la corriente.
0 .9 9 P ro b a b il id a d , q

0 .9 5 0 .9 0 0 .8 0

0 .6 0 0 .4 0

0 .2 0 0 .1 0 0 .0 5

0 .0 1 0 0 .5 1 .0 1 .5

c /

2 .0

Fig 3.1 Probabilidad de que una particula o piedra determinada no sea arrastrada por el flujo y forme parte de la coraza, segn Gessler

Al tratar con un cauce constituido por material no cohesivo y de granulometra extendida, uno de los problemas por resolver es el clculo del esfuerzo cortante mximo o esfuerzo cortante crtico que podra llegar a soportar dicho cauce. La coraza es funcin de la granulometra del material original, de las caractersticas del flujo y de las propiedades del fluido. Si el flujo no es muy intenso, ser capaz de arrastrar las partculas ms pequeas. Si esto ocurre, el dimetro medio de la granulometra del material permanece formando coraza, aunque un poco mayor, ser parecido a la del material original. A medida que la intensidad del flujo aumenta, la corriente arrastra partculas ms y ms grandes, y el dimetro medio de la coraza (Dmc) ir tambin aumentando hasta alcanzar un valor mximo. S el flujo an tiene capacidad para mover las partculas de dimetro mximo, stas a su vez sern arrastradas y dejarn al descubierto a las partculas finas que estaban protegiendo. Como resultado de la falta de proteccin, el dimetro medio de la coraza empezar a disminuir. Por tanto, el esfuerzo cortante mximo que puede resistir un material de granulometra no uniforme y extendida es aquel que corresponde a la condicin de flujo que deja en la superficie del lecho o coraza una distribucin granulomtrica tal cuyo dimetro medio es el mximo. Aunque durante el proceso de acorazamiento hay transporte de sedimentos, el acarreo de partculas no causa prcticamente ninguna erosin o abatimiento del fondo ni modifica la pendiente del cauce, se dice por ello que un cauce acorazado es estable. Sin embargo, si las caractersticas del flujo son tales que rebasan el esfuerzo cortante mximo que puede resistir el material que constituye el cauce, ste deja de ser estable, pues todo el material de la superficie del cauce estar en movimiento. Si esto acontece y

no hay continuidad en el acarreo de sedimentos, habr erosin continua del fondo hasta que el esfuerzo cortante del flujo se reduzca por disminucin de la pendiente o por aumento del rea hidrulica, y se iguale as con el que puede resistir el material del cauce; por el contrario, si hay continuidad en el transporte de sedimentos o gasto slido, no acontece ni erosin ni depsito, y slo se tendr un transporte continuo de sedimentos hasta que el tirante disminuya por reduccin del gasto lquido. Gessler desarroll un mtodo para estimar el esfuerzo cortante crtico o esfuerzo cortante mximo que puede resistir un material no cohesivo de granulometra no uniforme y extendida. 3.3.1 Mtodo de Gessler Gessler define como esfuerzo cortante crtico de un material no cohesivo de granulometra extendida, al esfuerzo que corresponde a la condicin de flujo que hace que el dimetro medio de la coraza del cauce alcance su valor mximo. El mtodo consiste, por consiguiente, en suponer diferentes condiciones de flujo o esfuerzos cortantes, 0, y en calcular, para cada esfuerzo cortante supuesto, el dimetro medio de la coraza o armadura, Dmc , que cada 0 deja en la superficie del lecho. El esfuerzo 0 con el cual se logra el mximo Dmc corresponde a la condicin crtica de arrastre, o sea que para dicha condicin de flujo 0 = c. El algoritmo de Gessler se presenta y detalla a continuacin en el que se explica el significado de las variables y parametros, as como la secuencia y la forma en que deben valuar. Po : distribucin de frecuencias relativas acumuladas (por ciento que pasa) del material original o caracterstico del cauce.

El eje de las ordenadas (porcentaje que pasa) de la curva granulomtrica caracterstica del material del cauce se divide en pequeos intervalos, Po. El tamao de estos se escoge de manera que la curva granulomtrica quede bien representada, por lo que la amplitud de los intervalos Po puede ser variable o constante. Sin embargo se recomienda dividirla en diez intervalos iguales, o sea en deciles: de 0 a 0.10, de 0.10 a 0.20,., de 0.90 a 1.00; de esta manera el tamao de los intervalos es constante e igual a 0.10, es decir Po = 0.10 po : Dn : marca de clase o punto medio del intervalo. dimetro o tamao correspondiente a la marca de clase del intervalo, en mm.

Si la granulometra caracterstica del cauce no se ajusta a ninguna distribucin terica, o no se cuenta con datos del cribado de las muestras, en tonces los dimetros caractersticos Dn se determinan directamente de la curva granulomtrica, leyendo lo ms preciso posible. Por el contrario, si la granulometra caracterstica sigue un modelo matemtico o distribucin terica, los dimetros Dn se generan con la ecuacin del modelo. En esta explicacin se asume que se conocen los dimetros tocantes a la marca de clase de cada intervaloi en que se dividi la curva granulomtrica. c : esfuerzo cortante crtico de las partculas, en N/m2

La resistencia del material del cauce se debe al efecto estabilizador de cada grano individual que permanece formando coraza. El esfuerzo cortante mximo que puede soportar la coraza antes de ser destruida o arrastrada por el flujo, se define por la suma de las contribuciones de los granos individuales que la constituyen. La contribucin de cada grano o esfuerzo cortante crtico de cada partcula se calcula empleando la ecuacin de Meyer-Peter y Mller c = 0.047 (s - ) Dn El calculo se efecta para cada dimetro representativo de cada intervalo, Dn, es decir, sustituyendo los dimetros en la ecuacin. o : esfuerzo cortante medio producido por el flujo, en N/m2.

Este esfuerzo es el que esta actuando sobre las partculas que constituyen el lecho; su valor o intensidad se supone y se mantiene constante hasta que se inicie una nueva aproximacin o iteracin. Si no se tiene idea con qu valor de o empezar el clculo, se recomienda, por sencillez, estimarlo tambin con el criterio de Meyer-Peter y Mller, pero utilizando el dimetro medio de la curva granulomtrica o el D65. c /o : relacin de esfuerzos o variable aleatoria. Z: variable aleatoria normalizada.

Como ya se indic, la fig 3.1 permite estimar la probabilidad q. Sin embargo, el clculo de dicha probabilidad puede hacerse con mayor precisin o analticamente, para lo cual hay que estandarizar o normalizar la variable aleatoria c /o. Q: probabilidad de que no sean arrastradas las partculas de tamao Dn.

Cada probabilidad q se estima con la ayuda de la tabla 1.4, segn el rea bajo la curva normal estndar que define cada valor referido de Z. Cuando el valor de Z no coincida con los que se dan en la tabla 1.4, la probabilidad q se calcula interpolando linealmente entre los valores ms prximos. Sin embargo, en lugar de utilizar la tabla 1.4 para valuar q, se sugiere emplear las ecuaciones: Si Z = 0 q = 1 1/2 exp{-1/2 Z2}( c1 t + c2 t2 + c3 t3 + c4 t4 + c5 t5) Si Z < 0 q = 1/2 exp{-1/2 Z2}( c1 t + c2 t2 + c3 t3 + c4 t4 + c5 t5) siendo t= 1 1 + 0.2316419Z (3.4) (3.3) (3.2)

c1 = 0.31938153 c2 = -0.356563782 c3 = 1.78147937 c4 = -1.821255978 c5 = 1.330274429 q Po :frecuencia del material original que permanece formando coraza Primero se multiplica la amplitud Po de cada intervalo por la probabilidad q vinculada al dimetro representativo del intervalo en consideracin. Una vez efectuados los productos, se suma y se anota el resultado. Pa : frecuencia relativa del material de la armaduara o porcentaje de partculas de tamao Dn que permanece formando coraza. Cada uno de estos porcentajes se valora dividiendo cada valor de q Po, entre la suma total de la misma Pa : distribucin de frecuencias relativas acumuladas (porcentaje acumulado que pasa) del material que permanece formando coraza.

Los porcentajes se calculan de la siguiente manera: en el primer rengln de la columna Pa se anota el cero y se suma al primer valor de la columna Pa, y el resultado se apunta en el tercer rengln de la columna Pa; luego este nmero se adiciona a la segunda cantidad que aparece en la columna Pa, y la suma se anota en el quinto rengln de la columna Pa, y as sucesivamente. Las cantidades en los renglones pares de la columna Pa son las marcas de clase de los intervalos y dibujando estos valores o porcentajes contra los respectivos diametros Dn, se determina la curva granulomtrica de la coraza. DnPa : porcentaje de partculas de tamao Dn que constituye la coraza o armadura. Las cantidades que forman esta columna se obtienen multiplicando los dimetros Dn por los correspondientes valores Pa. La suma de todos estos productos da el tamao del dimetro medio de la distribucin granulomtrica del material de la coraza. Dmc=

D ( Pa)dPa
i 0

(3.5)

Con el computo de Dmc se termina lo que sera una aproximacin; las siguientes se inician suponiendo otros esfuezos cortantes, o, se recomienda dibujar dichos valores contra los respectivos Dmc, ya que as se puede visualizar para qu intervalo de valores de o aumenta o disminuye el Dmc, fig 3.2. Una vez definida la condicin crtica de arrastre, Gessler recomienda para el diseo de canales sin revestimiento que el valor de c as calculado sea disminuido, dividindolo entre un factor de seguridad f, y el valor que l recomendo fue f = 1.3, esto es

cdiseo =

calculada f

(3.6)

Fig 3.2 Determinacin de la condicin crtica de arrastre

(1 q) : probabilidad que tienen las partculas Dn de ser arrastradas o de no formar parte de la coraza. (1 q) Po : frecuencia del material original que es erosionado o arrastrado. Pe : frecuencia relativa del material erosionado o porcentaje de partculas de tamao Dn que es arrastrado o lavado por el flujo. Pe : frecuencia relativa acumulada del material erosionado.

En la fig 3.3 se muestra la curva granulomtrica del material original o caracterstico del cauce, la del material que forma coraza y la del material erosionado. Si bien es cierto que el mtodo de Gessler es muy laborioso, es el ms racional que existe y es relativamente fcil de programar. Adems, cuando la granulometra del cauce sigue una distribucin terica, el mtodo de Gessler se simplifica mucho, como se explica a continuacin.

100% 90% 80% 70% 60% 50% 40% 30% 20% 10% 0% 100.0 CORAZA 10.0 Arrastre Dimetro mm Original 1.0 % que pasa

Fig 3.3 Curva granulomtrica del material original del cauce, coraza y del material original

3.1.2 Mtodo de Cruickshank-Garca Cruickshank y Garca Flores demostraron la forma de adimensionar el mtodo de Gessler y analizaron dos distribuciones de tamaos tericas, la log-normal (ec 1.34, cap 1) y la logartmica (ec 1.44, cap 1). Los resultados de su anlisis los presentaron en una serie de diagramas con parmetros adimensionales, y propusieron un procedimiento sencillo para calcular c : Si la granulometra del material de un cauce se ajusta a una distribucin log-normal o a una distribucin logartmica, el esfuerzo cortante crtico se determina utilizando la ecuacin de Meyer-Peter y Muller, esto es c = 0.047 (s - ) Dn (3.7)

en la ec 3.7, Dn es el dimetro efectivo del material del cauce. El porcentaje que pasa o valor de n, que define el tamao del dimetro efectivo se vala con ayuda de la fig 3.4, donde g es la desviacin estndar geomtrica de la curva granulomtrica caracterstica del material del cauce. Por ejemplo, si g = 2.4, de la fig 3.4 se obtiene que para una distribucin log-normal n = 88.466%, y para una distribucin logartmica n = 80.292%. Es decir, si la granulometra del material del cauce sigue una distribucin log-normal, el tamao del dimetro que habra que utilizar en la ec 3.7 sera el D88.466 ; pero si la granulometra se ajustara a una logartmica, el tamao del dimetro efectivo sera D80.292.

Fig 3.4 Relacin entre n% y g para determinar el porcentaje que define el dimetro efectivo, Dn, asociado al esfuerzo cortante crtico, segn Cruickshank y Garca Flores.

En la fig 3.4 se observa que los resultados para una distribucin logartmica son diferentes y, en general, proporciona los menores valores de n, o sea esfuerzos cortantes crticos tambin menores que los de una distribucin log-normal. Esto se debe a que una distribucin logartmica tiene un intervalo ms chico de variacin de tamaos en sus extremos, y el de la parte superior es precisamente el que se presta para el acorazamiento. Por ello es importante conocer y detallar el extremo de la curva granulomtrica que compete a su fraccin gruesa. 3.1.3 Contribuciones de Garca Flores para distribuciones de tamaos Con objeto de simplificar an ms el procedimiento de Gessler para los casos de granulometras tericas, Garca Flores extendi el estudio y concret las ecuaciones para evaluar el esfuerzo cortante crtico, c , y el mximo dimetro medio de la coraza, (Dmc)mx, de las siguientes distribuciones de tamaos tericas: log-normal, logartmica, log-log, normal y circular. A continuacin se presentan sus resultados. 3.1.3.1 Ecuaciones para la distribucin log-normal Cuando los logaritmos de los dimetros siguen una ley normal (una recta en papel log-normal), se dice que la distribucin de tamaos es log-normal (cap 1). Si la granulometra caracterstica del material del cauce se ajusta a una distribucin del tipo log-normal, entonces las ecuaciones para evaluar el esfuerzo cortante crtico y el mximo dimetro medio de la coraza son: Para c : Si 1.25 = g = 2.30

c = (0.7005 + 0.0427 g + 0.3578 g2) [0.047 (s - ) D50] Si 2.30 < g = 10 c = (-0.6877 + 1.1665 g + 0.1286 g2) [0.047 (s - ) D50] Para (Dmc)mx : Si 1.25 < g = 2.60 (Dmc)mx = (0.7953 - 0.2997 g + 0.4219 g2) D50 Si 2.60 < g = 10 (Dmc)mx = (-0.8666 + 0.8626 g + 0.2149 g2) D50

(3.8)

(3.9)

(3.10)

(3.11)

Para diseo, los valores de c , deducidos de las ecs 3.8 y 3.9, se deben reducir dividindolos entre un factor de seguridad, f = 1.5. 3.1.3.2 Ecuaciones para la distribucin logartmica Cuando los logaritmos de los dimetros siguen una ley lineal (una recta en papel semilogartmico), se dice que la distribucin de tamaos es logartmica (cap 1). Si la granulometra caracterstica del material del cauce se ajusta a una distribucin logartmica, las ecuaciones para determinar el esfuerzo cortante crtico y el mximo dimetro medio de la coraza son: Para c : Si 1.1 = g = 2.40 c = (1.0657 - 0.1391 g + 0.2524 g2) [0.047 (s - ) D50] Si 2.40 < g = 10 c = (-0.1189 + 0.8632 g + 0.0392 g2) [0.047 (s - ) D50] Para (Dmc)mx : Si 1.25 = g = 2.40 (3.14) (3.13) (3.12)

(Dmc)mx = (0.6220 + 0.0876 g + 0.2132 g2) D50 Si 2.40 < g = 10 (Dmc)mx = (-0.3846 + 0.9010 g + 0.0479 g2) D50

(3.15)

Para diseo, los valores de c , inferidos de las ecs 3.12 y 3.13, tienen que ser disminuidos, dividindolos entre un factor de seguridad, f = 1.5.

3.1.3.3 Ecuaciones para la distribucin log-log Cuando los logaritmos de los dimetros siguen una ley logartmica (una recta en un papel logartmico), se dice que la distribucin de tamaos es del tipo log-log (cap 1). Si la granulometra caracterstica del material del cauce se ajusta a una distribucin de tamaos del tipo log-log, las ecuaciones para estimar el esfuerzo cortante crtico y el mximo dimetro medio de la coraza son: Para c : Si 1.25 = g < 2.30 c = (0.7608 0.0282 g + 0.1993 g2) [0.047 (s - ) D50] (3.16) Si 2.30 = g = 10 c = (-0.1071 + 0.8054 g + 0.0225 g2) [0.047 (s - ) D50] (3.17) Para (Dmc)mx : Si 1.20 < g = 2.40 (Dmc)mx = (0.3738 + 0.3879 g + 0.1256 g2) D50 Si 2.40 < g = 10 (Dmc)mx = (-0.2233 + 0.8713 g + 0.0269 g2) D50 (3.19) (3.18)

Para diseo, los valores de c, que resultan de las ecs 3.16 y 3.17, se deben reducir dividindolos entre un factor de seguridad, f = 1.5. 3.1.3.4 Ecuaciones para la distribucin normal Cuando los dimetros siguen una ley normal o de Gauss (una recta en papel probabilidad), se dice que la distribucin de tamaos es normal o gaussiana. Si la granulometra caracterstica del material del cauce se ajusta a una distribucin de tamaos normal, las ecuaciones para calcular el esfuerzo cortante crtico y el mximo dimetro medio de la coraza son: Para c : Si 0.025 = g/D50 = 0.425 (3.20)

c = [1.2183 - 0.0163 (g/D50) + 0.6916 (g/D50)2] (0.047 (s - ) D50) Para (Dmc)mx : Si 0.025 = g/D50 = 0.425

(Dmc)mx = [0.9966 + 0.0731 (g/D50) + 0.9654 (g/D50)2] D50

(3.21)

Para diseo, el c obtenido con la ec 3.20 tiene que ser disminuido, para lo cual se divide entre un factor de seguridad, f = 1.5.

3.1.3.5 Ecuaciones para la distribucin circular Cuando los dimetros siguen una ley circular (una recta en papel para distribucin circular), se dice que la distribucin de tamaos es circular (cap 1). Si la granulometra caracterstica del material del cauce se ajusta a una distribucin de tamaos del tipo circular, entonces las ecuaciones para estimar el esfuerzo cortante crtico y el mximo dimetro medio de la coraza son: Para c : y dado que D82.6145 = 0.43654 Dmx = 0.95437 D84 = 3.2584 D50 la ec 3.22 puede expresarse en las formas siguientes c = 0.020517 (s - ) Dmx c = 0.04486 (s - ) D84 c = 0.15314 (s - ) D50 Para (Dmc)mx : (Dmc)mx = 3.2782 D50 = 0.9602 D84 = 0.4392 Dmx (3.27) (3.24) (3.25) (3.26) (3.23) c = 0.047 (s - ) D82.6145 (3.22)

Para diseo, el c calculado con las ecs 3.24 a 3.26 se debe reducir, dividindolo entre un factor de seguridad, como indica la ec 3.6. 3.1.4 Mtodo simplificado para calcular el fenmeno de acorazamiento La mayora de los cauces estn formados por material de granulometra extendida, es decir, en el lecho existen desde arenas finas hasta gravas. Cuando el esfuerzo cortante producido por un escurrimiento es menor que el resistido por los granos ms grandes, en el transporte de sedimentos predomina el material de menor tamao; esto hace que, en la superficie del fondo, se incremente el nmero de partculas gruesas no arrastradas por el flujo, y se forme una capa denominada coraza. Conforme la coraza va cubriendo el fondo del cauce, el transporte de sedimentos disminuye hasta llegar a una condicin de estabilidad en la cual el material de la superficie del fondo queda inmvil y el transporte de sedimentos tiende a ser nulo. El acorazamiento se presenta, por ejemplo, en el tramo aguas abajo de una presa que, por un lado, regula los escurrimientos y, por el otro, corta el transporte de sedimentos hacia aguas abajo de la cortina. Tambin hay acorazamiento en los canales artificiales no revestidos (como los de riego). En esta tesis se renen las contribuciones y experiencias de diferentes autores sobre el tema del acorazamiento. En particular, se presenta un mtodo simplificado para calcular cmo va evolucionando la coraza; se discuten tambin ideas acerca de la influencia del proceso de acorazamiento del fondo en la resistencia al flujo o friccin y en el transporte de sedimentos. El modelo se calibra con los datos experimentales de Little (1972).

3.1.4.1 Revisin de mtodos para calcular el acorazamiento Segn Graf (1971), el primero en publicar una descripcin cualitativa del acorazamiento fue Du Buat en 1796. Entre los primeros experimentos sobre acorazamiento sobresalen los de Harrison (1950), quien trabajo con 3 granulometras. De estos estudios se destacan las conclusiones siguientes: la acumulacin de partculas inmviles en la superficie del fondo causa un incremento en la rugosidad efectiva; una capa de partculas inmviles cuyo espesor sea igual a un dimetro de dichas partculas es suficiente para detener la erosin; las partculas inmviles se arreglan ellas mismas en una formacin con ondulaciones caractersticas; el fondo se degrada como una unidad por lo que, si las partculas inmviles quedan uniformemente distribuidas a lo largo del lecho, el perfil longitudinal del fondo al final del proceso tiende a ser paralelo al del inicio del mismo. A partir de 1960, surgi gran inters en el tema del acorazamiento. Las investigaciones realizadas pueden clasificarse en dos grupos: las basadas en la formulacin de Gessler, con la cual se obtiene la granulometra de un fondo acorazado, y las que plantean otras formulaciones para calcular, en forma dinmica, cmo las distintas fracciones del material del fondo se integran al transporte de sedimentos o se quedan en el fondo formando la coraza. Little (1972) realiz experimentos en laboratorio usando granulometras tipo logartmico normal (tambin denominada log-normal) y ampli el banco de datos disponible. Hizo adems una comparacin entre los dimetros medios de la coraza obtenidos en sus experimentos y los calculados con el mtodo de Gessler (1970), en la cual, los valores calculados resultaron hasta 29 por ciento menores que los medidos. Al aplicar el mtodo de Gessler, se obtiene la granulometra que habr despus de que un o acte por un tiempo suficientemente largo; esto es, cuando se llegue a la coraza estable. Diversos autores, utilizando los datos disponibles, han propuesto expresiones (logradas en general por regresin) para calcular dimetros caractersticos de la coraza; entre ellas destacan las de Little (1972) Y Shen y Lu (1983). Cruickshank y Garca (1976) presentan una grfica que permite obtener el esfuerzo cortante crtico mximo para el cual la coraza es estable, vlida cuando la distribucin granulomtrica original es del tipo log-normal o logartmica. Maza y Garca (1988 Y 1990) han propuesto diversas expresiones para calcular los dimetros caractersticos de la granulometra final de la coraza y el esfuerzo cortante mximo, vlidas para granulometra logartmica o lognormal. Entre las ltimas investigaciones experimentales est la de Masconi y Subhash (l986), quienes trabajaron con material de granulometra log-normal para estudiar la evolucin del proceso de acorazamiento y la manera en que se rompe la coraza. Ellos observaron que la coraza se forma lentamente, y se rompe de manera relativamente rpida. Cuando ocurre la rotura de la coraza, hay un transporte generalizado del material; por tanto, la granulometra del material arrastrado es prcticamente la misma que el fondo tena originalmente. Con el mtodo de Gessler se predice la distribucin final de tamaos del sedimento en la coraza; sin embargo, no es posible determinar la forma en que se va desarrollando la coraza, el tiempo que tarda el fondo en quedar acorazado, ni el cambio en la resistencia

al flujo conforme el lecho se va acorazando. Los mtodos que se resumen a continuacin tienen el propsito de calcular cmo evoluciona la coraza.

Fig. 3.5 Representacin esquemtica de la composicin de la capa activa Estos procedimientos utilizan el concepto de zona o capa de mezclado (Borah et al, 1982). Conforme los rizos y dunas se mueven lentamente hacia aguas abajo, existe un mezclado continuo de sedimentos. El rea de estas formas de fondo puede ser considerada como zona de mezclado, bajo la cual el material del fondo permanece inalterado. Cuando el material es de granulometra no uniforme, los granos gruesos tienden a alojarse en la base de la zona de mezclado, constituyendo una lente de granos grandes. Para algunas condiciones de flujo, el material grueso puede permanecer inmvil; en este caso, el flujo lava los finos formando una capa o coraza que protege al material original del fondo. La erosin del material fino ocurre en una capa de espesor pequeo, que Borah et al (1982) definen como la capa activa. Cuando el fondo queda inmvil, la capa activa es la coraza. Por el contrario, durante la sedimentacin se deposita material en diversas capas activas formndose una nueva zona de mezclado. Para modelar el proceso de erosin o formacin de la capa activa, Borah et al (1982) proponen dividir la zona de mezclado en varias capas; aquella en contacto con el flujo se considera como la capa activa. El espesor, porosidad y distribucin de tamaos de las partculas de esa capa puede variar durante el proceso, pero se supone que la capa es homognea en un instante dado. El espesor de la capa activa se define considerando los volmenes de cada fraccin de la granulometra como T= 100 * Dl 1 el (3.28)

Pi
i =1

donde, Pi es el porcentaje de la fraccin i; Dl y son el tamao y la porosidad de la fraccin l, respectivamente. La fraccin l es el tamao ms pequeo (Dl) de material que no puede trasportar el flujo. Algunas de las fracciones en la zona de mezclado no pueden ser movidas por el flujo y entonces contribuyen a la formacin de la coraza. A fin de tomar en cuenta la evolucin en el tiempo del espesor de la capa activa, Borah et al (1982) proponen calcular la composicin de tamaos de la misma segn evoluciona. Para ilustrar el procedimiento, considrese una capa compuesta por 3 tamaos de sedimento: D1< D2 < D3 (fig 3.5); el flujo en el cauce es capaz de trasportar los tres tamaos. Debido a su dimensin, las primeras partculas erosionadas son las de

dimetro D1 expuestas en la superficie; poco despus son removidas de la superficie las partculas de tamao D2, lo cual permitir que partculas ocultas de tamao D1 se muevan. Finalmente se movern las partculas de tamao D3 y, por lo tanto, las de los otros dos tamaos que estaban ocultas. As, una parte de las partculas de tamao D1 es removida de la superficie, y otra se remueve cada vez que algunos granos de tamao D2 son trasportados; cuando se mueven los granos de tamao D3, tambin lo hacen una parte de los de tamao D2 y dos partes de los de tamao D1. Generalizando este ordenamiento para una granulometra cualquiera, se forma una matriz cuyos elementos por renglones indican el nmero de veces que cada fraccin Dj,i j < i, est expuesta a ser erosionada una vez que el grano de la fraccin Di se movi. Los detalles del mtodo se pueden consultar en Borah et al (1982). En el desarrollo del modelo matemtico IALLUVIAL, investigadores de la Universidad de Iowa, EUA, presentan un mtodo relativamente sencillo, con el cual se calcula el proceso de acorazamiento la idea de Karim y Holly (1986) es ir obteniendo el porcentaje de la superficie que va quedando cubierta por sedimentos gruesos conforme se forma la coraza. La virtud de este enfoque es la de modelar el desarrollo del fenmeno desde su etapa inicial hasta la final, incluyendo la obtencin de cmo evoluciona la erosin, la extensin del acorazamiento, as como la duracin aproximada del proceso. Una de las ltimas contribuciones al tema es el desarrollo de un modelo matemtico de fondo mvil generalizado propuesto por Rahuel et al (1989), y Holly y Rahuel (1990). Para incluir el acorazamiento, estos autores toman en cuenta el concepto de capa de mezclado introducido por Borah et al (1982). El espesor de esta capa (ec 3.28) es una medida de la cantidad del sedimento que debe ser erosionada del fondo para formar la coraza con granos de tamao Dl o mayores. Se considera que, conforme el fondo se va acorazando, no existe intercambio de sedimento entre la capa de mezclado y el material original subyacente. Por tanto, una vez arrastradas todas las partculas que pueden ser removidas, se forma la coraza con partculas de dimetro Dl o mayor. El espesor de la capa de coraza tiende a Dl, que ser en caso extremo el espesor de la capa de mezclado. Como en casi todos los temas de ingeniera fluvial, solo existen por el momento mtodos aproximados para modelar el acorazamiento, aunque se ha avanzado mucho en el tema. Se recalca que este es un fenmeno complejo y, tal vez por ello, los mtodos son cada vez ms sofisticados (Borah, 1987). A continuacin se plantean las ideas de la Universidad de Iowa al respecto, para luego discutir la metodologa propuesta. 3.1.4.2 Bases del mtodo Cuando el material que constituye un cauce es de granulometra extendida, los procesos de erosin y acorazamiento son simultneos: para cierta socavacin de la cota del fondo, una parte del material es removida y trasportada por la corriente, mientras otra contribuye a la formacin de la coraza. As, si el material del cauce est constituido por un total de m fracciones de diferentes tamaos, de las cuales las fracciones ms grandes (de la l a la m) no pueden ser movidas por el esfuerzo cortante que produce el flujo, es posible modelar el proceso del acorazamiento calculando el rea de la superficie del fondo que va quedando cubierta por los granos de tamaos l a m. Para determinar esta rea acorazada, se obtiene primero el volumen del sedimento que queda inmvil conforme el fondo se degrada; se calcula despus el rea ocupada por dicho volumen,

suponiendo que el espesor de la coraza es equivalente al dimetro medio de los granos ms grandes. Esto se muestra de manera esquemtica en la fig 3.6, donde se nota que el rea de la superficie acorazada se incrementa con la degradacin del cauce y que tal proceso termina cuando se llega a formar una capa completa de material grueso sobre la superficie, la cual se conoce como coraza. El desarrollo presentado a continuacin se basa en el planteado originalmente por Karim et al (1983); la versin formal del mtodo puede verse en Karim y Holly (1986). Al hacer la revisin bibliogrfica se discuti que la coraza contiene sedimentos de todos los tamaos del material original del cauce, pero en el presente anlisis se supondr que en la superficie del lecho hay una parte de material constitutivo de la coraza y otra susceptible de moverse; esto significa que las particu1as finas que se encuentran entre otras ms grandes no contribuyen a la formacin de la coraza, ya que algunas de ellas podrn pasar a formar parte del trasporte de sedimento de un momento a otro y otras quedarn protegidas por las partculas de la coraza. Desde este punto de vista no es necesario determinar la distribucin de tamaos del material de la coraza al final del proceso; sin embargo, dado que durante el proceso del acorazamiento hay un cambio en la resistencia al flujo, se requiere estimar la granulometra de la coraza para tomar en cuenta el incremento del tamao de los granos en la resistencia al flujo. Por esto, posteriormente se describe un procedimiento para considerar el cambio en la granulometra durante el proceso de acorazamiento.

Fig. 3.6 Representacin esquemtica del acorazamiento de un fondo plano Deduccin de la expresin para calcular el porcentaje de rea acorazada Supngase que, durante un intervalo de tiempo t, la cota de fondo del cauce en un tramo del ro disminuye en una magnitud Az(t), segn se muestra en la parte superior de la fig 3.6, y se desea conocer el volumen de todos los slidos no movidos por el flujo y depositados en el fondo a causa de su tamao. Para una superficie de rea unitaria y espesor Az(t), el volumen total es Vm = A z(t) (3.29)

Si el material es de tamao uniforme, el volumen de los slidos en dicho volumen es Vs = /e Vm = (1-) Vm donde e Vv Vs : : : : porosidad de la muestra (Vv/Vm) relacin de vacos (Vv/Vs) volumen de vacos volumen de slidos (3.30)

Sustituyendo la ec 3.29 en la 3.30, y considerando que el rea A es unitaria, Vs = ( 1 - ) z(t) (3.31)

En caso de que la muestra estuviera compuesta de material no uniforme y se conociera su granulometra, el volumen de slidos por unidad de rea en planta se calcula como

Vs = (1 )Z (t ) Pk
k =1

(3.32)

donde Pk es el tamao del intervalo del porcentaje en peso en la curva granulomtrica con dimetro caracterstico Dmk. Sea Dl el dimetro caracterstico ms pequeo de las partculas no movidas por el flujo (determinado empleando, por ejemplo, el criterio de Shields). El volumen de los granos de tamaos caractersticos l a m de la curva granulomtrica que no pueden ser trasportados durante la erosin del tramo, y que por tanto forman parte de la coraza, ser entonces la expresin de la ec 3.32. El volumen de una partcula de la fraccin k, con dimetro caracterstico Dmk y factor de forma fs (definido como fs = c/(ab), donde c es el semieje menor y los otros dos semiejes se consideran iguales, a = b, resulta ser Vk = /6 fs Dmk3 (3.33)

El nmero de slidos de tamao Dmk que se agregan a la coraza debido a una disminucin de magnitud z, para una superficie de rea unitaria, es el cociente del volumen total de partculas de tamao caracterstico k entre el volumen de un grano del mismo tamao, o sea, Nk = ( 1 - ) z(t) Pk /6 fs Dmk3 (3.34)

Cada una de las Nk partculas ocupa una parte de la superficie, cuya rea es la proyeccin de la partcula sobre la superficie. Si se considera que esta proyeccin es un crculo, el rea Afk de cada grano de dimetro Dmk se expresa Afk = 1.5 ( 1 - ) z(t) Pk fs Dmk (3.35)

A fin de representar cada partcula con forma esfrica, se propone calcular el dimetro equivalente Dk; para ello, se iguala el volumen de la partcula de dimetro medio Dmk y factor de forma fs con el volumen de la partcula esfrica de dimetro equivalente Dk, por lo que la ec 3.35 queda Afk = CA ( 1 - ) z(t) Pk Dk (3.36)

donde CA = 3/2 fa-2/3. Un valor comn del factor de forma, fs, es 0.7; por tanto, CA es igual a 1.90. Para tomar en cuenta el carcter aleatorio del movimiento de las partculas, se multiplica la ec 3.36 por la probabilidad qk de que una partcula de tamao Dk se quede como parte de la coraza, determinada con el criterio de Gessler (1967), esto es, Afk = CA ( 1 - ) z(t) Pk qk Dk (3.37)

el porcentaje de rea acorazada Afk calculado con la ec 3.37 debe sumarse al rea del fondo previamente acorazada con granos del mismo tamao; es decir, Af (t, k) = Af (t -t, k) + CA ( 1 - ) z(t) Pk qk Dk (3.38a)

En esta expresin, Af (t,k) es el porcentaje acumulado del rea del fondo cubierta por granos de tamao Dk hasta el tiempo t. El rea acorazada total en el tiempo t es, finalmente, la suma de las contribuciones en ese instante de todas las fracciones de la curva granulomtrica, desde la l hasta la m,
A f (t ) =

k =l (t )

(t , k )

(3.38b)

donde, para las fracciones k < l, Af (t,k) es igual acero, puesto que estas fracciones no contribuyen a la formacin de la coraza.
3.1.4.3 Procedimiento simplificado para calcular el rea acorazada, At(t)

Las ecs 3.38 tienen carcter general y son vlidas para cualquier granulometra; Holly y Karim (1986) presentan una verificacin de la metodologa con datos del ro Misuri, EUA. Al analizar las ecs 3.38, puede verse que, al aplicar el mtodo, se requiere almacenar los valores de la granulometra de cada seccin, as como el rea acorazada acumulada para cada fraccin. Por otro lado, es necesario ir haciendo gran cantidad de operaciones, incluida la suma de la ec 3.38b; adems, este clculo debe efectuarse en cada seccin. Se realizaron varios clculos para determinar qu tan sensible es el resultado al nmero de fracciones de la curva granulomtrica en la suma de la ec 3.38b y se encontr que, en la mayora de los casos, era necesario hacer la suma con un nmero grande de valores. Esto, adems del requerimiento de memoria para almacenar los datos de granulometra de cada seccin, implica gran cantidad de sumas y operaciones, que deben efectuarse para cada tramo e intervalo de tiempo. Por ello, se

busc un mtodo alternativo con el cual se pudieran simplificar las sumas, as como reducir lo ms posible el almacenamiento de datos de la curva granulomtrica para cada tramo. En el procedimiento presentado a continuacin se supone que el material del fondo en el tramo del ro de inters puede caracterizarse empleando una distribucin granulomtrica de tipo log-normal, representada con nicamente dos parmetros (el dimetro D50 y la desviacin geomtrica estndar,g). Esta es una hiptesis comnmente aceptada en relacin con cauces arenosos; se aclara que es posible generalizar la metodologa aqu propuesta a otras distribuciones granulomtricas. La profundidad de erosin Z(t) es, en general, independiente de las fracciones de la curva granulomtrica y se obtiene tomando en cuenta las condiciones del flujo con la ecuacin de continuidad del sedimento (por ejemplo, Berezowsky y Jimnez, 1989); y las constantes de las ecs 3.38 aqu se consideran independientes de la fraccin k. Por tanto, el operador afecta solamente a los trminos qk, Pk y Dk. La idea es buscar una expresin que sustituya a la suma. En primer lugar, se plantea remplazar los trminos de la suma relacionados con la distribucin granulomtrica por la siguiente igualdad.
m Pl P = k DMR k =1 D k

(3.39)

Se define P1 como el porcentaje total, en peso, del material que se queda en ese instante formando parte de la coraza; DMR es un dimetro representativo o equivalente de ese material (fig 3.7)
Clculo de P1

Para obtener el tamao lmite del material que se queda formando la coraza, se igualan los esfuerzos cortantes crtico Ck (calculado con el criterio de Meyer-Peter y Muller) y actuante 0. Si se supone un nmero alto de Reynolds del grano, el dimetro D1 ser Dl = donde D1 h S Ss : : : : hS 0.047(Ss 1) (3.40)

dimetro del grano ms pequeo no arrastrado por el flujo tirante del escurrimiento pendiente hidrulica densidad relativa del material del fondo.

Los dimetros D1 se supone que siguen una ley de distribucin log-normal y se calculan con Di = D50 gZi donde Di zi g : : : tamao del material por debajo del cual queda el i por ciento de la muestra de suelo, en peso variable normal estndar desviacin estndar geomtrica (para distribucin log-normal g = (D84/D16) (3.41)

Obtenido Di con la ec 3.40, se despeja zi de la ec 3.41; la probabilidad P1 se calcula entonces utilizando la funcin de distribucin de probabilidad acumulada de Gauss; este clculo puede hacerse de manera relativamente sencilla empleando las funciones propuestas por Abramowitz y Stegun (1972). El valor de P1 es el complemento de P1, es decir, Pl = 1 - Pl (3.42)

Fig 3.7 Esquema para ubicar DMR en la curva granulomtrica Clculo del dimetro DMR

Para un flujo y una distribucin granulomtrica dados, despus de obtener P1, el valor del dimetro representativo del material que se queda formando la coraza se despeja a partir de la ec 3.39

DMR =

Pl P Dk k =1 k
m

(3.43)

La idea es obtener alguna expresin que permita calcular el valor de DMR en funcin de variables conocidas. Para ello, se determinaron gran cantidad de valores de DMR usando la ec 3.43 y cinco granulometras tipo log-normal con desviacin estndar geomtrica g de 3 a 7; para cada valor de DMR se calcul su porcentaje PDMR en la curva granulomtrica (fig 3.7). Se graficaron los valores de P1 contra PDMR; en la fig 3.8 se muestran las lneas extremas que corresponden a las distribuciones con desviacin estndar geomtrica 3 y 7. La variable PDMR puede relacionarse con las variables P1 y g mediante una funcin del tipo PDMR = Pl g

Fig 3.8 Relacin entre los porcentajes P1 y PDMR

Dado que el efecto de g es relativamente pequeo, se consider despreciable el parmetro ; usando todos los valores de P1 y PDMR obtenidos con las cinco granulometras se obtuvo la siguiente expresin por regresin lineal PDMR = 0.98804 Pl0.458 (3.44)

que tiene un coeficiente de correlacin lineal de 0.996. La grfica de la ec 3.44 se incluye tambin en la fig 3.8; una manera sencilla de obtener PDMR es PDMR P1. En resumen, con la ec 3.44 se calcula el valor de PDMR y con la ley de distribucin de probabilidad log-normal se determina el dimetro DMR que corresponde al porcentaje PDMR, donde se emplean las expresiones aproximadas referidas por Abramowitz y Stegun (1972). Se aclara que, como para calcular los dimetros en la ec 3.43 se emplea la ec 3.41, el resultado es independiente del valor del D50 utilizado; es decir, la expresin 3.44 es vlida para cualquier D50 y g de granulometra tipo log-normal.
Movimiento estocstico del material de P1

La expresin que toma en cuenta la naturaleza estocstica del movimiento del sedimento tiene la forma

Pk qk k =1 D k
m

(3.45)

donde qk es la probabilidad, segn Gessler, de que un grano de tamao k no se mueva. Para considerar este efecto, aqu se propone encontrar la probabilidad q(Dq) de un grano representativo de tamao Dg que, aplicada al cociente P1 / DMR, produzca resultados similares a la suma de la ec 3.45

Pk Pl qk = q( D g ) DMR k =1 D k
m

(3.46)

El dimetro Dg se determin por prueba y error. Al emplear el dimetro DMR se vio que los resultados de la suma eran, en general, 20 por ciento menores que los obtenidos con P1 q(Dg) /DMR, de lo cual se dedujo que un dimetro un poco menor que DMR producira resultados ms aproximados. Con base en esto, se propuso utilizar un dimetro Dg obtenido con el promedio geomtrico de los dimetros D1 y DMR, es decir, Dg = Dl DMR (3.47)

y calcular la probabilidad de Gessler con ese dimetro Dg, lo cual produce excelentes resultados. De esta manera, la expresin simplificada para calcular el porcentaje de rea acorazada en una seccin es Af (t) = Af (t -t) + 1.902 (1 - ) z(t) Pl q(Dg) DMR (3.48)

En esta expresin se han eliminado las sumas en k, hecho que permite determinar el porcentaje de rea acorazada en forma simple y rpida, sin necesidad de almacenar gran cantidad de informacin.
3.1.4.4 Lmites en el clculo del rea acorazada

Durante el clculo del rea acorazada es necesario tomar en cuenta lo siguiente: a) El rea acorazada se incrementa nicamente cuando hay erosin; en caso de que en un tramo del ro se tuviera depsito de sedimento, el rea acorazada debera disminuir si granos finos cubren a los de la coraza. Dado que el clculo de esta condicin es muy complicado, en el procedimiento aqu descrito se considera que el rea acorazada es constante cuando hay depsito, y se incrementa si hay erosin. b) Existe un esfuerzo cortante mximo que un fondo acorazado resiste; dicho esfuerzo se puede estimar usando el procedimiento de Cruickshank y Garca (1976) basado en el mtodo de Gessler (1967) o con las expresiones propuestas por Maza y Garca (1990). El rea acorazada se incrementa si el cortante actuante es menor que este cortante mximo. Por el contrario, si el cortante actuante en un tramo dado rebasa el mximo que resiste la coraza, esta se rompe; es decir, todos los granos se mueven, lo cual equivale a que el rea acorazada sea cero.
3.1.4.5 Modificacin de la granulometra

Como se mencion al principio, conforme el fondo se va acorazando, se modifica la distribucin de tamaos del material del fondo. La granulometra de la coraza final en un tramo dado depende del esfuerzo cortante mximo que se presente durante la formacin de la coraza. Desafortunadamente, el valor de este esfuerzo se desconoce antes del clculo, por lo que no es posible determinar dicha distribucin previamente. Sin embargo, al aplicar el procedimiento de Gessler para determinar la distribucin de

tamaos de la coraza para el cortante mximo, se encontr que, sobre todo para la distribucin log-normal, si los cortantes son cercanos al mximo, la distribucin de tamaos de la coraza cambia relativamente poco. Por ello, se decidi considerar como granulometra de la coraza a la obtenida para el esfuerzo cortante mximo. La distribucin de tamaos conforme el fondo se va acorazando se calcula para cada tramo con una funcin de la forma Di = Dio + (Dic Dio) f donde Dio Dic f : : : dimetro de la granulometra original del cauce dimetro de la granulometra de la coraza mxima funcin de calibracin, que depende del rea acorazada. (3.49)

Si la funcin de calibracin vara entre cero y uno, con la ec 3.49 se hace una interpolacin entre el dimetro de la granulometra original y el de la coraza mxima. As, si f es cero, Di es el dimetro original, y si f es igual a uno, Di es el de la coraza final. La variacin de la granulometra depende de la funcin de calibracin.
3.1.4.6 Efecto del acorazamiento en la resistencia al flujo

En ros, la resistencia al flujo se debe a dos efectos, la rugosidad de los granos y la rugosidad de las formas de fondo. Mientras que la estimacin de la rugosidad de los granos es relativamente simple, no se dispone de un criterio general para estimar la de las formas de fondo. Por ello, para calcular la resistencia al flujo en ros se recurre a mtodos semiempricos; entre los ms conocidos estn los propuestos por Einstein, Cruickshank-Maza, Engelund-Hansen. Los mtodos de clculo de resistencia al flujo fueron planteados originalmente para obtener la velocidad media U de un flujo uniforme. En general, se establece una funcin del tipo U = U (g, h, Di, S, etc) donde Di es algn dimetro caracterstico. Cada autor establece un enfoque terico y ajusta datos experimentales para incluir los efectos no medibles (como la rugosidad de las formas de fondo). Por ejemplo, en el mtodo de Cruickshank-Maza, la velocidad media se calcula con una expresin del tipo U/w50 = C1 ( h/D84)C2 (S/)C3 donde w50 velocidad de cada del dimetro D50 peso especfico sumergido relativo [=(S - )/] C1 a C3 constantes. En esta ecuacin se emplean como dimetros caractersticos D50 y D84. Para considerar el cambio en la rugosidad se plantea incluir este en la granulometra, formulando las siguientes hiptesis: a) al principio del fenmeno se tiene la distribucin de tamaos original.

b) existen inicialmente formas de fondo, por lo que es necesario emplear mtodos para calcular la resistencia al flujo que incluyan ese efecto. c) estos mtodos son vlidos para flujo no uniforme y no permanente. d) en las experiencias de laboratorio para estudiar acorazamiento, se parte en general de un fondo plano (sin formas de fondo) y cuando el fondo se acoraza se tiene tambin fondo plano; se sabe que la erosin es muy pequea (del orden de unos cuantos dimetros segn Gessler, 1971). Sin embargo, esto no es necesariamente lo que sucede en un ro, pues en el inicio se tienen formas de fondo y aunque el fondo tiende a hacerse plano, no siempre llegan a desaparecer tales formas. Este efecto se incluye en el cambio de granulometra calculado con la ec 3.49, pues conforme los granos tienden a ser ms gruesos, la rugosidad de estos tiende a ser ms importante.
3.1.4.7 Efecto del acorazamiento en el trasporte de sedimentos

El gasto slido se obtiene en general como funcin de la velocidad, las caractersticas del material del fondo, la resistencia al flujo, etc; es decir, Gs = Gs (U, , Di, S, etc.) De acuerdo con los trabajos experimentales del proceso de acorazamiento (por ejemplo, Gessler [1967] y Little [1972]), conforme el fondo se va acorazando disminuye el trasporte de sedimentos. Ellos informan que an despus de un largo perodo de iniciado el experimento, hay trasporte de sedimentos, aunque este tiende a ser prcticamente nulo, formado solo de granos muy finos. Dado que la granulometra y la resistencia al flujo se modifican durante el proceso del acorazamiento, tambin se modifica el trasporte de sedimentos, ya que su clculo depende de estas variables. Karim et al (1983), adems de modificar el gasto slido por efecto del acorazamiento, proponen considerar su disminucin en funcin del porcentaje de rea acorazada con la siguiente expresin Gs* = (1 - f) Gs (3.50)

donde f es una funcin de calibracin, que vara entre 0 y 1, y depende del porcentaje de rea acorazada. La idea es que, conforme el fondo se va acorazando, el gasto slido disminuye. Se hace notar que se emplea la misma funcin f que en el apartado 3.1.4.6.
3.1.4.8 Calibracin del modelo, obtencin de la funcin f

La calibracin del modelo consiste en encontrar la funcin f capaz de lograr que el modelo matemtico se ajuste a algn conjunto de datos medidos. Esta funcin depende del grado de acorazamiento en el cauce; para obtenerla, se emplean los experimentos realizados por Little (1972). Este autor trabaj en un canal rectangular con paredes de vidrio de 18.3 m de largo y 0.6 m de ancho, donde coloc una cama de material arenoso de 0.1 m de espesor en una longitud de 12.2 m. En el canal se recirculaba agua limpia a

20 C, con peso especifico de 998.8 kg/m3 y viscosidad cinemtica de 9.29 X 10-7 m2/s. El peso especfico relativo del material arenoso fue 2.65 y la porosidad del material se estim en 0.4. La pendiente inicial de la plantilla del canal en los experimentos fue de 0.002. Con un procedimiento muy detallado, Little elabor muestras de arena con granulometras tipo log-normal, todas con D50 = 1 mm y desviacin estndar geomtrica, g, entre 1.12 y 3.05. Cada experimento se realiz para un gasto constante y se procur tener un flujo uniforme; asimismo se registr la cantidad de material recolectado en la salida aguas abajo del canal a intervalos de tiempo fijos. Las pruebas se llevaron a cabo hasta que el material recolectado al final del canal fue muy pequeo. El objetivo principal de los experimentos de Little fue la obtencin de la granulometra de la coraza final y, por tanto, hay pocos datos relativos al proceso mismo; no obstante, para la calibracin del proceso de acorazamiento que aqu interesa pueden emplearse las mediciones de la variacin del gasto slido al final del canal. Para calibrar el procedimiento propuesto, se escogieron los datos del experimento 6-1. La granulometra al inicio de este ensayo es tipo log-normal con D50 = 0.001 m y g = 3.01. Aunque Little indica la granulometra al final del proceso del acorazamiento de sus experimentos, este dato no se emple para la calibracin, pues en general no se conoce; la granulometra de la coraza se obtuvo empleando el mtodo de Gessler para el esfuerzo cortante critico mximo del material. En este caso el esfuerzo cortante crtico mximo es de 0.31 kg/m2; los dimetros caractersticos empleados se dan en la tabla 3.1. El gasto lquido en este experimento es de 0.0127 m3/s. Los valores de la evolucin del gasto slido en el tiempo obtenidos a partir de las mediciones se muestran en la fig 3.9.

Fig 3.9 Comparacin entre gasto slido calculado y medido (experimento 6-1 de Little)

Para los clculos, el canal se dividi en 15 tramos de 1.2 m de longitud, es decir se emple un canal de 18 m de longitud. Esto se hizo para disminuir el efecto de la condicin de frontera de aguas abajo en el tramo de inters (los primeros 12.2 m). Se us un incremento de tiempo, t, de 20 s.

Tabla 3.1 Granulometras del experimento 6-1 Dimetro caracterstico D35 D50 D65 D84 D90 Granulometra inicial (m)
0.000703 0.001100 0.001423 0.002487 0.003234

Granulometra final (m)

Como criterio de resistencia al flujo se usa la formulacin de Engelund-Hansen. Para las condiciones iniciales, el tirante de flujo uniforme resulta ser de 6.2 cm., mientras que Little midi 5.6 cm. La diferencia se debe a que, para las condiciones hidrulicas del experimento, con el mtodo de Engelund se tienen dunas en el canal y, al haber mayor resistencia al flujo, el tirante de flujo uniforme es mayor, mientras que Little inicia sus experimentos a partir de un fondo plano. Este autor seala que, aunque no se aporta sedimento al canal aguas arriba, hay una tendencia a la formacin de dunas que conforme el fondo se va acorazando tienden a desaparecer. Empleando los datos de gasto slido de los primeros cien minutos del experimento, se calibr una ecuacin de gasto slido de la forma Gs = K Q2 S3/2 h-1/2 B g D35 2 (3.51)

donde B es el ancho del canal, la constante K obtenida en la calibracin vale 0.075. Engelund encontr es sus experimentos que esta constante vale K = 0.05. En cuanto a la calibracin de la funcin f en las ecs 3.49 y 3.50, la idea ms simple es que, dado que Af (t) vara entre cero y uno, la funcin de calibracin sea linealmente proporcional a Af (t), como suponen Karim et al (1983), por lo que la granulometra y el gasto slido se modifican linealmente conforme el fondo se va acorazando. Al comparar los clculos basados en esta hiptesis con los experimentos se observ que, aunque se sigue la tendencia de las mediciones, el gasto slido calculado disminuye ms rpido que el medido. Se probaron distintas funciones de calibracin; se encontr que cuando f se calcula con la funcin f(Af) = tanh(L1(1-Af) + L2Af) tanh(Ll) donde DEN : tanh (L2) tanh (L1) (3.52)

se obtienen los resultados que se muestran en la fig 3.9; esta funcin depende de la variable Af y de los lmites L1 y L2. En la fig 3.10 se muestra la grfica de esta funcin con limites L1 = -3 y L2 = 3, encontrados por prueba y error.

Fig 3.10 Funcin de calibracin, f (Af)

Volviendo a la fig 3.9, en ella puede verse el gasto slido calculado en la seccin que corresponde a la salida del canal comparado con el medido; la concordancia es muy buena. Se aclara que, debido a la pequeez del gasto slido, el eje vertical logartmico aumenta en este caso las diferencias: si se usa un eje lineal en la figura, la concordancia es mejor. Una vez calibrada la funcin f, se procedi a revisar el funcionamiento del modelo propuesto en otro caso; para ello, se emplearon los datos y resultados del experimento 3-4 de Little. El gasto lquido fue de 0.0162 m3/s y la pendiente 0.002. La granulometra inicial tiene g = 2.50 y la de la coraza final, calculada con el criterio de Gessler, se indica en la tabla 3.2. Para este caso el esfuerzo cortante crtico de la coraza mxima es de 0.23 kg/m2. La variacin con el tiempo del gasto slido en la ltima seccin del canal, calculado con el modelo matemtico y los resultados experimentales se muestran en la fig 3.11. La comparacin es bastante satisfactoria, pues aunque hay una diferencia al principio del proceso esta podra corregirse cambiando la constante K en la ec 3.52. Por tanto, se considera que la funcin de la ec 3.53 es la adecuada.

Fig 3.11 Comparacin entre gasto slido calculado y medido (experimento 3-4 de Little)