Determinación de mecanismos de deslizamientosa en áreas urbanizadas de la Sierra de las Apretaderas, término de Alcuzcuz, Benahavís, Málaga

La naturaleza intrínseca del material mantiene una estrecha relación con el tipo de inestabilidad. El trabajo recoge la experiencia en la determinación de mecanismo de deslizamientos en el área de la Sierra de Las Apretaderas, término de Alcuzcuz, en Benahavís, Málaga. La intensa urbanización de amplias áreas serranas conlleva la transformación del medio a través de la construcción de infraestructuras de comunicación y extensas áreas residenciales y la interacción de dichas construcciones con el medio provoca a menudo el desencadenamiento de siniestros relacionados con los deslizamientos de ladera en ocasiones poco sospechosos. El trabajo que se presenta recoge la experiencia de un siniestro acaecido bajo estas circunstancias y provocando numerosos daños materiales en edificaciones de lujo e infraestructuras anejas.

Estudio de las concentraciones y mecanismos de desplazamiento de radón y descendientes en recintos cerrados

La importancia de la medida de concentración de radón en recintos cerrados radica en que la dosis debida a su presencia en el ambiente representa aproximadamente el 50% de la dosis total de la radiación de origen natural recibida por la población, estando presente tanto en espacios abiertos, como en el interior de edificios. Se han desarrollado numerosos estudios encaminados a determinar el riesgo radiológico por unidad de exposición Para alcanzar este fin se han puesto a punto diversas técnicas que permiten analizar, tanto en una instalación de laboratorio como en espacios habitados, el factor de equilibrio, la fracción libre y el espectro dimensional de descendientes del radón en estado libre. Por otra parte, si tenemos en cuenta los mecanismos de desplazamiento del 222 Rn sabemos que éste se difunde a través de las grietas del terreno y de los materiales de construcción. Por lo tanto, tendremos que en recintos cerrados el mecanismo difusivo es el más significativo, siendo en este caso determinante la longitud de difusión desintegrarse.

Efficient numerical methods for global linear instability. Analysis and modeling of non-orthogonal swept attachment-line boundary layer flow

The aim of this thesis is to study the mechanisms of instability that occur in swept wings when the angle of attack increases. For this, a simplified model for the a simplified model for the non-orthogonal swept leading edge boundary layer has been used as well as different numerical techniques in order to solve the linear stability problem that describes the behavior of perturbations superposed upon this base flow. Two different approaches, matrix-free and matrix forming methods, have been validated using direct numerical simulations with spectral resolution. In this way, flow instability in the non-orthogonal swept attachment-line boundary layer is addressed in a linear analysis framework via the solution of the pertinent global (Bi-Global) PDE-based eigenvalue problem. Subsequently, a simple extension of the extended G¨ortler-H¨ammerlin ODEbased polynomial model proposed by Theofilis, Fedorov, Obrist & Dallmann (2003) for orthogonal flow, which includes previous models as particular cases and recovers global instability analysis results, is presented for non-orthogonal flow. Direct numerical simulations have been used to verify the stability results and unravel the limits of validity of the basic flow model analyzed. The effect of the angle of attack, AoA, on the critical conditions of the non-orthogonal problem has been documented; an increase of the angle of attack, from AoA = 0 (orthogonal flow) up to values close to _/2 which make the assumptions under which the basic flow is derived questionable, is found to systematically destabilize the flow. The critical conditions of non-orthogonal flows at 0 _ AoA _ _/2 are shown to be recoverable from those of orthogonal flow, via a simple analytical transformation involving AoA. These results can help to understand the mechanisms of destabilization that occurs in the attachment line of wings at finite angles of attack. Studies taking into account variations of the pressure field in the basic flow or the extension to compressible flows are issues that remain open. El objetivo de esta tesis es estudiar los mecanismos de la inestabilidad que se producen en ciertos dispositivos aerodinámicos cuando se aumenta el ángulo de ataque. Para ello se ha utilizado un modelo simplificado del flujo de base, así como diferentes técnicas numéricas, con el fin de resolver el problema de estabilidad lineal asociado que describe el comportamiento de las perturbaciones. Estos métodos; sin y con formación de matriz, se han validado utilizando simulaciones numéricas directas con resolución espectral. De esta manera, la inestabilidad del flujo de capa límite laminar oblicuo entorno a la línea de estancamiento se aborda en un marco de análisis lineal por medio del método Bi-Global de resolución del problema de valores propios en derivadas parciales. Posteriormente se propone una extensión simple para el flujo no-ortogonal del modelo polinomial de ecuaciones diferenciales ordinarias, G¨ortler-H¨ammerlin extendido, propuesto por Theofilis et al. (2003) para el flujo ortogonal, que incluye los modelos previos como casos particulares y recupera los resultados del analisis global de estabilidad lineal. Se han realizado simulaciones directas con el fin de verificar los resultados del análisis de estabilidad así como para investigar los límites de validez del modelo de flujo base utilizado. En este trabajo se ha documentado el efecto del ángulo de ataque AoA en las condiciones críticas del problema no ortogonal obteniendo que el incremento del ángulo de ataque, de AoA = 0 (flujo ortogonal) hasta valores próximos a _/2, en el cual las hipótesis sobre las que se basa el flujo base dejan de ser válidas, tiende sistemáticamente a desestabilizar el flujo. Las condiciones críticas del caso no ortogonal 0 _ AoA _ _/2 pueden recuperarse a partir del caso ortogonal mediante el uso de una transformación analítica simple que implica el ángulo de ataque AoA. Estos resultados pueden ayudar a comprender los mecanismos de desestabilización que se producen en el borde de ataque de las alas de los aviones a ángulos de ataque finitos. Como tareas pendientes quedaría realizar estudios que tengan en cuenta variaciones del campo de presión en el flujo base así como la extensión de éste al caso de flujos compresibles.

Flexibilidad discreta de mecanismos y estructuras

Cada vez con más frecuencia los procesos de diseño y fabricación de componentes mecánicos demandan herramientas de simulación que permitan el estudio dinámico de mecanismos con piezas deformables. Esto plantea dificultades numéricas debidas fundamentalmente al carácter no lineal de la dinámica de estos dispositivos; no linealidad fundamentalmente geométrica asociada a los grandes desplazamientos y/o grandes deformaciones. Existen en la literatura modelos de vigas no lineales capaces de resolver con precisión este tipo de problemas, aunque suelen ser computacionalmente costosos. Este trabajo pretende explorar la posibilidad de sustituir los elementos continuos deformables tipo viga por otros equivalentes formados por elementos rígidos interconectados con una flexibilidad discreta concentrada en los nodos. De esta manera, se pretende sustituir elementos que tienen una flexibilidad continua por otros equivalentes que la tienen discreta, concentrada en las articulaciones. Existen precedentes de estudios de discretización de vigas en barras rígidas unidas por muelles y resortes de flexión [3]. La novedad de este trabajo es que representa la rigidez a flexión y tracción mediante potenciales de penalización asociados a restricciones. Se estudiará una viga biapoyada sencilla en pequeñas deformaciones para establecer una equivalencia entre la viga continua y la viga discretizada con restricciones que representan la flexión. Para ello se obtendrá una relación entre la penalización a emplear y las características mecánicas de la viga y su discretización. El modelo resultante se implementará en un entorno de cálculo de C++, disponible en el grupo de investigación, para realizar diferentes cálculos que permitan evaluar los resultados.

Integración de la plataforma 4mobile con mecanismos de atención en redes sociales

En este documento se plasma el trabajo que se ha realizado para la incorporación de un mecanismo de atención en las redes sociales para la plataforma 4mobile. Inicialmente se ha realizado un estudio de mercado sobre como gestionan las empresas, en el caso de que tengan, el mecanismo de atención social en las redes sociales. Tras el estudio de otras plataformas similares hemos realizado el análisis de la plataforma 4mobile para la cual se necesita agregar el mecanismo de atención social. Se ha plasmado a lo largo del documento todas y cada una de las fases por las que se ha tenido que pasar para la realización del proyecto desde la explicación de cuáles son las funcionalidades que ofrece 4mobile hasta la elección de una herramienta que sea capaz de satisfacer las necesidades que demanda la plataforma. Con ello se ha pretendido que la plataforma 4mobile pueda adaptarse a las necesidades de negocio demandadas por los clientes y de esta forma poder proporcionar valor añadido en el mercado.

Inestabilidad de pórticos planos metálicos de edificación con inclusión de la flexibilidad de las uniones y elementos no prismáticos : un enfoque para la optimización

En esta tesis se aborda el problema de la modelización, análisis y optimización de pórticos metálicos planos de edificación frente a los estados límites último y de servicio. El objetivo general es presentar una técnica secuencial ordenada de optimización discreta para obtener el coste mínimo de pórticos metálicos planos de edificación, teniendo en cuenta las especificaciones del EC-3, incorporando las uniones semirrígidas y elementos no prismáticos en el proceso de diseño. Asimismo se persigue valorar su grado de influencia sobre el diseño final. El horizonte es extraer conclusiones prácticas que puedan ser de utilidad y aplicación simple para el proyecto de estructuras metálicas. La cantidad de publicaciones técnicas y científicas sobre la respuesta estructural de entramados metálicos es inmensa; por ello se ha hecho un esfuerzo intenso en recopilar el estado actual del conocimiento, sobre las líneas y necesidades actuales de investigación. Se ha recabado información sobre los métodos modernos de cálculo y diseño, sobre los factores que influyen sobre la respuesta estructural, sobre técnicas de modelización y de optimización, al amparo de las indicaciones que algunas normativas actuales ofrecen sobre el tema. En esta tesis se ha desarrollado un procedimiento de modelización apoyado en el método de los elementos finitos implementado en el entorno MatLab; se han incluido aspectos claves tales como el comportamiento de segundo orden, la comprobación ante inestabilidad y la búsqueda del óptimo del coste de la estructura frente a estados límites, teniendo en cuenta las especificaciones del EC-3. También se ha modelizado la flexibilidad de las uniones y se ha analizado su influencia en la respuesta de la estructura y en el peso y coste final de la misma. Se han ejecutado algunos ejemplos de aplicación y se ha contrastado la validez del modelo con resultados de algunas estructuras ya analizadas en referencias técnicas conocidas. Se han extraído conclusiones sobre el proceso de modelización y de análisis, sobre la repercusión de la flexibilidad de las uniones en la respuesta de la estructura. El propósito es extraer conclusiones útiles para la etapa de proyecto. Una de las principales aportaciones del trabajo en su enfoque de optimización es la incorporación de una formulación de elementos no prismáticos con uniones semirrígidas en sus extremos. Se ha deducido una matriz de rigidez elástica para dichos elementos. Se ha comprobado su validez para abordar el análisis no lineal; para ello se han comparado los resultados con otros obtenidos tras aplicar otra matriz deducida analíticamente existente en la literatura y también mediante el software comercial SAP2000. Otra de las aportaciones de esta tesis es el desarrollo de un método de optimización del coste de pórticos metálicos planos de edificación en el que se tienen en cuenta aspectos tales como las imperfecciones, la posibilidad de incorporar elementos no prismáticos y la caracterización de las uniones semirrígidas, valorando la influencia de su flexibilidad sobre la respuesta de la estructura. Así, se han realizado estudios paramétricos para valorar la sensibilidad y estabilidad de las soluciones obtenidas, así como rangos de validez de las conclusiones obtenidas. This thesis deals with the problems of modelling, analysis and optimization of plane steel frames with regard to ultimate and serviceability limit states. The objective of this work is to present an organized sequential technique of discrete optimization for achieving the minimum cost of plane steel frames, taking into consideration the EC-3 specifications as well as including effects of the semi-rigid joints and non-prismatic elements in the design process. Likewise, an estimate of their influence on the final design is an aim of this work. The final objective is to draw practical conclusions which can be handful and easily applicable for a steel-structure project. An enormous amount of technical and scientific publications regarding steel frames is currently available, thus making the achievement of a comprehensive and updated knowledge a considerably hard task. In this work, a large variety of information has been gathered and classified, especially that related to current research lines and needs. Thus, the literature collected encompasses references related to state-of-the-art design methods, factors influencing the structural response, modelling and optimization techniques, as well as calculation and updated guidelines of some steel Design Codes about the subject. In this work a modelling procedure based on the finite element implemented within the MatLab programming environment has been performed. Several keys aspects have been included, such as second order behaviour, the safety assessment against structural instability and the search for an optimal cost considering the limit states according to EC-3 specifications. The flexibility of joints has been taken into account in the procedure hereby presented; its effects on the structural response, on the optimum weight and on the final cost have also been analysed. In order to confirm the validity and adequacy of this procedure, some application examples have been carried out. The results obtained were compared with those available from other authors. Several conclusions about the procedure that comprises modelling, analysis and design stages, as well as the effect of the flexibility of connections on the structural response have been drawn. The purpose is to point out some guidelines for the early stages of a project. One of the contributions of this thesis is an attempt for optimizing plane steel frames in which both non-prismatic beam-column-type elements and semi-rigid connections have been considered. Thus, an elastic stiffness matrix has been derived. Its validity has been tested through comparing its accuracy with other analytically-obtained matrices available in the literature, and with results obtained by the commercial software SAP2000. Another achievement of this work is the development of a method for cost optimization of plane steel building frames in which some relevant aspects have been taken in consideration. These encompass geometric imperfections, non-prismatic beam elements and the numerical characterization of semi-rigid connections, evaluating the effect of its flexibility on the structural response. Hence, some parametric analyses have been performed in order to assess the sensitivity, the stability of the outcomes and their range of applicability as well.

eCOMMONtech: plataforma sofrware para monitorización del balance de Gases de Efecto Invernadero en el Marco de Mecanismos de Desarrollo Limpio Forestales y Proyectos REDD+

La monitorización de las condiciones que debe cumplir un área forestal en proyectos MDL o REDD de manera tradicional, es decir, mediante mediciones y controles in situ, conlleva unos costes difíciles de asumir. Por ello, se ha planteado el desarrollo de una metodología capaz de integrar tecnologías orientadas a la realización de inventarios de carbono en áreas forestales de países en vías de desarrollo, mediante la utilización de diferentes tecnologías (sensorización ambiental, teledetección espacial, técnicas forestales, internet, etc.) que permiten determinar aquellos procedimientos más eficaces desde el punto de vista de la calidad y fiabilidad de la información obtenida y del coste/beneficio; analizando, las mejoras que suponen frente a los métodos tradicionales. Para ello, se desarrollan algoritmos y métodos de análisis necesarios para extraer las variables e indicadores medioambientales con el fin de realizar la monitorización de los ciclos de carbono en ámbitos forestales atribuibles a proyectos de absorciones de CO2.El resultado es la creación de una plataforma web que permite la monitorización remota y en tiempo real de inventarios de carbono a través de la integración de datos provenientes de sistemas de sensorización, imágenes tratadas con tecnologías de observación de la tierra y datos de campo.

Influencia de las imperfecciones en la carga crítica de estructuras de entramados planos

La aparición de errores o pequeñas diferencias entre la situación ideal prevista y la real durante la ejecución o proyecto de una obra constituye un fenómeno inherente a la limitación humana. Las Normas Oficiales suelen recoger en su articulado tolerancias máximas admisibles para estos errores deducidos, muchas veces, mediante criterios empíricos. Un planteamiento acaso más racional será indicar los valores límites de las imperfecciones de la obra en función del nivel de seguridad adoptado. En este artículo, se estudia desde este punto de vista, la influencia de los errores de ejecución en el valor de la carga crítica de la estructura. Evidentemente, las imperfecciones no deben ser limitadas únicamente por criterios de estabilidad global de la estructura, puesto que existen otras causas de colapso de la misma, que pueden verse afectados más seriamente por las imperfecciones y por consiguiente ser más exigentes en los niveles de las tolerancias máximas. Las imperfecciones que se consideran aquí, corresponden a valores relativamente pequeños producidos por un conjunto de causas simultáneas y por lo tanto son susceptibles de un tratamiento estadístico. Se excluyen de este estudio las equivocaciones o errores groseros cuya descripción matemática mediante métodos probabilistas no es adecuada. Se utilizarán aquí la conjunción de dos técnicas de cálculo -un modelo estructural basado en la teoría de la inestabilidad elástica lineal y un modelo probabilista con distribución gaussiana o uniforme- que se desarrollan de un modo numérico mediante el procedimiento de simulación de Monte-Carlo. Se comprende que la extensión del procedimiento de Monte-Carlo al análisis de otros tipos de modelos estructurales más refinados o bien que consideren otros mecanismos de colapso, así como distintas imperfecciones, es directo a causa del carácter eminentemente numérico del método

Integración en la plataforma 4mobile de mecanismos de atención en redes sociales

El presente texto se ha desarrollado como Proyecto Fin de Grado en la Escuela Técnica Superior de Ingeniería y Sistemas de Telecomunicación de la Universidad Politécnica de Madrid en colaboración con la Consultora Tecnológica everis. El mismo tiene como objetivo realizar un estudio de los requisitos necesarios para poder desplegar un servicio que permita ofrecer a las Operadoras Móviles Virtuales (OMV) soporte técnico a sus clientes a través de las redes sociales. Este módulo establecerá una nueva vía de comunicación entre el consumidor y el cliente de la OMV con la propia OMV, el cual se caracterizará por una alta accesibilidad y una gran rapidez en el servicio permitiendo de este modo, una mayor satisfacción del cliente con su operador móvil y por tanto de una nueva vía para conseguir la fidelización del mismo y la captación de nuevos clientes. Para ello, este proyecto se lleva a cabo sobre el sistema de información 4mobile el cual es comercializado por la mencionada Consultora Tecnológica. Este sistema consiste en una plataforma web la cual permite cubrir todos los procesos de negocio comunes que un OMV necesita gestionar. Es por ello, que el mencionado estudio se centra en la evaluación de los aspectos necesarios para la integración de un módulo de estas características dentro de la plataforma de 4mobile. Este módulo, estará basado en una herramienta software que permitirá gestionar el ciclo de vida completo del comentario realizado por el cliente a través de una red social, desde que es publicado por el cliente, hasta que se haya respondido al mismo y su solución sea considerada como satisfactoria por el cliente. Por ello, y de cara a definir correctamente esta herramienta, será necesario un detallado análisis el cual recoja diversos aspectos y que a lo largo del texto, será denominado como Plan de Marketing de Medios Sociales (PMMS). Estos aspectos versarán tanto de las necesidades tecnológicas para su mencionada integración, como de la serie de características funcionales que una solución basada en servicio técnico a través de las redes sociales deberá poseer con el objetivo de ofrecer un servicio técnico de calidad. Finalmente, estas funcionalidades y necesidades tecnológicas se expondrán en forma de propuesta a everis para su integración en la plataforma 4mobile junto con un análisis de diseño a alto nivel software de la solución a desarrollar. ABSTRACT. This text has been developed as Final Degree Project in the Escuela Técnica Superior de Ingeniería y Sistemas de Telecomunicación de la Universidad Politécnica de Madrid in collaboration with the Technology Consultant everis, aims to conduct a study of how to provide technical support through social networks and the evaluation of the integration of a social support module within the system platform 4mobile, which cover all business processes that need to manage a Mobile Virtual Network Operator (MVNO), which is marketed by above-mentioned consultant . This module will establish a new communication channel between the consumer and the client with the OMV and itself, which is characterized by high accessibility and great fast service. Thus will allow a higher customer satisfaction with him service and thus a new way to get the same loyalty and attract new customers To this end, this project is performed on the information system 4mobile which is marketed by the consulting part thereof everis. This system is a platform to cover all business processes that need to manage an MVNO. Therefore, the mentioned study focuses on the evaluation of the elements necessary for integrating module these features within the platform 4mobile. This module will be based on a software tool for managing the entire lifecycle of the comment made by the customer via social networking, from the moment it is published, until it has been respond and their solution is considered satisfactory for the customer. So, in order to properly define this tool, a detailed analysis which will be necessary to collect various aspects and throughout the text, it will be referred to as Plan of Social Media Marketing (PMMS). These aspects will be address both: the technological needs for the mentioned integration, and the number of functional characteristics-based service through social networking solution, as this must in order to provide a quality service Finally, these goals and requirements will be discussed as everis offer for integration into the platform 4mobile analysis along with a high level of software design to develop the solution.

Análisis y simulación numérica de problemas electrohidrodinámicos relacionados con el electrospray

La presente tesis es un estudio analítico y numérico del electrospray. En la configuración más sencilla, un caudal constante del líquido a atomizar, que debe tener una cierta conductividad eléctrica, se inyecta en un medio dieléctrico (un gas u otro líquido inmiscible con el primero) a través de un tubo capilar metálico. Entre este tubo y un electrodo lejano se aplica un voltaje continuo que origina un campo eléctrico en el líquido conductor y en el espacio que lo rodea. El campo eléctrico induce una corriente eléctrica en el líquido, que acumula carga en su superficie, y da lugar a un esfuerzo eléctrico sobre la superficie, que tiende a alargarla en la dirección del campo eléctrico. El líquido forma un menisco en el extremo del tubo capilar cuando el campo eléctrico es suficientemente intenso y el caudal suficientemente pequeño. Las variaciones de presión y los esfuerzos viscosos asociados al movimiento del líquido son despreciables en la mayor parte de este menisco, siendo dominantes los esfuerzos eléctrico y de tensión superficial que actúan sobre la superficie del líquido. En el modo de funcionamiento llamado de conochorro, el balance de estos esfuerzos hace que el menisco adopte una forma cónica (el cono de Taylor) en una región intermedia entre el extremo del tubo y la punta del menisco. La velocidad del líquido aumenta al acercarse al vértice del cono, lo cual propicia que las variaciones de la presión en el líquido generadas por la inercia o por la viscosidad entren en juego, desequilibrando el balance de esfuerzos mencionado antes. Como consecuencia, del vértice del cono sale un delgado chorro de líquido, que transporta la carga eléctrica que se acumula en la superficie. La acción del campo eléctrico tangente a la superficie sobre esta carga origina una tracción eléctrica que tiende a alargar el chorro. Esta tracción no es relevante en el menisco, donde el campo eléctrico tangente a la superficie es muy pequeño, pero se hace importante en el chorro, donde es la causa del movimiento del líquido. Lejos del cono, el chorro puede o bien desarrollar una inestabilidad asimétrica que lo transforma en una espiral (whipping) o bien romperse en un spray de gotas prácticamente monodispersas cargadas eléctricamente. La corriente eléctrica transportada por el líquido es la suma de la corriente de conducción en el interior del líquido y la corriente debida a la convección de la carga acumulada en su superficie. La primera domina en el menisco y la segunda en el chorro lejano, mientras que las dos son comparables en una región intermedia de transferencia de corriente situada al comienzo del chorro aunque aguas abajo de la región de transición cono-chorro, en la que el menisco deja de ser un cono de Taylor. Para un campo exterior dado, la acumulación de carga eléctrica en la superficie del líquido reduce el campo eléctrico en el interior del mismo, que llega a anularse cuando la carga alcanza un estado final de equilibrio. El tiempo característico de este proceso es el tiempo de relajación dieléctrica, que es una propiedad del líquido. Cuando el tiempo de residencia del líquido en la región de transición cono-chorro (o en otra región del campo fluido) es grande frente al tiempo de relajación dieléctrica, la carga superficial sigue una sucesión de estados de equilibrio y apantalla al líquido del campo exterior. Cuando esta condición deja de cumplirse, aparecen efectos de relajación de carga, que se traducen en que el campo exterior penetra en el líquido, a no ser que su constante dieléctrica sea muy alta, en cuyo caso el campo inducido por la carga de polarización evita la entrada del campo exterior en el menisco y en una cierta región del chorro. La carga eléctrica en equilibrio en la superficie de un menisco cónico intensifica el campo eléctrico y determina su variación espacial hasta distancias aguas abajo del menisco del orden de su tamaño. Este campo, calculado por Taylor, es independiente del voltaje aplicado, por lo que las condiciones locales del flujo y el valor de la corriente eléctrica son también independientes del voltaje en tanto los tamaños de las regiones que determinan estas propiedades sean pequeños frente al tamaño del menisco. Los resultados experimentales publicados en la literatura muestran que existe un caudal mínimo para el que el modo cono-chorro que acabamos de describir deja de existir. El valor medio y la desviación típica de la distribución de tamaños de las gotas generadas por un electrospray son mínimos cuando se opera cerca del caudal mínimo. A pesar de que los mecanismos responsables del caudal mínimo han sido muy estudiados, no hay aún una teoría completa del mismo, si bien su existencia parece estar ligada a la aparición de efectos de relajación de carga en la región de transición cono-chorro. En esta tesis, se presentan estimaciones de orden de magnitud, algunas existentes y otras nuevas, que muestran los balances dominantes responsables de las distintas regiones de la estructura asintótica de la solución en varios casos de interés. Cuando la inercia del líquido juega un papel en la transición cono-chorro, los resultados muestran que la región de transferencia de corriente, donde la mayor parte de la corriente pasa a la superficie, está en el chorro aguas abajo de la región de transición cono-chorro. Los efectos de relajación de carga aparecen de forma simultánea en el chorro y la región de transición cuando el caudal se disminuye hasta valores de un cierto orden. Para caudales aún menores, los efectos de relajación de carga se notan en el menisco, en una región grande comparada con la de transición cono-chorro. Cuando el efecto de las fuerzas de viscosidad es dominante en la región de transición, la región de transferencia de corriente está en el chorro pero muy próxima a la región de transición cono-chorro. Al ir disminuyendo el caudal, los efectos de relajación de carga aparecen progresivamente en el chorro, en la región de transición y por último en el menisco. Cuando el caudal es mucho mayor que el mínimo del modo cono-chorro, el menisco deja de ser cónico. El campo eléctrico debido al voltaje aplicado domina en la región de transferencia de corriente, y tanto la corriente eléctrica como el tamaño de las diferentes regiones del problema pasan a depender del voltaje aplicado. Como resultado de esta dependencia, el plano caudal-voltaje se divide en diferentes regiones que se analizan separadamente. Para caudales suficientemente grandes, la inercia del líquido termina dominando frente a las fuerzas de la viscosidad. Estos resultados teóricos se han validado con simulaciones numéricas. Para ello se ha formulado un modelo simplificado del flujo, el campo eléctrico y el transporte de carga en el menisco y el chorro del electrospray. El movimiento del líquido se supone casi unidireccional y se describe usando la aproximación de Cosserat para un chorro esbelto. Esta aproximación, ampliamente usada en la literatura, permite simular con relativa facilidad múltiples casos y cubrir amplios rangos de valores de los parámetros reteniendo los efectos de la viscosidad y la inercia del líquido. Los campos eléctricos dentro y fuera del liquido están acoplados y se calculan sin simplificación alguna usando un método de elementos de contorno. La solución estacionaria del problema se calcula mediante un método iterativo. Para explorar el espacio de los parámetros, se comienza calculando una solución para valores fijos de las propiedades del líquido, el voltaje aplicado y el caudal. A continuación, se usa un método de continuación que permite delinear la frontera del dominio de existencia del modo cono-chorro, donde el método iterativo deja de converger. Cuando el efecto de la inercia del líquido domina en la región de transición cono-chorro, el caudal mínimo para el cual el método iterativo deja de converger es del orden del valor estimado del caudal para el que comienza a haber efectos de relajación de carga en el chorro y el cono. Aunque las simulaciones no convergen por debajo de dicho caudal, el valor de la corriente eléctrica para valores del caudal ligeramente mayores parece ajustarse a las estimaciones para caudales menores, reflejando un posible cambio en los balances aplicables. Por el contrario, cuando las fuerzas viscosas dominan en la región de transición, se pueden obtener soluciones estacionarias para caudales bastante menores que aquel para el que aparecen efectos de relajación de carga en la región de transición cono-chorro. Los resultados numéricos obtenidos para estos pequeños caudales se ajustan perfectamente a las estimaciones de orden de magnitud que se describen en la memoria. Por último, se incluyen como anexos dos estudios teóricos que han surgido de forma natural durante el desarrollo de la tesis. El primero hace referencia a la singularidad en el campo eléctrico que aparece en la línea de contacto entre el líquido y el tubo capilar en la mayoría de las simulaciones. Primero se estudia en qué situaciones el campo eléctrico tiende a infinito en la línea de contacto. Después, se comprueba que dicha singularidad no supone un fallo en la descripción del problema y que además no afecta a la solución lejos de la línea de contacto. También se analiza si los esfuerzos eléctricos infinitamente grandes a los que da lugar dicha singularidad pueden ser compensados por el resto de esfuerzos que actúan en la superficie del líquido. El segundo estudio busca determinar el tamaño de la región de apantallamiento en un chorro de líquido dieléctrico sin carga superficial. En esta región, el campo exterior es compensado parcialmente por el campo que induce la carga de polarización en la superficie del líquido, de forma que en el interior del líquido el campo eléctrico es mucho menor que en el exterior. Una región como ésta aparece en las estimaciones cuando los efectos de relajación de carga son importantes en la región de transferencia de corriente en el chorro. ABSTRACT This aim of this dissertation is a theoretical and numerical analysis of an electrospray. In its most simple configuration, a constant flow rate of the liquid to be atomized, which has to be an electrical conductor, is injected into a dielectric medium (a gas or another inmiscible fluid) through a metallic capillary tube. A constant voltage is applied between this tube and a distant electrode that produces an electric field in the liquid and the surrounding medium. This electric field induces an electric current in the liquid that accumulates charge at its surface and leads to electric stresses that stretch the surface in the direction of the electric field. A meniscus appears on the end of the capillary tube when the electric field is sufficiently high and the flow rate is small. Pressure variations and viscous stresses due to the motion of the liquid are negligible in most of the meniscus, where normal electric and surface tension stresses acting on the surface are dominant. In the so-called cone-jet mode, the balance of these stresses forces the surface to adopt a conical shape -Taylor cone- in a intermediate region between the end of the tube and the tip of the meniscus. When approaching the cone apex, the velocity of the liquid increases and leads to pressure variations that eventually disturb the balance of surfaces tension and electric stresses. A thin jet emerges then from the tip of the meniscus that transports the charge accumulated at its surface. The electric field tangent to the surface of the jet acts on this charge and continuously stretches the jet. This electric force is negligible in the meniscus, where the component of the electric field tangent to the surface is small, but becomes very important in the jet. Far from the cone, the jet can either develop an asymmetrical instability named “whipping”, whereby the jet winds into a spiral, or break into a spray of small, nearly monodisperse, charged droplets. The electric current transported by the liquid has two components, the conduction current in the bulk of the liquid and the convection current due to the transport of the surface charge by the flow. The first component dominates in the meniscus, the second one in the far jet, and both are comparable in a current transfer region located in the jet downstream of the cone-jet transition region where the meniscus ceases to be a Taylor cone. Given an external electric field, the charge that accumulates at the surface of the liquid reduces the electric field inside the liquid, until an equilibrium is reached in which the electric field induced by the surface charge counters the external electric field and shields the liquid from this field. The characteristic time of this process is the electric relaxation time, which is a property of the liquid. When the residence time of the liquid in the cone-jet transition region (or in other region of the flow) is greater than the electric relaxation time, the surface charge follows a succession of equilibrium states and continuously shield the liquid from the external field. When this condition is not satisfied, charge relaxation effects appear and the external field penetrates into the liquid unless the liquid permittivity is large. For very polar liquids, the field due to the polarization charge at the surface prevents the external field from entering the liquid in the cone and in certain region of the jet. The charge at the surface of a conical meniscus intensifies the electric field around the cone, determining its spatial variation up to distances downstream of the apex of the order of the size of the meniscus. This electric field, first computed by Taylor, is independent of the applied voltage. Therefore local flow characteristics and the electric current carried by the jet are also independent of the applied voltage provided the size of the regions that determine these magnitudes are small compared with the size of the meniscus. Many experiments in the literature show the existence of a minimum flow rate below which the cone-jet mode cannot be established. The mean value and the standard deviation of the electrospray droplet size distribution are minimum when the device is operated near the minimum flow rate. There is no complete explanation of the minimum flow rate, even though possible mechanisms have been extensively studied. The existence of a minimum flow rate seems to be connected with the appearance of charge relaxation effects in the transition region. In this dissertation, order of magnitude estimations are worked out that show the dominant balances in the different regions of the asymptotic structure of the solution for different conditions of interest. When the inertia of the liquid plays a role in the cone-jet transition region, the region where most of the electric current is transfered to the surface lies in the jet downstream the cone-jet transition region. When the flow rate decreases to a certain value, charge relaxation effects appear simultaneously in the jet and in the transition region. For smaller values of the flow rate, charge relaxation effects are important in a region of the meniscus larger than the transition region. When viscous forces dominate in the flow in the cone-jet transition region, the current transfer region is located in the jet immediately after the transition region. When flow rate is decreased, charge relaxation effects appears gradually, first in the jet, then in the transition region, and finally in the meniscus. When flow rate is much larger than the cone-jet mode minimum, the meniscus ceases to be a cone. The electric current and the structure of the solution begin to depend on the applied voltage. The flow rate-voltage plane splits into different regions that are analyzed separately. For sufficiently large flow rates, the effect of the inertia of the liquid always becomes greater than the effect of the viscous forces. A set of numerical simulations have been carried out in order to validate the theoretical results. A simplified model of the problem has been devised to compute the flow, the electric field and the surface charge in the meniscus and the jet of an electrospray. The motion of the liquid is assumed to be quasi-unidirectional and described by Cosserat’s approximation for a slender jet. This widely used approximation allows to easily compute multiple configurations and to explore wide ranges of values of the governing parameters, retaining the effects of the viscosity and the inertia of the liquid. Electric fields inside and outside the liquid are coupled and are computed without any simplification using a boundary elements method. The stationary solution of the problem is obtained by means of an iterative method. To explore the parameter space, a solution is first computed for a set of values of the liquid properties, the flow rate and the applied voltage, an then a continuation method is used to find the boundaries of the cone-jet mode domain of existence, where the iterative method ceases to converge. When the inertia of the liquid dominates in the cone-jet transition region, the iterative method ceases to converge for values of the flow rate for which order-of-magnitude estimates first predict charge relaxation effects to be important in the cone and the jet. The electric current computed for values of the flow rate slightly above the minimum for which convergence is obtained seems to agree with estimates worked out for lower flow rates. When viscous forces dominate in the transition region, stationary solutions can be obtained for flow rates significantly smaller than the one for which charge relaxation effects first appear in the transition region. Numerical results obtained for those small values of the flow rate agree with our order of magnitude estimates. Theoretical analyses of two issues that have arisen naturally during the thesis are summarized in two appendices. The first appendix contains a study of the singularity of the electric field that most of the simulations show at the contact line between the liquid and the capillary tube. The electric field near the contact line is analyzed to determine the ranges of geometrical configurations and liquid permittivity where a singularity appears. Further estimates show that this singularity does not entail a failure in the description of the problem and does not affect the solution far from the contact line. The infinite electric stresses that appear at the contact line can be effectively balanced by surface tension. The second appendix contains an analysis of the size and slenderness of the shielded region of a dielectric liquid in the absence of free surface charge. In this region, the external electric field is partially offset by the polarization charge so that the inner electric field is much lower than the outer one. A similar region appears in the estimates when charge relaxation effects are important in the current transfer region.

Calle e Interior de manzana. Mecanismos de transición en la ciudad consolidada: Manzana 13, Rue des Suisses, Londres-Villarroel

Frente a la colmatación indiscriminada de los territorios colindantes, emerge la posibilidad de realizar acciones concretas sobre el tejido urbano consolidado, recuperando espacios olvidados para dotarlos de nuevas cualidades. El objetivo del trabajo es estudiar el potencial de transformación de los interiores de manzana en lugares de relación entre los ciudadanos a través del estudio de tres casos de vivienda colectiva y descubrir los mecanismos proyectivos que impulsan las nuevas relaciones entre el ámbito público de la ciudad y el colectivo de la vivienda, cambiando así la percepción de que nada interesante ocurre en áreas de transición entre la calle, la vivienda y el interior de la manzana. ABSTRACT. Facing the indiscriminate fouling adjacent territories, it rises the possibility of acting on the consolidated urban fabric, recovering forgotten spaces to give them new qualities. The aim of this work is to study the potential of transforming the block interiors in places of relationship between citizens through the study of three cases of collective housing and find projective mechanisms that drive the new relationship between the public realm of the city and the collective sphere of housing, thus changing the perception that nothing interesting happens in areas of transition between the street, the housing and the interior of the block.

Análisis de mecanismos de integración de herramientas para Suministro Continuo (Continuous Delivery)

Las empresas actuales centran sus esfuerzos en conseguir una buena calidad del software que desarrollan, haciendo uso de metodologías, procedimientos y estándares adecuados para poder desarrollar y lanzar un producto software de calidad con la suficiente certeza de obtener un menor coste y riesgo finales. Es por ello que se han establecido nuevos paradigmas de organización que proponen el uso de nuevas tecnologías y procesos para permitir una puesta en producción ágil y menos arriesgada. Este trabajo se centra en la implementación de una de estas técnicas, el Suministro Continuo (Continuous Delivery - CD). La implantación de esta técnica se basa en llevar a cabo la automatización de todos los pasos entre el desarrollo y la puesta en producción de software. Dentro de las tareas necesarias para implementar esta técnica, este trabajo lleva a cabo la identificación de las herramientas básicas que permiten el Suministro Continuo y el estudio de mecanismos de integración de todas estas herramientas, automatizando todo el proceso de desarrollo software desde la construcción del proyecto hasta su puesta en producción. Posibilitando así un proceso de Suministro Continuo, en el que se lleven a cabo despliegues del producto software diarios que permitan recoger realimentación más rápidamente de los usuarios y poder reaccionar ágilmente a los cambios de necesidad de estos.

Mecanismos implicados en las etapas iniciales de infección en la cancrosis de los cítricos provocada por Xanthomonas citri subsp. citri

La cancrosis o chancro bacteriano de los cítricos (CBC) causada por Xanthomonas citri subsp. citri (Xcc) y X. fuscans subsp. aurantifolii, afecta a un gran número de especies dentro de la familia de las rutáceas, especialmente cítricos. Esta enfermedad produce graves pérdidas económicas allí donde está presente, principalmente porque la comercialización de cítricos desde las zonas afectadas hacía zonas libres de cancrosis, está sujeta a fuertes medidas cuarentenarias. La cancrosis se encuentra distribuida a nivel mundial pero no se ha localizado ni en la Unión Europea ni en ningún área del Mediterráneo. Se han descrito tres tipos de cancrosis en función de la gama de huésped y de las características fenotípicas y genotípicas de las bacterias que las producen. La más extendida es la cancrosis tipo A producida por Xcc, dentro de la cual se distinguen los subtipos Aw y A*, originarios de Florida y Sudeste Asiático, respectivamente, que de forma natural solo son capaces de producir enfermedad en lima mejicana. En este trabajo se presentan estudios sobre mecanismos implicados en las primeras etapas de la infección, como la quimiotaxis y formación de biopelículas, en la cancrosis de los cítricos. La quimiotaxis es el proceso por el cual las bacterias se dirigen hacia zonas favorables para su supervivencia y desarrollo. Los perfiles quimiotácticos obtenidos frente a distintas fuentes de carbono, así como los estudios en relación al contenido de proteínas aceptoras de grupos metilo (MCPs), permitieron agrupar a las cepas de Xanthomonas estudiadas en este trabajo, de acuerdo a la enfermedad producida y a su gama de huésped. Todas las cepas mostraron quimiotaxis positiva frente a extractos de hoja y apoplasto de diferentes especies, sin embargo, Xcc 306, X. alfalfae subsp. citrumelonis (Xac) y X. campestris pv. campestris (Xc) manifestaron respuestas más específicas frente a extractos de apoplasto de hojas de naranjo dulce, lima y col china, respectivamente. Dicho resultado nos permite asociar el mecanismo de quimiotaxis con la capacidad de las cepas de Xanthomonas para colonizar estos huéspedes de forma específica. Las cepas estudiadas fueron capaces de realizar movimiento tipo swimming, twitching y sliding en distintos medios, siendo el movimiento swimming el único en el que se encontraron diferencias entre las cepas de Xcc con distinta gama de huésped. En este trabajo se ha estudiado además la formación de biopelículas en superficies bióticas y abióticas, un mecanismo importante tanto para la supervivencia en superficie vegetal como para el desarrollo de la infección. Las cepas de Xanthomonas estudiadas fueron capaces de formar biopelículas in vitro, siendo mayor en un medio que simula el apoplasto y que contiene una baja concentración de nutrientes en comparación con medios que contenían alta concentración de nutrientes. La formación de biopelículas en superficie vegetal se encontró relacionada, en las cepas patógenas de cítricos, con la capacidad para infectar un tejido o huésped determinado. Se han caracterizado algunos de los componentes de la matriz extracelular producida por Xcc, que compone hasta un 90% de las bipoelículas. Entre ellos destaca el ADN extracelular, que tiene un papel como adhesina en las primeras etapas de formación de biopelículas y estructural en biopelículas maduras. Además, se han identificado el pilus tipo IV como componente importante en las biopelículas, que también participa en motilidad. Finalmente, se han realizado estudios sobre la expresión de genes implicados en motilidad bacteriana y formación de biopelículas que han confirmado las diferencias existentes entre cepas de Xcc de amplia y limitada gama de huésped, así como el papel que juegan elementos como el pilus tipo IV o el flagelo en estos procesos. ABSTRACT Xanthomonas citri subsp. citri (Xcc) and X. fuscans subsp. aurantifolii are the causal agents of Citrus Bacterial Canker (CBC) which is one of the most important citrus diseases. CBC affects all Citrus species as well as other species from Rutaceae family. CBC produces strong economic losses; furthermore the commercialization of plants and fruits is restricted from infested to citrus canker free areas. The disease is worldwide distributed in tropical and subtropical areas, however it is not present in the European Union. Three types of CBC have been described according to the host range and phenotypic and genotypic characteristics. CBC type A caused by Xcc is he widest distributed. Within CBC A type two subtypes Aw and A* were described from Florida and Iran respectively, both infecting only Mexican lime. Herein mechanisms connected to early events in the citrus bacterial canker disease such as chemotaxis and biofilm formation, were studied. Chemotaxis allows bacteria to move towards the more suitable environments for its survival, host colonization and infection. Studies performed on citrus pathogenic Xanthomonas and X. campestris pv. campestris (Xc), a crucifer pathogen, have shown different chemotactic profiles towards carbon compound as well as different MCPs profile, which clustered strains according to host range and disease caused. Every strain showed positive chemotaxis toward leaf extracts and apoplastic fluids from sweet orange, Mexican lime and Chinese cabbage leaves. However, a more specific response was found for strains Xcc 306, X. alfalfae subsp. citrumelonis and Xc towards sweet orange, Mexican lime and Chinese cabbage apoplastic fluids, respectively. These results relate chemotaxis with the higher ability of those strains to specifically colonize their proper host. Xanthomonas strains studied were able to perform swimming, sliding and twitching motilities. The ability to swim was variable among CBC strains and seemed related to host range. Biofilm formation is an important virulence factor for Xcc because it allows a better survival onto the plant surface as well as facilitates the infection process. The studied Xanthomonas strains were able to form biofilm in vitro, on both nutrient rich and apoplast mimicking media, furthermore the biofilm formation by all the strains was higher in the apoplast mimicking media. The ability to form biofilm in planta by Xcc and Xac strains was dependent of the host and the tissue colonized. The wide host range CBC strain was able to form biofilm onto several citrus leaves and fruits, however the limited host range CBC strain produced biofilm solely onto Mexican lime leaves and fruits. Furthermore Xac strain, which solely infects leaves of young plants, was not able to develop biofilms on fruits. Some components of the extracellular matrix produced by Xcc strains have been characterized. Extracellular DNA acted as an adhesin at the very early stages of biofilm formation and as structural component of mature biofilm for citrus pathogenic Xanthomonas. Furthermore type IV pilus has been identified as a component of the extracellular matrix in biofilm and motility. Transcriptional studies of genes related with biofilm formation and motility have confirmed the differential behavior found among wide and limited host range CBC strains as well as the role of type IV pili and flagellum on those processes.

Simulación eficiente de mecanismos flexibles mediante formulaciones topológicas

El objetivo de esta Tesis es presentar un método eficiente para la evaluación de sistemas multi-cuerpo con elementos flexibles con pequeñas deformaciones, basado en métodos topológicos para la simulación de sistemas tan complejos como los que se utilizan en la práctica y en tiempo real o próximo al real. Se ha puesto un especial énfasis en la resolución eficiente de aquellos aspectos que conllevan mayor coste computacional, tales como la evaluación de las ecuaciones dinámicas y el cálculo de los términos de inercia. Las ecuaciones dinámicas se establecen en función de las variables independientes del sistema, y la integración de las mismas se realiza mediante formulaciones implícitas de index-3. Esta Tesis se articula en seis Capítulos. En el Capítulo 1 se realiza una revisión bibliográfica de la simulación de sistemas flexibles y los métodos más relevantes de integración de las ecuaciones diferenciales del movimiento. Asimismo, se presentan los objetivos de esta Tesis. En el Capítulo 2 se presenta un método semi-recursivo para la evaluación de las ecuaciones de los sistemas multi-cuerpo con elementos flexibles basado en formulaciones topológicas y síntesis modal. Esta Tesis determina la posición de cada punto del cuerpo flexible en función de un sistema de referencia flotante que se mueve con dicho cuerpo y de las amplitudes de ciertos modos de deformación calculados a partir de un mallado obtenido mediante el Método de Elementos Finitos. Se presta especial atención en las condiciones de contorno que se han de tener en cuenta a la hora de establecer las variables que definen la deformación del cuerpo flexible. El Capítulo 3 se centra en la evaluación de los términos de inercia de los sistemas flexibles que generalmente conllevan un alto coste computacional. Se presenta un método que permite el cálculo de dichos términos basado en el uso de 24 matrices constantes que pueden ser calculadas previamente al proceso de integración. Estas matrices permiten evaluar la matriz de masas y el vector de fuerzas de inercia dependientes de la velocidad sin que sea necesario evaluar la posición deformada de todos los puntos del cuerpo flexible. Se realiza un análisis pormenorizado de dichas matrices con el objetivo de optimizar su cálculo estableciendo aproximaciones que permitan reducir el número de dichos términos y optimizar aún más su evaluación. Se analizan dos posibles simplificaciones: la primera utiliza una discretización no-consistente basada en elementos finitos en los que se definen únicamente los desplazamientos axiales de los nodos; en la segunda propuesta se hace uso de una matriz de masas concentradas (Lumped Mass). Basándose en la formulación presentada, el Capítulo 4 aborda la integración eficiente de las ecuaciones dinámicas. Se presenta un método iterativo para la integración con fórmulas de index-3 basado en la proyección de las ecuaciones dinámicas según las variables independientes del sistema multi-cuerpo. El cálculo del residuo del sistema de ecuaciones no lineales que se ha de resolver de modo iterativo se realiza mediante un proceso recursivo muy eficiente que aprovecha la estructura topológica del sistema. Se analizan tres formas de evaluar la matriz tangente del citado sistema no lineal: evaluación aproximada, numérica y recursiva. El método de integración presentado permite el uso de distintas fórmulas. En esta Tesis se analizan la Regla Trapezoidal, la fórmula BDF de segundo orden y un método híbrido TR-BDF2. Para este último caso se presenta un algoritmo de paso variable. En el Capítulo 5 plantea la implementación del método propuesto en un programa general de simulación de mecanismos que permita la resolución de cualquier sistema multi-cuerpo definiéndolo mediante un fichero de datos. La implementación de este programa se ha realizado tanto en C++ como en Java. Se muestran los resultados de las formulaciones presentadas en esta Tesis mediante la simulación de cuatro ejemplos de distinta complejidad. Mediante análisis concretos se comparan la formulación presentada con otras existentes. También se analiza el efecto del lenguaje de programación utilizado en la implementación y los efectos de las posibles simplificaciones planteadas. Por último, el Capítulo 6 resume las principales conclusiones alcanzadas en la Tesis y las futuras líneas de investigación que con ella se abren. ABSTRACT This Thesis presents an efficient method for solving the forward dynamics of a multi-body sys-tem formed by rigid and flexible bodies with small strains for real-time simulation of real-life models. It is based on topological formulations. The presented work focuses on the efficient solution of the most time-consuming tasks of the simulation process, such as the numerical integration of the motion differential equations and in particular the evaluation of the inertia terms corresponding to the flexible bodies. The dynamic equations are formulated in terms of independent variables of the muti-body system, and they are integrated by means of implicit index-3 formulae. The Thesis is arranged in six chapters. Chapter 1 presents a review of the most relevant and recent contributions related to the modelization of flexible multi-body systems and the integration of the corresponding dynamic equations. The main objectives of the Thesis are also presented in detail. Chapter 2 presents a semi-recursive method for solving the equations of a multi-body system with flexible bodies based on topological formulations and modal synthesis. This Thesis uses the floating frame approach and the modal amplitudes to define the position of any point at the flexible body. These modal deformed shapes are obtained by means of the Finite Element Method. Particular attention has been taken to the boundary conditions used to define the deformation of the flexible bodies. Chapter 3 focuses on the evaluation of the inertia terms, which is usually a very time-consuming task. A new method based on the use of 24 constant matrices is presented. These matrices are evaluated during the set-up step, before the integration process. They allow the calculation of the inertia terms in terms of the position and orientation of the local coordinate system and the deformation variables, and there is no need to evaluate the position and velocities of all the nodes of the FEM mesh. A deep analysis of the inertia terms is performed in order to optimize the evaluation process, reducing both the terms used and the number of arithmetic operations. Two possible simplifications are presented: the first one uses a non-consistent approach in order to define the inertia terms respect to the Cartesian coordinates of the FEM mesh, rejecting those corresponding to the angular rotations; the second approach makes use of lumped mass matrices. Based on the previously presented formulation, Chapter 4 is focused on the numerical integration of the motion differential equations. A new predictor-corrector method based on index-3 formulae and on the use of multi-body independent variables is presented. The evaluation of the dynamic equations in a new time step needs the solution of a set on nonlinear equations by a Newton-Raphson iterative process. The computation of the corresponding residual vector is performed efficiently by taking advantage of the system’s topological structure. Three methods to compute the tangent matrix are presented: an approximated evaluation that considers only the most relevant terms, a numerical approach based on finite differences and a recursive method that uses the topological structure. The method presented for integrating the dynamic equations can use a variety of integration formulae. This Thesis analyses the use of the trapezoidal rule, the 2nd order BDF formula and the hybrid TR-BDF2 method. A variable-time step strategy is presented for the last one. Chapter 5 describes the implementation of the proposed method in a general purpose pro-gram for solving any multibody defined by a data file. This program is implemented both in C++ and Java. Four examples are used to check the validity of the formulation and to compare this method with other methods commonly used to solve the dynamic equations of multi-body systems containing flexible bodies. The efficiency of the programming methodology used and the effect of the possible simplifications proposed are also analyzed. Chapter 6 summarizes the main Conclusions obtained in this Thesis and the new lines of research that have been opened.