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.

Microglial activation precedes acute neurodegeneration in Sandhoff disease and is suppressed by bone marrow transplantation

Sandhoff disease is a lysosomal storage disorder characterized by the absence of β-hexosaminidase and storage of GM2 ganglioside and related glycolipids in the central nervous system. The glycolipid storage causes severe neurodegeneration through a poorly understood pathogenic mechanism. In symptomatic Sandhoff disease mice, apoptotic neuronal cell death was prominent in the caudal regions of the brain. cDNA microarray analysis to monitor gene expression during neuronal cell death revealed an upregulation of genes related to an inflammatory process dominated by activated microglia. Activated microglial expansion, based on gene expression and histologic analysis, was found to precede massive neuronal death. Extensive microglia activation also was detected in a human case of Sandhoff disease. Bone marrow transplantation of Sandhoff disease mice suppressed both the explosive expansion of activated microglia and the neuronal cell death without detectable decreases in neuronal GM2 ganglioside storage. These results suggest a mechanism of neurodegeneration that includes a vigorous inflammatory response as an important component. Thus, this lysosomal storage disease has parallels to other neurodegenerative disorders, such as Alzheimer's and prion diseases, where inflammatory processes are believed to participate directly in neuronal cell death.

Evolution and divergence of sodium channel genes in vertebrates

Invertebrate species possess one or two Na+ channel genes, yet there are 10 in mammals. When did this explosive growth come about during vertebrate evolution? All mammalian Na+ channel genes reside on four chromosomes. It has been suggested that this came about by multiple duplications of an ancestral chromosome with a single Na+ channel gene followed by tandem duplications of Na+ channel genes on some of these chromosomes. Because a large-scale expansion of the vertebrate genome likely occurred before the divergence of teleosts and tetrapods, we tested this hypothesis by cloning Na+ channel genes in a teleost fish. Using an approach designed to clone all of the Na+ channel genes in a genome, we found six Na+ channel genes. Phylogenetic comparisons show that each teleost gene is orthologous to a Na+ channel gene or gene cluster on a different mammalian chromosome, supporting the hypothesis that four Na+ channel genes were present in the ancestors of teleosts and tetrapods. Further duplications occurred independently in the teleost and tetrapod lineages, with a greater number of duplications in tetrapods. This pattern has implications for the evolution of function and specialization of Na+ channel genes in vertebrates. Sodium channel genes also are linked to homeobox (Hox) gene clusters in mammals. Using our phylogeny of Na+ channel genes to independently test between two models of Hox gene evolution, we support the hypothesis that Hox gene clusters evolved as (AB) (CD) rather than {D[A(BC)]}.

Primordial emergence of the recombination activating gene 1 (RAG1): sequence of the complete shark gene indicates homology to microbial integrases.

The rearrangement of antibody and T-cell receptor gene segments is indispensable to the vertebrate immune response. All extant jawed vertebrates can rearrange these gene segments. This ability is conferred by the recombination activating genes I and II (RAG I and RAG II). To elucidate their origin and function, the cDNA encoding RAG I from a member of the most ancient class of extant gnathostomes, the Carcharhine sharks, was characterized. Homology domains identified within shark RAG I prompted sequence comparison analyses that suggested similarity of the RAG I and II genes, respectively, to the integrase family genes and integration host factor genes of the bacterial site-specific recombination system. Thus, the apparent explosive evolution (or "big bang") of the ancestral immune system may have been initiated by a transfer of microbial site-specific recombinases.

Population dynamics of flaviviruses revealed by molecular phylogenies.

The phylogeny of 123 complete envelope gene sequences was reconstructed in order to understand the evolution of tick- and mosquito-borne flaviviruses. An analysis of phylogenetic tree structure reveals a continual and asymmetric branching process in the tick-borne flaviviruses, compared with an explosive radiation in the last 200 years in viruses transmitted by mosquitoes. The distinction between these two viral groups probably reflects differences in modes of dispersal, propagation, and changes in the size of host populations. The most serious implication of this work is that growing human populations are being exposed to an expanding range of increasingly diverse viral strains.

Research in speech communication.

Advances in digital speech processing are now supporting application and deployment of a variety of speech technologies for human/machine communication. In fact, new businesses are rapidly forming about these technologies. But these capabilities are of little use unless society can afford them. Happily, explosive advances in microelectronics over the past two decades have assured affordable access to this sophistication as well as to the underlying computing technology. The research challenges in speech processing remain in the traditionally identified areas of recognition, synthesis, and coding. These three areas have typically been addressed individually, often with significant isolation among the efforts. But they are all facets of the same fundamental issue--how to represent and quantify the information in the speech signal. This implies deeper understanding of the physics of speech production, the constraints that the conventions of language impose, and the mechanism for information processing in the auditory system. In ongoing research, therefore, we seek more accurate models of speech generation, better computational formulations of language, and realistic perceptual guides for speech processing--along with ways to coalesce the fundamental issues of recognition, synthesis, and coding. Successful solution will yield the long-sought dictation machine, high-quality synthesis from text, and the ultimate in low bit-rate transmission of speech. It will also open the door to language-translating telephony, where the synthetic foreign translation can be in the voice of the originating talker.

Screen-printed electrode-based electrochemical detector coupled with in-situ ionic-liquid-assisted dispersive liquid–liquid microextraction for determination of 2,4,6-trinitrotoluene

A novel method is reported, whereby screen-printed electrodes (SPELs) are combined with dispersive liquid–liquid microextraction. In-situ ionic liquid (IL) formation was used as an extractant phase in the microextraction technique and proved to be a simple, fast and inexpensive analytical method. This approach uses miniaturized systems both in sample preparation and in the detection stage, helping to develop environmentally friendly analytical methods and portable devices to enable rapid and onsite measurement. The microextraction method is based on a simple metathesis reaction, in which a water-immiscible IL (1-hexyl-3-methylimidazolium bis[(trifluoromethyl)sulfonyl]imide, [Hmim][NTf2]) is formed from a water-miscible IL (1-hexyl-3-methylimidazolium chloride, [Hmim][Cl]) and an ion-exchange reagent (lithium bis[(trifluoromethyl)sulfonyl]imide, LiNTf2) in sample solutions. The explosive 2,4,6-trinitrotoluene (TNT) was used as a model analyte to develop the method. The electrochemical behavior of TNT in [Hmim][NTf2] has been studied in SPELs. The extraction method was first optimized by use of a two-step multivariate optimization strategy, using Plackett–Burman and central composite designs. The method was then evaluated under optimum conditions and a good level of linearity was obtained, with a correlation coefficient of 0.9990. Limits of detection and quantification were 7 μg L−1 and 9 μg L−1, respectively. The repeatability of the proposed method was evaluated at two different spiking levels (20 and 50 μg L−1), and coefficients of variation of 7 % and 5 % (n = 5) were obtained. Tap water and industrial wastewater were selected as real-world water samples to assess the applicability of the method.

An Electronic Reconfigurable Neural Architecture for Intrusion Detection

The explosive growth of the traffic in computer systems has made it clear that traditional control techniques are not adequate to provide the system users fast access to network resources and prevent unfair uses. In this paper, we present a reconfigurable digital hardware implementation of a specific neural model for intrusion detection. It uses a specific vector of characterization of the network packages (intrusion vector) which is starting from information obtained during the access intent. This vector will be treated by the system. Our approach is adaptative and to detecting these intrusions by using a complex artificial intelligence method known as multilayer perceptron. The implementation have been developed and tested into a reconfigurable hardware (FPGA) for embedded systems. Finally, the Intrusion detection system was tested in a real-world simulation to gauge its effectiveness and real-time response.

From Lockerbie to Umar Farouk AbdulMutallab (Assessing the Destruction of Civilian Aircrafts by Terrorism). Transatlantic Security Paper N.1–April 2010

Data published by the Federal Aviation Administration (F.A.A.) in its "Annual report on the criminal acts against civil aviation" indicates that in the year 1960, there have been fifteen attacks on board planes leaving 286 dead; 44 in the year 1970 with 650 dead (mostly hijackings); and 26 in the 1980s leaving 1207 dead. In the 1970s, the record is established by the year 1976 (168 dead).The three years - 1985 (390), 1988 (287), and 1989 (278) - were more deadly than the 1960s all together. These casualties were largely provoked by IEDs. Since the end of 1980 - the deadliest decade, with the exception of September 11, 2001 -, it is however a rare practice. The reasons for the decline in big and politically motivated hijackings were varied. One could have been the improvement of the effectiveness of the safety response by States, airports and companies. The improvised explosive device (I.E.D.) posed a serious threat to the civil aviation industry in the 1980s. Since the 1990s, the jihadi networks have regularly tried to target aircrafts using various types of IEDs.

Sistema de gestão de energia e ativação de uma espoleta eletrónica

Os incêndios florestais constituem uma realidade que pode ser minimizada com o auxílio Forças Armadas, através de homens e de meios. O Projeto FIREND®, no qual este trabalho está inserido, tem como objetivo projetar uma munição que ao libertar uma substância química sobre o fogo, o impeça de se expandir e se possível, o extinga. O projétil tem 155mm de calibre e o seu compartimento de carga terá cerca de 7,5

Estudo de balística interna

O conhecimento e a compreensão dos processos que ocorrem no interior do tubo de uma arma de fogo, desde a ignição da carga propulsora, até à saída do projétil, é importante no desenvolvimento do armamento, de uma forma económica e segura. Com a modelação matemática destes processos, que é uma abordagem muito utilizada atualmente, é possível otimizar a carga explosiva e os materiais utilizados. Este estudo foi realizado no âmbito do projeto FIREND, resultante da colaboração entre a Academia Militar e o Instituto Superior Técnico O modelo desenvolvido foi aplicado a um sistema, arma, carga propulsora e projétil, já existente, tendo os resultados obtidos sido validados com o software PRODAS e tabelas referentes ao presente sistema. Para melhor compreender a importância das variáveis envolvidas, é feita uma análise de sensibilidade a diversos parâmetros balísticos, mantendo constante as dimensões da arma em estudo. Na análise dos resultados verifica-se que a cinta de travamento é um elemento muito importante, alterando significativamente os valores da pressão e da velocidade. Verificou-se também que com cargas propulsoras de maior taxa de queima é possível obter valores mais elevados de pressão máxima e velocidade à saída. Palavras chave: Balística interna; Projétil de Artilharia; Pressão; Velocidade; Modelação matemática; Análise de sensibilidade.

Ash layers, hyaloclastite, and alteration of basaltic glass at DSDP Leg 65 Holes

Three types of tephra deposits were recovered on Leg 65 of the Deep Sea Drilling Project (DSDP) from three drill sites at the mouth of the Gulf of California: (1) a series of white ash layers at Sites 483, 484, and 485; (2) a layer of plagioclase- phyric sideromelane shards at Site 483; and (3) an indurated, cross-bedded hyaloclastite in Hole 483B. The ash layers in (1) are composed of colorless, fresh rhyolitic glass shards with minor dacitic and rare basaltic shards. These are thought to be derived from explosive volcanoes on the Mexican mainland. Most of the shards in (2) are fresh, but some show marginal to complete alteration to palagonite. The composition of the glass is that of a MORB-type tholeiite, low in Fe and moderately high in Ti, and possibly erupted from off-axis seamounts. Basaltic glass shards occurring in silt about 45 meters above the basement at Site 484 A in the Tamayo Fracture Zone show a distinctly alkalic composition similar to that of the single basement basalt specimen drilled at this site. The hyaloclastite in (3) is made up chiefly of angular sideromelane shards altered to smectite and zeolites (mainly phillipsite) and minor admixtures of terrigenous silt. A very high K and Ba content indicates significant uptake of at least these elements from seawater. Nevertheless, the unusual chemical composition of the underlying massive basalt flow is believed to be reflected in that of the hyaloclastite. This is a powerful argument for interpreting the massive basalt as a surface flow rather than an intrusion. Glass alteration is different in the glassy margins of flows than in thicker glassy pillow rinds. Also, it appears to proceed faster in coarse- than fine-grained sediments.

Petrology and geochemistry of sideromelane glass shards from Pleistocene as layers north and south of Gran Canaria

Major elements, S, F, Cl concentrations and relative proportions of S6+ to total S were analyzed with electron microprobe in sideromelane glass shards from Pleistocene volcaniclastic sediments drilled during ODP Leg 157. Glasses are moderately to strongly evolved and represent a spectrum from alkali basalt, basanite and nephelinite through hawaiite, mugearite and tephrite to phonolitic tephrite. Measured S6+/SumS (0.03±0.98) and calculated Fe2+/Fe3+ (2.5±5.8) ratios in the melt yield preeruptive redox conditions ranging from NNO-1.4 to NNO+2.1. The morphology of the glass shards, variations of S and Cl concentrations (0.010±0.127 wt% S, 0.018±0.129 wt% Cl), calculated preeruptive temperatures (1030±1200 °C) and oxygen fugacities suggest that glasses deposited even within the same ash layers have diverse origin and may have resulted from both submarine and subaerial eruptions. Most vesicle-free glasses are characterized by high concentrations of S and represent undegassed or slightly degassed submarine lavas, whereas vesiculated glasses with low concentrations of S and Cl are strongly degassed and can be ascribed to the eruptions in shallow water or on land. Sideromelane glass shards at Sites 953 are thought to have resulted from submarine eruptions northeast of Gran Canaria, glasses at Site 954 represent mostly volcaniclastic material of shallow water submarine and subaerial eruptions on Gran Canaria and Tenerife, and glasses deposited at Site 956 resulted from submarine or explosive eruptions on Tenerife.

(Table1) Sedimentation rates and age-control points form ODP Hole 152-919A

Planktic d18O and d13C records and point count records of biogenic, volcanic, and nonvolcanic terrigenous [ice-rafted debris (IRD)] sediment components from Hole 919A in the Irminger basin, northern North Atlantic provide a comprehensive dataset from which a paleoceanographic reconstruction for the last 630 kyr has been developed. The paleoceanographic evolution of the Irminger basin during this time contains both long-term patterns and significant developmental steps. One long-term pattern observed is the persistent deposition of hematite-stained ice-rafted debris. This record suggests that the modern and late Pleistocene discharges of icebergs from northern redbed regions to the Irminger Sea lie in the low end of the range observed over the last 630 kyr. In addition, Arctic front fluctuations appear to have been the main controlling factor on the long-term accumulation patterns of IRD and planktic biogenic groups. The Hole 919A sediment record also contains a long-term association between felsic volcanic ash abundances and light d18O excursions in both interglacial and glacial stages, which suggests a causal link between deglaciations and explosive Icelandic eruptions. A significant developmental step in the paleoceanographic reconstruction based on benthic evidence was for diminished supply of Denmark Strait Overflow Water (DSOW) beginning at ~380 ka, possibly initiated by the influx of meltwater from broad-scale iceberg discharges along the east Greenland coast. There is also planktic evidence of a two-step cooling of sea surface conditions in the Irminger basin, first at ~338-309 ka and later at ~211-190 ka, after which both glacials and interglacials were colder as the Arctic front migrated southeast of Site 919. In addition to offering these findings, this reconstruction provides a longer-term geologic context for the interpretation of more recent paleoceanographic events and patterns of deposition from this region.