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.

Contribución del tablero y de la subestructura frente al pandeo global de pilas esbeltas de hormigón armado : determinación simplificada de una "coacción equivalente"

Las pilas de los puentes son elementos habitualmente verticales que, generalmente, se encuentran sometidos a un estado de flexión compuesta. Su altura significativa en muchas ocasiones y la gran resistencia de los materiales constituyentes de estos elementos – hormigón y acero – hace que se encuentren pilas de cierta esbeltez en la que los problemas de inestabilidad asociados al cálculo en segundo orden debido a la no linealidad geométrica deben ser considerados. Además, la mayoría de las pilas de nuestros puentes y viaductos están hechas de hormigón armado por lo que se debe considerar la fisuración del hormigón en las zonas en que esté traccionado. Es decir, el estudio del pandeo de pilas esbeltas de puentes requiere también la consideración de un cálculo en segundo orden mecánico, y no solo geométrico. Por otra parte, una pila de un viaducto no es un elemento que pueda considerarse como aislado; al contrario, su conexión con el tablero hace que aparezca una interacción entre la propia pila y aquél que, en cierta medida, supone una cierta coacción al movimiento de la propia cabeza de pila. Esto hace que el estudio de la inestabilidad de una pila esbelta de un puente no puede ser resuelto con la “teoría del pandeo de la pieza aislada”. Se plantea, entonces, la cuestión de intentar definir un procedimiento que permita abordar el problema complicado del pandeo de pilas esbeltas de puentes pero empleando herramientas de cálculo no tan complejas como las que resuelven “el pandeo global de una estructura multibarra, teniendo en cuenta todas las no linealidades, incluidas las de las coacciones”. Es decir, se trata de encontrar un procedimiento, que resulta ser iterativo, que resuelva el problema planteado de forma aproximada, pero suficientemente ajustada al resultado real, pero empleando programas “convencionales” de cálculo que sean capaces de : - por una parte, en la estructura completa: o calcular en régimen elástico lineal una estructura plana o espacial multibarra compleja; - por otra, en un modelo de una sola barra aislada: o considerar las no linealidades geométricas y mecánicas a nivel tensodeformacional, o considerar la no linealidad producida por la fisuración del hormigón, o considerar una coacción “elástica” en el extremo de la pieza. El objeto de este trabajo es precisamente la definición de ese procedimiento iterativo aproximado, la justificación de su validez, mediante su aplicación a diversos casos paramétricos, y la presentación de sus condicionantes y limitaciones. Además, para conseguir estos objetivos se han elaborado unos ábacos de nueva creación que permiten estimar la reducción de rigidez que supone la fisuración del hormigón en secciones huecas monocajón de hormigón armado. También se han creado unos novedosos diagramas de interacción axil-flector válidos para este tipo de secciones en flexión biaxial. Por último, hay que reseñar que otro de los objetivos de este trabajo – que, además, le da título - era cuantificar el valor de la coacción que existe en la cabeza de una pila debido a que el tablero transmite las cargas de una pila al resto de los integrantes de la subestructura y ésta, por tanto, colabora a reducir los movimientos de la cabeza de pila en cuestión. Es decir, la cabeza de una pila no está exenta lo cual mejora su comportamiento frente al pandeo. El régimen de trabajo de esta coacción es claramente no lineal, ya que la rigidez de las pilas depende de su grado de fisuración. Además, también influye cómo las afecta la no linealidad geométrica que, para la misma carga, aumenta la flexión de segundo orden de cada pila. En este documento se define cuánto vale esta coacción, cómo hay que calcularla y se comprueba su ajuste a los resultados obtenidos en el l modelo no lineal completo. The piers of the bridges are vertical elements where axial loads and bending moments are to be considered. They are often high and also the strength of the materials they are made of (concrete and steel) is also high. This means that slender piers are very common and, so, the instabilities produced by the second order effects due to the geometrical non linear effects are to be considered. In addition to this, the piers are usually made of reinforced concrete and, so, the effects of the cracking of the concrete should also be evaluated. That is, the analysis of the instabilities of te piers of a bridge should consider both the mechanical and the geometrical non linearities. Additionally, the pier of a bridge is not a single element, but just the opposite; the connection of the pier to the deck of the bridge means that the movements of the top of the pier are reduced compared to the situation of having a free end at the top of the pier. The connection between the pier and the deck is the reason why the instability of the pier cannot be analysed using “the buckling of a compressed single element method”. So, the question of defining an approximate method for analysing the buckling of the slender piers of a bridge but using a software less complex than what it is needed for analysing the “ global buckling of a multibeam structure considering all t”, is arisen. Then, the goal should be trying to find a procedure for analysing the said complex problem of the buckling of the slender piers of a bridge using a simplified method. This method could be an iterative (step by step) procedure, being accurate enough, using “normal” software having the following capabilities: - Related to the calculation of the global structure o Ability for calculating a multibesam strucutre using elastic analysis. - Related to the calculation of a single beam strcuture:: o Ability for taking into account the geometrical and mechanical () non linearities o Ability for taking into account the cracking of the concrete. o Ability for using partial stiff constraints (elastic springs) at the end of the elements One of the objectives of this document is just defining this simplified methodology, justifying the accuracy of the proposed procedure by using it on some different bridges and presenting the exclusions and limitations of the propose method. In addition to this, some new charts have been created for calculating the reduction of the stiffness of hollow cross sections made of reinforced concrete. Also, new charts for calculating the reinforcing of hollow cross sections under biaxial bending moments are also included in the document. Finally, it is to be said that another aim of the document – as it is stated on the title on the document – is defining the value of the constraint on the top of the pier because of the connection of the pier to the deck .. and to the other piers. That is, the top of the pier is not a free end of a beam and so the buckling resistance of the pier is significantly improved. This constraint is a non-elastic constraint because the stiffness of each pier depends on the level of cracking. Additionally, the geometrical non linearity is to be considered as there is an amplification of the bending moments due to the increasing of the movements of the top of the pier. This document is defining how this constraints is to be calculated; also the accuracy of the calculations is evaluated comparing the final results with the results of the complete non linear calculations

Producción de Hidrógeno mediante energía solar térmica de alta temperatura.

Este proyecto versa sobre el estudio y diseño de una central termosolar de torre central, así como de su rendimiento y producción. La característica especial de esta planta es que tiene como fin suministrar calor a un proceso químico: cracking térmico del metano para la producción de hidrógeno. Debido a que el cracking térmico tiene lugar en el interior de un tanque de metal líquido (reactor químico) y al peso del mismo, resulta conveniente dejar el reactor a nivel de suelo. Así pues, los rayos solares tienen que descender desde la parte superior de la torre hasta el reactor. Los rayos solares se reflejan en los heliostatos para volver a ser reflejados en la parte superior de la torre, donde hay otro espejo que conduce los rayos solares hasta el reactor. La localización elegida para realizar el estudio ha sido Tabernas, cerca de la Plataforma Solar de Almería (PSA). De la estación SIAR (Sistema de Información Agroclimática para el Regadío del Ministerio de Agricultura, Alimentación y Medio Ambiente del Gobierno de España). de Tabernas precisamente se han obtenido los datos de radiación global, y mediante una correlación kd-kt, se ha obtenido la radiación directa. Para el estudio de la energía aportada por el campo solar, se ha elegido un campo norte de heliostatos, debido a que un campo circular además de tener menor rendimiento en este proyecto no sería factible. Los heliostatos considerados son cuadrados de 11 x 11 m. Estos heliostatos apuntan hacia el punto (0, 0, 100) m en el sistema de coordenadas absoluto considerado para el proyecto, que coincide con una altura a 100 m en el centro interior de la torre. Este punto es uno de los focos del elipsoide virtual, del cual forma parte el reflector (espejo situado en la parte superior) y cuyo otro foco se sitúa en la parte superior del receptor (reactor), cuyo fin es dirigir los rayos hacia el reactor, como se ha indicado. Una vez definido el proyecto, se lleva a cabo un dimensionado del campo solar, con el cual puede obtenerse un campo de heliostatos. Tras realizar la simulación, se obtienen datos instantáneos y medios. A modo de ejemplo, se incluye el rendimiento por bloqueos y sombras a las 9 h de la mañana en enero, donde puede observarse que la sombra de la torre tiene una alta importancia sobre los heliostatos situados en la zona oeste. Del mismo modo, se han obtenido los datos de incidencia de los rayos solares sobre el receptor, de forma que puede caracterizarse la incidencia del flujo térmico sobre el mismo mediante un mallado. A continuación se incluye una imagen del número de rayos solares que inciden sobre el receptor a las 12 h en junio. Las conclusiones que pueden extraerse del proyecto son:  El rendimiento anual de la planta es del 26 %.  La producción anual de hidrógeno sería de 163,7 t.  Para mayores potencias de plantas el rendimiento por sombras y bloqueos sería menor, al igual que el asociado al efecto coseno, a la dispersión y absorción atmosférica y al factor de interceptación.  Para reducir el efecto del tamaño en los dos primeros rendimientos mencionados en el párrafo anterior podría elevarse la altura de la torre.  Para aumentar el factor de interceptación podría estudiarse la colocación de espejos en el interior de la torre, de modo que reflejasen los rayos solares hasta el receptor.

Límites de ductilidad y de servicio al cálculo en carga última de estructuras de hormigón armado y de fábrica

En este trabajo se aborda una cuestión central en el diseño en carga última de estructuras de hormigón armado y de fábrica: la posibilidad efectiva de que las deformaciones plásticas necesarias para verificar un estado de rotura puedan ser alcanzadas por las regiones de la estructura que deban desarrollar su capacidad última para verificar tal estado. Así, se parte de las decisiones de diseño que mediante mera estática aseguran un equilibrio de la estructura para las cargas últimas que deba resistir, pero determinando directamente el valor de las deformaciones necesarias para llegar a tal estado. Por tanto, no se acude a los teoremas de rotura sin más, sino que se formula el problema desde un punto de vista elastoplástico. Es decir, no se obvia el recorrido que la estructura deba realizar en un proceso de carga incremental monótono, de modo que las regiones no plastificadas contribuyen a coaccionar las libres deformaciones plásticas que, en la teoría de rotura, se suponen. En términos de trabajo y energía, se introduce en el balance del trabajo de las fuerzas externas y en el de la energía de deformación, aquella parte del sistema que no ha plastificado. Establecido así el balance energético como potencial del sistema es cuando la condición de estacionariedad del mismo hace determinados los campos de desplazamientos y, por tanto, el de las deformaciones plásticas también. En definitiva, se trata de un modo de verificar si la ductilidad de los diseños previstos es suficiente, y en qué medida, para verificar el estado de rotura previsto, para unas determinadas cargas impuestas. Dentro del desarrollo teórico del problema, se encuentran ciertas precisiones importantes. Entre ellas, la verificación de que el estado de rotura a que se llega de manera determinada mediante el balance energético elasto-plástico satisface las condiciones de la solución de rotura que los teoremas de carga última predicen, asegurando, por tanto, que la solución determinada -unicidad del problema elásticocoincide con el teorema de unicidad de la carga de rotura, acotando además cuál es el sistema de equilibrio y cuál es la deformada de colapso, aspectos que los teoremas de rotura no pueden asegurar, sino sólo el valor de la carga última a verificar. Otra precisión se basa en la particularidad de los casos en que el sistema presenta una superficie de rotura plana, haciendo infinitas las posibilidades de equilibrio para una misma deformada de colapso determinada, lo que está en la base de, aparentemente, poder plastificar a antojo en vigas y arcos. Desde el planteamiento anterior, se encuentra entonces que existe una condición inherente a cualquier sistema, definidas unas leyes constitutivas internas, que permite al mismo llegar al inicio del estado de rotura sin demandar deformación plástica alguna, produciéndose la plastificación simultánea de todas las regiones que hayan llegado a su solicitación de rotura. En cierto modo, se daría un colapso de apariencia frágil. En tal caso, el sistema conserva plenamente hasta el final su capacidad dúctil y tal estado actúa como representante canónico de cualquier otra solución de equilibrio que con idéntico criterio de diseño interno se prevea para tal estructura. En la medida que el diseño se acerque o aleje de la solución canónica, la demanda de ductilidad del sistema para verificar la carga última será menor o mayor. Las soluciones que se aparten en exceso de la solución canónica, no verificarán el estado de rotura previsto por falta de ductilidad: la demanda de deformación plástica de alguna región plastificada estará más allá de la capacidad de la misma, revelándose una carga de rotura por falta de ductilidad menor que la que se preveía por mero equilibrio. Para la determinación de las deformaciones plásticas de las rótulas, se ha tomado un modelo formulado mediante el Método de los Elementos de Contorno, que proporciona un campo continuo de desplazamientos -y, por ende, de deformaciones y de tensiones- incluso en presencia de fisuras en el contorno. Importante cuestión es que se formula la diferencia, nada desdeñable, de la capacidad de rotación plástica de las secciones de hormigón armado en presencia de cortante y en su ausencia. Para las rótulas de fábrica, la diferencia se establece para las condiciones de la excentricidad -asociadas al valor relativo de la compresión-, donde las diferencias entres las regiones plastificadas con esfuerzo normal relativo alto o bajo son reseñables. Por otro lado, si bien de manera un tanto secundaria, las condiciones de servicio también imponen un límite al diseño previo en carga última deseado. La plastificación lleva asociadas deformaciones considerables, sean locales como globales. Tal cosa impone que, en estado de servicio, si la plastificación de alguna región lleva asociadas fisuraciones excesivas para el ambiente del entorno, la solución sea inviable por ello. Asimismo, las deformaciones de las estructuras suponen un límite severo a las posibilidades de su diseño. Especialmente en edificación, las deformaciones activas son un factor crítico a la hora de decidirse por una u otra solución. Por tanto, al límite que se impone por razón de ductilidad, se debe añadir el que se imponga por razón de las condiciones de servicio. Del modo anterior, considerando las condiciones de ductilidad y de servicio en cada caso, se puede tasar cada decisión de diseño con la previsión de cuáles serán las consecuencias en su estado de carga última y de servicio. Es decir, conocidos los límites, podemos acotar cuáles son los diseños a priori que podrán satisfacer seguro las condiciones de ductilidad y de servicio previstas, y en qué medida. Y, en caso de no poderse satisfacer, qué correcciones debieran realizarse sobre el diseño previo para poderlas cumplir. Por último, de las consecuencias que se extraen de lo estudiado, se proponen ciertas líneas de estudio y de experimentación para poder llegar a completar o expandir de manera práctica los resultados obtenidos. ABSTRACT This work deals with a main issue for the ultimate load design in reinforced concrete and masonry structures: the actual possibility that needed yield strains to reach a ultimate state could be reached by yielded regions on the structure that should develop their ultimate capacity to fulfill such a state. Thus, some statically determined design decisions are posed as a start for prescribed ultimate loads to be counteracted, but finding out the determined value of the strains needed to reach the ultimate load state. Therefore, ultimate load theorems are not taken as they are, but a full elasto-plastic formulation point of view is used. As a result, the path the structure must develop in a monotonus increasing loading procedure is not neglected, leading to the fact that non yielded regions will restrict the supposed totally free yield strains under a pure ultimate load theory. In work and energy terms, in the overall account of external forces work and internal strain energy, those domains in the body not reaching their ultimate state are considered. Once thus established the energy balance of the system as its potential, by imposing on it the stationary condition, both displacements and yield strains appear as determined values. Consequently, what proposed is a means for verifying whether the ductility of prescribed designs is enough and the extent to which they are so, for known imposed loads. On the way for the theoretical development of the proposal, some important aspects have been found. Among these, the verification that the conditions for the ultimate state reached under the elastoplastic energy balance fulfills the conditions prescribed for the ultimate load state predicted through the ultimate load theorems, assuring, therefore, that the determinate solution -unicity of the elastic problemcoincides with the unicity ultimate load theorem, determining as well which equilibrium system and which collapse shape are linked to it, being these two last aspects unaffordable by the ultimate load theorems, that make sure only which is the value of the ultimate load leading to collapse. Another aspect is based on the particular case in which the yield surface of the system is flat -i.e. expressed under a linear expression-, turning out infinite the equilibrium possibilities for one determined collapse shape, which is the basis of, apparently, deciding at own free will the yield distribution in beams and arches. From the foresaid approach, is then found that there is an inherent condition in any system, once defined internal constitutive laws, which allows it arrive at the beginning of the ultimate state or collapse without any yield strain demand, reaching the collapse simultaneously for all regions that have come to their ultimate strength. In a certain way, it would appear to be a fragile collapse. In such a case case, the system fully keeps until the end its ductility, and such a state acts as a canonical representative of any other statically determined solution having the same internal design criteria that could be posed for the that same structure. The extent to which a design is closer to or farther from the canonical solution, the ductility demand of the system to verify the ultimate load will be higher or lower. The solutions being far in excess from the canonical solution, will not verify the ultimate state due to lack of ductility: the demand for yield strains of any yielded region will be beyond its capacity, and a shortcoming ultimate load by lack of ductility will appear, lower than the expected by mere equilibrium. For determining the yield strains of plastic hinges, a Boundary Element Method based model has been used, leading to a continuous displacement field -therefore, for strains and stresses as well- even if cracks on the boundary are present. An important aspect is that a remarkable difference is found in the rotation capacity between plastic hinges in reinforced concrete with or without shear. For masonry hinges, such difference appears when dealing with the eccentricity of axial forces -related to their relative value of compression- on the section, where differences between yield regions under high or low relative compressions are remarkable. On the other hand, although in a certain secondary manner, serviceability conditions impose limits to the previous ultimate load stated wanted too. Yield means always big strains and deformations, locally and globally. Such a thing imposes, for serviceability states, that if a yielded region is associated with too large cracking for the environmental conditions, the predicted design will be unsuitable due to this. Furthermore, displacements must be restricted under certain severe limits that restrain the possibilities for a free design. Especially in building structures, active displacements are a critical factor when chosing one or another solution. Then, to the limits due to ductility reasons, other limits dealing with serviceability conditions shoud be added. In the foresaid way, both considering ductility and serviceability conditions in every case, the results for ultimate load and serviceability to which every design decision will lead can be bounded. This means that, once the limits are known, it is possible to bound which a priori designs will fulfill for sure the prescribed ductility and serviceability conditions, and the extent to wich they will be fulfilled, And, in case they were not, which corrections must be performed in the previous design so that it will. Finally, from the consequences derived through what studied, several study and experimental fields are proposed, in order to achieve a completeness and practical expansion of the obtained results.

Avaliação da integridade estrutural de elementos de concreto armado a partir das propriedades modais obtidas por técnicas de excitações aleatórias e transientes

A análise dinâmica experimental tem sido amplamente pesquisada como uma ferramenta de avaliação de integridade de estruturas de concreto armado. Existem técnicas de identificação de danos baseadas em propriedades modais como frequências de ressonâncias, deformadas modais, curvaturas modais e amortecimento. Há também técnicas baseadas na não linearidade da resposta dinâmica, que apesar do grande potencial na detecção de danos, têm sido pouco exploradas nos últimos anos. Este trabalho tem por objetivo avaliar a integridade estrutural de vigas de concreto armado através do comportamento da resposta dinâmica. Foram realizados ensaios dinâmicos em duas vigas de concreto armado com 3,5 m de comprimento, 25 cm de largura, 35 cm de altura e idênticas taxas de armaduras, mas configuradas com barras de aço de diferentes diâmetros, 2 ϕ 16 mm e 8 ϕ 8 mm, respectivamente. Tais vigas, inicialmente íntegras, foram submetidas a ciclos de carregamento e descarregamento com intensidades crescentes até atingir a ruptura do elemento. Após cada ciclo, as propriedades dinâmicas foram avaliadas experimentalmente, com o emprego de técnicas de excitação por sinais do tipo aleatório e tipo transiente, respectivamente, visando determinar parâmetros que indiquem a deterioração gradativa do elemento. Nesses ensaios dinâmicos aplicaram-se diferentes amplitudes da força de excitação. Verificou-se que o aumento da amplitude da força dinâmica de excitação provocou reduções nos valores das frequências de ressonância de 1,1% e 2,4%, associadas, respectivamente, às excitações aleatórias e transientes; e um comportamento não linear dos índices de amortecimento, associados às excitações aleatórias, mantendo um crescimento linear com as excitações transientes. Constatou-se, ainda, que os valores das frequências de ressonância decrescem com a redução de rigidez mecânica, diminuída com o aumento do nível de fissuração induzido nos modelos. Já os valores dos índices de amortecimento, após cada ciclo, se comportaram de forma não linear e assumiram diferentes valores, conforme a técnica de excitação empregada. Acredita-se que esta não linearidade está relacionada aos danos provocados no elemento pela solicitação estrutural e, por consequência, ao processo de como a dissipação de energia é empregada no processo de instauração, configuração e propagação das fissuras nos elementos de concreto armado.

Influência do tempo de imersão em solução aquosa contendo H2S sobre a tenacidade de tubo API 5L X65 sour avaliada a partir de ensaio Charpy

Com o decorrer dos anos o consumo de petróleo e seus derivados aumentou significativamente e com isso houve a necessidade de se investir em pesquisas para descobertas de novas jazidas de petróleo como o pré-sal. Porém, não apenas a localização dessas jazidas deve ser estudada, mas, também, sua forma de exploração. Essa exploração e extração, na maioria das vezes, se dão em ambientes altamente corrosivos e o transporte do produto extraído é realizado através de tubulações de aço de alta resistência e baixa liga (ARBL). Aços ARBL expostos a ambientes contendo H2S e CO2 (sour gas) sofrem corrosão generalizada que promovem a entrada de hidrogênio atômico no metal, podendo diminuir sua tenacidade e causar falha induzida pela presença de hidrogênio (Hydrogen Induced Cracking HIC), gerando falhas graves no material. Tais falhas podem ser desastrosas para o meio ambiente e para a sociedade. O objetivo deste trabalho é estudar a tenacidade, utilizando ensaio Charpy, de um tubo API 5L X65 sour após diferentes tempos de imersão em uma solução saturada com H2S. O eletrólito empregado foi a solução A (ácido acético contendo cloreto de sódio) da norma NACE TM0284 (2011), fazendo-se desaeração com injeção de N2, seguida de injeções de H2S. Os materiais foram submetidos a: ensaios de resistência a HIC segundo a norma NACE TM0284 (2011) e exames em microscópio óptico e eletrônico de varredura para caracterização microestrutural, de inclusões e trincas. As amostras foram submetidas a imersão em solução A durante 96h e 360h, sendo que, após doze dias do término da imersão, foram realizados os ensaios Charpy e exames fractográficos. Foram aplicados dois métodos: o de energia absorvida e o da expansão lateral, conforme recomendações da norma ASTM E23 (2012). As curvas obtidas, em função da temperatura de impacto, foram ajustadas pelo método da tangente hiperbólica. Esses procedimentos foram realizados nas duas seções do tubo (transversal e longitudinal) e permitiram a obtenção dos seguintes parâmetros: energias absorvidas e expansão lateral nos patamares superior e inferior e temperaturas de transição dúctil-frágil (TTDF) em suas diferentes definições, ou seja, TTDFEA, TTDFEA-DN, TTDFEA-FN, TTDFEL, TTDFEL-DN e TTDFEL-FN (identificação no item Lista de Abreviaturas e Siglas). No exame fractográfico observou-se que o material comportou-se conforme o previsto, ou seja, em temperaturas mais altas ocorreu fratura dúctil, em temperaturas próximas a TTDF obteve-se fratura mista e nas temperaturas mais baixas observou-se o aparecimento de fratura frágil. Os resultados mostraram que quanto maior o tempo de imersão na solução A, menor é a energia absorvida e a expansão lateral no patamar superior, o que pode ser explicado pelo (esperado) aumento do teor de hidrogênio em solução sólida com o tempo de imersão. Por sua vez, os resultados mostraram que há tendência à diminuição da temperatura de transição dúctil-frágil com o aumento do tempo de imersão, particularmente, as TTDFEA-DN e TTDFEL-DN das duas seções do tubo (longitudinal e transversal). Esse comportamento controverso, que pode ser denominado de tenacificação com o decorrer do tempo de imersão na solução A, foi explicado pelo aparecimento de trincas secundárias durante o impacto (Charpy). Isso indica uma limitação do ensaio Charpy para a avaliação precisa de materiais hidrogenados.

Considerações sobre os métodos de dimensionamento estrutural de pavimentos com camadas recicladas a frio com espuma de asfalto.

A reciclagem profunda com espuma de asfalto tem sido uma alternativa de sucesso para a restauração de pavimentos degradados. Em relação às soluções tradicionais de reabilitação, como os recapeamentos, tem a vantagem de proporcionar a correção de defeitos em camadas inferiores, com a manutenção ou pequena elevação do greide da pista, além de ganhos ambientais, como um menor consumo de materiais virgens da natureza e redução do volume de material descartado. Entretanto, no Brasil não há método para dimensionamento estrutural para esta tecnologia, o que dificulta seu emprego. Para o desenvolvimento de um procedimento de dimensionamento que contemple este tipo de solução, foram estudados métodos presentes na bibliografia internacional: guia da AASHTO de 1993 e Caltrans, dos EUA, TRL386 e TRL611, da Inglaterra, as duas versões do guia sul-africano TG2 e os métodos oriundos do Austroads, tanto o procedimento interino de 2011 como adaptações de órgãos da Austrália e Nova Zelândia. Observou-se divergência de opiniões quanto ao comportamento do material reciclado com espuma de asfalto. Alguns órgãos e autores consideram o comportamento do mesmo mais próximo às misturas asfálticas, sendo o mecanismo de falha o trincamento, e outros o definem como semelhante a um material granular modificado com alta coesão e ruptura devido às deformações permanentes. Correlaciona-se tal associação ao teor de espuma usualmente utilizado nas obras rodoviárias. Outros aspectos que se destacam para este tipo de base são o ganho de resistência ao longo do tempo devido à cura, mesmo com início da operação da rodovia e a importância da infraestrutura remanescente no dimensionamento. Tais fatos foram corroborados pelos estudos de caso e resultados do trecho experimental construído na Rodovia Ayrton Senna - SP 070, monitorado por meio de ensaios deflectométricos com FWD durante um ano. Como resultado do trabalho, foi proposto um procedimento para o dimensionamento estrutural de pavimentos com camadas recicladas a frio com espuma de asfalto utilizando dados deflectométricos que atende o método do Manual de Pavimentação do Departamento Nacional de Infraestrutura de Transportes (DNIT) e incorpora diferentes aspectos na análise mecanicista. Outras conclusões são a viabilidade técnica a longo prazo da solução mencionada e a importância do controle tecnológico, com ênfase para o monitoramento deflectométrico nos primeiros meses de operação do pavimento para averiguar a evolução da cura do material.

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.

Evidence of Intracrystalline Mesostructured Porosity in Zeolites by Advanced Gas Sorption, Electron Tomography and Rotation Electron Diffraction

The small size of micropores (typically <1 nm) in zeolites causes slow diffusion of reactant and product molecules in and out of the pores and negatively impacts the product selectivity of zeolite based catalysts, for example, fluid catalytic cracking (FCC) catalysts. Size-tailored mesoporosity was introduced into commercial zeolite Y crystals by a simple surfactant-templating post-synthetic mesostructuring process. The resulting mesoporous zeolite Y showed significantly improved product selectivity in both laboratory testing and refinery trials. Advanced characterization techniques such as electron tomography, three-dimensional rotation electron diffraction, and high resolution gas adsorption coupled with hysteresis scanning and density functional theory, unambiguously revealed the intracystalline nature and connectivity of the introduced mesopores. They can be considered as molecular highways that help reactant and product molecules diffuse quickly to and away from the catalytically active sites within the zeolite crystals and, thus, shift the selectivity to favor the production of more of the valuable liquid fuels at reduced yields of coke and unconverted feed.

Ethanol dehydration via azeotropic distillation with gasoline fraction mixtures as entrainers: A pilot-scale study with industrially produced bioethanol and naphta

Various hydrocarbons (n-hexane, cyclohexane, toluene, isooctane) and mixtures of them (binary, ternary or quaternary), as well as two different types of industrially produced naphtha (one obtained by direct distillation and the other from a catalytic cracking process), have been tested as candidate entrainers to dehydrate ethanol. The tests were carried out in an azeotropic distillation column on a semi pilot plant. The results show that it is possible to dehydrate bioethanol using naphtha as entrainer, obtaining as a result a fuel blend with negligible water content and ready for immediate use in motor vehicles.

Mesostructured zeolites: bridging the gap between zeolites and MCM-41

Surfactant-templating is one of the most versatile and useful techniques to implement mesoporous systems into solid materials. Various strategies based on various interactions between surfactants and solid precursors have been explored to produce new structures. Zeolites are invaluable as size- and shape-selective solid acid catalysts. Nevertheless, their micropores impose limitations on the mass transport of bulky feed and/or product molecules. Many studies have attempted to address this by utilizing surfactant-assisting technology to alleviate the diffusion constraints. However, most efforts have failed due to micro/mesopore phase separation. Recently, a new technique combining the uses of cationic surfactants and mild basic solutions was introduced to synthesise mesostructured zeolites. These materials sustain the unique characteristics of zeolites (i.e., strong acidity, crystallinity, microporosity, and hydrothermal stability), including tunable mesopore sizes and degrees of mesoporosity. The mesostructured zeolites are now commercially available through Rive Technology, and show superior performance in VGO cracking. This feature article provides an overview of recent explorations in the introduction of mesoporosity into zeolites using surfactant-templating techniques. Various porous materials, preparation methods, physical and catalytic properties of mesostructured zeolites will be discussed.

Concreto com agregados reciclados de borracha de pneu: resistência à compressão a altas temperaturas

A sustentabilidade é uma preocupação para a indústria da construção civil, uma vez que é responsável pelo consumo de uma grande quantidade de recursos naturais e por impactos ambientais. Assim, a utilização de agregados reciclados em substituição dos agregados naturais mostra-se benéfica ao minimizar os impactos ambientais, o consumo de recursos naturais e na redução de alguns problemas urbanos associados à acumulação de lixo. Neste contexto, o trabalho de investigação desenvolvido teve como objetivo o estudo da utilização de agregados reciclados de borracha de pneu na composição do concreto, contribuindo com uma alternativa sustentável para o problema do depósito de pneus após a sua vida útil. O trabalho laboratorial realizado compreendeu ensaios de resistência à compressão do concreto. Estudaram-se três composições de concreto, uma composição de referência, uma composição com uma taxa de substituição de 15% e outra com taxa de substituição de 30% de agregados naturais por agregados reciclados de borracha de pneu usados. Os diferentes provetes de concreto foram submetidos a vários níveis de carregamento (0,15fcd; 0,3fcd e 0,7fcd) e a diferentes níveis de temperatura (20, 300, 500 e 700ºC). Os resultados obtidos permitiram verificar que o aumento da percentagem de agregados de borracha reciclados de pneu inserido no concreto conduz a um aumento do controlo de fendilhação e minimiza o surgimento de fissuração de origem térmica.

Turkey votes: Part II - Erdoğan reloaded. EPC Commentary, 3 April 2014

On 30 March, Turkey’s ruling Justice and Development Party (AKP) scooped a significant victory in local elections, taking almost 44 percent of the vote despite accusations of corruption, undermining the rule of law, fundamental rights and freedoms. While there have been claims of election fraud and the main opposition party, the Republican People’s Party (CHP), has demanded recounts in several cities including Istanbul and Ankara, it is clear that even allowing for some level of fraud the win was substantial and more than most people expected. Prime Minister Recep Tayyip Erdoğan has reached a juncture. He has two choices: return to the path of democracy after a period of democratic back-sliding which included passing several controversial reforms such as a new internet law which led to the recent banning of Twitter and Youtube; or alternatively he can forge ahead with his much talked of revenge campaign against those he has accused of creating a “parallel state” and conspiring to remove him from power. Given that Erdoğan viewed this election as a referendum on his popularity and leadership there is a serious risk that he will do the latter; using the significant mandate given to him to do whatever he wants, including further cracking down on democracy.

Mean grain size and structure of ODP Hole 140-504B samples

Microstructural investigations of ocean crust samples provide a complementary approach to both marine surveys and laboratory experiments. The recovery of relatively undeformed diabases from Deep Sea Drilling Project (DSDP)/Ocean Drilling Program (ODP) Hole 504B provides a first opportunity to examine a reference section of microstructural features that influence strain localization at depths of 2 km in the ocean crust. Syn- and post-crystallization features in plagioclase and augite crystals have been examined by optical microscopy and secondary and backscattered electron imaging. These features show a strong influence of modal composition and primary textures on early sites of strain localization. Thermal cracking and subsequent alteration intensities and distribution are strongly phase dependent. A consistently higher intragranular fracture density is observed in augite crystals relative to plagioclase. The impact of alteration on the mechanical response of diabases is likely to depend on the primary textural characteristics. Even where extensive augite alteration occurs, the rock remains supported by a framework of weakly altered plagioclase crystals. The Hole 504B diabases from Leg 140 provide a valuable comparison for future studies of more deformed sections likely to be encountered at depth. Advances in constraining the detailed rheology of the ocean crust at spreading centers would benefit from experimental deformation of texturally diverse diabase and gabbro samples.

C4-C7 alkenones in deep sea sediments

It is believed that C4 to C7 hydrocarbons in petroleum are formed by the cracking of organic matter at depths generally exceeding 1,000 m at temperatures in excess of 50 °C (Cordel, 1972; Dow, 1974; Tissot et al., 1974)). Also, none of the alkanes in the butane-heptane range are formed biologically as far as is known at present. Consequently, it is thought that they do not occur in shallow, Recent sediments. In 1962, I analysed 22 samples of Recent sediments from 7 different environments and verified that these hydrocarbons were not present at the p.p.m. level (Dunton and Hunt, 1962) although traces of a few hydrocarbons such as butane, isobutane, isopentane and n-heptane have been found (Sokolov, 1957; Veber and Turkeltaub, 1958; Erdman et al., 1958; Emery and Hoggan, 1958). No identification of individual hexanes or heptanes has been reported except when there has been clear evidence of seepage from deeper source sediments (McIver, 1973).