Optimización de cimentaciones directas de medianería y esquina mediante modelos de elementos finitos

Existe un amplio catálogo de posibles soluciones para resolver la problemática de las zapatas de medianería así como, por extensión, las zapatas de esquina como caso particular de las anteriores. De ellas, las más habitualmente empleadas en estructuras de edificación son, por un lado, la utilización de una viga centradora que conecta la zapata de medianería con la zapata del pilar interior más próximo y, por otro, la colaboración de la viga de la primera planta trabajando como tirante. En la primera solución planteada, el equilibrio de la zapata de medianería y el centrado de la respuesta del terreno se consigue gracias a la colaboración del pilar interior con su cimentación y al trabajo a flexión de la viga centradora. La modelización clásica considera que se logra un centrado total de la reacción del terreno, con distribución uniforme de las tensiones de contacto bajo ambas zapatas. Este planteamiento presupone, por tanto, que la viga centradora logra evitar cualquier giro de la zapata de medianería y que el pilar puede, por ello, considerarse perfectamente empotrado en la cimentación. En este primer modelo, el protagonismo fundamental recae en la viga centradora, cuyo trabajo a flexión conduce frecuentemente a unas escuadrías y a unas cuantías de armado considerables. La segunda solución, plantea la colaboración de la viga de la primera planta, trabajando como tirante. De nuevo, los métodos convencionales suponen un éxito total en el mecanismo estabilizador del tirante, que logra evitar cualquier giro de la zapata de medianería, dando lugar a una distribución de tensiones también uniforme. Los modelos convencionales existentes para el cálculo de este tipo de cimentaciones presentan, por tanto, una serie de simplificaciones que permiten el cálculo de las mismas, por medios manuales, en un tiempo razonable, pero presentan el inconveniente de su posible alejamiento del comportamiento real de la cimentación, con las consecuencias negativas que ello puede suponer en el dimensionamiento de estos elementos estructurales. La presente tesis doctoral desarrolla un contraste de los modelos convencionales de cálculo de cimentaciones de medianería y esquina, mediante un análisis alternativo con modelos de elementos finitos, con el objetivo de poner de manifiesto las diferencias entre los resultados obtenidos con ambos tipos de modelización, analizar cuáles son las variables que más influyen en el comportamiento real de este tipo de cimentaciones y proponer un nuevo modelo de cálculo, de tipo convencional, más ajustado a la realidad. El proceso de investigación se desarrolla mediante una etapa experimental virtual que utiliza como modelo un pórtico tipo de edificación, ortogonal, de hormigón armado, con dos vanos y número variable de plantas. Tras identificar el posible giro de la cimentación como elemento clave en el comportamiento de las zapatas de medianería y de esquina, se adoptan como variables de estudio aquellas que mayor influencia puedan tener sobre el citado giro de las zapatas y sobre la rigidez del conjunto del elemento estructural. Así, se han estudiado luces de 3 m a 7 m, diferente número de plantas desde baja+1 hasta baja+4, resistencias del terreno desde 100 kN/m2 hasta 300 kN/m2, relaciones de forma de la zapata de medianería de 1,5 : 1 y 2 : 1, aumento y reducción de la cuantía de armado de la viga centradora y variación del canto de la viga centradora desde el mínimo canto compatible con el anclaje de la armadura de los pilares hasta un incremento del 75% respecto del citado canto mínimo. El conjunto de pórticos generados al aplicar las variables indicadas, se ha calculado tanto por métodos convencionales como por el método de los elementos finitos. Los resultados obtenidos ponen de manifiesto importantes discrepancias entre ambos métodos que conducen a importantes diferencias en el dimensionamiento de este tipo de cimentaciones. El empleo de los métodos tradicionales da lugar, por un lado, a un sobredimensionamiento de la armadura de la viga centradora y, por otro, a un infradimensionamiento, tanto del canto de la viga centradora, como del tamaño de la zapata de medianería y del armado de la viga de la primera planta. Finalizado el análisis y discusión de resultados, la tesis propone un nuevo método alternativo, de carácter convencional y, por tanto, aplicable a un cálculo manual en un tiempo razonable, que permite obtener los parámetros clave que regulan el comportamiento de las zapatas de medianería y esquina, conduciendo a un dimensionamiento más ajustado a las necesidades reales de este tipo de cimentación. There is a wide catalogue of possible solutions to solve the problem of party shoes and, by extension, corner shoes as a special case of the above. From all of them, the most commonly used in building structures are, on one hand, the use of a centering beam that connects the party shoe with the shoe of the nearest interior pillar and, on the other hand, the collaboration of the beam of the first floor working as a tie rod. In the first proposed solution, the balance of the party shoe and the centering of the ground response is achieved thanks to the collaboration of the interior pillar with his foundation along with the bending work of the centering beam. Classical modeling considers that a whole centering of the ground reaction is achieved, with uniform contact stress distribution under both shoes. This approach to the issue presupposes that the centering beam manages to avoid any rotation of the party shoe, so the pillar can be considered perfectly embedded in the foundation. In this first model, the leading role lies in the centering beam, whose bending work usually leads to important section sizes and high amounts of reinforced. The second solution, consideres the collaboration of the beam of the first floor, working as tie rod. Again, conventional methods involve a total success in the stabilizing mechanism of the tie rod, that manages to avoid any rotation of the party shoe, resulting in a stress distribution also uniform. Existing conventional models for calculating such foundations show, therefore, a series of simplifications which allow calculation of the same, by manual means, in a reasonable time, but have the disadvantage of the possible distance from the real behavior of the foundation, with the negative consequences this could bring in the dimensioning of these structural elements. The present thesis develops a contrast of conventional models of calculation of party and corner foundations by an alternative analysis with finite element models with the aim of bring to light the differences between the results obtained with both types of modeling, analysis which are the variables that influence the real behavior of this type of foundations and propose a new calculation model, conventional type, more adjusted to reality. The research process is developed through a virtual experimental stage using as a model a typical building frame, orthogonal, made of reinforced concrete, with two openings and variable number of floors. After identifying the possible spin of the foundation as the key element in the behavior of the party and corner shoes, it has been adopted as study variables, those that may have greater influence on the spin of the shoes and on the rigidity of the whole structural element. So, it have been studied lights from 3 m to 7 m, different number of floors from lower floor + 1 to lower floor + 4, máximum ground stresses from 100 kN/m2 300 kN/m2, shape relationships of party shoe 1,5:1 and 2:1, increase and decrease of the amount of reinforced of the centering beam and variation of the height of the centering beam from the minimum compatible with the anchoring of the reinforcement of pillars to an increase of 75% from the minimum quoted height. The set of frames generated by applying the indicated variables, is calculated both by conventional methods such as by the finite element method. The results show significant discrepancies between the two methods that lead to significant differences in the dimensioning of this type of foundation. The use of traditional methods results, on one hand, to an overdimensioning of the reinforced of the centering beam and, on the other hand, to an underdimensioning, both the height of the centering beam, such as the size of the party shoe and the reinforced of the beam of the first floor. After the analysis and discussion of results, the thesis proposes a new alternative method, conventional type and, therefore, applicable to a manual calculation in a reasonable time, that allows to obtain the key parameters that govern the behavior of party and corner shoes, leading to a dimensioning more adjusted to the real needings of this type of foundation.

Optimización de la agitación de un digestor anaerobio mediante mecánica de fluidos computacional

Debido al gran interés existente en el ahorro y recuperación de energía, y en el deseo de obtener productos que permitan usos beneficiosos del fango procedente de la depuración del agua residual, la digestión anaerobia es el proceso de estabilización de uso más extendido. El tiempo de retención de sólidos es un factor clave en el proceso de digestión anaerobia. En base al tiempo de retención de sólidos, se dimensiona el volumen de los digestores anaerobios para así obtener la reducción de materia orgánica deseada, con la correspondiente producción de biogás. La geometría del digestor y su sistema de agitación deben ser adecuados para alcanzar el tiempo de retención de sólidos de diseño. Los primeros trabajos sobre la agitación de los digestores realizaban únicamente experimentos con trazadores y otros métodos de medición. En otros casos, la mezcla era evaluada mediante la producción de biogás. Estas técnicas tenían el gran hándicap de no conocer lo que sucedía realmente dentro del digestor y sólo daban una idea aproximada de su funcionamiento. Mediante aplicación de la mecánica de fluidos computacional (CFD) es posible conocer con detalle las características del fluido objeto de estudio y, por lo tanto, simular perfectamente el movimiento del fango de un digestor anaerobio. En esta tesis se han simulado mediante CFD diferentes digestores a escala real (unos 2000 m3 de volumen) agitados con bomba/s de recirculación para alcanzar los siguientes objetivos: establecer la influencia de la relación entre el diámetro y la altura, de la pendiente de la solera, del número de bombas y del caudal de recirculación en dichos digestores, definir el campo de velocidades en la masa de fango y realizar un análisis energético y económico. Así, es posible conocer mejor cómo funciona el sistema de agitación de un digestor anaerobio a escala real equipado con bomba/s de recirculación. Los resultados obtenidos muestran que una relación diámetro/altura del digestor por encima de 1 empeora la agitación del mismo y que la pendiente en la solera del digestor favorece que la masa de fango esté mejor mezclada, siendo más determinante la esbeltez del tanque que la pendiente de su solera. No obstante, también es necesario elegir adecuadamente los parámetros de diseño del sistema de agitación, en este caso el caudal de recirculación de fango, para obtener una agitación completa sin apenas zonas muertas. En el caso de un digestor con una geometría inadecuada es posible mejorar su agitación aumentando el número de bombas de recirculación y el caudal de las mismas, pero no se llegará a alcanzar una agitación total de la masa de fango debido a su mal diseño original. Anaerobic digestion is the process for waste water treatment sludge stabilization of more widespread use due to the huge interest in saving and recovering energy and the wish to obtain products that allow beneficial uses for the sludge. The solids retention time is a key factor in the anaerobic digestion. Based on the solids retention time, volume anaerobic digester is sized to obtain the desired reduction in organic matter, with the corresponding production of biogas. The geometry of the digester and the stirring system should be adequate to achieve the design solid retention time. Early works on digesters stirring just performed tracer experiments and other measurement methods. In other cases, mixing was evaluated by biogas production. These techniques had the great handicap of not knowing what really happened inside the digester and they only gave a rough idea of its operation. By application of computational fluid dynamics (CFD), it is possible to know in detail the characteristics of the fluid under study and, therefore, simulate perfectly the sludge movement of an anaerobic digester. Different full-scale digesters (about 2000 m3 of volume) agitated with pump/s recirculation have been simulated by CFD in this thesis to achieve the following objectives: to establish the influence of the relationship between the diameter and height, the slope of the bottom, the number of pumps and the recirculation flow in such digesters, to define the velocity field in the mass of sludge and carry out an energy and economic analysis. Thus, it is possible to understand better how the agitation system of a full-scale anaerobic digester equipped with pump/s recirculation works. The results achieved show that a diameter/height ratio of the digester above 1 worsens its stirring and that the slope of the digester bottom favors that the mass of sludge is better mixed, being more decisive the tank slenderness than the slope of its bottom. However, it is also necessary to select properly the design parameters of the agitation system, in this case the sludge recirculation flow rate, for a complete agitation with little dead zones. In the case of a digester with inadequate geometry, its agitation can be improved by increasing the number of recirculation pumps and flow of them, but it will not reach a full agitation of the mass of sludge because of the poor original design.

Modelos simplificados de elementos finitos para el estudio del colapso progresivo en estructuras de hormigón armado sometido a explosiones

La mayoría de las estructuras construidas actualmente están constituidas por hormigón armado, el cual, a pesar de poseer una gran resistencia, está expuesto a múltiples amenazas. Una de ellas, y de las que más preocupa a la sociedad actualmente, son los posibles ataques terroristas como impactos o colocación de explosivos en las proximidades de las estructuras o en el interior de las mismas. Por ello, cada vez más se está mostrando un gran interés en el comportamiento de la estructura y sus diferentes componentes para poder prevenir cualquier tipo de daño. Con este fin, el presente proyecto aborda el estudio de un caso de colapso de estructura sometida a una carga explosiva. El estudio y la modelización del hormigón armado presenta un gran reto debido a los múltiples factores que definen su comportamiento, marcadamente no lineal cuando se producen en él grandes deformaciones. Se ha procedido a la realización del estudio de una estructura porticada de dos plantas, utilizando diferentes técnicas de simulación mediante el programa de elementos finitos LS-DYNA. En concreto se han empleado tres modelos diferentes, un modelo de elementos finitos sólidos tridimensionales con materiales segregados (hormigón y acero), un modelo de elementos finitos estructurales (viga, lámina) con material homogeneizado, y un modelo de elementos finitos estructurales con materiales segregados, comparando la evolución del colapso de la estructura tras la demolición de uno de los pilares centrales. Se ha analizado también el desplazamiento y la velocidad de colapso para un punto concreto. Estos resultados han sido comparados con los datos obtenidos de la fase de ensayo del Proyecto SEGTRANS sobre el mismo punto de la estructura. Con todo ello se ha podido determinar la viabilidad de los modelos para el estudio del colapso de la estructura porticada empleada en este trabajo, señalando las ventajas de cada uno de ellos, y la posibilidad de la aplicación de estos modelos para el estudio de la estructura de un edificio completo teniendo en cuenta las limitaciones de los mismos y de los recursos existentes actualmente.

Evaluación del tamaño óptimo de mallado para la modelización mediante elementos finitos de losas de hormigón frente a cargas explosivas

Los elementos estructurales empleados en construcción no han sido en general diseñados para soportar acciones impulsivas, como la detonación de artefactos explosivos. Desde el siglo pasado el mundo ha sufrido ataques terroristas en los que en muchos casos se han producido explosiones que causaron víctimas, heridos y la destrucción de las construcciones próximas. Debido a este hecho, instituciones públicas y privadas comenzaron a mostrar interés por el comportamiento de los elementos estructurales que componen sus instalaciones. El hormigón armado es uno de los principales materiales utilizados en las estructuras de obras debido a sus buenas características, cuyo análisis y modelización en deformaciones dinámicas supone un campo de desafíos al que se está prestando gran atención en los últimos años. LS-DYNA® es un programa basado en elementos finitos capaz de simular problemas reales complejos en el que se han desarrollado distintos modelos de hormigón. Tres de esos modelos (K&C, RHT y CSCM) son evaluados con losas de distintos tamaños de mallado de elementos finitos frente a la detonación de 2 kg de TNT situados a 1 m de distancia. Dichos modelos son simulados y se obtienen los valores de las aceleraciones máximas en unos determinados puntos de las losas. Los valores son sometidos a la aplicación del Método GCI (Grid Convergence Index) para una relación de refinamiento de mallado no constante, cuyos resultados se comparan con aquellos registrados en los acelerómetros empleados durante la primera fase de ensayos del Proyecto SEGTRANS. Mediante el análisis de los resultados obtenidos se determina cual es el modelo de material y tamaño de mallado más adecuado que pueda emplearse en un futuro para poder modelar estructuras más complejas y con niveles de explosivo más elevados.

Apuntes de Mecánica de Fluidos

Esta publicacion contiene unos apuntes para la parte de Mecánica de Fluidos de las asignaturas de la ETSAM, junto con cierta cantidad de material adicional que excede el ámbito del curso. Aunque está lejos de ser completa y de estar libre de errores, el autor espera sea útil.

El método de los elementos de contorno en los problemas de mecánica del suelo

La preparación de estas notas ha llevado, al más veterano de los autores, a rememorar sus primeros tanteos con los métodos numéricos. Tratando de desarrollar su tesis doctoral sobre efectos dinámicos en puentes de ferrocarril, descubrió, en 1968, en la biblioteca del Laboratorio de Transporte (donde el profesor ]iménez Salas era Director) las Actas de la reunión ASTM en las que Quilan y Sung proponían la asimilación del comportamiento dinámico del semiespacio elástico a un sistema con un grado de libertad. Además de incorporar estos resultados a un modelo de puente para tener en cuenta los fenómenos de interacción dinámica terreno-estructura dicho autor entró en contacto con algunos miembros del equipo de investigación del Prof. ]iménez Salas que, por entonces, estaba explorando la posibilidad de aplicación del ordenador y los métodos numéricos necesarios para tratar los problemas más difíciles de Mecánica de los Medios Continuos. De hecho fue ese grupo quien contribuyó a introducir en España el método de los elementos finitos en la ingeniería civil, pero además, y en relación directa con el título de este artículo fue el propio ]iménez Salas quién inició la línea de trabajo de lo que mas tarde se ha llamado Método Indirecto de Elementos de Contorno que luego fue seguida por otros miembros de su grupo. En aquélla época poco podía sospechar el autor precitado que iba a dedicar una parte sustancial de su vida al desarrollo de ese método numérico en su versión Directa y mucho menos que gran parte de la motivación vendría del problema de interacción dinámica terreno-estructura, una de cuyas primeras soluciones había obtenido en la mencionada visita al Laboratorio de Transporte. En efecto los autores trataban en 1975 de encontrar un procedimiento que les permitiera afrontar el estudio de la interacción en túneles sometidos a carga sísmica y tropezaron, al utilizar el método de elementos finitos, con el problema de las reflexiones de ondas en los contornos artificiales creados al truncar la malla de cálculo. Deseando evitar el uso contornos absorbentes y otros recursos similares se exploró la posibilidad de soluciones fundamentales que incorporasen el comportamiento en el infmito y, fruto de ello, fueron los primeros trabajos que introdujeron el Método Directo de los Elementos de Contorno en España en problemas estáticos. La extensión a teoría del potencial, dinámica en el dominio de la frecuencia, plasticidad, etc tuvo lugar inmediatamente siendo en la mayoría de los casos los problemas típicos de mecánica del suelo los que motivaron y justifican el esfuerzo realizado. Un campo apasionante, el de la poroelasticidad ha dado lugar a nuevas contribuciones y también se han escrito libros de diverso calado que describen las posibilidades del método para dar contestación a preguntas de gran importancia técnica. Los autores quieren poner de manifiesto que la redacción de este trabajo, debe considerarse no solo como la muestra de algunos resultados de aplicación a problemas prácticos, sino también como un homenaje y reconocimiento explícito a la labor precursora del Prof. ]iménez Salas y a su espíritu de permanente curiosidad por el conocimiento científico.

Numerical Error Prediction and its applications in CFD using tau-estimation

Nowadays, Computational Fluid Dynamics (CFD) solvers are widely used within the industry to model fluid flow phenomenons. Several fluid flow model equations have been employed in the last decades to simulate and predict forces acting, for example, on different aircraft configurations. Computational time and accuracy are strongly dependent on the fluid flow model equation and the spatial dimension of the problem considered. While simple models based on perfect flows, like panel methods or potential flow models can be very fast to solve, they usually suffer from a poor accuracy in order to simulate real flows (transonic, viscous). On the other hand, more complex models such as the full Navier- Stokes equations provide high fidelity predictions but at a much higher computational cost. Thus, a good compromise between accuracy and computational time has to be fixed for engineering applications. A discretisation technique widely used within the industry is the so-called Finite Volume approach on unstructured meshes. This technique spatially discretises the flow motion equations onto a set of elements which form a mesh, a discrete representation of the continuous domain. Using this approach, for a given flow model equation, the accuracy and computational time mainly depend on the distribution of nodes forming the mesh. Therefore, a good compromise between accuracy and computational time might be obtained by carefully defining the mesh. However, defining an optimal mesh for complex flows and geometries requires a very high level expertize in fluid mechanics and numerical analysis, and in most cases a simple guess of regions of the computational domain which might affect the most the accuracy is impossible. Thus, it is desirable to have an automatized remeshing tool, which is more flexible with unstructured meshes than its structured counterpart. However, adaptive methods currently in use still have an opened question: how to efficiently drive the adaptation ? Pioneering sensors based on flow features generally suffer from a lack of reliability, so in the last decade more effort has been made in developing numerical error-based sensors, like for instance the adjoint-based adaptation sensors. While very efficient at adapting meshes for a given functional output, the latter method is very expensive as it requires to solve a dual set of equations and computes the sensor on an embedded mesh. Therefore, it would be desirable to develop a more affordable numerical error estimation method. The current work aims at estimating the truncation error, which arises when discretising a partial differential equation. These are the higher order terms neglected in the construction of the numerical scheme. The truncation error provides very useful information as it is strongly related to the flow model equation and its discretisation. On one hand, it is a very reliable measure of the quality of the mesh, therefore very useful in order to drive a mesh adaptation procedure. On the other hand, it is strongly linked to the flow model equation, so that a careful estimation actually gives information on how well a given equation is solved, which may be useful in the context of _ -extrapolation or zonal modelling. The following work is organized as follows: Chap. 1 contains a short review of mesh adaptation techniques as well as numerical error prediction. In the first section, Sec. 1.1, the basic refinement strategies are reviewed and the main contribution to structured and unstructured mesh adaptation are presented. Sec. 1.2 introduces the definitions of errors encountered when solving Computational Fluid Dynamics problems and reviews the most common approaches to predict them. Chap. 2 is devoted to the mathematical formulation of truncation error estimation in the context of finite volume methodology, as well as a complete verification procedure. Several features are studied, such as the influence of grid non-uniformities, non-linearity, boundary conditions and non-converged numerical solutions. This verification part has been submitted and accepted for publication in the Journal of Computational Physics. Chap. 3 presents a mesh adaptation algorithm based on truncation error estimates and compares the results to a feature-based and an adjoint-based sensor (in collaboration with Jorge Ponsín, INTA). Two- and three-dimensional cases relevant for validation in the aeronautical industry are considered. This part has been submitted and accepted in the AIAA Journal. An extension to Reynolds Averaged Navier- Stokes equations is also included, where _ -estimation-based mesh adaptation and _ -extrapolation are applied to viscous wing profiles. The latter has been submitted in the Proceedings of the Institution of Mechanical Engineers, Part G: Journal of Aerospace Engineering. Keywords: mesh adaptation, numerical error prediction, finite volume Hoy en día, la Dinámica de Fluidos Computacional (CFD) es ampliamente utilizada dentro de la industria para obtener información sobre fenómenos fluidos. La Dinámica de Fluidos Computacional considera distintas modelizaciones de las ecuaciones fluidas (Potencial, Euler, Navier-Stokes, etc) para simular y predecir las fuerzas que actúan, por ejemplo, sobre una configuración de aeronave. El tiempo de cálculo y la precisión en la solución depende en gran medida de los modelos utilizados, así como de la dimensión espacial del problema considerado. Mientras que modelos simples basados en flujos perfectos, como modelos de flujos potenciales, se pueden resolver rápidamente, por lo general aducen de una baja precisión a la hora de simular flujos reales (viscosos, transónicos, etc). Por otro lado, modelos más complejos tales como el conjunto de ecuaciones de Navier-Stokes proporcionan predicciones de alta fidelidad, a expensas de un coste computacional mucho más elevado. Por lo tanto, en términos de aplicaciones de ingeniería se debe fijar un buen compromiso entre precisión y tiempo de cálculo. Una técnica de discretización ampliamente utilizada en la industria es el método de los Volúmenes Finitos en mallas no estructuradas. Esta técnica discretiza espacialmente las ecuaciones del movimiento del flujo sobre un conjunto de elementos que forman una malla, una representación discreta del dominio continuo. Utilizando este enfoque, para una ecuación de flujo dado, la precisión y el tiempo computacional dependen principalmente de la distribución de los nodos que forman la malla. Por consiguiente, un buen compromiso entre precisión y tiempo de cálculo se podría obtener definiendo cuidadosamente la malla, concentrando sus elementos en aquellas zonas donde sea estrictamente necesario. Sin embargo, la definición de una malla óptima para corrientes y geometrías complejas requiere un nivel muy alto de experiencia en la mecánica de fluidos y el análisis numérico, así como un conocimiento previo de la solución. Aspecto que en la mayoría de los casos no está disponible. Por tanto, es deseable tener una herramienta que permita adaptar los elementos de malla de forma automática, acorde a la solución fluida (remallado). Esta herramienta es generalmente más flexible en mallas no estructuradas que con su homóloga estructurada. No obstante, los métodos de adaptación actualmente en uso todavía dejan una pregunta abierta: cómo conducir de manera eficiente la adaptación. Sensores pioneros basados en las características del flujo en general, adolecen de una falta de fiabilidad, por lo que en la última década se han realizado grandes esfuerzos en el desarrollo numérico de sensores basados en el error, como por ejemplo los sensores basados en el adjunto. A pesar de ser muy eficientes en la adaptación de mallas para un determinado funcional, este último método resulta muy costoso, pues requiere resolver un doble conjunto de ecuaciones: la solución y su adjunta. Por tanto, es deseable desarrollar un método numérico de estimación de error más asequible. El presente trabajo tiene como objetivo estimar el error local de truncación, que aparece cuando se discretiza una ecuación en derivadas parciales. Estos son los términos de orden superior olvidados en la construcción del esquema numérico. El error de truncación proporciona una información muy útil sobre la solución: es una medida muy fiable de la calidad de la malla, obteniendo información que permite llevar a cabo un procedimiento de adaptación de malla. Está fuertemente relacionado al modelo matemático fluido, de modo que una estimación precisa garantiza la idoneidad de dicho modelo en un campo fluido, lo que puede ser útil en el contexto de modelado zonal. Por último, permite mejorar la precisión de la solución resolviendo un nuevo sistema donde el error local actúa como término fuente (_ -extrapolación). El presenta trabajo se organiza de la siguiente manera: Cap. 1 contiene una breve reseña de las técnicas de adaptación de malla, así como de los métodos de predicción de los errores numéricos. En la primera sección, Sec. 1.1, se examinan las estrategias básicas de refinamiento y se presenta la principal contribución a la adaptación de malla estructurada y no estructurada. Sec 1.2 introduce las definiciones de los errores encontrados en la resolución de problemas de Dinámica Computacional de Fluidos y se examinan los enfoques más comunes para predecirlos. Cap. 2 está dedicado a la formulación matemática de la estimación del error de truncación en el contexto de la metodología de Volúmenes Finitos, así como a un procedimiento de verificación completo. Se estudian varias características que influyen en su estimación: la influencia de la falta de uniformidad de la malla, el efecto de las no linealidades del modelo matemático, diferentes condiciones de contorno y soluciones numéricas no convergidas. Esta parte de verificación ha sido presentada y aceptada para su publicación en el Journal of Computational Physics. Cap. 3 presenta un algoritmo de adaptación de malla basado en la estimación del error de truncación y compara los resultados con sensores de featured-based y adjointbased (en colaboración con Jorge Ponsín del INTA). Se consideran casos en dos y tres dimensiones, relevantes para la validación en la industria aeronáutica. Este trabajo ha sido presentado y aceptado en el AIAA Journal. También se incluye una extensión de estos métodos a las ecuaciones RANS (Reynolds Average Navier- Stokes), en donde adaptación de malla basada en _ y _ -extrapolación son aplicados a perfiles con viscosidad de alas. Este último trabajo se ha presentado en los Actas de la Institución de Ingenieros Mecánicos, Parte G: Journal of Aerospace Engineering. Palabras clave: adaptación de malla, predicción del error numérico, volúmenes finitos

Cargas explosivas sobre pilares y forjados modelados con elementos tipo viga y lámina

Los fenómenos de impacto y explosión sobre estructuras de hormigón tienen efectos en muchos casos catastróficos a pesar de su reducida probabilidad. Las estructuras de hormigón no suelen estar diseñadas para resistir este tipo de solicitaciones dinámicas. El análisis numérico mediante elementos finitos con integración explícita permite una aproximación suficiente a los efectos de la onda explosiva sobre pilares y forjados de estructuras reticuladas de hormigón. Los materiales recientemente implementados en LS-Dyna para hormigón como el CSCM [1], para elementos de continuo 3D, y la formulación que proporciona la debida compatibilidad con los elementos viga de acero dispuestos de forma segregada, permite estudiar de forma realista modelos detallados de pilares y forjados. Pero las limitaciones computacionales hacen inviable emplear estos métodos en estructuras completas. Como alternativa es posible usar modelos de elementos estructurales de vigas y láminas para el análisis de estas estructuras. Sin embargo es necesario un adecuado ajuste de parámetros y propiedades en estos modelos. Este trabajo muestra un método con en el que obtener modelos de elementos estructurales, elementos viga y lámina, usando modelos de material [2] adecuados para ellos, junto a un procedimiento para incluir la armadura de forma adecuada. Utilizando este método es posible representar con suficiente aproximación el comportamiento de modelos detallados realistas de forjados y pilares de estructuras reticuladas de hormigón frente a acciones explosivas, posibilitando el análisis de una estructura completa frente a explosión.

Estudio dinámico de puentes de ferrocarril sometidos a nuevos modelos de carga en el marco europeo

Más allá de la alta velocidad, y debido en parte a su éxito, en los últimos años han aparecido trenes de pasajeros capaces de circular a velocidades más elevadas y ser más eficientes (como, por ejemplo, los de tipo automotor). Además, algunos de los antiguos modelos recogidos en las normativas parecen no cubrir determinados casos, lo que ha llevado a la reciente modificación de éstas. A todo esto se le añade la necesidad de crear una normativa común europea que establezca las bases para lograr la interoperabilidad. Con este proyecto se pretende estudiar los efectos dinámicos en puentes de estas nuevas tipologías de trenes, haciendo especial hincapié en los trenes automotores. El trabajo se centra principalmente en puentes isostáticos de luces cortas, por ser éstos especialmente sensibles a este tipo de acciones. El parámetro crítico, en estos casos, son las aceleraciones en centro de vano. Se realizara el cálculo dinámico por diferentes métodos: elementos finitos, análisis modal, así como programas de integración numérica en Matlab de los modos de vibración, analizando y comparando la efectividad de estos. Las diferencias entre la consideración o no de la distribución de las cargas transmitidas al tablero de la estructura de la vía a serán discutidas en el proyecto, así como la interacción dinámica entre estructura y vehículo. Finalmente se aplicarán los conocimientos obtenidos a varios casos reales de puentes incluidos en la red ferroviaria española, como por ejemplo, el viaducto de alta velocidad del Jalón. Los resultados de estos estudios tienen una repercusión en diversas normativas europeas actualmente en revisión y en las cuales está involucrado el Grupo de Mecánica Computacional: EN15528 European Committee for Standarization (2008), EN1991-2 European Committee for Standarization (2003) y las Especificaciones Técnicas de Interoperabilidad Ferroviaria (TSI).

Fenómenos de resonancia asociados al transporte de mercancías líquidas por ferrocarril

Los vagones de mercancías se unen entre sí mediante acoplamientos UIC estándar. Aunque cada plataforma está equipada con un tensor de enganche en ambos extremos, sólo uno de ellos se utiliza, mientras que el otro descansa apoyado en un gancho. Debido a los efectos de oleaje del líquido contenido en las cisternas, este tensor en reposo puede entrar en resonancia con la frecuencia natural de algunos fluidos, llegando a oscilar de forma apreciable. En ocasiones llega a saltar fuera de su alojamiento, provocando ciertas molestias durante la circulación, aunque sin entrañar ningún riesgo para la seguridad de la marcha. Con el propósito de simular la influencia de diversos factores sobre el comportamiento dinámico asociado a estos fenómenos de oleaje, se ha desarrollado un modelo de la interacción fluido-vehículo, empleado técnicas de modelado de sistemas multicuerpo y de elementos finitos.

Estudios de confort de vehículos ferroviarios mediante técnicas combinadas de simulación FEM y MBS

La elasticidad de la caja de viajeros de los vehículos ferroviarios tiene una gran influencia sobre el confort. Por esta razón, cuando se desea simular el comportamiento dinámico del vehículo para estudios de confort, resulta conveniente construir un modelo elástico de caja, a fin de obtener resultados más precisos. La construcción de este tipo de modelos pasa por el desarrollo de dos etapas fundamentales, que comprenden la generación de un modelo de caja mediante la técnica de los elementos finitos (FEM) para su caracterización dinámica como cuerpo elástico y la definición de un modelo de sistema multicuerpo (MBS) que englobe los restantes componentes del vehículo. En este artículo se presentan los resultados obtenidos en un estudio comparativo llevado a cabo con modelos de caja rígida y elástica, en los que se ha valorado el nivel de confort obtenido con ambas configuraciones. Para ello, se ha simulado el comportamiento del vehículo a dos velocidades distintas, de 70 y 110km/h, y con dos niveles de irregularidades. Se han analizado las aceleraciones de la caja, que se han procesado de acuerdo a las especificaciones de la norma EN12299, a fin de obtener el índice de comodidad. Este parámetro se ha utilizado para comparar el nivel de confort obtenido con ambos modelos, habiéndose encontrado una gran diferencia en los índices calculados con caja rígida y con caja elástica, lo que confirma la gran influencia de la elasticidad de la caja en los estudios de confort llevados a cabo mediante técnicas de simulación dinámica.

Desarrollo en la definición geométrica y en el modelado de la degradación mecánica de las paredes arteriales en el cayado aórtico

Desarrollo de técnicas para la reconstrucción de imágenes médicas y mallas de elementos finitos. Desarrollo de modelos de comportamiento mecánico para reproducir la degradación de la capacidad portante en la pared arterial.

Estabilidad del tráfico en viaductos de carretera y ferrocarril

Cuando un automóvil o un tren circula sobre un viaducto se producen efectos dinámicos que, además de aumentar los esfuerzos que ha de soportar el puente, afectan a los propios vehículos. Los fenómenos de vibración lateral en los viaductos de ferrocarril están contemplados de manera muy sucinta en la normativa y pueden afectar a la seguridad y comodidad del tráfico. En el caso de las carreteras, esa dinámica también puede contribuir al riesgo de accidente y al aumento de la sensación de incomodidad, sobre todo si se considera la presencia de vientos laterales. La preocupación por estos fenómenos queda reflejada en el número de artículos publicados en revistas científicas en los últimos años. Para este trabajo se han desarrollado modelos de interacción dinámica que representan el comportamiento de las estructuras, mediante el método de los elementos finitos, y el de los vehículos, mediante sistemas multicuerpo.

Análisis estructural 3D de los invernaderos multitunel del sudeste de España. Elementos verticales

Este trabajo analiza la capacidad estructural de los elementos verticales que componen el pórtico resistente más limitativo de la variante de invernadero multitúnel más utilizada en el sudeste de España. Se estudia la influencia de la longitud y anchura del invernadero en la distribución de tensiones que se produce en dichos elementos. Se ha realizado un modelo tridimensional de elementos finitos con el programa ROBOT MILLENNIUM en su versión 20.1 en el que se ha representado la estructura de barras que conforma el invernadero. Se han analizado los resultados de tensiones normales máximas de 16 estructuras con geometría de túnel similar, pero de diferente anchura y longitud, de tal forma que la superficie en planta oscila entre los 600 m 2 y los 38.400 m 2 . El cálculo de las acciones existentes sobre los invernaderos se ha realizado según la norma UNE-EN 13031-1. Los resultados numéricos obtenidos han sido estudiados para determinar la influencia de las dimensiones del invernadero en la distribución de tensiones. Aunque se observan significativas reducciones en las tensiones con la adición de túneles, no se puede generalizar la influencia del aumento de las dimensiones del invernadero en los cambios tensionales del mismo. Los resultados obtenidos se consideran de gran utilidad para la optimización de las estructuras de los invernaderos multitúnel así como para garantizar la seguridad mecánica de las mismas

Análisis estructural 3D de los invernaderos multitunel del sudeste de España. Elementos horizontales

Este trabajo analiza la capacidad estructural de los Este trabajo analiza la capacidad estructural de los elementos horizontales que componen el pórtico resistente más limitativo de la variante de invernadero multitúnel más utilizada en el sudeste de España. Se analiza la influencia de la anchura y longitud del invernadero en la distribución de tensiones que se producen en dichos elementos. Se ha realizado un modelo tridimensional de elementos finitos con el programa ROBOT MILLENNIUM en su versión 20.1 en el que se ha representado la estructura de barras que conforma el invernadero. Se han analizado los resultados de tensiones normales máximas de 16 estructuras con geometría de túnel similar, pero de diferente anchura y longitud, de tal forma que la superficie en planta oscila entre los 600 m 2 y los 38.400 m 2 .El cálculo de las acciones existentes sobre los invernaderos se ha realizado según la norma UNE-EN 13031-1. El aumento de las dimensiones del invernadero conlleva importantes variaciones en las tensiones máximas registradas por canalones. Sin embargo, en arcos y tirantes de cercha, las tensiones permanecen constantes. Los resultados obtenidos se consideran de gran utilidad para la optimización de las estructuras de los invernaderos multitúnel así como para garantizar la seguridad mecánica de las mismas.