Modelos de elementos finitos explícitos para explosiones en estructuras reticuladas de hormigón armado : aplicaciones al estudio del colapso de edificios

La frecuencia con la que se producen explosiones sobre edificios, ya sean accidentales o intencionadas, es reducida, pero sus efectos pueden ser catastróficos. Es deseable poder predecir de forma suficientemente precisa las consecuencias de estas acciones dinámicas sobre edificaciones civiles, entre las cuales las estructuras reticuladas de hormigón armado son una tipología habitual. En esta tesis doctoral se exploran distintas opciones prácticas para el modelado y cálculo numérico por ordenador de estructuras de hormigón armado sometidas a explosiones. Se emplean modelos numéricos de elementos finitos con integración explícita en el tiempo, que demuestran su capacidad efectiva para simular los fenómenos físicos y estructurales de dinámica rápida y altamente no lineales que suceden, pudiendo predecir los daños ocasionados tanto por la propia explosión como por el posible colapso progresivo de la estructura. El trabajo se ha llevado a cabo empleando el código comercial de elementos finitos LS-DYNA (Hallquist, 2006), desarrollando en el mismo distintos tipos de modelos de cálculo que se pueden clasificar en dos tipos principales: 1) modelos basados en elementos finitos de continuo, en los que se discretiza directamente el medio continuo mediante grados de libertad nodales de desplazamientos; 2) modelos basados en elementos finitos estructurales, mediante vigas y láminas, que incluyen hipótesis cinemáticas para elementos lineales o superficiales. Estos modelos se desarrollan y discuten a varios niveles distintos: 1) a nivel del comportamiento de los materiales, 2) a nivel de la respuesta de elementos estructurales tales como columnas, vigas o losas, y 3) a nivel de la respuesta de edificios completos o de partes significativas de los mismos. Se desarrollan modelos de elementos finitos de continuo 3D muy detallados que modelizan el hormigón en masa y el acero de armado de forma segregada. El hormigón se representa con un modelo constitutivo del hormigón CSCM (Murray et al., 2007), que tiene un comportamiento inelástico, con diferente respuesta a tracción y compresión, endurecimiento, daño por fisuración y compresión, y rotura. El acero se representa con un modelo constitutivo elastoplástico bilineal con rotura. Se modeliza la geometría precisa del hormigón mediante elementos finitos de continuo 3D y cada una de las barras de armado mediante elementos finitos tipo viga, con su posición exacta dentro de la masa de hormigón. La malla del modelo se construye mediante la superposición de los elementos de continuo de hormigón y los elementos tipo viga de las armaduras segregadas, que son obligadas a seguir la deformación del sólido en cada punto mediante un algoritmo de penalización, simulando así el comportamiento del hormigón armado. En este trabajo se denominarán a estos modelos simplificadamente como modelos de EF de continuo. Con estos modelos de EF de continuo se analiza la respuesta estructural de elementos constructivos (columnas, losas y pórticos) frente a acciones explosivas. Asimismo se han comparado con resultados experimentales, de ensayos sobre vigas y losas con distintas cargas de explosivo, verificándose una coincidencia aceptable y permitiendo una calibración de los parámetros de cálculo. Sin embargo estos modelos tan detallados no son recomendables para analizar edificios completos, ya que el elevado número de elementos finitos que serían necesarios eleva su coste computacional hasta hacerlos inviables para los recursos de cálculo actuales. Adicionalmente, se desarrollan modelos de elementos finitos estructurales (vigas y láminas) que, con un coste computacional reducido, son capaces de reproducir el comportamiento global de la estructura con una precisión similar. Se modelizan igualmente el hormigón en masa y el acero de armado de forma segregada. El hormigón se representa con el modelo constitutivo del hormigón EC2 (Hallquist et al., 2013), que también presenta un comportamiento inelástico, con diferente respuesta a tracción y compresión, endurecimiento, daño por fisuración y compresión, y rotura, y se usa en elementos finitos tipo lámina. El acero se representa de nuevo con un modelo constitutivo elastoplástico bilineal con rotura, usando elementos finitos tipo viga. Se modeliza una geometría equivalente del hormigón y del armado, y se tiene en cuenta la posición relativa del acero dentro de la masa de hormigón. Las mallas de ambos se unen mediante nodos comunes, produciendo una respuesta conjunta. En este trabajo se denominarán a estos modelos simplificadamente como modelos de EF estructurales. Con estos modelos de EF estructurales se simulan los mismos elementos constructivos que con los modelos de EF de continuo, y comparando sus respuestas estructurales frente a explosión se realiza la calibración de los primeros, de forma que se obtiene un comportamiento estructural similar con un coste computacional reducido. Se comprueba que estos mismos modelos, tanto los modelos de EF de continuo como los modelos de EF estructurales, son precisos también para el análisis del fenómeno de colapso progresivo en una estructura, y que se pueden utilizar para el estudio simultáneo de los daños de una explosión y el posterior colapso. Para ello se incluyen formulaciones que permiten considerar las fuerzas debidas al peso propio, sobrecargas y los contactos de unas partes de la estructura sobre otras. Se validan ambos modelos con un ensayo a escala real en el que un módulo con seis columnas y dos plantas colapsa al eliminar una de sus columnas. El coste computacional del modelo de EF de continuo para la simulación de este ensayo es mucho mayor que el del modelo de EF estructurales, lo cual hace inviable su aplicación en edificios completos, mientras que el modelo de EF estructurales presenta una respuesta global suficientemente precisa con un coste asumible. Por último se utilizan los modelos de EF estructurales para analizar explosiones sobre edificios de varias plantas, y se simulan dos escenarios con cargas explosivas para un edificio completo, con un coste computacional moderado. The frequency of explosions on buildings whether they are intended or accidental is small, but they can have catastrophic effects. Being able to predict in a accurate enough manner the consequences of these dynamic actions on civil buildings, among which frame-type reinforced concrete buildings are a frequent typology is desirable. In this doctoral thesis different practical options for the modeling and computer assisted numerical calculation of reinforced concrete structures submitted to explosions are explored. Numerical finite elements models with explicit time-based integration are employed, demonstrating their effective capacity in the simulation of the occurring fast dynamic and highly nonlinear physical and structural phenomena, allowing to predict the damage caused by the explosion itself as well as by the possible progressive collapse of the structure. The work has been carried out with the commercial finite elements code LS-DYNA (Hallquist, 2006), developing several types of calculation model classified in two main types: 1) Models based in continuum finite elements in which the continuous medium is discretized directly by means of nodal displacement degrees of freedom; 2) Models based on structural finite elements, with beams and shells, including kinematic hypothesis for linear and superficial elements. These models are developed and discussed at different levels: 1) material behaviour, 2) response of structural elements such as columns, beams and slabs, and 3) response of complete buildings or significative parts of them. Very detailed 3D continuum finite element models are developed, modeling mass concrete and reinforcement steel in a segregated manner. Concrete is represented with a constitutive concrete model CSCM (Murray et al., 2007), that has an inelastic behaviour, with different tension and compression response, hardening, cracking and compression damage and failure. The steel is represented with an elastic-plastic bilinear model with failure. The actual geometry of the concrete is modeled with 3D continuum finite elements and every and each of the reinforcing bars with beam-type finite elements, with their exact position in the concrete mass. The mesh of the model is generated by the superposition of the concrete continuum elements and the beam-type elements of the segregated reinforcement, which are made to follow the deformation of the solid in each point by means of a penalty algorithm, reproducing the behaviour of reinforced concrete. In this work these models will be called continuum FE models as a simplification. With these continuum FE models the response of construction elements (columns, slabs and frames) under explosive actions are analysed. They have also been compared with experimental results of tests on beams and slabs with various explosive charges, verifying an acceptable coincidence and allowing a calibration of the calculation parameters. These detailed models are however not advised for the analysis of complete buildings, as the high number of finite elements necessary raises its computational cost, making them unreliable for the current calculation resources. In addition to that, structural finite elements (beams and shells) models are developed, which, while having a reduced computational cost, are able to reproduce the global behaviour of the structure with a similar accuracy. Mass concrete and reinforcing steel are also modeled segregated. Concrete is represented with the concrete constitutive model EC2 (Hallquist et al., 2013), which also presents an inelastic behaviour, with a different tension and compression response, hardening, compression and cracking damage and failure, and is used in shell-type finite elements. Steel is represented once again with an elastic-plastic bilineal with failure constitutive model, using beam-type finite elements. An equivalent geometry of the concrete and the steel is modeled, considering the relative position of the steel inside the concrete mass. The meshes of both sets of elements are bound with common nodes, therefore producing a joint response. These models will be called structural FE models as a simplification. With these structural FE models the same construction elements as with the continuum FE models are simulated, and by comparing their response under explosive actions a calibration of the former is carried out, resulting in a similar response with a reduced computational cost. It is verified that both the continuum FE models and the structural FE models are also accurate for the analysis of the phenomenon of progressive collapse of a structure, and that they can be employed for the simultaneous study of an explosion damage and the resulting collapse. Both models are validated with an experimental full-scale test in which a six column, two floors module collapses after the removal of one of its columns. The computational cost of the continuum FE model for the simulation of this test is a lot higher than that of the structural FE model, making it non-viable for its application to full buildings, while the structural FE model presents a global response accurate enough with an admissible cost. Finally, structural FE models are used to analyze explosions on several story buildings, and two scenarios are simulated with explosive charges for a full building, with a moderate computational cost.

Estudio teórico-experimental para la mejora de la resistencia estructural de autocares basado en el diseño de uniones y asientos

La presente tesis analiza la mejora de la resistencia estructural ante vuelco de autocares enfocando dos vías de actuación: análisis y propuestas de requisitos reglamentarios a nivel europeo y la generación de herramientas que ayuden al diseño y a la verificación de estos requisitos. Los requisitos reglamentarios de resistencia estructural a vuelco contemplan la superestructura de los vehículos pero no para los asientos y sistemas de retención. La influencia de los pasajeros retenidos es superior a la incluida en reglamentación (Reg. 66.01) debiendo considerarse unida al vehículo un porcentaje de la masa de los pasajeros del 91% para cinturón de tres puntos y del 52% para cinturón subabdominal frente al 50% reglamentario para todos los casos. Se ha determinado la cinemática y dinámica del vuelco normativo en sus diferentes fases, formulando las energías en las fases iniciales (hasta el impacto contra el suelo) y determinando la fase final de deformación a través del análisis secuencial de ensayos de módulos reales. Se han determinado los esfuerzos para los asientos que se dividen en dos fases diferenciadas temporalmente: una primera debida a la deformación estructural y una segunda debida al esfuerzo del pasajero retenido que se produce en sentido opuesto (con una deceleración del pasajero en torno a 3.3 g). Se ha caracterizado a través de ensayos cuasi.estáticos el comportamiento de perfiles a flexión y de las uniones estructurales de las principales zonas del vehículo (piso, ventana y techo) verificándose la validez del comportamiento plástico teórico Kecman.García para perfiles de hasta 4 mm de espesor y caracterizando la resistencia y rigidez en la zona elástica de las uniones en función del tipo de refuerzo, materiales y perfiles (análisis de más de 180 probetas). Se ha definido un método de ensayo cuasi.estático para asientos ante esfuerzos de vuelco, ensayándose 19 butacas y determinándose que son resistentes (salvo las uniones a vehículo con pinzas), que son capaces de absorber hasta más de un 17% de la energía absorbida, aunque algunos necesitan optimización para llegar a contribuir en el mecanismo de deformación estructural. Se han generado modelos simplificados para introducir en los modelos barra.rótula plástica: un modelo combinado unión+rótula plástica (que incluye la zona de rigidez determinada en función del tipo de unión) para la superestructura y un modelo simplificado de muelles no.lineales para los asientos. Igualmente se ha generado la metodología de diseño a través de ensayos virtuales con modelos de detalle de elementos finitos tanto de las uniones como de los asientos. Se ha propuesto una metodología de diseño basada en obtener el “mecanismo óptimo de deformación estructural” (elevando la zona de deformación lateral a nivel de ventana y en pilar o en costilla en techo). Para ello se abren dos vías: diseño de la superestructura (selección de perfiles y generación de uniones resistentes) o combinación con asientos (que en lugar de solo resistir las cargas pueden llegar a modificar el mecanismo de deformación). Se ha propuesto una metodología de verificación alternativa al vuelco de vehículo completo que contempla el cálculo cuasi.estático con modelos simplificados barra.rótula plástica más el ensayo de una sección representativa con asientos y utillajes antropomórficos retenidos que permite validar el diseño de las uniones, determinar el porcentaje de energía que debe absorberse por deformación estructural (factor C) y verificar el propio asiento como sistema de retención. ABSTRACT This research analyzes the improvement of the structural strength of buses and coaches under rollover from two perspectives: regulatory requirements at European level and generation of tools that will help to the design and to the verification of requirements. European Regulations about rollover structural strength includes requirements for the superstructure of the vehicles but not about seats, anchorages and restraint systems. The influence of the retained passengers is higher than the one included currently in the Regulations (Reg. 66.01), being needed to consider a 91% of the passenger mass as rigidly joint to the vehicle (for 3 points’ belt, a 52% for 2 points’ belt) instead of the 50% included in the Regulation. Kinematic and dynamic of the normative rollover has been determined from testing of different sections, formulating the energies of the first phases (up to the first impact with the ground) and determining the last deformation phase through sequential analysis of movements and deformations. The efforts due to rollover over the seats have been established, being divided in two different temporal phases: a first one due to the structural deformation of the vehicle and a second one due to the effort of the restrained passenger being this second one in opposite sense (with a passenger deceleration around 3.3 g). From quasi.static testing, the behavior of the structural tubes under flexural loads, including the principal joints in the vehicle (floor, window and roof), the validity of the theoretical plastic behavior according Kecman.García theories have been verified up to 4 mm of thickness. Strength of the joints as well as the stiffness of the elastic zone has been determined in function of main parameters: type of reinforcement, materials and section of the tubes (more than 180 test specimens). It has been defined a quasi.static testing methodology to characterize the seats and restrain system behavior under rollover, testing 19 double seats and concluding that they are resistant (excepting clamping joints), that they can absorb more than a 17 of the absorbed energy, and that some of them need optimization to contribute in the structural deformation mechanism. It has been generated simplified MEF models, to analyze in a beam.plastic hinge model: a combined model joint+plastic hinge (including the stiffness depending on the type of joint) for the superstructure and a simplified model with non.lineal springs to represent the seats. It has been detailed methodologies for detailed design of joints and seats from virtual testing (MEF models). A design methodology based in the “optimized structural deformation mechanism” (increasing the height of deformation of the lateral up to window level) is proposed. Two possibilities are analyzed: design of the superstructure based on the selection of profiles and design of strength joints (were seats only resist the efforts and contribute in the energy absorption) or combination structure.seats, were seats contributes in the deformation mechanism. An alternative methodology to the rollover of a vehicle that includes the quasi.static calculation with simplified models “beam.joint+plastic hinge” plus the testing of a representative section of the vehicle including seats and anthropomorphic ballast restrained by the safety belts is presented. The test of the section allows validate the design of the joints, determine the percentage of energy to be absorbed by structural deformation (factor C) and verify the seat as a retention system.

Métodos diferenciales e interferométricos para la evaluación de deformaciones de la corteza terrestre mediante técnicas GNSS e InSAR = Differential and interferometric methods for evaluation of crustal deformation using GNSS and InSAR techniques

La actividad volcánica interviene en multitud de facetas de la propia actividad humana, no siempre negativas. Sin embargo, son más los motivos de peligrosidad y riesgo que incitan al estudio de la actividad volcánica. Existen razones de seguridad que inciden en el mantenimiento del seguimiento y monitorización de la actividad volcánica para garantizar la vida y la seguridad de los asentamientos antrópicos en las proximidades de los edificios volcánicos. En esta tesis se define e implementa un sistema de monitorización de movimientos de la corteza en las islas de Tenerife y La Palma, donde el impacto social que representa un aumento o variación de la actividad volcánica en las islas es muy severo. Aparte de la alta densidad demográfica del Archipiélago, esta población aumenta significativamente, en diferentes periodos a lo largo del año, debido a la actividad turística que representa la mayor fuente de ingresos de las islas. La población y los centros turísticos se diseminan predominantemente a lo largo de las costas y también a lo largo de los flancos de los edificios volcánicos. Quizá el mantenimiento de estas estructuras sociales y socio-económicas son los motivos más importantes que justifican una monitorización de la actividad volcánica en las Islas Canarias. Recientemente se ha venido trabajando cada vez más en el intento de predecir la actividad volcánica utilizando los nuevos sistemas de monitorización geodésica, puesto que la actividad volcánica se manifiesta anteriormente por deformación de la corteza terrestre y cambios en la fuerza de la gravedad en la zona donde más tarde se registran eventos volcánicos. Los nuevos dispositivos y sensores que se han desarrollado en los últimos años en materias como la geodesia, la observación de la Tierra desde el espacio y el posicionamiento por satélite, han permitido observar y medir tanto la deformación producida en el terreno como los cambios de la fuerza de la gravedad antes, durante y posteriormente a los eventos volcánicos que se producen. Estos nuevos dispositivos y sensores han cambiado las técnicas o metodologías geodésicas que se venían utilizando hasta la aparición de los mismos, renovando métodos clásicos y desarrollando otros nuevos que ya se están afianzando como metodologías probadas y reconocidas para ser usadas en la monitorización volcánica. Desde finales de la década de los noventa del siglo pasado se han venido desarrollando en las Islas Canarias varios proyectos que han tenido como objetivos principales el desarrollo de nuevas técnicas de observación y monitorización por un lado y el diseño de una metodología de monitorización volcánica adecuada, por otro. Se presenta aquí el estudio y desarrollo de técnicas GNSS para la monitorización de deformaciones corticales y su campo de velocidades para las islas de Tenerife y La Palma. En su implementación, se ha tenido en cuenta el uso de la infraestructura geodésica y de monitorización existente en el archipiélago a fin de optimizar costes, además de complementarla con nuevas estaciones para dar una cobertura total a las dos islas. Los resultados obtenidos en los proyectos, que se describen en esta memoria, han dado nuevas perspectivas en la monitorización geodésica de la actividad volcánica y nuevas zonas de interés que anteriormente no se conocían en el entorno de las Islas Canarias. Se ha tenido especial cuidado en el tratamiento y propagación de los errores durante todo el proceso de observación, medida y proceso de los datos registrados, todo ello en aras de cuantificar el grado de fiabilidad de los resultados obtenidos. También en este sentido, los resultados obtenidos han sido verificados con otros procedentes de sistemas de observación radar de satélite, incorporando además a este estudio las implicaciones que el uso conjunto de tecnologías radar y GNSS tendrán en un futuro en la monitorización de deformaciones de la corteza terrestre. ABSTRACT Volcanic activity occurs in many aspects of human activity, and not always in a negative manner. Nonetheless, research into volcanic activity is more likely to be motivated by its danger and risk. There are security reasons that influence the monitoring of volcanic activity in order to guarantee the life and safety of human settlements near volcanic edifices. This thesis defines and implements a monitoring system of movements in the Earth’s crust in the islands of Tenerife and La Palma, where the social impact of an increase (or variation) of volcanic activity is very severe. Aside from the high demographic density of the archipelago, the population increases significantly in different periods throughout the year due to tourism, which represents a major source of revenue for the islands. The population and the tourist centres are mainly spread along the coasts and also along the flanks of the volcanic edifices. Perhaps the preservation of these social and socio-economic structures is the most important reason that justifies monitoring volcanic activity in the Canary Islands. Recently more and more work has been done with the intention of predicting volcanic activity, using new geodesic monitoring systems, since volcanic activity is evident prior to eruption because of a deformation of the Earth’s crust and changes in the force of gravity in the zone where volcanic events will later be recorded. The new devices and sensors that have been developed in recent years in areas such as geodesy, the observation of the Earth from space, and satellite positioning have allowed us to observe and measure the deformation produced in the Earth as well as the changes in the force of gravity before, during, and after the volcanic events occur. The new devices and sensors have changed the geodetic techniques and methodologies that were used previously. The classic methods have been renovated and other newer ones developed that are now vouched for as proven recognised methodologies to be used for volcanic monitoring. Since the end of the 1990s, in the Canary Islands various projects have been developed whose principal aim has been the development of new observation and monitoring techniques on the one hand, and the design of an appropriate volcanic monitoring methodology on the other. The study and development of GNSS techniques for the monitoring of crustal deformations and their velocity field is presented here. To carry out the study, the use of geodetic infrastructure and existing monitoring in the archipelago have been taken into account in order to optimise costs, besides complementing it with new stations for total coverage on both islands. The results obtained in the projects, which are described below, have produced new perspectives in the geodetic monitoring of volcanic activity and new zones of interest which previously were unknown in the environment of the Canary Islands. Special care has been taken with the treatment and propagation of errors during the entire process of observing, measuring, and processing the recorded data. All of this was done in order to quantify the degree of trustworthiness of the results obtained. Also in this sense, the results obtained have been verified with others from satellite radar observation systems, incorporating as well in this study the implications that the joint use of radar technologies and GNSS will have for the future of monitoring deformations in the Earth’s crust.

Nepfix: on-demand cloud elastic simulation of networks of evolutionary processors

Una Red de Procesadores Evolutivos o NEP (por sus siglas en ingles), es un modelo computacional inspirado por el modelo evolutivo de las celulas, específicamente por las reglas de multiplicación de las mismas. Esta inspiración hace que el modelo sea una abstracción sintactica de la manipulation de information de las celulas. En particu¬lar, una NEP define una maquina de cómputo teorica capaz de resolver problemas NP completos de manera eficiente en tóerminos de tiempo. En la praóctica, se espera que las NEP simuladas en móaquinas computacionales convencionales puedan resolver prob¬lemas reales complejos (que requieran ser altamente escalables) a cambio de una alta complejidad espacial. En el modelo NEP, las cóelulas estóan representadas por palabras que codifican sus secuencias de ADN. Informalmente, en cualquier momento de cómputo del sistema, su estado evolutivo se describe como un coleccion de palabras, donde cada una de ellas representa una celula. Estos momentos fijos de evolucion se denominan configuraciones. De manera similar al modelo biologico, las palabras (celulas) mutan y se dividen en base a bio-operaciones sencillas, pero solo aquellas palabras aptas (como ocurre de forma parecida en proceso de selection natural) seran conservadas para la siguiente configuracióon. Una NEP como herramienta de computation, define una arquitectura paralela y distribuida de procesamiento simbolico, en otras palabras, una red de procesadores de lenguajes. Desde el momento en que el modelo fue propuesto a la comunidad científica en el año 2001, múltiples variantes se han desarrollado y sus propiedades respecto a la completitud computacional, eficiencia y universalidad han sido ampliamente estudiadas y demostradas. En la actualidad, por tanto, podemos considerar que el modelo teórico NEP se encuentra en el estadio de la madurez. La motivación principal de este Proyecto de Fin de Grado, es proponer una aproxi-mación práctica que permita dar un salto del modelo teórico NEP a una implantación real que permita su ejecucion en plataformas computacionales de alto rendimiento, con el fin de solucionar problemas complejos que demanda la sociedad actual. Hasta el momento, las herramientas desarrolladas para la simulation del modelo NEP, si bien correctas y con resultados satisfactorios, normalmente estón atadas a su entorno de ejecucion, ya sea el uso de hardware específico o implementaciones particulares de un problema. En este contexto, el propósito fundamental de este trabajo es el desarrollo de Nepfix, una herramienta generica y extensible para la ejecucion de cualquier algo¬ritmo de un modelo NEP (o alguna de sus variantes), ya sea de forma local, como una aplicación tradicional, o distribuida utilizando los servicios de la nube. Nepfix es una aplicacion software desarrollada durante 7 meses y que actualmente se encuentra en su segunda iteration, una vez abandonada la fase de prototipo. Nepfix ha sido disenada como una aplicacion modular escrita en Java 8 y autocontenida, es decir, no requiere de un entorno de ejecucion específico (cualquier maquina virtual de Java es un contenedor vólido). Nepfix contiene dos componentes o móodulos. El primer móodulo corresponde a la ejecución de una NEP y es por lo tanto, el simulador. Para su desarrollo, se ha tenido en cuenta el estado actual del modelo, es decir, las definiciones de los procesadores y filtros mas comunes que conforman la familia del modelo NEP. Adicionalmente, este componente ofrece flexibilidad en la ejecucion, pudiendo ampliar las capacidades del simulador sin modificar Nepfix, usando para ello un lenguaje de scripting. Dentro del desarrollo de este componente, tambióen se ha definido un estóandar de representacióon del modelo NEP basado en el formato JSON y se propone una forma de representation y codificación de las palabras, necesaria para la comunicación entre servidores. Adicional-mente, una característica importante de este componente, es que se puede considerar una aplicacion aislada y por tanto, la estrategia de distribution y ejecución son total-mente independientes. El segundo moódulo, corresponde a la distribucióon de Nepfix en la nube. Este de-sarrollo es el resultado de un proceso de i+D, que tiene una componente científica considerable. Vale la pena resaltar el desarrollo de este modulo no solo por los resul-tados prócticos esperados, sino por el proceso de investigation que se se debe abordar con esta nueva perspectiva para la ejecución de sistemas de computación natural. La principal característica de las aplicaciones que se ejecutan en la nube es que son gestionadas por la plataforma y normalmente se encapsulan en un contenedor. En el caso de Nepfix, este contenedor es una aplicacion Spring que utiliza el protocolo HTTP o AMQP para comunicarse con el resto de instancias. Como valor añadido, Nepfix aborda dos perspectivas de implementation distintas (que han sido desarrolladas en dos iteraciones diferentes) del modelo de distribution y ejecucion, que tienen un impacto muy significativo en las capacidades y restricciones del simulador. En concreto, la primera iteration utiliza un modelo de ejecucion asincrono. En esta perspectiva asincrona, los componentes de la red NEP (procesadores y filtros) son considerados como elementos reactivos a la necesidad de procesar una palabra. Esta implementation es una optimization de una topologia comun en el modelo NEP que permite utilizar herramientas de la nube para lograr un escalado transparente (en lo ref¬erente al balance de carga entre procesadores) pero produce efectos no deseados como indeterminacion en el orden de los resultados o imposibilidad de distribuir eficiente-mente redes fuertemente interconectadas. Por otro lado, la segunda iteration corresponde al modelo de ejecucion sincrono. Los elementos de una red NEP siguen un ciclo inicio-computo-sincronizacion hasta que el problema se ha resuelto. Esta perspectiva sincrona representa fielmente al modelo teórico NEP pero el proceso de sincronizacion es costoso y requiere de infraestructura adicional. En concreto, se requiere un servidor de colas de mensajes RabbitMQ. Sin embargo, en esta perspectiva los beneficios para problemas suficientemente grandes superan a los inconvenientes, ya que la distribuciín es inmediata (no hay restricciones), aunque el proceso de escalado no es trivial. En definitiva, el concepto de Nepfix como marco computacional se puede considerar satisfactorio: la tecnología es viable y los primeros resultados confirman que las carac-terísticas que se buscaban originalmente se han conseguido. Muchos frentes quedan abiertos para futuras investigaciones. En este documento se proponen algunas aproxi-maciones a la solucion de los problemas identificados como la recuperacion de errores y la division dinamica de una NEP en diferentes subdominios. Por otra parte, otros prob-lemas, lejos del alcance de este proyecto, quedan abiertos a un futuro desarrollo como por ejemplo, la estandarización de la representación de las palabras y optimizaciones en la ejecucion del modelo síncrono. Finalmente, algunos resultados preliminares de este Proyecto de Fin de Grado han sido presentados recientemente en formato de artículo científico en la "International Work-Conference on Artificial Neural Networks (IWANN)-2015" y publicados en "Ad-vances in Computational Intelligence" volumen 9094 de "Lecture Notes in Computer Science" de Springer International Publishing. Lo anterior, es una confirmation de que este trabajo mas que un Proyecto de Fin de Grado, es solo el inicio de un trabajo que puede tener mayor repercusion en la comunidad científica. Abstract Network of Evolutionary Processors -NEP is a computational model inspired by the evolution of cell populations, which might model some properties of evolving cell communities at the syntactical level. NEP defines theoretical computing devices able to solve NP complete problems in an efficient manner. In this model, cells are represented by words which encode their DNA sequences. Informally, at any moment of time, the evolutionary system is described by a collection of words, where each word represents one cell. Cells belong to species and their community evolves according to mutations and division which are defined by operations on words. Only those cells are accepted as surviving (correct) ones which are represented by a word in a given set of words, called the genotype space of the species. This feature is analogous with the natural process of evolution. Formally, NEP is based on an architecture for parallel and distributed processing, in other words, a network of language processors. Since the date when NEP was pro¬posed, several extensions and variants have appeared engendering a new set of models named Networks of Bio-inspired Processors (NBP). During this time, several works have proved the computational power of NBP. Specifically, their efficiency, universality, and computational completeness have been thoroughly investigated. Therefore, we can say that the NEP model has reached its maturity. The main motivation for this End of Grade project (EOG project in short) is to propose a practical approximation that allows to close the gap between theoretical NEP model and a practical implementation in high performing computational platforms in order to solve some of high the high complexity problems society requires today. Up until now tools developed to simulate NEPs, while correct and successful, are usu¬ally tightly coupled to the execution environment, using specific software frameworks (Hadoop) or direct hardware usage (GPUs). Within this context the main purpose of this work is the development of Nepfix, a generic and extensible tool that aims to execute algorithms based on NEP model and compatible variants in a local way, similar to a traditional application or in a distributed cloud environment. Nepfix as an application was developed during a 7 month cycle and is undergoing its second iteration once the prototype period was abandoned. Nepfix is designed as a modular self-contained application written in Java 8, that is, no additional external dependencies are required and it does not rely on an specific execution environment, any JVM is a valid container. Nepfix is made of two components or modules. The first module corresponds to the NEP execution and therefore simulation. During the development the current state of the theoretical model was used as a reference including most common filters and processors. Additionally extensibility is provided by the use of Python as a scripting language to run custom logic. Along with the simulation a definition language for NEP has been defined based on JSON as well as a mechanisms to represent words and their possible manipulations. NEP simulator is isolated from distribution and as mentioned before different applications that include it as a dependency are possible, the distribution of NEPs is an example of this. The second module corresponds to executing Nepfix in the cloud. The development carried a heavy R&D process since this front was not explored by other research groups until now. It's important to point out that the development of this module is not focused on results at this point in time, instead we focus on feasibility and discovery of this new perspective to execute natural computing systems and NEPs specifically. The main properties of cloud applications is that they are managed by the platform and are encapsulated in a container. For Nepfix a Spring application becomes the container and the HTTP or AMQP protocols are used for communication with the rest of the instances. Different execution perspectives were studied, namely asynchronous and synchronous models were developed for solving different kind of problems using NEPs. Different limitations and restrictions manifest in both models and are explored in detail in the respective chapters. In conclusion we can consider that Nepfix as a computational framework is suc-cessful: Cloud technology is ready for the challenge and the first results reassure that the properties Nepfix project pursued were met. Many investigation branches are left open for future investigations. In this EOG implementation guidelines are proposed for some of them like error recovery or dynamic NEP splitting. On the other hand other interesting problems that were not in the scope of this project were identified during development like word representation standardization or NEP model optimizations. As a confirmation that the results of this work can be useful to the scientific com-munity a preliminary version of this project was published in The International Work- Conference on Artificial Neural Networks (IWANN) in May 2015. Development has not stopped since that point and while Nepfix in it's current state can not be consid¬ered a final product the most relevant ideas, possible problems and solutions that were produced during the seven months development cycle are worthy to be gathered and presented giving a meaning to this EOG work.

The folded plate model in structural analysis. Elastic, plastic and dynamic formulations for nonprismatic structures

A specific numerical procedure for the analysis of arbitrary nonprismatic folded plate structures is presented. An elastic model is studied and compared with a harmonic solution for a prismatic structure. An extension to the plastic analysis is developed, and the influence of the structural geometry and loading pattern is analyzed. Nonprismatic practical cases, with arbitrary geometry and loading are shown, as well in the elastic range as in the plastic one. Finally, a dynamic formulation is outlined

High Resolution Detection of Mechanical Forces Exerted by Locomoting Fibroblasts on the Substrate

We have developed a new approach to detect mechanical forces exerted by locomoting fibroblasts on the substrate. Cells were cultured on elastic, collagen-coated polyacrylamide sheets embedded with 0.2-μm fluorescent beads. Forces exerted by the cell cause deformation of the substrate and displacement of the beads. By recording the position of beads during cell locomotion and after cell removal, we discovered that most forces were radially distributed, switching direction in the anterior region. Deformations near the leading edge were strong, transient, and variable in magnitude, consistent with active local contractions, whereas those in the posterior region were weaker, more stable, and more uniform, consistent with passive resistance. Treatment of cells with cytochalasin D or myosin II inhibitors caused relaxation of the forces, suggesting that they are generated primarily via actin–myosin II interactions; treatment with nocodazole caused no immediate effect on forces. Immunofluorescence indicated that the frontal region of strong deformation contained many vinculin plaques but no apparent concentration of actin or myosin II filaments. Strong mechanical forces in the anterior region, generated by locally activated myosin II and transmitted through vinculin-rich structures, likely play a major role in cell locomotion and in mechanical signaling with the surrounding environment.

Characterization of an In Vitro Model of Elastic Fiber Assembly

Elastic fibers consist of two morphologically distinct components: elastin and 10-nm fibrillin-containing microfibrils. During development, the microfibrils form bundles that appear to act as a scaffold for the deposition, orientation, and assembly of tropoelastin monomers into an insoluble elastic fiber. Although microfibrils can assemble independent of elastin, tropoelastin monomers do not assemble without the presence of microfibrils. In the present study, immortalized ciliary body pigmented epithelial (PE) cells were investigated for their potential to serve as a cell culture model for elastic fiber assembly. Northern analysis showed that the PE cells express microfibril proteins but do not express tropoelastin. Immunofluorescence staining and electron microscopy confirmed that the microfibril proteins produced by the PE cells assemble into intact microfibrils. When the PE cells were transfected with a mammalian expression vector containing a bovine tropoelastin cDNA, the cells were found to express and secrete tropoelastin. Immunofluorescence and electron microscopic examination of the transfected PE cells showed the presence of elastic fibers in the matrix. Biochemical analysis of this matrix showed the presence of cross-links that are unique to mature insoluble elastin. Together, these results indicate that the PE cells provide a unique, stable in vitro system in which to study elastic fiber assembly.

Climatic and halokinetic controls on alluvial–lacustrine sedimentation during compressional deformation, Andean forearc, northern Chile

Peer reviewed

Integral elastic, electronic-state, ionization, and total cross sections for electron scattering with furfural

We report absolute experimental integral cross sections (ICSs) for electron impact excitation of bands of electronic-states in furfural, for incident electron energies in the range 20-250 eV. Wherever possible, those results are compared to corresponding excitation cross sections in the structurally similar species furan, as previously reported by da Costa et al. [Phys. Rev. A 85, 062706 (2012)] and Regeta and Allan [Phys. Rev. A 91, 012707 (2015)]. Generally, very good agreement is found. In addition, ICSs calculated with our independent atom model (IAM) with screening corrected additivity rule (SCAR) formalism, extended to account for interference (I) terms that arise due to the multi-centre nature of the scattering problem, are also reported. The sum of those ICSs gives the IAM-SCAR+I total cross section for electron-furfural scattering. Where possible, those calculated IAM-SCAR+I ICS results are compared against corresponding results from the present measurements with an acceptable level of accord being obtained. Similarly, but only for the band I and band II excited electronic states, we also present results from our Schwinger multichannel method with pseudopotentials calculations. Those results are found to be in good qualitative accord with the present experimental ICSs. Finally, with a view to assembling a complete cross section data base for furfural, some binary-encounter-Bethe-level total ionization cross sections for this collision system are presented. (C) 2016 AIP Publishing LLC.

Comportamento em desgaste por erosão cavitação, erosão - corrosão e em ensaios de microesclerometria linear instrumentada de um aço inoxidável martensítico AISI 410 nitretado a plasma em baixa temperatura, utilizando a tecnologia de tela ativa.

Amostras de um aço inoxidável martensítico AISI 410 temperado e revenido foram nitretadas a plasma em baixa temperatura usando o tratamento de nitretação plasma DC e a nitretação a plasma com tela ativa. Ambos os tratamentos foram realizados a 400 °C, utilizando mistura gasosa de 75 % de nitrogênio e 25 % de hidrogênio durante 20 horas e 400 Pa de pressão. As amostras de aço AISI 410 temperado e revenido foram caracterizadas antes e depois dos tratamentos termoquímicos, usando as técnicas de microscopia óptica, microscopia eletrônica de varredura, medidas de microdureza, difração de raios X e medidas de teor de nitrogênio em função da distância à superfície por espectrometria WDSX de raios X. A resistência à erosão por cavitação do aço AISI 410 nitretado DC e com tela ativa foi avaliada segundo a norma ASTM G32 (1998). Os ensaios de erosão, de erosão - corrosão e de esclerometria linear instrumentada segundo norma ASTM C1624 (2005) somente foram realizados no aço AISI 410 nitretado com tela ativa. Ensaios de nanoindentação instrumentada forma utilizados para medir a dureza (H) e o módulo de elasticidade reduzido (E*) e calcular as relações H/E* e H3/E*2 e a recuperação elástica (We), utilizando o método proposto por Oliver e Pharr. Ambos os tratamentos produziram camadas nitretadas de espessura homogênea constituídas por martensita expandida supersaturada em nitrogênio e nitretos de ferro com durezas superiores a 1200 HV, porém, a nitretação DC produziu maior quantidade de nitretos de ferro do que o tratamento de tela ativa. Os resultados de erosão por cavitação do aço nitretado DC mostraram que a precipitação de nitretos de ferro é prejudicial para a resistência à cavitação já que reduziu drasticamente o período de incubação e aumentou a taxa de perda de massa nos estágios iniciais do ensaio; entretanto, depois da remoção desses nitretos de ferro, a camada nitretada formada somente por martensita expandida resistiu bem ao dano por cavitação. Já no caso do aço nitretado com tela ativa, a resistência à erosão por cavitação aumentou 27 vezes quando comparada com o aço AISI 410 sem nitretar, fato atribuído à pequena fração volumétrica e ao menor tamanho dos nitretos de ferro presente na camada nitretada, às maiores relações H/E* e H3/E*2 e à alta recuperação elástica da martensita expandida. A remoção de massa ocorreu, principalmente, pela formação de crateras e de destacamento de material da superfície dos grãos por fratura frágil sem evidente deformação plástica. As perdas de massa acumulada mostradas pelo aço nitretado foram menores do que aquelas do aço AISI 410 nos ensaios de erosão e de erosão corrosão. O aço nitretado apresentou uma diminuição nas taxas de desgaste em ambos os ensaios de aproximadamente 50 % quando comparadas com o aço AISI 410. O mecanismo de remoção de material foi predominantemente dúctil, mesmo com o grande aumento na dureza. Os resultados de esclerometria linear instrumentada mostraram que a formação de martensita expandida possibilitou uma diminuição considerável do coeficiente de atrito em relação ao observado no caso do aço AISI 410 sem nitretar. O valor de carga crítica de falha foi de 14 N. O mecanismo de falha operante no aço nitretado foi trincamento por tensão.

Resistência à fadiga de tubo API 5L X65 cladeado e soldado circunferencialmente com eletrodos de Inconel® 625

As recentes descobertas de petróleo e gás na camada do Pré-sal representam um enorme potencial exploratório no Brasil, entretanto, os desafios tecnológicos para a exploração desses recursos minerais são imensos e, consequentemente, têm motivado o desenvolvimento de estudos voltados a métodos e materiais eficientes para suas produções. Os tubos condutores de petróleo e gás são denominados de elevadores catenários ou do inglês \"risers\", e são elementos que necessariamente são soldados e possuem fundamental importância nessa cadeia produtiva, pois transportam petróleo e gás natural do fundo do mar à plataforma, estando sujeitos a carregamentos dinâmicos (fadiga) durante sua operação. Adicionalmente, um dos problemas centrais à produção de óleo e gás das reservas do Pré-Sal está diretamente associado a meios altamente corrosivos, tais como H2S e CO2. Uma forma mais barata de proteção dos tubos é a aplicação de uma camada de um material metálico resistente à corrosão na parte interna desses tubos (clad). Assim, a união entre esses tubos para formação dos \"risers\" deve ser realizada pelo emprego de soldas circunferenciais de ligas igualmente resistentes à corrosão. Nesse contexto, como os elementos soldados são considerados possuir defeitos do tipo trinca, para a garantia de sua integridade estrutural quando submetidos a carregamentos cíclicos, é necessário o conhecimento das taxas de propagação de trinca por fadiga da solda circunferencial. Assim, neste trabalho, foram realizados ensaios de propagação de trinca por fadiga na região da solda circunferencial de Inconel® 625 realizada em tubo de aço API 5L X65 cladeado, utilizando corpos de prova do tipo SEN(B) (Single Edge Notch Bending) com relações entre espessura e largura (B/W) iguais a 0,5, 1 e 2. O propósito central deste trabalho foi de obter a curva da taxa de propagação de trinca por fadiga (da/dN) versus a variação do fator de intensidade de tensão (ΔK) para o metal de solda por meio de ensaios normatizados, utilizando diferentes técnicas de acompanhamento e medição da trinca. A monitoração de crescimento da trinca foi feita por três técnicas: variação da flexibilidade elástica (VFE), queda de potencial elétrico (QPE) e análise de imagem (Ai). Os resultados mostraram que as diferentes relações B/W utilizadas no estudo não alteraram significantemente as taxas de propagação de trinca por fadiga, respeitado que a propagação aconteceu em condições de escoamento em pequena escala na frente da trinca. Os resultados de propagação de trinca por fadiga permitiram a obtenção das regiões I e II da curva da/dN versus ΔK para o metal de solda. O valor de ΔKlim obtido para o mesmo foi em torno de 11,8 MPa.m1/2 e os valores encontrados das constantes experimentais C e m da equação de Paris-Erdogan foram respectivamente iguais a 1,55 x10-10 [(mm/ciclo)/(MPa.m1/2)m] e 4,15. A propagação de trinca no metal de solda deu-se por deformação plástica, com a formação de estrias de fadiga.

Análise da abrangência de modelo modificado para mancais curtos com deformação.

Este trabalho apresenta uma discussão sobre o estudo dos efeitos térmicos e elásticos decorrentes da pressão de sustentação presentes nos mancais. Para tanto, propõe-se um modelo matemático baseado nas equações para mancais curtos considerando a região de cavitação e utilizando o princípio da continuidade de massa. Com isto, deduzem-se as equações para o mancal a partir das equações de Reynolds e da energia, aplicando uma solução modificada para a solução de Ocvirk, sendo as equações resolvidas numericamente pelo Método das Diferenças Finitas. Somado o tratamento de mecânica dos fluidos, o trabalho discute dois modelos térmicos de previsão de temperatura média do fluido e sua influência no campo de pressão, apresentando gráficos representativos do campo de pressão e de temperatura, assim como as diferenças e implicações das diferenças. Para o cálculo de deformação da estrutura, utiliza-se um Modelo de Elementos Finitos para uma dada geometria, fazendo-se uma avaliação da variação do campo de pressão e o quanto essa diferença afeta as demais propriedades do fluido. Por fim, com o modelo completo, calcula-se o quanto esse modelamento para mancais curtos se aproxima de soluções para mancais finitos, com base em resultados da literatura, chegando a desvios quase oito vezes menores que os previstos pela literatura. Além disso, pode-se estabelecer a abrangência do modelo, ou seja, prever as condições em que suas propriedades são válidas e podem ser utilizadas para estudos iniciais.

Relationship between static and dynamic elastic modulus of calcarenite heated at different temperatures: the San Julián's stone

The San Julián’s stone is the main material used to build the most important historical buildings in Alicante city (Spain). This paper describes the analysis developed to obtain the relationship between the static and the dynamic modulus of this sedimentary rock heated at different temperatures. The rock specimens have been subjected to heating processes at different temperatures to produce different levels of weathering on 24 specimens. The static and dynamic modulus has been measured for every specimen by means of the ISRM standard and ultrasonic tests, respectively. Finally, two analytic formulas are proposed for the relationship between the static and the dynamic modulus for this stone. The results have been compared with some relationships proposed by different researchers for other types of rock. The expressions presented in this paper can be useful for the analysis, using non-destructive techniques, of the integrity level of historical constructions built with San Julián’s stone affected by fires.

Revisiting the effect of foundation embedment on elastic settlement: A new approach

The effect of foundation embedment on settlement calculation is a widely researched topic in which there is no scientific consensus regarding the magnitude of settlement reduction. In this paper, a non-linear three dimensional Finite Element analysis has been performed with the aim of evaluating the aforementioned effect. For this purpose, 1800 models were run considering different variables, such as the depth and dimensions of the foundation and the Young’s modulus and Poisson’s ratio of the soil. The settlements from models with foundations at surface level and at depth were then compared and the relationship between them established. The statistical analysis of this data allowed two new expressions, with a mean maximum error of 1.80%, for the embedment influence factor of a foundation to be proposed and these to be compared with commonly used corrections. The proposed equations were validated by comparing the settlements calculated with the proposed influence factors and the true settlements measured in several real foundations. From the comprehensive study of all modelled cases, an improved approach, when compared to those proposed by other authors, for the calculation of the true elastic settlements of an embedded foundation is proposed.