Presa del arroyo Manina (Madrid)

Este proyecto tiene por objeto la construcción de la presa de Arroyo Manina, para poder hacer frente al aumento de población que va a sufrir el municipio de Torrelodones en los próximos años. En los alrededores existen pequeñas presas pero no albergan un volumen suficiente para que sirvan de abastecimiento o su cota queda excesivamente baja. Además, las poblaciones se han acercado a la cota máxima de estos embalses, por lo que un posible recrecimiento implicaría el desplazamiento de gente. Por tanto, se ha visto la necesidad de construir una presa desde cero.También hay que añadir, que tanto el municipio de Torrelodones como el de Hoyo de Manzanares se encuentran situados dentro del Parque Regional de la Cuenca Alta del Manzanares. Este parque tiene zonas altamente protegidas y otras que se pretende se utilicen como zonas recreativas.

Aliviaderos tipo laberinto y en tecla de piano

Las conferencias presentadas durante la Jornada Técnica sobre Avances en investigación aplicada en seguridad hidráulica de presas organizada por el Centro de Estudios Hidrográficos del CEDEX y la Universidad Politécnica de Madrid el mes de junio de 2013 versan sobre la mejora de la seguridad hidráulica de las presas, una materia de especial importancia en nuestro país (habida cuenta del elevado número de infraestructuras de regulación en servicio y de su antigüedad media) y constituye una tarea en la que tanto el CEDEX como la Universidad Politécnica de Madrid han trabajado activamente en los últimos años. El Departamento de presas de la ETSICCP, el Centro de Estudios Hidrográficos del CEDEX el Centro Internacional de Modelación Numérica (CIMNE), la empresa VEMSA y la compañía de ingeniería JESÚS GRANELL Ingenieros Consultores viene acometiendo un proyecto de investigación en materia de vertederos tipo laberinto. El objetivo principal de la investigación es la obtención de conocimientos y herramientas que permitan un diseño razonado de este tipo de estructuras hidráulicas. Además, se estudian otros efectos en su comportamiento, como son las condiciones hidráulicas de aproximación a la embocadura. La investigación se basa en el ensayo hidráulico en modelos físicos complementado por modelos numéricos mediante aplicaciones de última generación. Descripción del proceso de diseño y construcción del modelo.

Modelación física y numérica de aliviaderos con cajeros altamente convergentes

Presentación de las conclusiones del proyecto de investigación ALCON. La publicación del Reglamento Técnico sobre seguridad de Presas y Embalses en el año 1996 ha generado la necesidad de adaptar muchas de las presas en explotación en España para cumplir esta nueva normativa, debiendo ampliar gran parte de ellas su capacidad de desagüe y también aumentar los resguardos. Muchas de las presas españolas se encuentran aún en proceso de adaptación a los nuevos criterios debido en ocasiones a que la ampliación del aliviadero existente es muy costosa o plantea dificultades técnicas importantes. Por esta razón, en los últimos años se han desarrollado nuevas tipologías de aliviaderos, capaces de satisfacer las necesidades planteadas. Una de estas soluciones innovadoras son los aliviaderos con cajeros altamente convergentes. Esta tipología, adecuada para presas de gravedad de hormigón, tiene la particularidad de que la longitud del vertedero es sensiblemente mayor que el ancho del cuenco amortiguador. Para conducir el caudal vertido al cuenco, se disponen unos cajeros con fuerte convergencia. Si bien este tipo de aliviaderos puede plantearse en presas de nueva construcción, la ampliación de capacidad de presas existentes mediante esta tipología resulta sencilla.Desde el punto de vista de la explotación hidroeléctrica, la ampliación de capacidad conseguida con esta tipología permite aumentar la cota del nivel máximo normal (NMN) gracias a la reducción de la altura de lámina vertiente, manteniendo los resguardos exigidos por la normativa. Ello se traduce en un incremento considerable de la capacidad de generación hidroeléctrica, aumentando tanto la altura del salto como el volumen embalsado. Los aliviaderos con cajeros altamente convergentes ya han sido utilizados en varias presas tanto en España como en el resto del mundo. Sin embargo no existen referencias con criterios de diseño de carácter general para esta tipología, por lo que se hace necesario a la hora de proyectar un aliviadero de este tipo realizar ensayos en modelo físico específicos en cada caso. Con el objetivo desarrollar una metodología para el cálculo hidráulico y proyecto de este tipo de aliviaderos se está llevando a cabo el proyecto ALCON, incluido en el Plan Nacional de Investigación Científica del Ministerio de Economía y Competitividad, en el que participan ALATEC, UPM Y CIMNE. En el proyecto se combinan la modelación física y numérica para alcanzar los objetivos planteados.

Estudio experimental sobre la reacción álcali-sílice en el hormigón producida por áridos de reacción lenta

La finalidad de esta investigación se enmarca dentro de los estudios sobre hormigones de presas llevados a cabo en el Laboratorio Central de Estructuras y Materiales del CEDEX. En España se han diagnosticado 18 obras afectadas tanto por la reacción álcali-sílice de tipo rápido como por la de tipo lento. Dos de de las obras, fabricadas con áridos graníticos, no presentan signos de deterioro pero en laboratorio se han hallado productos expansivos. En las 16 presas españolas restantes, el hormigón de 10 de ellas estaba fabricado con áridos graníticos, de las cuales, 7 están afectadas por la reacción de tipo lento. Sin embargo, en las clasificaciones internacionales de las rocas potencialmente reactivas se establece habitualmente que las rocas graníticas son inocuas o de reactividad muy baja. En los casos puntuales encontrados en la literatura, la reactividad de este tipo de áridos se encuentra frecuentemente asociada al cuarzo microcristalino que contienen y ocasionalmente se asocia esta reactividad también al cuarzo deformado y/o microfisurado. En la tesis doctoral de Víctor Daniel Lanza Fernández, también realizada en el CEDEX, se han tratado los áridos de reacción rápida, dando una descripción detallada de los componentes que intervienen en este tipo de reacción y presentando un ensayo eficaz para su detección en el laboratorio. La investigación desarrollada en la presente tesis doctoral se ha centrado en los áridos de reacción lenta, mucho menos estudiados. A partir del estudio bibliográfico realizado sobre este tipo de reacción se ha detectado ciertas lagunas, basadas principalmente en la falta de un método de detección en el laboratorio para los áridos de reacción lenta. No se ha encontrado un procedimiento en el estudio petrográfico que sea fiable y cuantificable para este tipo de áridos. El ensayo acelerado de barras de mortero, que para los áridos de reacción rápida supone un método rápido y fiable, falla en la detección de los áridos de reacción lenta al aplicar los límites normalizados. Si bien en publicaciones recientes se ha propuesto la posibilidad de ampliar el tiempo de tratamiento hasta los 90 días, la duración del ensayo es más larga de lo deseado. Para resolver estas lagunas, se han tomado áridos de obras reales afectadas por la reacción álcali-sílice (en algún caso también de las canteras de donde se extrajeron los áridos para la ejecución), con lo que la reactividad de estos áridos queda demostrada por su comportamiento en obra. Al objeto de aumentar la muestra de ensayo y utilizar también áridos inocuos, se han tomado muestras de canteras, tanto en uso como abandonadas. Sobre estos áridos se ha realizado una caracterización completa con los ensayos actualmente disponibles para el estudio de su reactividad con los álcalis del cemento: estudio petrográfico, ensayo acelerado de barras de mortero y Gel-Pat modificado. En la investigación se ha desarrollado una metodología de tinción de geles álcali-sílice para la cuantificación de los mismos en el interior de las barras de mortero. Se ha relacionado el volumen de gel generado en la reacción con la expansión producida, tanto para áridos lentos, como rápidos, estableciendo analogías y diferencias en el comportamiento de ambos. La cuantificación de este tipo de compuestos realizada en las barras de mortero, abre la posibilidad de estudiar la capacidad expansiva que presentan en testigos de hormigón de obras afectadas. Este dato podría ser una herramienta importante para el desarrollo de modelos matemáticos que puedan predecir el futuro comportamiento de las estructuras afectadas por la reacción álcali-sílice. La tinción ha sido asimismo utilizada para establecer un método de detección de áridos reactivos basado en el ensayo de Gel-Pat modificado. La metodología de detección propuesta para estos áridos es cuantitativa y proporciona los resultados en 14 días, suponiendo una ventaja importante sobre los métodos actuales, que permiten detectarlos en 90 días (de forma orientativa en 56 días). This research is part of the studies on durability of concrete dams carried out in the Laboratorio Central de Estructuras y Materiales of CEDEX during the last years. At the present time, 18 public works affected by alkali-silica reaction have been diagnosed in Spain. In 12 of them the concrete was mixed with granitic aggregates, 7 of which were damaged by the slow alkali-silica reaction, 3 by the rapid alkali-silica reaction and in the other 2 cases, although there was no visual evidence of the reaction in field, the laboratory tests showed the presence of expansive products inside the concrete. However, in the international classifications of potentially reactive aggregates granitic rocks are pointed out as innocuous or as low reactivity and usually their reactivity is attributed to the presence of microcrystalline quartz and occasionally with the strained and microcracked quartz. The rapid reactive aggregates were deeply studied in a previous research (Víctor Daniel Lanza’s Thesis) also carried out in CEDEX. In this research, a detailed description of the components involved in the rapid alkali-silica reaction was given. Also, a fast method to detect them in the laboratory was developed. The research of the present PhD Thesis is focused in the slow-reactive aggregates, much less studied. The state of the art has shown some gaps in the knowledge about this reaction, mainly the absence of a reliable method to detect slow-reactive aggregates. On one side, there is no systematic and quantitative petrographic test for these aggregates. On the other side, the accelerated mortar bar test has been proved to be effective to detect rapid reactive aggregates but the standard limits fail in the detection of the slowreactive aggregates. In recent investigations it has been proposed extending the test until 90 days, but this period is considered too long. In this research, aggregates taken from affected structures have been used in the test, so their reactivity has been shown in practice. With the aim of increasing the number of samples aggregates from quarries have been also used. A method to stain alkali-silica gels has been applied to mortar bars, thereby quantifying the volume of gel generated inside. This volume has been related to the expansion of the bars, both for rapid for slow reactive aggregates, establishing similarities and differences. The quantification of the gel volume, applied to concrete cores extracted from affected structures, could be an important tool in mathematical models to predict the future behaviour of the structures affected by the alkali-silica reaction. The staining method has been also used to develop a method to detect slow reactive aggregates, based on the optimised Gel-Pat test, obtaining results in 14 days. This represents a major improvement compared to the actual methods (accelerated mortar bar test) which give an indicative value of reactivity at 56 days, and a final result at 90 days.

Aplicación de método de los elementos finitos a los problemas de amplificación dinámica de suelos

En el análisis de estructuras situadas en un emplazamiento de riesgo sísmico, es necesario conocer las característiias de un posible movimiento sísmico actuante, en la cota de cimentación de la estructura. La importancia del papel que juega el suelo circundante a la estructura en la modificación de las características del movimientos sísmico (principalmente en cuanto a amplitudes y contenido frecuencial), es bien conocida, denominándose amplificación sísmica a este fenómeno. En el caso de estructuras importantes, como Centrales Nucleares, presas, etc., o terrenos muy blandos, la presencia de la propia estructura también afecta de forma significativa al movimiento, dando lu gar al fenómeno de la interacción suelo-estructura. En esta comunicación se tratará únicamente del problema de la amplificación de las ondas sísmicas causantes del movimiento del suelo, y del análisis de la respuesta del mismo en ausencia de la estructura. Este movimiento del suelo denominado "movimiento del campo libre"(free field motion) puede, bien ser utilizado directamente en el análisis de la estructura, despreciando los efectos de la interacción, o emplearse como "input" en ciertas técnicas de análisis de interacción suelo-estructura.

Modelización de la rotura y fluidificación en geomateriales y geoestructuras

La rotura de las geoestructuras puede ser causada por cambios en las tensiones efectivas debidos bien por cargas externas (terremotos, por ejemplo), bien por variación de las presiones intersticiales (lluvia), o cambios en la geometría (erosión), así como por una disminución de las propiedades resistentes de los materiales (meteorización, ataque químico, etc). El caso particular de los deslizamientos es interesante, existiendo diversas clasificaciones que tienen en cuenta su forma, velocidad de propagación, etc. Dos de estos casos son los deslizamientos propiamente dichos y los flujos. En el primer caso, la deformación se concentra en zonas de pequeño espesor, que se idealiza como una superficie (superficie de rotura). La cinemática de esta rotura se puede considerar como el movimiento relativo de dos masas cuyas deformaciones no son grandes. Este mecanismo está usualmente asociado a materiales sobreconsolidados que presentan reblandecimiento. Los flujos se producen generalmente en materiales de baja densidad y estructura metaestable, con tendencia a compactar, de forma que se generan presiones intersticiales que aumentan el ángulo de rozamiento movilizado, pudiéndose llegar en algunos casos a la licuefacción. Este mecanismo de rotura se conoce como rotura difusa, y no ha sido tan estudiado como el de localización a pesar de se trata frecuentemente de roturas de tipo catastrófico. De hecho, el suelo pasa de un estado sólido a un estado fluidificado, con una gran movilidad. Es importante para el ingeniero predecir tanto el comportamiento de las geoestructuras bajo las cargas de cálculo como las condiciones en las que se producirá la rotura. De esta manera, en algunos casos, se podrán reforzar las zonas más débiles aumentando así su seguridad. En otros casos, no se podrá realizar este refuerzo (grandes deslizamientos como avalanchas, lahares, etc), pero sí se podrán conocer las consecuencias de la rotura: velocidad de propagación, alcance, espesores, etc. La modelización de estos problemas es compleja, ya que aparecen dificultades en los modelos matemáticos, constitutivos o reológicos y numéricos. Dado que en los geomateriales aparece una interacción fuerte entre el esqueleto sólido y los fluidos intersticiales, esto debe ser tenido en cuenta en los modelos matemáticos. En este trabajo se describirán, pues, el desarrollo y aplicación de técnicas avanzadas de modelización; matemática, constitutiva/reológica y numérica. Se dedicará especial atención a los problemas de transición entre suelos y suelos fluidificados, que hoy en día se estudian en una gran mayoría de los casos con modelos diferentes. Así por ejemplo, se emplean modelos constitutivos para el comportamiento previo a la rotura, y reológicos para los materiales fluidificados. En lo que respecta a los modelos matemáticos, existen formulaciones nuevas en velocidades (o desplazamientos), tensiones, y presiones de aire y agua intersticial, de los que se pueden obtener modelos simplificados integrados en profundidad para deslizamientos rápidos. Respecto de los modelos constitutivos, es interesante la teoría de la Plasticidad Generalizada (modelo básico de Pastor-Zienkiewicz y extensiones a suelos no saturados). Se estudiará la extensión de estos modelos elastoplásticos a la viscoplasticidad (Perzyna y consistente), explorando la posibilidad de emplearlos tanto antes como después de la rotura. Finalmente, en lo que a modelos numéricos se refiere, se describirá la implementación de los modelos matemáticos y constitutivos desarrollados en (i) modelos clásicos de elementos finitos, como el GeHoMadrid desarrollado en los grupo investigador M2i al que pertenece el autor de este trabajo, (ii) Métodos de tipo Taylor Galerkin, y (iii) métodos sin malla como el SPH y el Material Point Model. Estos modelos se aplicarán, principalmente a (i) Licuefacción de estructuras cimentadas en el fondo marino (ii) presas de residuos mineros (iii) deslizamientos rápidos de laderas.

Evaluación de efectos de biocidas contenidos en recubrimientos “antifouling “(AF coatings) en ecosistemas marinos

De entre los diversos grupos de contaminantes que pueden ser dañinos para ecosistemas acuáticos, sobresalen en los últimos años los biocidas utilizados como principio activo en recubrimientos antifouling ó pinturas anti-incrustantes para cascos de barcos y todo tipo de equipamiento sumergido ó en contacto con agua. Estos recubrimientos se aplican como sistema de protección para combatir la formación y asentamiento de comunidades bioincrustantes (fouling) frente a superficies expuestas al agua, también tienen como finalidad proteger frente a la corrosión de tipo químico y biológico. Normalmente estas pinturas anti-incrustantes son aplicadas en embarcaciones comerciales y de recreo, plataformas petrolíferas, tuberías submarinas, compuertas de presas, instalaciones destinadas a acuicultura, entre otro equipamiento. La utilización de biocidas en la formulación de pinturas anti-incrustantes para barcos ha sido propuesta por muchos investigadores como la aplicación de más impacto para los ecosistemas marinos, debido al intenso tráfico marítimo mundial que provoca la difusión de biocidas contaminantes en los mares, sobre todo en zonas costeras, bahías y puertos que es donde se magnifica el problema. Los llamados recubrimientos antifouling de segunda generación fueron los primeros que emplearon el tributilestaño (TBT) como principio activo en su formulación, a pesar de su gran eficacia y amplia utilización, al cabo del tiempo se ha demostrado su alta toxicidad y persistencia; por lo que existe una gran actividad investigadora en la búsqueda de alternativas asimismo eficientes pero más respetuosas con el medio marino. En este trabajo se pretende comparar los efectos para ecosistemas marinos del TBT y cinco biocidas como son dibutilestaño dicloruro, diuron (diclorofenil dimetil urea), piritionato de Zn ó Zn omadine, óxido de cobre (I) y DCOIT (Dicloro-2-n-octil-4-isotiazol-3-ona). Estos biocidas se han seleccionado en función de su naturaleza química diferenciada, distintas solubilidades en agua, eficacia, ecotoxicidad, persistencia y bioacumulación fundamentalmente. Para proceder a la clasificación del riesgo para ecosistema marino de los biocidas mencionados, nos valdremos de dos metodologías, una será la evaluación de unos índices de riesgo de biocidas para ecosistemas acuáticos con el fin de realizar una clasificación prospectiva de los mismos, basada en criterios PBT (Persistencia, Bioacumulación y Toxicidad para organismos acuáticos) y otra será una evaluación de la razón de la EXPOSICIÓN, valorada como PEC (Predictive Environmental Concentration) con relación a los EFECTOS originados, valorados como PNEC (Predicted No-Effect Concentration). Para cuantificar PEC nos valdremos de modelizaciones para ecosistemas marinos como MAMPEC (Marine Antifoulant Model to Predict Environmental Concentrations). Para cuantificar PNEC se hace uso de bioensayos a corto o largo plazo en organismos acuáticos. De esta manera podremos agrupar los biocidas en diversas categorías de riesgo y así poder decidir cuando su impacto medioambiental es asumible ó no asumible, alternativas posibles y en todo caso que decisiones se deben tomar al respecto.

Flood frequency analysis by a bivariate model based on copulas

La adecuada estimación de avenidas de diseño asociadas a altos periodos de retorno es necesaria para el diseño y gestión de estructuras hidráulicas como presas. En la práctica, la estimación de estos cuantiles se realiza normalmente a través de análisis de frecuencia univariados, basados en su mayoría en el estudio de caudales punta. Sin embargo, la naturaleza de las avenidas es multivariada, siendo esencial tener en cuenta características representativas de las avenidas, tales como caudal punta, volumen y duración del hidrograma, con el fin de llevar a cabo un análisis apropiado; especialmente cuando el caudal de entrada se transforma en un caudal de salida diferente durante el proceso de laminación en un embalse o llanura de inundación. Los análisis de frecuencia de avenidas multivariados han sido tradicionalmente llevados a cabo mediante el uso de distribuciones bivariadas estándar con el fin de modelar variables correlacionadas. Sin embargo, su uso conlleva limitaciones como la necesidad de usar el mismo tipo de distribuciones marginales para todas las variables y la existencia de una relación de dependencia lineal entre ellas. Recientemente, el uso de cópulas se ha extendido en hidrología debido a sus beneficios en relación al contexto multivariado, permitiendo superar los inconvenientes de las técnicas tradicionales. Una copula es una función que representa la estructura de dependencia de las variables de estudio, y permite obtener la distribución de frecuencia multivariada de dichas variables mediante sus distribuciones marginales, sin importar el tipo de distribución marginal utilizada. La estimación de periodos de retorno multivariados, y por lo tanto, de cuantiles multivariados, también se facilita debido a la manera en la que las cópulas están formuladas. La presente tesis doctoral busca proporcionar metodologías que mejoren las técnicas tradicionales usadas por profesionales para estimar cuantiles de avenida más adecuados para el diseño y la gestión de presas, así como para la evaluación del riesgo de avenida, mediante análisis de frecuencia de avenidas bivariados basados en cópulas. Las variables consideradas para ello son el caudal punta y el volumen del hidrograma. Con el objetivo de llevar a cabo un estudio completo, la presente investigación abarca: (i) el análisis de frecuencia de avenidas local bivariado centrado en examinar y comparar los periodos de retorno teóricos basados en la probabilidad natural de ocurrencia de una avenida, con el periodo de retorno asociado al riesgo de sobrevertido de la presa bajo análisis, con el fin de proporcionar cuantiles en una estación de aforo determinada; (ii) la extensión del enfoque local al regional, proporcionando un procedimiento completo para llevar a cabo un análisis de frecuencia de avenidas regional bivariado para proporcionar cuantiles en estaciones sin aforar o para mejorar la estimación de dichos cuantiles en estaciones aforadas; (iii) el uso de cópulas para investigar tendencias bivariadas en avenidas debido al aumento de los niveles de urbanización en una cuenca; y (iv) la extensión de series de avenida observadas mediante la combinación de los beneficios de un modelo basado en cópulas y de un modelo hidrometeorológico. Accurate design flood estimates associated with high return periods are necessary to design and manage hydraulic structures such as dams. In practice, the estimate of such quantiles is usually done via univariate flood frequency analyses, mostly based on the study of peak flows. Nevertheless, the nature of floods is multivariate, being essential to consider representative flood characteristics, such as flood peak, hydrograph volume and hydrograph duration to carry out an appropriate analysis; especially when the inflow peak is transformed into a different outflow peak during the routing process in a reservoir or floodplain. Multivariate flood frequency analyses have been traditionally performed by using standard bivariate distributions to model correlated variables, yet they entail some shortcomings such as the need of using the same kind of marginal distribution for all variables and the assumption of a linear dependence relation between them. Recently, the use of copulas has been extended in hydrology because of their benefits regarding dealing with the multivariate context, as they overcome the drawbacks of the traditional approach. A copula is a function that represents the dependence structure of the studied variables, and allows obtaining the multivariate frequency distribution of them by using their marginal distributions, regardless of the kind of marginal distributions considered. The estimate of multivariate return periods, and therefore multivariate quantiles, is also facilitated by the way in which copulas are formulated. The present doctoral thesis seeks to provide methodologies that improve traditional techniques used by practitioners, in order to estimate more appropriate flood quantiles for dam design, dam management and flood risk assessment, through bivariate flood frequency analyses based on the copula approach. The flood variables considered for that goal are peak flow and hydrograph volume. In order to accomplish a complete study, the present research addresses: (i) a bivariate local flood frequency analysis focused on examining and comparing theoretical return periods based on the natural probability of occurrence of a flood, with the return period associated with the risk of dam overtopping, to estimate quantiles at a given gauged site; (ii) the extension of the local to the regional approach, supplying a complete procedure for performing a bivariate regional flood frequency analysis to either estimate quantiles at ungauged sites or improve at-site estimates at gauged sites; (iii) the use of copulas to investigate bivariate flood trends due to increasing urbanisation levels in a catchment; and (iv) the extension of observed flood series by combining the benefits of a copula-based model and a hydro-meteorological model.

Evaluación de la seguridad al deslizamiento de una presa de fábrica en un cimiento con una familia de discontinuidades y con un criterio de rotura con ley de fluencia no asociada

En España existen del orden de 1,300 grandes presas, de las cuales un 20% fueron construidas antes de los años 60. El hecho de que existan actualmente una gran cantidad de presas antiguas aún en operación, ha producido un creciente interés en reevaluar su seguridad empleando herramientas nuevas o modificadas que incorporan modelos de fallo teóricos más completos, conceptos geotécnicos más complejos y nuevas técnicas de evaluación de la seguridad. Una manera muy común de abordar el análisis de estabilidad de presas de gravedad es, por ejemplo, considerar el deslizamiento a través de la interfase presa-cimiento empleando el criterio de rotura lineal de Mohr-Coulomb, en donde la cohesión y el ángulo de rozamiento son los parámetros que definen la resistencia al corte de la superficie de contacto. Sin embargo la influencia de aspectos como la presencia de planos de debilidad en el macizo rocoso de cimentación; la influencia de otros criterios de rotura para la junta y para el macizo rocoso (ej. el criterio de rotura de Hoek-Brown); las deformaciones volumétricas que ocurren durante la deformación plástica en el fallo del macizo rocoso (i.e., influencia de la dilatancia) no son usualmente consideradas durante el diseño original de la presa. En este contexto, en la presente tesis doctoral se propone una metodología analítica para el análisis de la estabilidad al deslizamiento de presas de hormigón, considerando un mecanismo de fallo en la cimentación caracterizado por la presencia de una familia de discontinuidades. En particular, se considera la posibilidad de que exista una junta sub-horizontal, preexistente y persistente en el macizo rocoso de la cimentación, con una superficie potencial de fallo que se extiende a través del macizo rocoso. El coeficiente de seguridad es entonces estimado usando una combinación de las resistencias a lo largo de los planos de rotura, cuyas resistencias son evaluadas empleando los criterios de rotura no lineales de Barton y Choubey (1977) y Barton y Bandis (1990), a lo largo del plano de deslizamiento de la junta; y el criterio de rotura de Hoek y Brown (1980) en su versión generalizada (Hoek et al. 2002), a lo largo del macizo rocoso. La metodología propuesta también considera la influencia del comportamiento del macizo rocoso cuando este sigue una ley de flujo no asociada con ángulo de dilatancia constante (Hoek y Brown 1997). La nueva metodología analítica propuesta es usada para evaluar las condiciones de estabilidad empleando dos modelos: un modelo determinista y un modelo probabilista, cuyos resultados son el valor del coeficiente de seguridad y la probabilidad de fallo al deslizamiento, respectivamente. El modelo determinista, implementado en MATLAB, es validado usando soluciones numéricas calculadas mediante el método de las diferencias finitas, empleando el código FLAC 6.0. El modelo propuesto proporciona resultados que son bastante similares a aquellos calculados con FLAC; sin embargo, los costos computacionales de la formulación propuesta son significativamente menores, facilitando el análisis de sensibilidad de la influencia de los diferentes parámetros de entrada sobre la seguridad de la presa, de cuyos resultados se obtienen los parámetros que más peso tienen en la estabilidad al deslizamiento de la estructura, manifestándose además la influencia de la ley de flujo en la rotura del macizo rocoso. La probabilidad de fallo es obtenida empleando el método de fiabilidad de primer orden (First Order Reliability Method; FORM), y los resultados de FORM son posteriormente validados mediante simulaciones de Monte Carlo. Los resultados obtenidos mediante ambas metodologías demuestran que, para el caso no asociado, los valores de probabilidad de fallo se ajustan de manera satisfactoria a los obtenidos mediante las simulaciones de Monte Carlo. Los resultados del caso asociado no son tan buenos, ya que producen resultados con errores del 0.7% al 66%, en los que no obstante se obtiene una buena concordancia cuando los casos se encuentran en, o cerca de, la situación de equilibrio límite. La eficiencia computacional es la principal ventaja que ofrece el método FORM para el análisis de la estabilidad de presas de hormigón, a diferencia de las simulaciones de Monte Carlo (que requiere de al menos 4 horas por cada ejecución) FORM requiere tan solo de 1 a 3 minutos en cada ejecución. There are 1,300 large dams in Spain, 20% of which were built before 1960. The fact that there are still many old dams in operation has produced an interest of reevaluate their safety using new or updated tools that incorporate state-of-the-art failure modes, geotechnical concepts and new safety assessment techniques. For instance, for gravity dams one common design approach considers the sliding through the dam-foundation interface, using a simple linear Mohr-Coulomb failure criterion with constant friction angle and cohesion parameters. But the influence of aspects such as the persistence of joint sets in the rock mass below the dam foundation; of the influence of others failure criteria proposed for rock joint and rock masses (e.g. the Hoek-Brown criterion); or the volumetric strains that occur during plastic failure of rock masses (i.e., the influence of dilatancy) are often no considered during the original dam design. In this context, an analytical methodology is proposed herein to assess the sliding stability of concrete dams, considering an extended failure mechanism in its rock foundation, which is characterized by the presence of an inclined, and impersistent joint set. In particular, the possibility of a preexisting sub-horizontal and impersistent joint set is considered, with a potential failure surface that could extend through the rock mass; the safety factor is therefore computed using a combination of strength along the rock joint (using the nonlinear Barton and Choubey (1977) and Barton and Bandis (1990) failure criteria) and along the rock mass (using the nonlinear failure criterion of Hoek and Brown (1980) in its generalized expression from Hoek et al. (2002)). The proposed methodology also considers the influence of a non-associative flow rule that has been incorporated using a (constant) dilation angle (Hoek and Brown 1997). The newly proposed analytical methodology is used to assess the dam stability conditions, employing for this purpose the deterministic and probabilistic models, resulting in the sliding safety factor and the probability of failure respectively. The deterministic model, implemented in MATLAB, is validated using numerical solution computed with the finite difference code FLAC 6.0. The proposed deterministic model provides results that are very similar to those computed with FLAC; however, since the new formulation can be implemented in a spreadsheet, the computational cost of the proposed model is significantly smaller, hence allowing to more easily conduct parametric analyses of the influence of the different input parameters on the dam’s safety. Once the model is validated, parametric analyses are conducting using the main parameters that describe the dam’s foundation. From this study, the impact of the more influential parameters on the sliding stability analysis is obtained and the error of considering the flow rule is assessed. The probability of failure is obtained employing the First Order Reliability Method (FORM). The probabilistic model is then validated using the Monte Carlo simulation method. Results obtained using both methodologies show good agreement for cases in which the rock mass has a nonassociate flow rule. For cases with an associated flow rule errors between 0.70% and 66% are obtained, so that the better adjustments are obtained for cases with, or close to, limit equilibrium conditions. The main advantage of FORM on sliding stability analyses of gravity dams is its computational efficiency, so that Monte Carlo simulations require at least 4 hours on each execution, whereas FORM requires only 1 to 3 minutes on each execution.

Modelos hidráulicos bidimensionales de lluvia directa aplicados al cálculo de eventos hidrológicos

Una de las principales preocupaciones a la hora de realizar cualquier proyecto de ingeniería es la adecuada determinación de los hidrogramas de avenida de los cursos de agua que lo afecten. Hoy en día es posible la simulación de modelos hidrológicos complejos de lluvia directa sobre malla 2D de elementos finitos, empleando los avances en cálculo vectorial que ofrece la tecnología CUDA (Arquitectura Unificada de Dispositivos de Cálculo), sin la necesidad de determinar subcuencas y tránsitos (como en HEC-HMS). Con los datos que ofrecen los Modelos Digitales del Terreno GRID del PNOA español (Plan Nacional de Ortofotografía Aérea) con precisión de 5 x 5 m y las Coberturas del Terreno SIG CORINE (Coordinación de Información del Medio Ambiente), que permite la evaluación de la rugosidad del suelo, se tiene la información necesaria para construir fácilmente este tipo de modelos. En este documento se presenta un procedimiento para delimitar fácilmente los cursos de agua principales y disponer en estas zonas un tamaño más pequeño de los elementos de la malla. El método propuesto permite una mejor definición de estas áreas después de unas pocas iteraciones, acelerando por tanto el proceso de construcción del modelo sin comprometer la calidad de los resultados. Para demostrar el método, se modelizaron dos cuencas completas mediante lluvia directa sobre una malla de elementos 2D. Se ensayaron diferentes niveles de precisión (variando el número de elementos) y se incluyeron dos presas, usando los datos del Modelo Digital del Terreno del PNOA. También se completó el cálculo hidrológico de ambas mediante HEC-HMS. En primer lugar se estudió la cuenca del río Zapardiel (cuenca del Duero, con superficie total de 1.450 km2, no presenta regulación artificial, pero sí algunas áreas de almacenamiento naturales). La segunda cuenca ensayada correspondió al río Zújar (cuenca del Guadiana, con superficie total de 8.500 km2, que incluye dos grandes presas: La Serena, con una capacidad de 3.200 hm3 y el Embalse del Zújar, situado aguas abajo de la primera, con un volumen útil de 300 hm3). Los resultados confirmaron que la forma de modelizar clásica con enfoque de subcuencas y tránsitos (tipo HEC-HMS) puede ser reemplazada por simulaciones de lluvia directa sobre malla 2D, por lo tanto reduciendo la incertidumbre inherente en la estimación de tiempos de concentración de subcuencas y parámetros de tránsito. Gracias a la disponibilidad de datos públicos de GRID PNOA y cobertura del terreno CORINE, actualmente es posible construir modelos de mallas de elementos 2D para simular lluvia directa sobre cuenca sin necesidad de costosos levantamientos topográficos. Con el hardware disponible hoy en día en ordenadores de uso general, se puede establecer un límite superior razonable para mallas 2D de 4 millones de elementos, lo que permite la correcta simulación de la lluvia directa en cuencas de hasta 10.000 km2. Para cuencas más grandes, este proceso se podría aplicar de forma repetida a varias zonas más pequeñas o en sucesivas secciones de la cuenca.

Los molinos del tramo medio del río Azuer. Ciudad Real.

Se trata de un estudio sobre los molinos hidráulicos del río Azuer, afluente del margen izquierdo (aguas abajo) del Guadiana, que discurre entre el Campo de Montiel y la Mancha Baja, hasta las Tablas de Daimiel, en su tramo medio, con un total de 33 molinos en todo el cauce y 25 en el tramo estudiado. Debido al espacio limitado de la comunicación solamente desarrollaré un estudio grafico, con planimetrías descriptivas y fotográfico, de tres molinos, uno restaurado, El Rezuelo de Membrilla, y dos en un estado aceptable, los Moros en San Carlos del Valle y El Grande de Manzanares. Concretaré el estudio de los edificios, desde un nivel arquitectónico, en cuanto a distribución, organización de las plantas, tamaño, programa y técnicas constructivas, además de describir los elementos propios de un molino, presas, ladrones, cárcavos, socaz, etc. Por último comentaré las diversas posibilidades de reutilización y de restauración

Revisión de la norma sismorresistente PDS-1. Sugerencias de algunas posibles modificaciones

Las mayores dificultades que aparecen en una futura reedición de la Norma Sismorresistente son probablemente las derivadas de una carencia casi absoluta de datos básicos recogidos en el territorio nacional, y en particular los referentes a registros de acelerogramas acaecidos en los últimos años. Por otra parte, si bien la información histórica de los seísmos descritos en términos cualitativos - intensidades- es suficientemente extensa, no existe hasta la fecha un mapa adecuado de riesgo sísmico que comprenda todo el país. La incertidumbre representada por las anteriores limitaciones debe de ser subsanada en la redacción de una nueva norma, mediante un enérgico programa de trabajo, que incluya como objetivo prioritario, la aportación de los anteriores datos básicos referidos a España. Esta nota hace referencia, de un modo específico, a los criterios que deberían ser considerados en la modificación del espectro básico de respuesta al tener en cuenta la influencia del terreno bajo la cimentación (amplificación del suelo). Aquí no se incluyen los efectos del tipo de cimentación, que deberían ser considerados de un modo separado así como la interacción dinámica suelo-estructura, que constituye una situación más compleja y cuyo estudio exige al menos el conocimiento de la amplificación del suelo y que se suele llevar a cabo en el caso de estructuras importantes (centrales nucleares, presas de tierra, etc.) que no deben ser objetivo prioritario de una norma sismorresistente de aplicación más general a edificación normal

Contribución al conocimiento de los primeros geotécnicos del siglo XX

Leyendo distintos artículos en la Revista de Obras Públicas (Jiménez Salas, 1945) uno recuerda a las grandes figuras como Coulomb (1773), Poncelet (1840), Rankine (1856), Culmann (1866), Mohr (1871), Boussinesq (1876) y otros muchos, que construyeron la base de un conocimiento que poco a poco irían facilitando la complicada tarea que suponía la construcción. Pero sus avances eran aproximaciones que presentaban notables diferencias frente al comportamiento de la naturaleza. Esas discrepancias con la naturaleza llegó un momento que se hicieron demasiado patentes. Importantes asientos en la construcción de los modernos edificios, rotura de presas de materiales sueltos y grandes corrimientos de tierras, por ejemplo durante la construcción del canal de Panamá, llevaron a la Sociedad Americana de Ingenieros Civiles (ASCE) a crear un comité que analizase las prácticas de la construcción de la época. Hechos similares se producían en Europa, por ejemplo en desmontes para ferrocarriles, que en el caso de Suecia supusieron unas cuantiosas perdidas materiales y humanas. El ingeniero austriaco-americano Karl Terzaghi (1883) había podido comprobar, en su práctica profesional, la carencia de conocimientos para afrontar muchos de los retos que la naturaleza ofrecía. Inicialmente buscó la respuesta en la geología pero encontró que ésta carecía de la definición necesaria para la práctica de la ingeniería, por lo que se lanzó a una denodada tarea investigadora basada en el método experimental. Comenzó en 1917 con escasos medios, pero pronto llegó a desarrollar algunos ensayos que le permitieron establecer los primeros conceptos de una nueva ciencia, la Mecánica de Suelos. Ciencia que ve la luz en 1925 con la publicación de su libro Erdbaumechanik auf bodenphysikalischer Grundlage. Rápidamente otras figuras empezaron a hacer sus contribuciones científicas y de divulgación, como es el caso del ingeniero austriaco-americano Arthur Casagrande (1902), cuya iniciativa de organizar el primer Congreso Internacional de Mecánica de Suelos e Ingeniería de Cimentaciones proporcionó el altavoz que necesitaba esa nueva ciencia para su difusión. Al mismo tiempo, más figuras internacionales se fueron uniendo a este período de grandes avances e innovadores puntos de vista. Figuras como Alec Skempton (1914) en el Reino Unido, Ralph Peck (1912) en los Estados Unidos o Laurits Bjerrum (1918) en Noruega sobresalieron entre los grandes de la época. Esta tesis investiga las vidas de estos geotécnicos, artífices de múltiples avances científicos de la nueva ciencia denominada Mecánica de Suelos. Todas estas grandes figuras de la geotecnia fueron presidentes, en distintos periodos, de la Sociedad Internacional de Mecánica de Suelos e Ingeniería de Cimentaciones. Se deja constancia de ello en las biografías que han sido elaboradas a partir de fuentes de variada procedencia y de los datos cruzados encontrados sobre estos extraordinarios geotécnicos. Así, las biografías de Terzaghi, Casagrande, Skempton, Peck y Bjerrum contribuyen no solo a su conocimiento individual sino que constituyen conjuntamente un punto de vista privilegiado para la comprensión de los acontecimientos vividos por la Mecánica de Suelos en el segundo tercio del siglo XX, extendiéndose en algunos casos hasta los albores del siglo XXI. Las aportaciones científicas de estos geotécnicos encuentran también su lugar en la parte técnica de esta tesis, en la que sus contribuciones individuales iniciales que configuran los distintos capítulos conservan sus puntos de vista originales, lo que permite tener una visión de los principios de la Mecánica de Suelos desde su mismo origen. On reading several articles in the journal, Revista de Obras Públicas (Jiménez Salas, 1945), one recalls such leading figures as Coulomb (1773), Poncelet (1840), Rankine (1856), Culmann (1866), Mohr (1871) and Boussinesq (1876) among many others, who created the basis of scientific knowledge that would make the complicated task of construction progressively easier. However, their advances were approximations which suffered considerable discrepancies when faced with the behaviour of the forces of nature. There came a time when such discrepancies became all too evident. Substantial soil settlements when constructing modern buildings, embankment dam failures and grave landslides, during the construction of the Panama Canal for example, led the American Society of Civil Engineers (ASCE) to form a committee in order to analyse construction practices of the time. Similar incidents had taken place in Europe, for example with railway slides, which in the case of Sweden, had resulted in heavy losses in both materials and human lives. During the practice of his career, the Austrian-American engineer Karl Terzaghi (1883) had encountered the many challenges posed by the forces of nature and the lack of knowledge at his disposal with which to overcome them. Terzaghi first sought a solution in geology only to discover that this lacked the necessary accuracy for the practice of engineering. He therefore threw himself into tireless research based on the experimental method. He began in 1917 on limited means but soon managed to develop several tests, which would allow him to establish the basic fundamentals of a new science; Soil Mechanics, a science which first saw the light of day on the publication of Terzaghi’s book, Erdbaumechanik auf bodenphysikalischer Grundlage. Other figures were quick to make their own scientific contributions. Such was the case of Austrian-American engineer, Arthur Casagrande (1902), whose initiative to organize the first International Congress of Soil Mechanics and Foundation Engineering provided the springboard that this science needed. At the same time, other international figures were becoming involved in this period of great advances and innovative concepts. Figures including the likes of Alec Skempton (1914) in the United Kingdom, Ralph Peck (1912) in the United States, and Laurits Bjerrum (1918) in Norway stood out amongst the greatest of their time. This thesis investigates the lives of these geotechnical engineers to whom we are indebted for a great many scientific advances in this new science known as Soil Mechanics. Moreover, each of these eminent figures held the presidency of the International Society of Soil Mechanics and Foundation Engineering, record of which can be found in their biographies, drawn from diverse sources, and by crosschecking and referencing all the available information on these extraordinary geotechnical engineers. Thus, the biographies of Terzaghi, Casagrande, Skempton, Peck and Bjerrum not only serve to provide knowledge on the individual, but moreover, as a collective, they present us with an exceptional insight into the important developments which took place in Soil Mechanics in the second third of the 20th century, and indeed, in some cases, up to the dawn of the 21st. The scientific contributions of these geotechnical engineers also find their place in the technical part of this thesis in which the initial individual contributions which make up several chapters retain their original approaches allowing us a view of the principles of Soil Mechanics from its very beginnings.

De la utopía a la realidad: Casto Fernández-Shaw en la cuenca del Guadalquivir

El objetivo de esta ponencia es estudiar la aportación de Casto Fernández-Shaw a las importantes obras hidráulicas que se acometieron en Andalucía durante el siglo XX. Este arquitecto diseñó y construyó cuatro presas y centrales hidroeléctricas en la cuenca del río Guadalquivir y sus afluentes, entre los años 1920 y 1931, para la empresa Mengemor. Estudiaremos los casos de El Carpio, Jándula, Encinarejo y Alcalá del Río. Las examinaremos comparándolas además no sólo con el resto de su repertorio (como su proyecto de Monumento al Triunfo de la Civilización o sus arquitecturas futuristas), sino también con otras construcciones hidráulicas similares de otros arquitectos como Joaquín Vaquero Palacios (Presa de Grandas de Salime), Miguel Fisac (Hidroeléctrica de Canfranc), o Ignacio Álvarez Castelao (Conjuntos de Silvón y Arbón); sin olvidar tampoco otras importantes presas españolas como Aldeadávila, Belesar, Susqueda... Todas estos ejemplos mencionados tienen en común la imbricación en su conjunto de las artes estéticas, paisajísticas y funcionales propias de una obra arquitectónica pero aplicadas a una fábrica hidráulica, y es interesante estudiar cómo los arquitectos, acostumbrados a trabajar en enclaves urbanos, se enfrentan a una obra industrial enmarcada en un entorno natural y a una escala de paisaje.

Arquitecturas del Agua: las centrales hidroeléctricas de Antonio Palacios, entre el clasicismo y la modernidad

Desde los dibujos industriales futuristas de Sant Elia (1913-14) hasta la central hidroeléctrica de aires deconstructivistas de Becker Architekten (Kempten, Alemania; 2011), los arquitectos se han interesado continuamente por el diseño y construcción de presas y centrales. Las enormes paredes de hormigón o piedra que cierran los valles para embalsar el agua, que después caerá como un torrente a través de los aliviaderos constituyen, por su verticalidad, su potencia visual y sus evocaciones sonoras, uno de los elementos más interesantes a los que se tienen que enfrentar los técnicos a la hora de la construcción de este tipo de ingenios hidráulicos. Uno de los casos más interesantes de colaboración entre arquitectura e ingeniería es el de Antonio Palacios (1874-1945), relevante arquitecto de principios del s.XX, responsable de emblemáticos edificios de la capital (Círculo de Bellas Artes, Palacio de Correos...), y que realizó dos centrales hidroeléctricas, en Mengíbar (Jaén), y en Tambre (A Coruña). En esta comunicación analizaremos estos dos interesantes ejemplos de arquitectura hidroeléctrica comparándolos no sólo entre sí sino también con el resto de la obra industrial y civil de su autor.