SteelConverter and Caltech VirtualShaker: Rapid Nonlinear Cloud-Based Structural Model Conversion and Analysis

STEEL, the Caltech created nonlinear large displacement analysis software, is currently used by a large number of researchers at Caltech. However, due to its complexity, lack of visualization tools (such as pre- and post-processing capabilities) rapid creation and analysis of models using this software was difficult. SteelConverter was created as a means to facilitate model creation through the use of the industry standard finite element solver ETABS. This software allows users to create models in ETABS and intelligently convert model information such as geometry, loading, releases, fixity, etc., into a format that STEEL understands. Models that would take several days to create and verify now take several hours or less. The productivity of the researcher as well as the level of confidence in the model being analyzed is greatly increased.

It has always been a major goal of Caltech to spread the knowledge created here to other universities. However, due to the complexity of STEEL it was difficult for researchers or engineers from other universities to conduct analyses. While SteelConverter did help researchers at Caltech improve their research, sending SteelConverter and its documentation to other universities was less than ideal. Issues of version control, individual computer requirements, and the difficulty of releasing updates made a more centralized solution preferred. This is where the idea for Caltech VirtualShaker was born. Through the creation of a centralized website where users could log in, submit, analyze, and process models in the cloud, all of the major concerns associated with the utilization of SteelConverter were eliminated. Caltech VirtualShaker allows users to create profiles where defaults associated with their most commonly run models are saved, and allows them to submit multiple jobs to an online virtual server to be analyzed and post-processed. The creation of this website not only allowed for more rapid distribution of this tool, but also created a means for engineers and researchers with no access to powerful computer clusters to run computationally intensive analyses without the excessive cost of building and maintaining a computer cluster.

In order to increase confidence in the use of STEEL as an analysis system, as well as verify the conversion tools, a series of comparisons were done between STEEL and ETABS. Six models of increasing complexity, ranging from a cantilever column to a twenty-story moment frame, were analyzed to determine the ability of STEEL to accurately calculate basic model properties such as elastic stiffness and damping through a free vibration analysis as well as more complex structural properties such as overall structural capacity through a pushover analysis. These analyses showed a very strong agreement between the two softwares on every aspect of each analysis. However, these analyses also showed the ability of the STEEL analysis algorithm to converge at significantly larger drifts than ETABS when using the more computationally expensive and structurally realistic fiber hinges. Following the ETABS analysis, it was decided to repeat the comparisons in a software more capable of conducting highly nonlinear analysis, called Perform. These analyses again showed a very strong agreement between the two softwares in every aspect of each analysis through instability. However, due to some limitations in Perform, free vibration analyses for the three story one bay chevron brace frame, two bay chevron brace frame, and twenty story moment frame could not be conducted. With the current trend towards ultimate capacity analysis, the ability to use fiber based models allows engineers to gain a better understanding of a building’s behavior under these extreme load scenarios.

Following this, a final study was done on Hall’s U20 structure [1] where the structure was analyzed in all three softwares and their results compared. The pushover curves from each software were compared and the differences caused by variations in software implementation explained. From this, conclusions can be drawn on the effectiveness of each analysis tool when attempting to analyze structures through the point of geometric instability. The analyses show that while ETABS was capable of accurately determining the elastic stiffness of the model, following the onset of inelastic behavior the analysis tool failed to converge. However, for the small number of time steps the ETABS analysis was converging, its results exactly matched those of STEEL, leading to the conclusion that ETABS is not an appropriate analysis package for analyzing a structure through the point of collapse when using fiber elements throughout the model. The analyses also showed that while Perform was capable of calculating the response of the structure accurately, restrictions in the material model resulted in a pushover curve that did not match that of STEEL exactly, particularly post collapse. However, such problems could be alleviated by choosing a more simplistic material model.

Aerodynamics of vertical-axis wind turbines in full-scale and laboratory-scale experiments

Within a wind farm, multiple turbine wakes can interact and have a substantial effect on the overall power production. This makes an understanding of the wake recovery process critically important to optimizing wind farm efficiency. Vertical-axis wind turbines (VAWTs) exhibit features that are amenable to dramatically improving this efficiency. However, the physics of the flow around VAWTs is not well understood, especially as it pertains to wake interactions, and it is the goal of this thesis to partially fill this void. This objective is approached from two broadly different perspectives: a low-order view of wind farm aerodynamics, and a detailed experimental analysis of the VAWT wake.

One of the contributions of this thesis is the development of a semi-empirical model of wind farm aerodynamics, known as the LRB model, that is able to predict turbine array configurations to leading order accuracy. Another contribution is the characterization of the VAWT wake as a function of turbine solidity. It was found that three distinct regions of flow exist in the VAWT wake: (1) the near wake, where periodic blade shedding of vorticity dominates; (2) a transition region, where growth of a shear-layer instability occurs; (3) the far wake, where bluff-body oscillations dominate. The wake transition can be predicted using a new parameter, the dynamic solidity, which establishes a quantitative connection between the wake of a VAWT and that of a circular cylinder. The results provide insight into the mechanism of the VAWT wake recovery and the potential means to control it.

A hybrid-parallel framework for the nonlinear seismic analysis of very tall buildings

FRAME3D, a program for the nonlinear seismic analysis of steel structures, has previously been used to study the collapse mechanisms of steel buildings up to 20 stories tall. The present thesis is inspired by the need to conduct similar analysis for much taller structures. It improves FRAME3D in two primary ways.

First, FRAME3D is revised to address specific nonlinear situations involving large displacement/rotation increments, the backup-subdivide algorithm, element failure, and extremely narrow joint hysteresis. The revisions result in superior convergence capabilities when modeling earthquake-induced collapse. The material model of a steel fiber is also modified to allow for post-rupture compressive strength.

Second, a parallel FRAME3D (PFRAME3D) is developed. The serial code is optimized and then parallelized. A distributed-memory divide-and-conquer approach is used for both the global direct solver and element-state updates. The result is an implicit finite-element hybrid-parallel program that takes advantage of the narrow-band nature of very tall buildings and uses nearest-neighbor-only communication patterns.

Using three structures of varied sized, PFRAME3D is shown to compute reproducible results that agree with that of the optimized 1-core version (displacement time-history response root-mean-squared errors are ~〖10〗^(-5) m) with much less wall time (e.g., a dynamic time-history collapse simulation of a 60-story building is computed in 5.69 hrs with 128 cores—a speedup of 14.7 vs. the optimized 1-core version). The maximum speedups attained are shown to increase with building height (as the total number of cores used also increases), and the parallel framework can be expected to be suitable for buildings taller than the ones presented here.

PFRAME3D is used to analyze a hypothetical 60-story steel moment-frame tube building (fundamental period of 6.16 sec) designed according to the 1994 Uniform Building Code. Dynamic pushover and time-history analyses are conducted. Multi-story shear-band collapse mechanisms are observed around mid-height of the building. The use of closely-spaced columns and deep beams is found to contribute to the building's “somewhat brittle” behavior (ductility ratio ~2.0). Overall building strength is observed to be sensitive to whether a model is fracture-capable.

Petrologia dos diabásios da região de Angra dos Reis, RJ

O tema abordado nesta dissertação de mestrado é a petrogênese dos diabásios do Eocretáceo que integram o Enxame de Diques da Serra do Mar na região de Angra dos Reis (RJ). Estes diques intrudem rochas do Terreno Ocidental e do Arco Magmático Rio Negro, inserindo-se no contexto geodinâmico de fragmentação do Gondwana no Cretáceo Inferior. Essas intrusões são subverticais e orientam-se preferencialmente a N039E. Intrusões retilíneas são as mais observadas em campo. Os diabásios são compostos, essencialmente, por plagioclásio e clinopiroxênios (augita e pigeonita). Quartzo, minerais opacos, apatita e olivina ocorrem como minerais acessórios. A composição mineral secundária é marcada pela presença de biotita, uralita e saussurita. Os diabásios estudados compreendem uma série transicional com afinidade toleítica. Os estudos petrológicos revelam a existência de cinco suítes de alto-TiO2, tendo sido as mesmas discriminadas com base nas razões La/Yb de amostras com o mesmo teor de MgO. Destas, apenas a suíte 2 têm quantidade mínima de amostras para avaliação de processos evolutivos. O resultado da modelagem petrogenética aponta que o processo evolutivo mais provável para a geração dos diabásios da suíte 2 é cristalização fracionada sem mudança na assembleia fracionante. As amostras representativas dos líquidos parentais das cinco suítes apresentam razões La/Nb e La/Yb maiores que a unidade, típicas de derivação de fonte enriquecida, que, em várias províncias basálticas continentais, tem sido associada à fusão exclusiva do manto litosférico subcontinental ou a sua participação em processos de mistura de fontes. A discriminação de fontes mantélicas foi feita exclusivamente com base nos teores de elementos traços de amostras representativas de líquidos parentais das cinco suítes discriminadas na área de estudo. Os modelos indicam que todas as cinco suítes estudadas estão relacionadas a fontes harzburgíticas com granada residual. A suíte 5 da área de estudo não pode ser associada à mesma fonte mantélica geradora das suítes 1, 2, 3 e 4 que, por sua vez, se originaram por diferentes quantidades de fusão parcial (respectivamente, 44%, 35%, 31% e 25%) de uma mesma fonte. Cabe ressaltar que a ocorrência de duas suítes, com razões La/Yb muito distintas, porém geograficamente muito próximas, sugere que processos de heterogeneidade mantélica, tipicamente laterais, devem também ocorrer verticalmente.

Avaliação numérica de estabilidade lateral de vigas casteladas

Restrições de espaço e altura são frequentemente impostas às edificações residenciais, comerciais, industriais, depósitos e galpões com um ou diversos pavimentos em função de aspectos de regulamentos regionais, técnicos, econômicos ou ainda de natureza estética. A fim de proporcionar a passagem de tubulações e dutos de grande diâmetro sob vigas de aço, grandes alturas são normalmente requeridas, demandando por vezes, magnitudes de altura inviáveis entre pavimentos de edificações. Diversas soluções estruturais podem ser utilizadas para equacionar tais obstáculos, onde dentre outras, pode-se citar as vigas com inércia variável, stub-girders, treliças mistas, vigas misuladas e vigas com uma ou múltiplas aberturas na alma com geometrias variadas. No que tange às vigas casteladas, solução estrutural pautada neste estudo, a estabilidade é sempre um motivo de preocupação tipicamente durante a construção quando os contraventamentos laterais ainda não estão instalados. De qualquer forma, o comprimento destravado em geral alcançado pelos vãos destas vigas, são longos o suficiente para que a instabilidade ocorra. Todavia, o acréscimo substancial da resistência à flexão de tais membros devido ao aumento da altura oriundo de seu processo fabril em relação ao perfil matriz, aliada a economia de material e utilidade fim de serviço, garante a atratividade no aproveitamento destas, para grandes vãos junto aos projetistas. Não obstante, este aumento proporcional no comprimento dos vãos faz com que a instabilidade lateral ganhe importância especial. Neste contexto, o presente trabalho tem por objetivo desenvolver um modelo numérico que permita a realização de uma avaliação paramétrica a partir da calibração do modelo com resultados experimentais, efetuar a análise do comportamento de vigas casteladas e verificar seus mecanismos de falha, considerando comportamento elasto-plástico, além das não-linearidades geométricas. Também é objetivo deste trabalho, avaliar, quantificar e determinar a influência das diferenças geométricas características das vigas casteladas em relação às vigas maciças com as mesmas dimensões, analisando e descrevendo o comportamento estrutural destas vigas de aço para diversos comprimentos de vãos. A metodologia empregada para tal estudo baseou-se em uma análise paramétrica com o auxílio do método numérico dos elementos finitos.

Charge and spin transport in graphene devices

170 p.

Physiological modules for generating discrete and rhythmic movements: Component analysis of EMG signals

A central question in Neuroscience is that of how the nervous system generates the spatiotemporal commands needed to realize complex gestures, such as handwriting. A key postulate is that the central nervous system (CNS) builds up complex movements from a set of simpler motor primitives or control modules. In this study we examined the control modules underlying the generation of muscle activations when performing different types of movement: discrete, point-to-point movements in eight different directions and continuous figure-eight movements in both the normal, upright orientation and rotated 90 degrees. To test for the effects of biomechanical constraints, movements were performed in the frontal-parallel or sagittal planes, corresponding to two different nominal flexion/abduction postures of the shoulder. In all cases we measured limb kinematics and surface electromyographic activity (EMB) signals for seven different muscles acting around the shoulder. We first performed principal component analysis (PCA) of the EMG signals on a movement-by-movement basis. We found a surprisingly consistent pattern of muscle groupings across movement types and movement planes, although we could detect systematic differences between the PCs derived from movements performed in each sholder posture and between the principal components associated with the different orientations of the figure. Unexpectedly we found no systematic differences between the figute eights and the point-to-point movements. The first three principal components could be associated with a general co-contraction of all seven muscles plus two patterns of reciprocal activatoin. From these results, we surmise that both "discrete-rhythmic movements" such as the figure eight, and discrete point-to-point movement may be constructed from three different fundamental modules, one regulating the impedance of the limb over the time span of the movement and two others operating to generate movement, one aligned with the vertical and the other aligned with the horizontal.

A-SCALA: an age-structured statistical catch-at-length analysis for assessing tuna stocks in the eastern tropical Pacific Ocean

English: We describe an age-structured statistical catch-at-length analysis (A-SCALA) based on the MULTIFAN-CL model of Fournier et al. (1998). The analysis is applied independently to both the yellowfin and the bigeye tuna populations of the eastern Pacific Ocean (EPO). We model the populations from 1975 to 1999, based on quarterly time steps. Only a single stock for each species is assumed for each analysis, but multiple fisheries that are spatially separate are modeled to allow for spatial differences in catchability and selectivity. The analysis allows for error in the effort-fishing mortality relationship, temporal trends in catchability, temporal variation in recruitment, relationships between the environment and recruitment and between the environment and catchability, and differences in selectivity and catchability among fisheries. The model is fit to total catch data and proportional catch-at-length data conditioned on effort. The A-SCALA method is a statistical approach, and therefore recognizes that the data collected from the fishery do not perfectly represent the population. Also, there is uncertainty in our knowledge about the dynamics of the system and uncertainty about how the observed data relate to the real population. The use of likelihood functions allow us to model the uncertainty in the data collected from the population, and the inclusion of estimable process error allows us to model the uncertainties in the dynamics of the system. The statistical approach allows for the calculation of confidence intervals and the testing of hypotheses. We use a Bayesian version of the maximum likelihood framework that includes distributional constraints on temporal variation in recruitment, the effort-fishing mortality relationship, and catchability. Curvature penalties for selectivity parameters and penalties on extreme fishing mortality rates are also included in the objective function. The mode of the joint posterior distribution is used as an estimate of the model parameters. Confidence intervals are calculated using the normal approximation method. It should be noted that the estimation method includes constraints and priors and therefore the confidence intervals are different from traditionally calculated confidence intervals. Management reference points are calculated, and forward projections are carried out to provide advice for making management decisions for the yellowfin and bigeye populations. Spanish: Describimos un análisis estadístico de captura a talla estructurado por edad, A-SCALA (del inglés age-structured statistical catch-at-length analysis), basado en el modelo MULTIFAN- CL de Fournier et al. (1998). Se aplica el análisis independientemente a las poblaciones de atunes aleta amarilla y patudo del Océano Pacífico oriental (OPO). Modelamos las poblaciones de 1975 a 1999, en pasos trimestrales. Se supone solamente una sola población para cada especie para cada análisis, pero se modelan pesquerías múltiples espacialmente separadas para tomar en cuenta diferencias espaciales en la capturabilidad y selectividad. El análisis toma en cuenta error en la relación esfuerzo-mortalidad por pesca, tendencias temporales en la capturabilidad, variación temporal en el reclutamiento, relaciones entre el medio ambiente y el reclutamiento y entre el medio ambiente y la capturabilidad, y diferencias en selectividad y capturabilidad entre pesquerías. Se ajusta el modelo a datos de captura total y a datos de captura a talla proporcional condicionados sobre esfuerzo. El método A-SCALA es un enfoque estadístico, y reconoce por lo tanto que los datos obtenidos de la pesca no representan la población perfectamente. Además, hay incertidumbre en nuestros conocimientos de la dinámica del sistema e incertidumbre sobre la relación entre los datos observados y la población real. El uso de funciones de verosimilitud nos permite modelar la incertidumbre en los datos obtenidos de la población, y la inclusión de un error de proceso estimable nos permite modelar las incertidumbres en la dinámica del sistema. El enfoque estadístico permite calcular intervalos de confianza y comprobar hipótesis. Usamos una versión bayesiana del marco de verosimilitud máxima que incluye constreñimientos distribucionales sobre la variación temporal en el reclutamiento, la relación esfuerzo-mortalidad por pesca, y la capturabilidad. Se incluyen también en la función objetivo penalidades por curvatura para los parámetros de selectividad y penalidades por tasas extremas de mortalidad por pesca. Se usa la moda de la distribución posterior conjunta como estimación de los parámetros del modelo. Se calculan los intervalos de confianza usando el método de aproximación normal. Cabe destacar que el método de estimación incluye constreñimientos y distribuciones previas y por lo tanto los intervalos de confianza son diferentes de los intervalos de confianza calculados de forma tradicional. Se calculan puntos de referencia para el ordenamiento, y se realizan proyecciones a futuro para asesorar la toma de decisiones para el ordenamiento de las poblaciones de aleta amarilla y patudo.

Reservatório de gás não convencional: história de caso na Bacia do Recôncavo - Brasil

O principal objetivo deste trabalho foi à realização de um estudo geológico sobre os folhelhos fraturados da Formação Candeias que ocorrem em diferentes partes da Bacia do Recôncavo, sendo enfatizadas suas principais relações com a acumulação de petróleo. Este trabalho contempla, ainda, uma revisão bibliográfica acerca da evolução da campanha exploratória nesta bacia, desde seus primórdios, na década de trinta, além de uma análise sob as perspectivas na área de exploração e produção. A Bacia do Recôncavo, situada na região Nordeste do Brasil, está inserida numa grande bacia tafrogênica (Rifte abortado Recôncavo-Tucano-Jatobá), originada nos momentos iniciais de ruptura do paleocontinente Gondwana ocidental, no Eocretáceo, sob os campos de tensões que produziram o Oceano Atlântico. Em termos estratigráficos, a bacia guarda um importante registro geológico, compreendido por depósitos continentais (lacustres, fluvio-deltáicos, neríticos, eólicos e leques aluviais), sendo subdividido em seis seqüências sedimentares, referentes aos estágios pré-rifte e rifte. As principais acumulações de petróleo na bacia podem ser agrupadas em três grandes sistemas: pré-rifte, rifte-Candeias e rifte-Ilhas, representando plays exploratórios distintos. As perspectivas futuras para a Bacia do Recôncavo devem considerar a importância da identificação de novos plays exploratórios, através da caracterização de novos intervalos geradores e reservatórios, como é o caso dos folhelhos fraturados e do shale gas. Não menos importante é o desenvolvimento de novas técnicas de produção, aumentando o fator de recuperação de campos produtores.

Rapid quantitative detection of Yersinia pestis by lateral-flow immunoassay and up-converting phosphor technology-based biosensor

Up-converting phosphor technology (UPT)-based lateral-flow immunoassay has been developed for quantitative detection of Yersinia pestis rapidly and specifically. In this assay, 400 nm up-converting phosphor particles were used as the reporter. A sandwich immumoassay was employed by using a polyclonal antibody against F1 antigen of Y. pestis immobilized on the nitrocellulose membrane and the same antibody conjugated to the UPT particles. The signal detection of the strips was performed by the UPT-based biosensor that could provide a 980 nm IR laser to excite the phosphor particles, then collect the visible luminescence emitted by the UPT particles and finally convert it to the voltage as a signal. V-T and V-c stand for the multiplied voltage units for the test and the control line, respectively, and the ratio V-T/V-C is directly proportional to the number of Y pestis in a sample. We observed a good linearity between the ratio and log CFU/ml of Y pestis above the detection limit, which was approximately 10(4) CFU/mI. The precision of the intra- and inter-assay was below 15% (coefficient of variation, CV). Cross-reactivity with related Gram-negative enteric bacteria was not found. The UPT-LF immunoassay system presented here takes less than 30 min to perform from the sample treatment to the data analysis. The current paper includes only preliminary data concerning the biomedical aspects of the assay, but is more concentrated on the technical details of establishing a rapid manual assay using a state-of-the-art label chemistry. (c) 2006 Elsevier B.V. All rights reserved.

Variation des bio-volumes des zooplancton du plateau continental entre le Cap des Palmes et le Cap des Trois Pointes

Three years of weekly sampling from a coastal station and 29 monthly cruises over the whole continental shelf were studied for zooplankton quantitative variation. Settled volumes were preferred to displacement volumes. At the coastal station, near Abidjan, a negative correlation was found between the log2 of zooplankton volume and the preceding fortnight temperature. On the whole shelf, the differences between the 6 considered areas were tested by the variance analysis. There were significative differences in shallow waters only (20 m). During the main cold season, the upwelling of Tabou causes a very important enrichment 30 to 60 nautical miles to the east. Eastwards the plankton drifts and decreases in abundance. The zooplankton maximum is not always inshore, but often in the middle of the shelf and sometimes over the slope. During the little cold season the enrichments caused by coastal upwelling are less abundant and restricted to smaller areas. During the warm season, the waters are uniformly poor. During the cold season, over the 60m depths, the zooplankton maximum lies between 10 and 20 m and seems to sink in deeper waters. In warm season the vertical repartition is rather homogeneous in the first 40 meters. The diel vertical migrations show a very consistent rhythm, varying with the season.

On some new mathematical models for infective diseases: analysis, equilibrium, positivity and vaccination controls

196 p.

Nuevos restos viarios romanos en el norte de Burgos. ¿Vía "Flaviobriga-luliobriga"? Intento de reconstrucción del trazado por medio de sistemas de información geográficos

Este artículo forma parte de la investigación realizada en el marco de la tesis de licenciatura de su autora: "La red viaria romana en el norte de Burgos. Valles de Mena, Losa y Sotoscueva. Vía Pisoraca- Flauiobriga. Via Flauiobriga-Iuliobriga. Vías secundarias", leída en el Departamento de Estudios Clásicos de la Facultad de Filología, Geografía e Historia de la UPV/EHU el 9 de octubre de 1996.

Postrelease survival, vertical and horizontal movements, and thermal habitats of five species of pelagic sharks in the central Pacific Ocean

From 2001 to 2006, 71 pop-up satellite archival tags (PSATs) were deployed on five species of pelagic shark (blue shark [Prionace glauca]; shortfin mako [Isurus oxyrinchus]; silky shark [Carcharhinus falciformis]; oceanic whitetip shark [C. longimanus]; and bigeye thresher [Alopias superciliosus]) in the central Pacific Ocean to determine species-specific movement patterns and survival rates after release from longline fishing gear. Only a single postrelease mortality could be unequivocally documented: a male blue shark which succumbed seven days after release. Meta-analysis of published reports and the current study (n=78 reporting PSATs) indicated that the summary effect of postrelease mortality for blue sharks was 15% (95% CI, 8.5–25.1%) and suggested that catch-and-release in longline fisheries can be a viable management tool to protect parental biomass in shark populations. Pelagic sharks displayed species-specific depth and temperature ranges, although with significant individual temporal and spatial variability in vertical movement patterns, which were also punctuated by stochastic events (e.g., El Niño-Southern Oscillation). Pelagic species can be separated into three broad groups based on daytime temperature preferences by using the unweighted pair-group method with arithmetic averaging clustering on a Kolmogorov-Smirnov Dmax distance matrix: 1) epipelagic species (silky and oceanic whitetip sharks), which spent >95% of their time at temperatures within 2°C of sea surface temperature; 2) mesopelagic-I species (blue sharks and shortfin makos, which spent 95% of their time at temperatures from 9.7° to 26.9°C and from 9.4° to 25.0°C, respectively; and 3) mesopelagic-II species (bigeye threshers), which spent 95% of their time at temperatures from 6.7° to 21.2°C. Distinct thermal niche partitioning based on body size and latitude was also evident within epipelagic species.