Simulación de la calidad del aire en la isla de Gran Canaria mediante el método de los elementos finitos y su validación con datos experimentales

[EN]Este estudio muestra una simulación 3D de la calidad del aire en la isla de Gran Canaria. Se simula el transporte difusión y reacción de los gases emitidos en dos centrales térmicas de la isla. Se presenta una estrategia de modelización y simulación mediante elementos finitos. Los resultados se comparan con datos experimentales. Se usa el método del Meccano para mallar la isla de Gran Canaria, a partir de una triangulación de la topografia de la isla. A continuación, con los datos de las estaciones de viento, se calcula un campo de viento para la zona de estudio…

Desarrollo de un método ensemble para la predicción del viento a escala local usando elementos finitos

[EN]Los autores han desarrollado un modelo de viento de masa consistente especialmente diseñado para su aplicación en la escala local y en zonas de orografía compleja. Se ha dotado a este modelo de carácter predictivo usando como entrada resultados de HARMONIE. El HARMONIE es un modelo meteorológico predictivo de escala regional usado en la AEMET. Por otra parte, en los últimos años los métodos ensemble se han consolidado en la predicción meteorológica a escala regional…

Adaptación del método Demit al contexto mexicano : Identificación de elementos a reformular para su aplicación eficaz

El método DEMIT (Desarrollo Energético por Modelización e Inteligencia Territorial, véase figura 1) fue implementado para dar respuesta a las solicitudes de los ciudadanos de Quebec, Canadá, de transparencia y de participación en la toma de decisión para la puesta en marcha de parques eólicos. DEMIT articula dos métodos multicriterio: la Ayuda a la Decisión Multicriterio (ADMC) y los Sistemas de Información Geográfica (SIG) colaborativos y participativos. Dichas herramientas multicriterio proveen a los tomadores de decisión de un conjunto de condiciones que permitirán una toma de decisión transparente y participativa durante la implantación de un parque eólico. Este método fue objeto de un caso de estudio, entre 2010 y 2011, de un parque eólico en operación en Quebec. En 2012, gracias a un financiamiento para un proyecto bilateral Quebec-México, se realizó una primera y somera adaptación al contexto mexicano, en donde se identificaron algunos elementos que deben reformularse para su aplicación eficaz en México. Dichos elementos se presentan en esta comunicación

Adaptación del método Demit al contexto mexicano : Identificación de elementos a reformular para su aplicación eficaz

El método DEMIT (Desarrollo Energético por Modelización e Inteligencia Territorial, véase figura 1) fue implementado para dar respuesta a las solicitudes de los ciudadanos de Quebec, Canadá, de transparencia y de participación en la toma de decisión para la puesta en marcha de parques eólicos. DEMIT articula dos métodos multicriterio: la Ayuda a la Decisión Multicriterio (ADMC) y los Sistemas de Información Geográfica (SIG) colaborativos y participativos. Dichas herramientas multicriterio proveen a los tomadores de decisión de un conjunto de condiciones que permitirán una toma de decisión transparente y participativa durante la implantación de un parque eólico. Este método fue objeto de un caso de estudio, entre 2010 y 2011, de un parque eólico en operación en Quebec. En 2012, gracias a un financiamiento para un proyecto bilateral Quebec-México, se realizó una primera y somera adaptación al contexto mexicano, en donde se identificaron algunos elementos que deben reformularse para su aplicación eficaz en México. Dichos elementos se presentan en esta comunicación

Adaptación del método Demit al contexto mexicano : Identificación de elementos a reformular para su aplicación eficaz

El método DEMIT (Desarrollo Energético por Modelización e Inteligencia Territorial, véase figura 1) fue implementado para dar respuesta a las solicitudes de los ciudadanos de Quebec, Canadá, de transparencia y de participación en la toma de decisión para la puesta en marcha de parques eólicos. DEMIT articula dos métodos multicriterio: la Ayuda a la Decisión Multicriterio (ADMC) y los Sistemas de Información Geográfica (SIG) colaborativos y participativos. Dichas herramientas multicriterio proveen a los tomadores de decisión de un conjunto de condiciones que permitirán una toma de decisión transparente y participativa durante la implantación de un parque eólico. Este método fue objeto de un caso de estudio, entre 2010 y 2011, de un parque eólico en operación en Quebec. En 2012, gracias a un financiamiento para un proyecto bilateral Quebec-México, se realizó una primera y somera adaptación al contexto mexicano, en donde se identificaron algunos elementos que deben reformularse para su aplicación eficaz en México. Dichos elementos se presentan en esta comunicación

Modelación numérica de deslizamientos de ladera en embalses mediante el método de partículas y elementos finitos (PFEM)

Se presentan resultados obtenidos mediante la aplicación del MétododePartículasyElementosFinitos (PFEM) en la simulación dedeslizamientosdeladeraenembalses. Es un fenómeno complejo, por la interacción entre el material deslizado, la masa de agua del embalse, yel conjunto formado por el vaso y la presa. PFEM es un esquema numérico original con el que se ha afrontado con éxito la resolución de problemas de interacción fluido-estructura. Combina un enfoque Lagrangiano con la resolución de las ecuaciones deelementosfinitosmediante la generación de una malla, que se actualiza en cada paso de tiempo. Se presentan resultados de casos de validación en los que se han reproducido ensayos en laboratorio existentes en la bibliografía técnica. Se muestran también otros cálculos más complejos, sobre la cartografía a escala real de un embalse, donde se aprecia el fenómeno de generación de la ola, su propagación por el embalse y la afección a la presa. Por último, se ha modelado el deslizamiento ocurrido en 1958 en la bahía de Lituya (Alaska), enel que la caída de 90 millones de toneladas de roca produjo una ola que alcanzó una sobreelevación máxima de 524 m en la ladera opuesta. Los resultados permiten afirmar que PFEM puede ser una herramienta útil enel análisis de riesgos frente a este tipo de fenómenos, ofreciendo una buena aproximación de las afecciones potenciales

Determinación del coeficiente de rozamiento grano-pared en silos de pared corrugada mediante el método de los elementos discretos

El comportamiento de los materiales granulares almacenados en silos se ve afectado por varios parámetros, tanto aquellos característicos del material como de la geometría del silo. La determinación del coeficiente de rozamiento pared-partícula es uno de los parámetros de mayor importancia, siendo habituales en su determinación el uso de ensayos de corte directo. En el presente trabajo se estudia dicho coeficiente para el supuesto de una pared corrugada, el cual, teóricamente, debe representar a un valor efectivo que depende tanto del coeficiente de rozamiento grano-pared (para el caso de una pared lisa) y el ángulo de rozamiento interno del material. La determinación del rozamiento efectivo se ha realizado a través de la simulación por elementos discretos de un ensayo de corte sobre una pared corrugada. Los valores obtenidos han sido comparados con las prescripciones expuestas en la normativa vigente. La potencialidad del método de los elementos discretos permite el estudio de diversas configuraciones geométricas de la pared corrugada de los silos sin necesidad de realizar ensayos de laboratorio. Esto permitirá, en trabajos sucesivos, investigar la influencia de muy diversos factores en el valor del rozamiento efectivo grano-pared en este tipo de paredes. The mechanical behaviour of granular materials stored in silos is affected by numerous parameters, some of them being related to the characteristic of the stored materials and others to the geometry of the silo. The determination of the particle-wall friction coefficient (or wall friction) is of great importance and direct shear tests are usually conducted in order to obtain its value. In this work this variable is analysed for the case of a corrugated wall. This value is expected to be an effective value between the particle-wall friction coefficient (obtained for a flat wall) and the internal friction coefficient of the material under study. The effective wall friction determination has been carried out by using a discrete element model to simulate a direct shear test on a corrugated wall. The values obtained have been compared with prescriptions given in the current standards. The potential of the discrete element method allows different geometries of the silo corrugated wall to be considered without the necessity of developing laboratory tests. In future works it will allow the influence of numerous parameters on the effective wall friction in corrugated walls to be studied.

Análisis del estado de tensiones en uniones carpinteras de empalme de llave por el método de los elementos finitos

El objeto del presente trabajo es la determinación de la distribución de tensiones en la unión carpintera de empalme de llave mediante el método de los elementos finitos y su comparación con los valores obtenidos mediante la teoría de Resistencia de Materiales. Se analizan las zonas donde se produce concentración de tensiones y se estudia la influencia del refino de la malla sobre los resultados con el fin de determinar el tamaño que proporciona los valores de tensión más ajustados a la teoría. En las zonas donde la concentración de tensiones es menor, distintos tamaños de la malla muestran valores de tensión similares. Los resultados muestran una simetría central de la distribución de isotensiones donde el centro de simetría se corresponde con el centro geométrico de la unión. En las zonas donde se produce concentración de tensiones, los valores de las mismas aumentan considerablemente con el refino de la malla. La comparación de los niveles de tensiones normales obtenidos por el MEF y la teoría clásica muestran diferencias reducidas, excepto en los puntos de concentración de tensiones. The purpose of this study is to determine the stress distribution in the joint carpenter of halved and tabled joint with the finite element method and its comparison with the values obtained using the theory of Strength of Materials. The stress concentration areas where analyzed and the influence of mesh refinement was studied on the results in order to determine the mesh size that provides the stress values more consistent with the theory. In areas where stress concentration is lower, different mesh sizes show similar stress values. The results show a central symmetry of the isobar lines distribution where the centre of symmetry corresponds to the geometric centre of the joint. In areas where stress concentration occurs, the same values increase considerably with the refinement of the mesh. Comparison of normal stress levels obtained by the FEM and the classical theory shows small differences, except at points of stress concentration.

Propuesta de Configuración y de Método de Inspección de Uniones Mixtas Mediante Pernos Conectores

La utilización conjunta de los soportes de acero con los forjados de hormigón armado aúna las ventajas constructivas y arquitectónicas de ambos sistemas de ejecución. El principal inconveniente que se presenta es la complejidad de sus nudos. Las uniones estructurales viga-pilar mediante pernos tipo Nelson o Köco aparecen sólo ocasionalmente en el mundo de la edificación. Los casos en los que esto ocurre normalmente corresponden a obras de especial singularidad, a uniones secundarias o a anclajes de reducida responsabilidad estructural. La presente Tesis Doctoral tiene el objetivo principal de aportar una alternativa a las uniones convencionales entre losas o vigas de hormigón armado con pilares metálicos, utilizando para ello los pernos conectores. En este trabajo se establece un campo de aplicación concreto de las estructuras de edificación en el que la utilización de uniones mediante pernos conectores resulta óptima, dada la posibilidad de industrialización de estos nudos, su sencillez y la alta resistencia de los pernos en la transmisión del cortante. Se ha procedido al estudio teórico-analítico de la resistencia de estas uniones y a su posterior comparación experimental. Se presentan las diferencias obtenidas entre los resultados experimentales y los previstos por las formulaciones teóricas, debidas principalmente a las particularidades de las configuraciones de la unión y al efecto del rozamiento. El excelente comportamiento mecánico obtenido en los resultados experimentales induce a presagiar que la configuración de unión que se propone, junto con el método de cálculo aportado, puede llegar a sustituir a las soluciones convencionales utilizadas hasta el momento en las estructuras mixtas de edificación. Los pernos conectores aportan grandes ventajas al campo de la construcción aunque sus soldaduras, pese a la gran seguridad que las caracteriza, son difícilmente comprobables mediante ensayos no destructivos. Los ensayos destructivos no son en absoluto recomendables para las configuraciones de unión estudiadas en esta Tesis. Esto se debe principalmente al reducido número de pernos que configuran las uniones y a la gran responsabilidad que cada perno asume respecto a la seguridad global de la estructura. Los métodos convencionales de inspección presentan grandes dificultades e inconvenientes para comprobar las soldaduras de los pernos, por lo que se han investigado nuevas posibilidades de inspección y control. Esta Tesis plantea un método alternativo para la inspección de las soldaduras de los pernos conectores. El método desarrollado se basa en la acústica y consiste en el cálculo analítico de los modos propios de vibración del perno y su posterior comparación con las frecuencias naturales obtenidas a partir de una medición experimental. La investigación analítica comprende numerosas simulaciones por elementos finitos en las que se varían las condiciones de contorno para simular los defectos de soldadura. La investigación experimental consiste en la medición acústica de los pernos conectores y su posterior ensayo mediante métodos destructivos y no destructivos. Finalmente se han comparado los resultados analíticos con los obtenidos en el laboratorio y se ha obtenido una clara correlación entre las frecuencias calculadas y la calidad de las soldaduras. De este modo se puede averiguar el estado interno de una soldadura a partir de medición de su frecuencia natural de vibración. Se propone un método de inspección acústica que permite comprobar la calidad de las soldaduras de modo sencillo, económico y de fácil ejecución en obra. En conclusión, esta Tesis presenta una nueva configuración de unión mixta mediante pernos conectores y un completo método de calculo y comprobación. Finalmente se ha desarrollado un método no destructivo para la inspección de los pernos que permite garantizar la calidad de sus soldaduras en obra.

Notas sobre el método de los elementos finitos

Ideas básicas: reducción del medio contínuo a discreto (elementos, nodos y funciones de forma); definición de un estado "elemental" relacionado con otro "general", mediante el formalismo "matriz de rigidez" y vectores consistentes; método de síntesis de propiedades generales: método directo de la rigidez. Posibilidades de estudio: derivación de la teoría, modelización de problemas reales con F.E.M., desarrollo de programas de ordenador, estudio de métodos numéricos (sistemas ecuacionales, autovalores, convergencia, integración etc.). Contrariamente a lo que se diga el F.E.M. no elimina la necesidad de profundizar en el estudio de la elasticidad, matemáticas, etc. Se trata tan sólo de un instrumento potentísimo que, en principio, permite abordar cualquier problema.

Análisis del estado de tensiones en uniones carpinteras de empalme de llave por el método de los elementos finitos

El objeto del presente trabajo es la determinación de la distribución de tensiones en la unión carpintera de empalme de llave mediante el método de los elementos finitos y su comparación con los valores obtenidos mediante la teoría de Resistencia de Materiales. Se analizan las zonas donde se produce concentración de tensiones y se estudia la influencia del refino de la malla sobre los resultados con el fin de determinar el tamaño que proporciona los valores de tensión más ajustados a la teoría. En las zonas donde la concentración de tensiones es menor, distintos tamaños de la malla muestran valores de tensión similares. Los resultados muestran una simetría central de la distribución de isotensiones donde el centro de simetría se corresponde con el centro geométrico de la unión. En las zonas donde se produce concentración de tensiones, los valores de las mismas aumentan considerablemente con el refino de la malla. La comparación de los niveles de tensiones normales obtenidos por el MEF y la teoría clásica muestran diferencias reducidas, excepto en los puntos de concentración de tensiones.

Análisis experimental y por el método de los elementos finitos del estado de tensiones en uniones carpinteras de empalme de llave

In this research, the halved and tabled traditional timber scarf joint is analyzed. This joint consists in two end joint pieces usually subjected to tension. Initially, the study is discussed from an experimental point of view. In this way, 3 critical cross-sections are established (section of the notch, section of the horizontal plane and reduced section) and mechanical tests are performed to achieve the failure on each of critical sections by changing the geometry of the joint. The study is completed by developing a finite element model which allows verify experimental results and extend the analysis to other geometries. This model has to simulate the real behavior of the material which is being studied, so mechanical tests are performed to obtain the elastic constants and the coefficients of friction of the material. In the reduced section, an abrupt decrease of the effective cross-section takes place, and this effect is also experimentally analyzed. These tests indicate that a crack is initiated before the bending-tension failure occurs in the reduced section. The test material consists of wood of Pinus sylvestris L. coming from the “Valsaín´s Sawmill” (Segovia) with “premium quality” according to the nonstructural wood visual classification of sawmill. It is observed that initiation of a crack, in the mortise (bottom of reduced section), and shear stress concentration, at the initial part of the heel (beginning of horizontal plane), completely determine the mechanical behaviour of the joint, resulting in 3 failure modes: local compression failure in the section of the notch, shear failure in the horizontal plane, and failure of stresses concentration, mainly perpendicular to the grain tension, at the bottom of reduced section. The geometric optimization is obtained for halved and tabled traditional scarf joint, when the joint has made with similar properties of wood than tested specimens, for any height and width of the cross-section. It is considered the failure due to the initiation of a crack in reduced section, by applying a correction coefficient into the usual equation used to design the members subjected to both tension and bending. Therefore, it is possible to obtain, analytically, the design conditions to be met of the 3 critical cross-sections. According to the theoretical optimization, the tension strength of complete cross-section is reduced until 14%, when using this type of joint. The experimental optimization indicates even a greater reduction, until 6%. En el presente trabajo de investigación se analiza el comportamiento mecánico de las uniones tradicionales de empalme de llave, que consisten en dos piezas unidas por sus testas transmitiéndose entre ellas principalmente un esfuerzo de tracción. Inicialmente, el estudio se aborda desde un punto de vista experimental. De este modo, se establecen las 3 secciones críticas o de estudio (sección del encaje, sección rasante del cogote y sección reducida) y se realizan ensayos mecánicos, variando la geometría de la unión, para alcanzar la rotura en cada una de ellas. Se completa el estudio mediante la elaboración de un modelo por elementos finitos que permite verificar los resultados experimentales y ampliar el análisis a otras geometrías. Este modelo debe simular el comportamiento real del material objeto de estudio, por lo que se realizan ensayos para obtener las constantes elásticas y los coeficientes de rozamiento del mismo. También se analiza, experimentalmente, el efecto entalladura que reduce bruscamente la sección completa del tirante, estableciendo que el fallo por flexotracción en la sección reducida de la pieza, no llega a producirse por el inicio previo de una grieta. El material de ensayo consiste en madera de Pinus sylvestris L. (pino silvestre) procedente del Aserradero de Valsaín (Segovia) y de calidad “Extra” o “Primera” según la clasificación visual no estructural del aserradero. Se observa que el inicio de una grieta en la mortaja del rediente y la concentración de tensiones tangenciales en la parte inicial del cogote, determinan completamente el comportamiento mecánico de la unión, dando lugar a 3 modos distintos de rotura: fallo por compresión en la sección del encaje, fallo por cortante en la sección rasante y fallo por concentración de tensiones, principalmente tracciones perpendiculares, en el rebaje de la sección reducida. Se consigue optimizar geométricamente cualquier empalme de llave confeccionado con madera de características similares a la ensayada, para cualquier valor de la altura y de la anchura de la sección. Se considera el agotamiento en la sección reducida causado por el inicio de grieta, mediante la aplicación de un coeficiente corrector en la expresión habitual de agotamiento por flexotracción, en consecuencia, finalmente es posible obtener, de modo analítico, un valor del índice de agotamiento en cada una de las 3 secciones de estudio. La optimización teórica del empalme de llave indica que la capacidad resistente del tirante bruto se reduce al 14%, cuando se coloca este tipo de unión tradicional. Experimentalmente se obtiene, que, para la sección ensayada, la capacidad resistente del tirante bruto se reduce todavía más, llegando al 6%.

Análisis del proceso de fisuración en estructuras de hormigón por el Método de los Elementos Finitos mediante un modelo plástico

El Método de los Elementos Finitos ha demostrado ser una poderosa herramienta para el análisis de problemas complejos en Ingeniería Estructural, tales como la resolución de estructuras cuyos materiales son susceptibles de plastificación, por ejemplo. En cualquier caso, tal éxito sólo es posible cuando el modelo elegido de comportamiento corresponde a la realidad experimental; siendo éste uno de los problemas aún no resueltos para el hormigón por no existir modelos suficientemente generales que predigan su comportamiento dentro de las teorías de la Mecánica del Continuo. En este trabajo se presenta un modelo para el estudio de la fisuración en el hormigón dentro de un contexto matemático análogo al utilizado en la Teoría de la Plasticidad; de manera que los esquemas utilizados en el M.E.F. para la implementación de ésta siguen siendo válidos. El resultado es un método para el análisis de estructuras de hormigón que recoge el complejo fenómeno de la fisuración con relativa sencillez. También son presentados algunos ejemplos resueltos con un programa de ordenador que implementa dicho modelo dentro de las técnicas del M.E.F.

Implementación del método simplet en un programa de elementos finitos basados en volumenes de control

En este trabajo se presenta la implantación del método SIMPLET en un programa de elementos finitos basados en volúmenes de control. Este método toma en cuanta las variaciones de temperatura en la determinación del campo de presiones en problemas de convección libre, con el próposito de acelerar la convergencia y disminuir los tiempos de computación. El método SIMPLET, que originalmente ha sido propuesto en base al método de diferencias finitas con mallas desplazadas, se implementó en un programa de elementos finitos basados en volúmenes de control con mallas no desplazadas. Se resolvió un problema de convección libre en una cavidad cuadrada y los resultados obtenidos en términos del número de iteracciones y tiempo de computación se compararon se compararon con los resultados del método SIMPLE. Los resultados muestran que en este tipo de problemas el método SIMPLET es más rápido que le SIMPLE cuando el número de Rayleigh es bajo, mientras que para números de Rayleigh elevados, el desempeño de ambos métodos es similar.

Interpolación de Birkhoff a partir de una formulación variacional asociada al método de los elementos finitos

El problema de la interpolación polinómica de Birkhoff, como es conocido, no siempre admite solución. Se presenta un método como alternativa a este tipo de problemas para la obtención de interpolantes con unas condiciones de continuidad determinadas y en el cual el criterio de aproximación es el de minimización de un cierto funcional real, utilizando el método de los elementos finitos. Se analiza la eficiencia del método en diversos ejemplos de aplicación 1-0 y 2-D.