Métodos simplificados para cálculo de la respuesta sísmica, no lineal, de puentes

Tradicionalmente, el cálculo dinámico no-lineal de estructuras ha sido una tarea ardua para el ingeniero ya que comporta cierta especialización en el ingeniero habituado al cálculo estático. El necesario conocimiento de técnicas de modelado en dinámica, de tiempos de integración, etc., y la dificultad de reconocer la influencia en el comportamiento de ciertas modificaciones de diseño, hace útil la existencia de un método más simplificado y transparente que permita análisis lineal, imprescindible por ejemplo en procesos de reacondicionamiento sísmico. En este artículo se presenta un método simplificado para analizar dinámicamente puentes a través de acelerogramas y cuyas principales ventajas son su sencillez, rápidez y bajo coste, permitiendo recoger el comportamiento dinámico y determinar el grado de plastificación, esfuerzos y desplazamientos con suficiente aproximación en puentes de altas complejidades.

Los plásticos en el diseño de la maquinaria agrícola. Parte IV procesos de transformación

Existen técnicas muy variadas de moldeo: extrusión, inyección, moldeo, prensado, etc. Si bien, parte de estos procesos pueden emplearse indistintamente en la elaboración de elementos con materiales termoplásticos o termoestables, otros únicamente están indicados para la fabricación de elementos con uno de estos tipos de plástico

Intereses exógenos y resistencia ciudadana en la transformación físicosocial de un asentamiento auto-producido de Dakar (2005-2011)

La presión sobre el acceso a la tierra en la península de Dakar, Senegal, ha generado un proceso de transformación físico-social del mayor asentamiento auto-producido de la periferia urbana de Dakar (Pikine Irregular Sur) en el periodo 2005-2011. Esta investigación se propone analizar en qué medida las acciones llevadas a cabo: i) satisfacen o no el sistema integral de necesidades de los afectados (Demanda); ii) favorecen la implicación de los afectados o, por el contrario, fortalecen las estructuras de poder dominantes (Oferta), y iii) promueven o no la autonomía económica, la sostenibilidad ambiental, la equidad social y la legitimidad de los actores que participan (Entorno). La investigación se ha llevado a cabo por el método de estudio de casos, mediante un proceso de recogida de información participativo, articulado de abajo hacia arriba, en el que han participado 196 personas; la información recogida se ha procesado mediante la utilización de dos herramientas de análisis de realidades complejas: flujo de programa y sistema funcional. Los resultados del estudio de los casos 1 y 2 muestran que las acciones se han llevado a cabo en beneficio de los intereses de los mercados y de los actores que toman las decisiones y gestionan los recursos, incrementando la dependencia económica, empobreciendo a los directamente afectados y fragmentando su tejido social, promoviendo un desarrollo urbano desequilibrado, disperso y segregado, y deslegitimando a las instituciones involucradas. En el estudio del caso 3, en cambio, se verifica que la acción se ha llevado a cabo en beneficio de los intereses de los directamente afectados, utilizando recursos endógenos, fortaleciendo su tejido social, promoviendo un desarrollo urbano concentrado e inclusivo, e incrementando la legitimidad de las instituciones involucradas. El análisis cruzado de los resultados muestra que la pobreza, la exclusión social y la precariedad habitacional, en el contexto específico de la periferia de Dakar, no son un problema de falta de recursos, sino de falta de ética por parte de quienes gestionan y financian los procesos de transformación físico-social, y de falta de voluntad por parte de las autoridades de implicar a los ciudadanos, y en particular a los directamente afectados, en los procesos de toma de decisiones de las acciones y los programas que les afectan. Finalmente se concluye que no se trata de luchar contra la pobreza, sino de luchar contra los sistemas injustos que generan la pobreza y la exclusión. ABSTRACT The pressure on the access to land in the area of Dakar, Senegal, has generated a physical/social transformation process in the largest self-produced settlement in the urban periphery of Dakar (South Irregular Pikine), in the period 2005-2011. This research aims to analyze to what extent the actions carried out: i) satisfy the integral needs of the affected population (Demand), ii) encourage the involvement of those affected or otherwise strengthen dominant power structures (Supply), and iii) promote or inhibit economic autonomy, environmental sustainability, social equity and legitimacy of the actors involved (Environment). The research was conducted by the case study method. 196 people participated in a bottom-up participatory process of information gathering. The information collected was processed using two custom tools aimed to analyse complex systems: the program flow and the functional system. The results of case studies 1 and 2 show that the actions were carried out in the best interest of the markets and the actors who made decisions and managed resources, increasing economic dependence, impoverishing those directly affected and fragmenting their social network by promoting an unbalanced, dispersed and segregated urban development, and de-legitimizing the institutions involved in the process. On the contrary, in the case study 3 it is verified that the action has been carried out on behalf of the interests of those directly affected, using endogenous resources, strengthening the social network, promoting a focused and inclusive urban development, and increasing legitimacy of the institutions involved. The comparative analysis of the results shows that poverty, social exclusion and precarious housing in South Irregular Pikine are not a problem of lack of resources but of unethical behaviour of those who manage and finance physical and social transformation processes, as well as unwillingness on the part of the authorities to involve citizens in the decision making processes for the actions and programs that affect them. Finally we conclude that the enemy to fight against is not poverty but unfair systems that generate poverty and exclusion. RÉSUMÉ La pression sur le foncier dans la presqu'île de Dakar, au Sénégal, a généré un processus de transformation physico-sociale du plus grand établissement auto-produit de la périphérie urbaine de Dakar (Pikine Irrégulier Sud), durant la période 2005-2011. Cette recherche vise à analyser dans quelle mesure les actions menées: i) répondent ou non à l’ensemble des besoins intégrales des personnes concernés (Demande), ii) encouragent l'implication des concernées ou au contraire renforcent les structures de pouvoir dominant (Offre), et iii) favorisent ou non l'autonomie économique, la durabilité environnementale, l'équité sociale et la légitimité des acteurs concernés (Contexte). La recherche a été menée par la méthode des études de cas, par un processus de collecte d'information participative articulée du bas vers le haut dans lequel ont participé 196 personnes. Les informations recueillies ont été traitées à l'aide de deux outils d'analyse de réalités complexes: le système fonctionnel et le flux de programme. Les résultats de l'étude de cas 1 et 2 montrent que les actions ont été menées au profit des intérêts des marchés et des acteurs qui prennent les décisions et gèrent les ressources, augmentant la dépendance économique, appauvrissant ceux qui sont directement concernés et fragmentant leur tissu social, encourageant un développement urbain déséquilibré, dispersé et ségrégé, et délégitimant les institutions concernées. Dans l'étude de cas 3, cependant, il est vérifié que l'action a été menée au profit des intérêts de ceux qui sont directement concernés, en utilisant des ressources endogènes, en renforçant son tissu social, en encourageant un développement urbain concentré et inclusif, et en augmentant la légitimité des institutions concernées. L'analyse comparée des résultats montre que la pauvreté, l'exclusion sociale et le logement précaire, dans le contexte spécifique de la banlieue de Dakar, ne sont pas un problème de manque de ressources, mais de manque d'éthique de la part de ceux qui gèrent et financent les processus de transformation physico-sociale, et de manque de volonté de la part des autorités d'impliquer les citoyens, en particulier ceux qui sont directement concernés, dans le processus de prise de décision des actions et des programmes qui les impliquent. Finalement, on conclut qu'il ne s'agit pas de lutter contre la pauvreté, mais de lutter contre les systèmes injustes qui produisent la pauvreté et l'exclusion.

Polynomial approximation using particle swarm optimization of lineal enhanced neural networks with no hidden layers.

This paper presents some ideas about a new neural network architecture that can be compared to a Taylor analysis when dealing with patterns. Such architecture is based on lineal activation functions with an axo-axonic architecture. A biological axo-axonic connection between two neurons is defined as the weight in a connection in given by the output of another third neuron. This idea can be implemented in the so called Enhanced Neural Networks in which two Multilayer Perceptrons are used; the first one will output the weights that the second MLP uses to computed the desired output. This kind of neural network has universal approximation properties even with lineal activation functions. There exists a clear difference between cooperative and competitive strategies. The former ones are based on the swarm colonies, in which all individuals share its knowledge about the goal in order to pass such information to other individuals to get optimum solution. The latter ones are based on genetic models, that is, individuals can die and new individuals are created combining information of alive one; or are based on molecular/celular behaviour passing information from one structure to another. A swarm-based model is applied to obtain the Neural Network, training the net with a Particle Swarm algorithm.

Implementación de una práctica virtual de biotecnología: la transformación bacteriana

En este documento se realiza la implementación de la fase "Transformación Bacteriana" en el "Laboratorio Virtual de Biotecnología Agroforestal" de la Universidad Politécnica de Madrid. Esta fase representa una continuación de un trabajo previo, en el que se implementó el laboratorio virtual y se diseñó una arquitectura para el desarrollo de las fases posteriores. La Transformación Bacteriana es la tercera fase del proceso de modificación genética de un chopo para dotarle de resistencia frente a ciertos hongos y tiene como objetivo la introduccion de un plásmido modificado genéticamente en la bacteria Agrobacterium tumefaciens. Para el desarrollo de esta fase se determinaron las acciones que debían de ser agregadas al tutor automático y, además, se programaron los scripts de ciertos objetos 3D ya existentes y se modelaron en 3D nuevos instrumentos y maquinarias necesarios para la realización de la práctica. Por otra parte, luego de la etapa de pruebas de las dos primeras fases, se concluyó que se debían de resolver algunos problemas de usabilidad del visor Firestorm. Gracias a que este visor es un proyecto de código abierto, fue posible corregir los problemas identificados. La resolución de estos problemas se explica como parte del presente trabajo.

Adaptive control system based on lineal control theory for the path-following problem of a car-like mobile robot

The objective of this paper is to design a path following control system for a car-like mobile robot using classical linear control techniques, so that it adapts on-line to varying conditions during the trajectory following task. The main advantages of the proposed control structure is that well known linear control theory can be applied in calculating the PID controllers to full control requirements, while at the same time it is exible to be applied in non-linear changing conditions of the path following task. For this purpose the Frenet frame kinematic model of the robot is linearised at a varying working point that is calculated as a function of the actual velocity, the path curvature and kinematic parameters of the robot, yielding a transfer function that varies during the trajectory. The proposed controller is formed by a combination of an adaptive PID and a feed-forward controller, which varies accordingly with the working conditions and compensates the non-linearity of the system. The good features and exibility of the proposed control structure have been demonstrated through realistic simulations that include both kinematics and dynamics of the car-like robot.

Criptoanálisis de un criptosistema de dos canales basado en una función no lineal caótica

En este artículo se realiza el criptoanálisis de un criptosistema caótico basado en el sistema de Lorenz, que utiliza dos canales de comunicación y una función no lineal caótica. Se demuestra que el sistema propuesto es inseguro ya que los valores de los parámetros del sistema pueden determinarse con gran precisión utilizando un receptor intruso.

La silvicultura como primera operación de transformación de la madera Forestry as the first operation of wood transformation

Se analiza la influencia de las principales actuaciones selvicolas en la calidad de la madera. Turno, Espaciamiento, Poda, Abonado, Riego y Tratamiento contra plagas.

Sensor magneto-óptico de corriente con supresión de birrefringencia lineal mediante espejo de Faraday

Fiber optic sensors have some advantages in subjects related with electrical current and magnetic field measurement. In spite of the optical fiber utilization advantages we have to take into account undesirable effects, which are present in real non-ideal optical fibers. In telecommunication and sensor application fields the presence of inherent and induced birefringence is crucial. The presence of birefringence may cause an undesirable change in the polarization state. In order to compensate the linear birefringence a promising method has been chosen. This method employs orthogonal polarization conjugation in the back propagation direction of the light wave in the fiber. A study and a simulation of an experimental setup are realized with the advantage of a significant sensitivity improvement.

Nuevo algoritmo de generalización lineal y ortogonalización de construcciones en un entorno urbano.

Resumen: La generalización cartográfica es el proceso que permite la obtención de cartografía derivada a partir de la generada mediante el registro de información geográfica. Normalmente, la generalización implica una reducción de escala, aunque podría ser una reducción de la cantidad de información sin reducir la escala, debido a un cambio de objetivo de la cartografía. La dificultad de la generalización aumenta cuanto menor es la escala destino y mayor la diferencia entre ésta y la escala origen. Estas dificultades se concentran en aquellas entidades que representan objetos modificados por el hombre, es decir, principalmente en los núcleos urbanos, ya que deben preservar ciertas propiedades en su geometría una vez generalizados. En este artículo, se propone un algoritmo para la generalización lineal de núcleos urbanos, basándose en las características constructivas y relacionales de las edificaciones que los conforman. Abstract: Mapping generalization is the process which derived maps are obtained from others generated by recording geographic information. Usually, generalization involves a reduction of scale, although it could be a reduction in the amount of information without reducing the scale, due to a change in target mapping. The smaller target scale, the more difficult is the generalization. Also, the more difference between origin scale and target scale, the more difficult is the generalization. The difficulties are concentrated in entities that represent objects modified by man, mainly in urban areas, as they must preserve certain geometry properties once generalized. In this paper, an algorithm is proposed for linear generalization of urban areas, based on constructive and relational characteristics of the buildings that shape them.

Caos en evolución no lineal de ondas de Alfven usando un modelo resistivo

Interacciones no lineales de ondas de Alfvén existen tanto para plasmas en el espacio como en laboratorios, con efectos que van desde calentamiento hasta conducción de corriente. Un ejemplo de emisión de ondas de Alfvén en ingeniería aeroespacial aparece en amarras espaciales. Estos dispositivos emiten ondas en estructuras denominadas “Alas de Alfvén”. La ecuación derivada no lineal de Schrödinger (DNLS) posee la capacidad de describir la propagación de ondas de Alfvén de amplitud finita circularmente polarizadas en un plasma frío. En esta investigación, dicha ecuación es truncada con el objetivo de explorar el acoplamiento coherente, débilmente no lineal y cúbico de tres ondas cerca de resonancia. De las tres ondas, una es linealmente inestable y las otras dos son amortiguadas. Por medio de la utilización de este modelo se genera un flujo 4D formado por tres amplitudes y una fase relativa. En trabajos anteriores se analizó la transición dura hacia caos en flujos 3D.2005). Se presenta en este artículo un análisis teórico y numérico del comportamiento del sistema cuando la tasa de crecimiento de la onda inestable es muy próxima a cero y considerando amortiguamiento resistivo, es decir se satisface una relación cuadrática entre amortiguamientos y números de onda. Al igual que en los trabajos anteriores, se ha encontrado que sin importar cuan pequeña es la tasa de crecimiento de la onda inestable existe un dominio paramétrico, en el espacio de fase, donde aparecen oscilaciones caóticas que están ausentes para una tasa de crecimiento nula. Sin embargo diagramas de bifurcación y dominios de estabilidad presentan diferencias con respecto a lo estudiado anteriormente.

Modelo lineal para la predicción del daño producido por impacto mecánico en pera blanquilla.

Frutos de pera cv. "Blanquilla" de distintos grados de madurez se ensayaron a impacto mecánico con un indentador esférico de 50.6 g de masa desde tres alturas de caída ( 4,6 y 10 cm) y se determinó un modelo lineal para cada uno de los parámetros que definían la magulladura producida (profundidad y diámetro) y en cada uno de los dos grupos de frutos estudiados: no senescentes y senescentes en función de las variables definitorias del fenómeno del impacto y del estado del fruto que se consideraron mas representativas. Se observa que los modelos lineales determinados explican la magulladura en menos del 57% siendo la energía máxima de impacto la que adquiere mayor importancia en ambos grupos. En los senescentes, el estado fisiológico (expresado por el índice de madurez E/A y la acidez titulable) posee una relevancia mayor que el resto de variables y las deformaciones (indicativas de un comportamiento plástico del producto) toman gran importancia.

Medida de la conductividad térmica con el método de la aguja térmica, basado en la fuente lineal de calor transitorio, para su aplicación en los cerramientos de adobes y bloques de tierra comprimida

Existe un creciente interés internacional por el ahorro energético y la sostenibilidad en la edificación con importantes repercusiones en la Arquitectura. La inercia térmica es un parámetro fundamental para poder valorar energéticamente un edificio en condiciones reales. Para ello es necesario cambiar el enfoque tradicional de transmisión de calor en régimen estacionario por otro en régimen dinámico en el que se analizan las ondas térmicas y el flujo de calor oscilante que atraviesan los cerramientos. Los parámetros que definen la inercia térmica son: el espesor, la difusividad y el ciclo térmico. A su vez la difusividad está determinada por la conductividad térmica, la densidad y el calor específico del material. De estos parámetros la conductividad es el más complejo, variable y difícil de medir, especialmente en los cerramientos de tierra debido a su heterogeneidad y complejidad higrotérmica. En general, los métodos de medida de la conductividad o transmitancias en los paramentos presentan inconvenientes a la hora de medir un edificio construido con tierra: dificultades de implementación, el elevado coste o la fiabilidad de los resultados, principalmente. El Método de la Aguja Térmica (MAT) se basa en el principio de la evolución en el tiempo del calor emitido por una fuente lineal al insertarse en el seno de un material. Se ha escogido este método porque resulta práctico, de bajo coste y de fácil aplicación a gran escala pero tiene serios problemas de fiabilidad y exactitud. En esta tesis se desarrolla un método de medida de la conductividad térmica para Piezas de Albañilería de Tierra Cruda en laboratorio basado en el MAT, se mejora su fiabilidad, se analiza su incertidumbre, se compara con otros métodos de referencia y se aplica en adobes, Bloques de Tierra Comprimida y probetas de tierra estabilizada con distintas proporciones de paja. Este método servirá de base a una posterior aplicación in situ. Finalmente se proponen modelos matemáticos para mejorar la exactitud del dispositivo utilizado y para la estimación de la conductividad de cerramientos de tierra en función de su densidad. Con los resultados obtenidos se analizan las posibilidades de amortiguación y retardo de las ondas térmicas y capacidad de almacenaje de energía de los cerramientos en función de su densidad y humedad. There is growing international interest in energy saving and sustainability in buildings with significant impact on Architecture. Thermal inertia is a key parameter to assess energy in buildings in real conditions. This requires changing the traditional approach to heat transfer in steady state by another in dynamic regime which analyzes the thermal waves and oscillating heat flux passing through the external walls. The parameters defining the thermal inertia are: the thickness, the diffusivity and the thermal cycle. In turn, the diffusivity is determined by the thermal conductivity, density and specific heat of the material. Of these parameters, thermal conductivity is the most complex, variable and difficult to measure, especially in earth walls due to their heterogeneity and hygrothermal complexity. In general, the methods of measurement of conductivity and transmittance in walls have drawbacks when measuring a building with earth: implementation difficulties, high cost, or reliability of the results, mainly. The Thermal Needle Procedure (TNP) is based on the principle of evolution in time of heat from a line source when inserted within a material. This method was chosen because it is a practical, low cost and easy to implement on a large scale but has serious problems of reliability and accuracy. This thesis develops a laboratory method for measuring the thermal conductivity of Masonry Units Unfire Earth-based based on TNP, its uncertainty is analyzed, compared to other reference methods and applies in adobes, Compressed Earth Blocks and stabilized soil specimens with different proportions of straw. This method will form the basis of a subsequent application in situ. Finally, mathematical models are proposed to improve the accuracy of the device used, and to estimate the conductivity of earth enclosures depending on its density. With the results obtained earth enclosures are analyzed to estimate their possibilities of delay and buffer of termal waves and energy storage capacity according to their density and moisture.

La bóveda tabicada en España en el siglo XIX: la transformación de un sistema constructivo

Las bóvedas tabicadas experimentan una importante evolución, casi una transformación, en España a lo largo del s. XIX. Al inicio, son un sistema constructivo tradicional y codificado durante al menos 3 siglos, con usos, formas y tamaños limitados. 100 años después se emplean de manera generalizada en nuevos tipos edificatorios, con soluciones constructivas novedosas, tamaños mucho mayores y una gran libertad formal. Esta tesis doctoral indaga en el porqué de esta transformación, desde el punto de vista de la técnica, así como en sus consecuencias, que llegan hasta la actualidad. El texto agrupa en tres partes: La primera, a la que he llamado «La construcción tabicada tradicional» sirve de punto de partida. En ella se ponen en contexto las bóvedas tabicadas antes del s. XIX, como aparecen en España y cual es su desarrollo y distribución territorial hasta finalizar el s. XVIII. En esta parte se ha estudiado: -la aparición de las bóvedas tabicadas en España, en los últimos años del s. XIV, y como se convierten en un sistema constructivo cada vez más frecuente y consolidado, empezando por el territorio de la Corona de Aragón (actuales Valencia, Cataluña y Aragón) y extendiéndose desde ahí a Castilla y Andalucía. -la difusión del sistema tabicado en los tratados de arquitectura españoles de los siglos XVII y XVIII: Fray Lorenzo de San Nicolás, Tosca, Briguz y Brú, García Berruguilla, Plo y Camín; también a través de documentos más prácticos: contratos de obra y pleitos. -la expansión de la construcción tabicada por Francia desde España. Como allí es aceptada por académicos y tratadistas; como vuelve a España influenciada por esos tratadistas, los más importantes de la época. Y como los tratadistas españoles reproducen las ideas francesas sobre la construcción tabicada en sus escritos. Este aspecto es importante ya que, a pesar de que la construcción tabicada tenía una tradición mucho mayor aquí, las ideas francesas sobre su comportamiento se imponen a las españolas incluso en España y esto influye en su desarrollo posterior. La segunda parte indaga en los aspectos constructivos del desarrollo de las bóvedas tabicadas a lo largo del siglo XIX. Dentro de ella se describe: -como se mantiene la construcción tabicada tradicional a lo largo del siglo XIX, conviviendo con una línea que avanza en la búsqueda de una nueva forma de construir, separándose cada vez más de la inicial. Para mostrar esta pervivencia, se estudian en esta parte tres documentos que se mueven en la línea tradicional, uno de principio de siglo, otro de los años centrales y otro de final del XIX. -la aparición de dos nuevos materiales: el cemento y el hierro, que en pocos años pasan de tener un uso accesorio en la construcción a fabricarse y utilizarse masivamente. La unión de estos dos materiales con los tradicionales en la construcción de bóvedas tabicadas, yeso y ladrillo, será decisiva en su evolución. -la creación de dos tipos arquitectónicos nuevos: los edificios industriales y la vivienda masiva. En ambos tipos tiene aplicación la construcción tabicada. Además, estos nuevos tipos son muy importantes en su evolución, ya que sus condicionantes específicos: rapidez de ejecución en ambos casos y búsqueda de luces cada vez mayores y de estructuras incombustibles en el segundo fuerzan al uso de los nuevos materiales disponibles (cemento, hierro) y al proyecto de nuevas formas y soluciones. En este aspecto me he centrado en Cataluña, la primera región industrializada de España y en la que se aplica la bóveda tabicada a los nuevos tipos, especialmente a la arquitectura industrial. Y la tercera parte estudia los aspectos estructurales del desarrollo de las bóvedas tabicadas a lo largo del siglo XIX. Se estudia aquí: -la aplicación a las estructuras tabicadas de las teorías de arcos y bóvedas que se desarrollan desde el siglo XVIII. Las bóvedas tabicadas, cuyas dimensiones se obtenian mediante reglas de proporción, se calculan sucesivamente mediante lineas de empuje, análisis de membrana y también se les intenta aplicar la teoría elástica, sin llegar a conseguir resultados válidos. Siguiendo el marco teórico del moderno análisis límite de fábrica, se concluye que es el más adecuado para estudiar estas bóvedas, como cualquier otras construidas con fábrica. -un estudio exhaustivo de los ensayos de resistencia que se realizan sobre estructuras tabicadas a lo largo de los siglos XIX y la primera mitad del XX. La razón de estos ensayos es la nueva necesidad, aparecida con la mentalidad científica, de ofrecer garantías técnicas y numéricas sobre los sistemas constructivos y estructurales. La tesis concluye que las bóvedas tabicadas aparecen en España en el siglo XIV en la zona de Levante; que se extienden por casi todo el territorio siendo un sistema codificado en los últimos años del siglo XVI y que empiezan a ser recogidas en los tratados ya en el siglo XVII; hasta la mitad del siglo XIX la forma en que se construyen y los edificios para los que se utiliza permanecen invariables, sin más novedad que la exportación del sistema constructivo a Francia y la reintroducción en España a través de los influyentes tratados franceses. En el siglo XIX confluyen varios motivos por los que las bóvedas tabicadas se transforman: la generalización del uso del hierro y el cemento; la construcción de grandes edificios industriales y la aparición de la mentalidad científica, con el desarrollo de la teoría de estructuras. Los nuevos edificios industriales requieren grandes luces y nuevas formas; los nuevos materiales (especialmente el hierro) permiten estos avances. La teoría de estructuras se aplicará a estas nuevas bóvedas; en conjunto con los datos recogidos en los ensayos realizados dará seguridad a arquitectos y constructores para construir mayores tamaños y formas novedosas. A pesar de esta transformación, la mejor manera de enfrentarse al cálculo de una bóveda tabicada, bien sea para diseñarla en un edificio nuevo como para comprobar su validez en uno ya construido, es englobarla, como cualquier otra fábrica, dentro del análisis límite. Es el método más seguro y a la vez, el más sencillo de todos los disponibles actualmente.

Simulación en tiempo real de vehículos industriales con modelos multicuerpo de gran complejidad

El objetivo de esta Tesis ha sido la consecución de simulaciones en tiempo real de vehículos industriales modelizados como sistemas multicuerpo complejos formados por sólidos rígidos. Para el desarrollo de un programa de simulación deben considerarse cuatro aspectos fundamentales: la modelización del sistema multicuerpo (tipos de coordenadas, pares ideales o impuestos mediante fuerzas), la formulación a utilizar para plantear las ecuaciones diferenciales del movimiento (coordenadas dependientes o independientes, métodos globales o topológicos, forma de imponer las ecuaciones de restricción), el método de integración numérica para resolver estas ecuaciones en el tiempo (integradores explícitos o implícitos) y finalmente los detalles de la implementación realizada (lenguaje de programación, librerías matemáticas, técnicas de paralelización). Estas cuatro etapas están interrelacionadas entre sí y todas han formado parte de este trabajo. Desde la generación de modelos de una furgoneta y de camión con semirremolque, el uso de tres formulaciones dinámicas diferentes, la integración de las ecuaciones diferenciales del movimiento mediante métodos explícitos e implícitos, hasta el uso de funciones BLAS, de técnicas de matrices sparse y la introducción de paralelización para utilizar los distintos núcleos del procesador. El trabajo presentado en esta Tesis ha sido organizado en 8 capítulos, dedicándose el primero de ellos a la Introducción. En el Capítulo 2 se presentan dos formulaciones semirrecursivas diferentes, de las cuales la primera está basada en una doble transformación de velocidades, obteniéndose las ecuaciones diferenciales del movimiento en función de las aceleraciones relativas independientes. La integración numérica de estas ecuaciones se ha realizado con el método de Runge-Kutta explícito de cuarto orden. La segunda formulación está basada en coordenadas relativas dependientes, imponiendo las restricciones por medio de penalizadores en posición y corrigiendo las velocidades y aceleraciones mediante métodos de proyección. En este segundo caso la integración de las ecuaciones del movimiento se ha llevado a cabo mediante el integrador implícito HHT (Hilber, Hughes and Taylor), perteneciente a la familia de integradores estructurales de Newmark. En el Capítulo 3 se introduce la tercera formulación utilizada en esta Tesis. En este caso las uniones entre los sólidos del sistema se ha realizado mediante uniones flexibles, lo que obliga a imponer los pares por medio de fuerzas. Este tipo de uniones impide trabajar con coordenadas relativas, por lo que la posición del sistema y el planteamiento de las ecuaciones del movimiento se ha realizado utilizando coordenadas Cartesianas y parámetros de Euler. En esta formulación global se introducen las restricciones mediante fuerzas (con un planteamiento similar al de los penalizadores) y la estabilización del proceso de integración numérica se realiza también mediante proyecciones de velocidades y aceleraciones. En el Capítulo 4 se presenta una revisión de las principales herramientas y estrategias utilizadas para aumentar la eficiencia de las implementaciones de los distintos algoritmos. En primer lugar se incluye una serie de consideraciones básicas para aumentar la eficiencia numérica de las implementaciones. A continuación se mencionan las principales características de los analizadores de códigos utilizados y también las librerías matemáticas utilizadas para resolver los problemas de álgebra lineal tanto con matrices densas como sparse. Por último se desarrolla con un cierto detalle el tema de la paralelización en los actuales procesadores de varios núcleos, describiendo para ello el patrón empleado y las características más importantes de las dos herramientas propuestas, OpenMP y las TBB de Intel. Hay que señalar que las características de los sistemas multicuerpo problemas de pequeño tamaño, frecuente uso de la recursividad, y repetición intensiva en el tiempo de los cálculos con fuerte dependencia de los resultados anteriores dificultan extraordinariamente el uso de técnicas de paralelización frente a otras áreas de la mecánica computacional, tales como por ejemplo el cálculo por elementos finitos. Basándose en los conceptos mencionados en el Capítulo 4, el Capítulo 5 está dividido en tres secciones, una para cada formulación propuesta en esta Tesis. En cada una de estas secciones se describen los detalles de cómo se han realizado las distintas implementaciones propuestas para cada algoritmo y qué herramientas se han utilizado para ello. En la primera sección se muestra el uso de librerías numéricas para matrices densas y sparse en la formulación topológica semirrecursiva basada en la doble transformación de velocidades. En la segunda se describe la utilización de paralelización mediante OpenMP y TBB en la formulación semirrecursiva con penalizadores y proyecciones. Por último, se describe el uso de técnicas de matrices sparse y paralelización en la formulación global con uniones flexibles y parámetros de Euler. El Capítulo 6 describe los resultados alcanzados mediante las formulaciones e implementaciones descritas previamente. Este capítulo comienza con una descripción de la modelización y topología de los dos vehículos estudiados. El primer modelo es un vehículo de dos ejes del tipo chasis-cabina o furgoneta, perteneciente a la gama de vehículos de carga medianos. El segundo es un vehículo de cinco ejes que responde al modelo de un camión o cabina con semirremolque, perteneciente a la categoría de vehículos industriales pesados. En este capítulo además se realiza un estudio comparativo entre las simulaciones de estos vehículos con cada una de las formulaciones utilizadas y se presentan de modo cuantitativo los efectos de las mejoras alcanzadas con las distintas estrategias propuestas en esta Tesis. Con objeto de extraer conclusiones más fácilmente y para evaluar de un modo más objetivo las mejoras introducidas en la Tesis, todos los resultados de este capítulo se han obtenido con el mismo computador, que era el top de la gama Intel Xeon en 2007, pero que hoy día está ya algo obsoleto. Por último los Capítulos 7 y 8 están dedicados a las conclusiones finales y las futuras líneas de investigación que pueden derivar del trabajo realizado en esta Tesis. Los objetivos de realizar simulaciones en tiempo real de vehículos industriales de gran complejidad han sido alcanzados con varias de las formulaciones e implementaciones desarrolladas. ABSTRACT The objective of this Dissertation has been the achievement of real time simulations of industrial vehicles modeled as complex multibody systems made up by rigid bodies. For the development of a simulation program, four main aspects must be considered: the modeling of the multibody system (types of coordinates, ideal joints or imposed by means of forces), the formulation to be used to set the differential equations of motion (dependent or independent coordinates, global or topological methods, ways to impose constraints equations), the method of numerical integration to solve these equations in time (explicit or implicit integrators) and the details of the implementation carried out (programming language, mathematical libraries, parallelization techniques). These four stages are interrelated and all of them are part of this work. They involve the generation of models for a van and a semitrailer truck, the use of three different dynamic formulations, the integration of differential equations of motion through explicit and implicit methods, the use of BLAS functions and sparse matrix techniques, and the introduction of parallelization to use the different processor cores. The work presented in this Dissertation has been structured in eight chapters, the first of them being the Introduction. In Chapter 2, two different semi-recursive formulations are shown, of which the first one is based on a double velocity transformation, thus getting the differential equations of motion as a function of the independent relative accelerations. The numerical integration of these equations has been made with the Runge-Kutta explicit method of fourth order. The second formulation is based on dependent relative coordinates, imposing the constraints by means of position penalty coefficients and correcting the velocities and accelerations by projection methods. In this second case, the integration of the motion equations has been carried out by means of the HHT implicit integrator (Hilber, Hughes and Taylor), which belongs to the Newmark structural integrators family. In Chapter 3, the third formulation used in this Dissertation is presented. In this case, the joints between the bodies of the system have been considered as flexible joints, with forces used to impose the joint conditions. This kind of union hinders to work with relative coordinates, so the position of the system bodies and the setting of the equations of motion have been carried out using Cartesian coordinates and Euler parameters. In this global formulation, constraints are introduced through forces (with a similar approach to the penalty coefficients) are presented. The stabilization of the numerical integration is carried out also by velocity and accelerations projections. In Chapter 4, a revision of the main computer tools and strategies used to increase the efficiency of the implementations of the algorithms is presented. First of all, some basic considerations to increase the numerical efficiency of the implementations are included. Then the main characteristics of the code’ analyzers used and also the mathematical libraries used to solve linear algebra problems (both with dense and sparse matrices) are mentioned. Finally, the topic of parallelization in current multicore processors is developed thoroughly. For that, the pattern used and the most important characteristics of the tools proposed, OpenMP and Intel TBB, are described. It needs to be highlighted that the characteristics of multibody systems small size problems, frequent recursion use and intensive repetition along the time of the calculation with high dependencies of the previous results complicate extraordinarily the use of parallelization techniques against other computational mechanics areas, as the finite elements computation. Based on the concepts mentioned in Chapter 4, Chapter 5 is divided into three sections, one for each formulation proposed in this Dissertation. In each one of these sections, the details of how these different proposed implementations have been made for each algorithm and which tools have been used are described. In the first section, it is shown the use of numerical libraries for dense and sparse matrices in the semirecursive topological formulation based in the double velocity transformation. In the second one, the use of parallelization by means OpenMP and TBB is depicted in the semi-recursive formulation with penalization and projections. Lastly, the use of sparse matrices and parallelization techniques is described in the global formulation with flexible joints and Euler parameters. Chapter 6 depicts the achieved results through the formulations and implementations previously described. This chapter starts with a description of the modeling and topology of the two vehicles studied. The first model is a two-axle chassis-cabin or van like vehicle, which belongs to the range of medium charge vehicles. The second one is a five-axle vehicle belonging to the truck or cabin semi-trailer model, belonging to the heavy industrial vehicles category. In this chapter, a comparative study is done between the simulations of these vehicles with each one of the formulations used and the improvements achieved are presented in a quantitative way with the different strategies proposed in this Dissertation. With the aim of deducing the conclusions more easily and to evaluate in a more objective way the improvements introduced in the Dissertation, all the results of this chapter have been obtained with the same computer, which was the top one among the Intel Xeon range in 2007, but which is rather obsolete today. Finally, Chapters 7 and 8 are dedicated to the final conclusions and the future research projects that can be derived from the work presented in this Dissertation. The objectives of doing real time simulations in high complex industrial vehicles have been achieved with the formulations and implementations developed.