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

Modelos de elementos finitos con aplicación tridimensional de la analogía fluido-fibra, para predecir el comportamiento mecánico de piezas estructurales de madera bajo la presencia de los principales tipos de nudos

Los nudos son los defectos que más disminuyen la resistencia de piezas de madera en la escala estructural, al ocasionar no solo una discontinuidad material, sino también la desviación de las fibras que se encuentran a su alrededor. En la década de los 80 se introdujo la teoría de la analogía fluido-fibra, como un método que aproximaba adecuadamente todas estas desviaciones. Sin embargo en su aplicación tridimensional, nunca se consideraron las diferencias geométricas en el sentido perpendicular al eje longitudinal de las piezas estructurales, lo cual imposibilitaba la simulación numérica de algunos de los principales tipos de nudos, y disminuía la precisión obtenida en aquellos nudos en los que la modelización sí era viable. En este trabajo se propone un modelo programado en lenguaje paramétrico de un software de elementos finitos que, bajo una formulación en tres dimensiones más general, permitirá estudiar de forma automatizada el comportamiento estructural de la madera bajo la influencia de los principales tipos de nudos, a partir de la geometría visible de los mismos y la posición de la médula en la pieza, y el cual ha sido contrastado experimentalmente, simulando de forma muy precisa el comportamiento mecánico de vigas sometidas a ensayos de flexión a cuatro puntos. Knots are the defects that most reduce the strength of lumber at the structural level, by causing not only a material discontinuity but also the deviation of the fibers that surround them. In the 80's it was introduced the theory of the flow-grain analogy as a method to approximating adequately these deviations. However, in three-dimensional applications, geometrical differences in the direction perpendicular to the longitudinal axis of the structural specimens were never considered before, which prevented the numerical simulation of some of the main types of knots, and decreased the achieved precision in those kind of knots where modeling itself was possible. This paper purposes a parametric model programmed in a finite element software, in the way that with a more general three-dimensional formulation, an automated study of the structural behavior of timber under the influence of the main types of knots is allowed by only knowing the visible geometry of such defects, and the position of the pith. Furthermore that has been confirmed experimentally obtaining very accurately simulations of the mechanical behavior of beams under four points bending test.

Análisis estructural de algunas obras de Félix Candela mediante modelos de Elementos Finitos

Por su facilidad constructiva y la relativa sencillez de su análisis estructural, el paraboloide hiperbólico de hormigón armado es una forma profusamente utilizada por el arquitecto Félix Candela (1910-1997). En este artículo se analiza el comportamiento de varias de esas obras. Concretamente dos paralelogramos alabeados [1], el restaurante Los Manantiales (México) y la Iglesia de Guadalupe (Madrid). Se explica la geometría del paraboloide hiperbólico y se describe brevemente el análisis estructural presentado por Félix Candela. Los resultados presentados posteriormente confirman la distribución de esfuerzos descrita por Candela. En todos los casos estudiados las tensiones resultan ser admisibles para el hormigón. Se ha verificado también el comportamiento de membrana de las estructuras mediante un factor cuyo valor oscila entre 0 (flexión pura) y 1 (membrana pura). El comportamiento se aleja o acerca al de membrana pura en función de las condiciones de apoyo como ya recalcó el arquitecto madrileño

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.

Amplificación sísmica: una aproximación por Elementos Finitos

En las páginas que siguen se presenta el estudio de uno de los más típicos problemas de dinámica estructural, cual es la obtención de la respuesta de una estructura excitada por un movimiento de la base. Este es un caso muy frecuente en ingeniería sísmica, donde el objeto del estudio puede ser el edificio (sometido a un movimiento en la cimentación) o un estrato de terreno sobre fondo rígido. Al objeto de facilitar un soporte intuitivo a la exposicióri, ésta se organiza en base al segundo de los casos citados (estrate en base rígida). La aproximación escogida, elementos finitos, pone de relieve una vez más la potencia y generalidad del método en lo que respecta a la formulación del sistema de equilibrio. La discusión se centra en un aspecto concreto del método: la elección de funciones de forma.

Consolidación elastoplástica con deformaciones finitas. Implementación con elementos finitos y ejemplos numéricos

En un programa de elementos finitos se implementa un modelo matemático para la consolidación elastoplástica con deformaciones finitas en un medio representativo de un suelo totalmente saturado. El tratamiento algorítmico de la elasticidad en deformaciones finitas para la fase sólida está basado en una descomposición multiplicativa y acoplado con el algoritmo de flujo del fluido mediante la presión intersticial de Kirchho. Se utiliza una formulación mixta de elementos finitos con dos campos en que los desplazamientos nodales del sólido y las presiones nodales de agua en los poros están acoplados mediante las ecuaciones de equilibrio de masa y cantidad de movimiento. La ley de comportamiento de la fase sólida se representa mediante una teoría de tipo Cam-Clay modificada, formulada en el espacio de las tensiones principales de Kirchho, y se utiliza una aplicación de retorno que se lleva a cabo en el espacio de deformaciones definido por los invariantes de los alargamientos elásticos logarítmicos principales. El comportamiento de la fase fluida se representa mediante una ley de Darcy generalizada formulada respecto a la configuración actual. El modelo de elementos finitos es completamente linealizable con exactitud. Se presentan varios ejemplos numéricos con y sin efectos de deformaciones finitas para demostrar el impacto de la no linealidad geométrica en las correspondientes respuestas. El artículo finaliza con un estudio del comportamiento del modelo de elementos finitos en relación con la precisión y la estabilidad numérica.

Simulación de excavaciones mediante un método de acoplamiento elementos finitos – elementos de contorno

Cuando se modelan sistemas físicos no lineales de extensión infinita, como las excavaciones, se hace necesario simular adecuadamente tanto la solución en el infinito como la no linealidad. El método de elementos finitos es una herramienta efectiva para representar la no linealidad. Sin embargo, el tratamiento del campo infinito truncando el dominio es bastante cuestionable. Por otro lado, el método de elementos de contorno es adecuado para simular el comportamiento en el infinito sin truncamientos. Por combinación de ambos métodos, se puede obtener un uso adecuado de las ventajas de cada uno. En este trabajo se proponen diversas posibilidades de acoplamiento entre los dos métodos. Se desarrollan algoritmos de acoplamiento basados en una descomposición de dominios y se comparan con los esquemas más tradicionales de acoplamiento.

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.

Elementos Finitos no Lineales en Mecánica de Sólidos y Estructuras. Introducción básica y formulación del Método de los Elementos Finitos : Curso

Con el método de elementos finitos se ha desarrollado un procedimiento para el tratamiento de ecuaciones diferenciales y su resolución numérica. En este capítulo se va a realizar la aplicación a problemas de potencial regidos por la ecuación de Laplace. Se trata de una situación relativamente sencilla pero que permite comentar la mayoría de las características del método de elementos finitos en desplazamientos y sirve como introducción a situaciones más complicadas. Además, esta ecuación rige problemas de interés como la filtración, la torsión de cilindros, la distribución de termperaturas en régimen estacionario, etc. En lo que sigue la exposición se limita a problemas planos. Tras un primer ejemplo representativo del tema se realiza la formulación débil del problema de potencial. Acto seguido se introduce el concepto de discretización con toda su parafernalia de elementos, nudos, variables nodales, funciones de interpolación etc. El capítulo finaliza con la mención de algunos temas más especializados que el lector estudioso deberá proseguir en otro tipo de textos.

Fallo mecánico debido a grietas en Niobato de Litio bajo irradiación, simulación mediante elementos finitos

Ion-induced nano-track in LiNbO3. Motivation. From macro to nanoscale. Finite element method for nano-structured materials. Simulations of X-cut and Z-cut in LiNbO3. Experiments versus Simulations. Conclusions

Segmentation of RV in 4D Cardiac MR Volumes using region-merging graph cuts

Non-invasive quantitative assessment of the right ventricular anatomical and functional parameters is a challenging task. We present a semi-automatic approach for right ventricle (RV) segmentation from 4D MR images in two variants, which differ in the amount of user interaction. The method consists of three main phases: First, foreground and background markers are generated from the user input. Next, an over-segmented region image is obtained applying a watershed transform. Finally, these regions are merged using 4D graph-cuts with an intensity based boundary term. For the first variant the user outlines the inside of the RV wall in a few end-diastole slices, for the second two marker pixels serve as starting point for a statistical atlas application. Results were obtained by blind evaluation on 16 testing 4D MR volumes. They prove our method to be robust against markers location and place it favourably in the ranks of existing approaches.

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.

Contribución al desarrollo de técnicas CAD para el diseño de antenas impresas y dispositivos pasivos de microondas basadas en el método de elementos finitos

El objetivo principal de esta tesis es el desarrollo de herramientas numéricas basadas en técnicas de onda completa para el diseño asistido por ordenador (Computer-Aided Design,‘CAD’) de dispositivos de microondas. En este contexto, se desarrolla una herramienta numérica basada en el método de los elementos finitos para el diseño y análisis de antenas impresas mediante algoritmos de optimización. Esta técnica consiste en dividir el análisis de una antena en dos partes. Una parte de análisis 3D que se realiza sólo una vez en cada punto de frecuencia de la banda de funcionamiento donde se sustituye una superficie que contiene la metalización del parche por puertas artificiales. En una segunda parte se inserta entre las puertas artificiales en la estructura 3D la superficie soportando una metalización y se procede un análisis 2D para caracterizar el comportamiento de la antena. La técnica propuesta en esta tesis se puede implementar en un algoritmo de optimización para definir el perfil de la antena que permite conseguir los objetivos del diseño. Se valida experimentalmente dicha técnica empleándola en el diseño de antenas impresas de banda ancha para diferentes aplicaciones mediante la optimización del perfil de los parches. También, se desarrolla en esta tesis un procedimiento basado en el método de descomposición de dominio y el método de los elementos finitos para el diseño de dispositivos pasivos de microonda. Se utiliza este procedimiento en particular para el diseño y sintonía de filtros de microondas. En la primera etapa de su aplicación se divide la estructura que se quiere analizar en subdominios aplicando el método de descomposición de dominio, este proceso permite analizar cada segmento por separado utilizando el método de análisis adecuado dado que suele haber subdominios que se pueden analizar mediante métodos analíticos por lo que el tiempo de análisis es más reducido. Se utilizan métodos numéricos para analizar los subdominios que no se pueden analizar mediante métodos analíticos. En esta tesis, se utiliza el método de los elementos finitos para llevar a cabo el análisis. Además de la descomposición de dominio, se aplica un proceso de barrido en frecuencia para reducir los tiempos del análisis. Como método de orden reducido se utiliza la técnica de bases reducidas. Se ha utilizado este procedimiento para diseñar y sintonizar varios ejemplos de filtros con el fin de comprobar la validez de dicho procedimiento. Los resultados obtenidos demuestran la utilidad de este procedimiento y confirman su rigurosidad, precisión y eficiencia en el diseño de filtros de microondas. ABSTRACT The main objective of this thesis is the development of numerical tools based on full-wave techniques for computer-aided design ‘CAD’ of microwave devices. In this context, a numerical technique based on the finite element method ‘FEM’ for the design and analysis of printed antennas using optimization algorithms has been developed. The proposed technique consists in dividing the analysis of the antenna in two stages. In the first stage, the regions of the antenna which do not need to be modified during the CAD process are initially characterized only once from their corresponding matrix transfer function (Generalized Admittance matrix, ‘GAM’). The regions which will be modified are defined as artificial ports, precisely the regions which will contain the conducting surfaces of the printed antenna. In a second stage, the contour shape of the conducting surfaces of the printed antenna is iteratively modified in order to achieve a desired electromagnetic performance of the antenna. In this way, a new GAM of the radiating device which takes into account each printed antenna shape is computed after each iteration. The proposed technique can be implemented with a genetic algorithm to achieve the design objectives. This technique is validated experimentally and applied to the design of wideband printed antennas for different applications by optimizing the shape of the radiating device. In addition, a procedure based on the domain decomposition method and the finite element method has been developed for the design of microwave passive devices. In particular, this procedure can be applied to the design and tune of microwave filters. In the first stage of its implementation, the structure to be analyzed is divided into subdomains using the domain decomposition method; this process allows each subdomains can be analyzed separately using suitable analysis method, since there is usually subdomains that can be analyzed by analytical methods so that the time of analysis is reduced. For analyzing the subdomains that cannot be analyzed by analytical methods, we use the numerical methods. In this thesis, the FEM is used to carry out the analysis. Furthermore the decomposition of the domain, a frequency sweep process is applied to reduce analysis times. The reduced order model as the reduced basis technique is used in this procedure. This procedure is applied to the design and tune of several examples of microwave filters in order to check its validity. The obtained results allow concluding the usefulness of this procedure and confirming their thoroughness, accuracy and efficiency for the design of microwave filters.