Métodos numéricos para a integração das equações da cinética pontual de um reator nuclear

As equações da cinétiica pontual de um reator nuclear térmico são integradas numericamente, utilizando um método matricial de continuação analitica. Essas equações são essencialmente não-negativas e possuem um autovalor dominante vinculado à reatividade do sistema. Também, descrevem-se os métodos de Hansen e Porsching.

Diagramas momento-curvatura para elementos estruturais de concreto armado submetidos a carregamento monotonico ou cíclico

Este trabalho tem por objetivo desenvolver e implementar, computacionalmente, procedimentos numéricos eficientes, aplicados à determinacão do diagrama momento-curvatura, correspondentes à: -uma seção tipica, em vigas de concreto armado, submetida à carga monotônica ou cíclica de curta duração; - um ponto genérico da superficie média em placas de concreto armado, submetidas à carga monotônica de curta duração. Ainda à luz dos resultados obtidos, visa também propôr um modelo simplificado em termos de resultantes de tensões e deformações generalizadas. Inicialmente, é descrito um modelo laminar para vigas, no qual a carga é aplicada de forma incremental sendo que para cada etapa, as equaçães de equilibrio não-lineares são resolvidas de maneira iterativa. Como consequência é proposta uma relação momento-curvatura em termos de resultantes. A fim de verificar a validade e aplicabilidade dos métodos e dos algorítmos estudados e comparar-se os resultados com dados experimentais e respostas obtidas por outros pesquisadores, é apresentada uma série de exemplos numéricos. A continuação, é aplicado o procedimento anterior para modelos de laje, livres de solicitações de membrana. Finalmente através de um estudo paramétrico dos diversos fatores que afetam o diagrama momento-curvatura, propõe-se uma relação simplificada.

Inversão numérica da transformada de Laplace por polinômios trigonométricos e de Laguerre

Neste trabalho são desenvolvidos métodos numéricos para inversão da transformada de Laplace, fazendo-se uso de polinômios trigonométricos e de Laguerre. Sua utilização é ilustrada num problema de fronteira móvel da área de engenharia nuclear, através do algoritmo computacional ALG-619. Uma revisão dos aspectos analíticos básicos da transformada de Laplace e sua utilização na resolução de equações diferenciais parciais é apresentada de maneira suscinta.

Métodos recursivos para o cálculo da integral de convolução

O objetivo principal deste trabalho é apresentar um método recursivo para a determinação da resposta forçada de sistema de segunda ordem na forma de uma íntegra de concolução, proveniente da utilização de propriedades de transição da resposta impulso de tais sistemas. Descrevem-se também diversos métodos analíticos e numéricos desenvolvidos para o cálculo da resposta forçada, bem como as limitações de cada método. As vantagens do método recursivo proposto são notáveis já que não é requerido o cálculo de autovalores das matrizes nem a redução à primeira ordem, e nem o uso de hipóteses adicionais sobre natureza dos coeficientes matriciais do sistema. Como aplicação do método proposto, considera-se o cálculo da resposta dinâmica de estruturas flexíveis sujeitas a excitações arbitrárias tais como terremotos.

Simulação numérica de escoamentos hipersônicos em não-equilíbrio termo-químico através do método dos elementos finitos

o presente trabalho tem por objetivo desenvolver um código computacional destinado à simulaçáo numérica tridimensional e de escoamentos compressíveis hipersônicos em especial para analisar o processo de reentrada de veículos espaciais na atmosfera terrestre. Devido ás temperaturas elevadas encontradas neste tipo de escoamento, efeitos de altas temperaturas sáo incluídos, como excitação vibracional das moléculas e reações químícas. Utiliza-se o modelo de 2 temperaturas (translacional e vibracional);5 espécíes (N2, O2, NO, N, O) e 17 reações químicas de Park [1986 para o modelamento do fenômeno termo-químico. Introduz-se uma formulação baseada no método dos elementos finitos, através do esquema de Taylor-Calerkin, para a discretização do tempo e do espaço e para a solução do sistema de equações. Alguns exemplos e aplicações de problemas hipersônicos bi e tridimensionais, difusivos e não difusivos, são apresentados, incluindo comparações entre os resultados obtidos no presente trabalho e na literatura.

Análise elastoplástica de grelhas planas de concreto armado

O objetivo deste trabalho é apresentar um programa computacional para analisar problemas elastoplásticos com pequenas deformações em grelhas planas de concreto armado, utilizando o método dos deslocamentos. O material é definido através de um critério de escoamento generalizado que leva em consideração o efeito dos esforços de corte, momento fletor e momento torsor nas seções críticas e apresenta endurecimento isótropo. Uma solução incremental empregando o método da matriz de rigidez tangente foi adotado. Algumas comparações foram feitas com resultados experimentais, numéricos e teóricos.

Processamento e interpretação de dados aerogeofísicos magnetométricos e gamaespectrométricos no estudo do controle das mineralizações auríferas na região do Complexo Bossoroca, São Sepé, RS

Sete ocorrências auriferas estão inseridas na Seqüência Campestre do Complexo Bossoroca, no Escudo Sul-rio-grandense, municípios de São Sepé, no Estado do Rio Grande do Sul. Nessa região, foi realizado um aerolevantamento geofisico no qual foram coletados dados magnetométricos e aerogamaespectrométricos (K, U, Th e contagem total), em convênio finnado entre órgãos federais brasileiros (CPRM, DNPM e CNEN) e a empresa canadense Texas Instruments no período 1972-1973. Nessa dissertação são apresentados as interpretações e resultados dos dados e a metodologia utilizada para a interpretação das principais estruturas magnéticas e radiométricas visando aplicação na exploração mineral de ouro. Esse tipo de processamento, correlacionado aos dados geológicos, é uma prática usual de campanhas de prospecção mineral desde 1940 e é considerada como uma ferramenta de grande potencial para esse propósito. O processamento de dados magnetométricos forneceu mapas de primeira e segunda derivadas verticais, sinal analítico, gradiente horizontal, integral vertical, continuação para cima (1000, 2000 e 3000 metros) e seus principais lineamentos do campo magnético residual. Os dados foram tratados e posterionnente transfonnados em mapas de contorno e imagens pseudocoloridos e em tons de cinza com relevo sombreado (iluminação N, NE e SE) para realçar os principais lineamentos. Os mapas produzidos a partir dos dados aerogamaespectrométricos foram os de contagem total, K, U e Th, suas razões Th/K, U/K e UITh; mapas temários nos padrões RGB e CMY, potássio e urânio anômalos e parfunetro F, também foram transfonnados em mapas de contorno, pseudocolridos e tons de cinza com relevo sombreado (iluminação a N, NE e SE) e um mapa geológico interpretativo da região do Complexo Bossoroca As imagens magnetométricas mostraram-se importantes na caracterização estrutural regional, possibilitando a identificação de lineamentos NE-SW e NW-SE, não mapeados anterionnente e que podem estar relacionadas diretamente às ocorrências auriferas da região do Complexo Bossoroca. A aplicação das técnicas de processamento e interpretação de estruturas magnéticas por meio dos filtros como gradiente horizontal, primeira e segunda derivada, sinal analítico, continuação para cima, possibilitaram a determinação de descontinnidades fisicas que podem ser traduzidas em estruturas geológicas. O processamento de dados aerogamaespectrométricos possibilitou a caracterização de unidades geológicas regionais, a identificação de possíveis trends e zonas anômalas em potássio, relacionadas com as principais ocorrências auriferas da região. As imagens ternárias nos padrões RGB e CMY pennitiram a distinção e a caracterização das grandes unidades geológicas. Pode-se individualizar três regiões principais dentro da Seqüência Campestre, não mapeadas anterionnente, onde se encontram as ocorrências auriferas.

Danos ambientais causados por empresas privadas no Vale do Taquari : principais características e conseqüências jurídicas

O presente trabalho identifica e analisa danos ambientais causados por empresas privadas no Vale do Taquari, RS, Brasil. O foco principal do mesmo recai sobre os danos que foram objeto de investigação pelo Ministério Público, através de Inquéritos Civis, nas Comarcas de Estrela e Lajeado. O período examinado foi de janeiro de 1999 a dezembro de 2000. Foram constatados 73 Inquéritos Civis e nestes, foram identificados os danos causados, o agente causador, as medidas tomadas e a atual posição da investigação. Buscou-se verificar se as denúncias feitas e a instauração dos Inquéritos Civis interferem em questões relevantes em nível gerencial, como a continuação da atividade da empresa e manutenção dos empregos. Demonstra este trabalho que as principais conseqüências dos Inquéritos Civis por danos ambientais, no período avaliado, foram a reparação do dano e a adoção de medidas necessárias para a proteção ambiental.

Estudo da distribuição e dos territórios das artérias cerebrais rostral, média e caudal e cerebelares rostral, média e caudal na superfície do encéfalo em Didelphis albiventris (Gambá).

Estudaram-se as artérias cerebrais e cerebelares rostral, média e caudal na superfície do encéfalo de 30 Didelphis albiventris, 14 machos e 16 fêmeas, com o sistema arterial preenchido com Látex 603 pigmentado, e fixados em formol 20%. Sistematizaram-se as artérias e seus principais ramos, determinando a área territorial de irrigação, e as principais variações. A artéria cerebral rostral apresentou-se, à direita em 100% das peças e à esquerda esteve presente em 96,7% e ausente em 3,3% dos casos. Dividiu-se em ramo medial, tronco hemisférico rostral, e ramo lateral, artéria lateral do bulbo olfatório, presentes em 100% das amostras, bilateralmente. O tronco hemisférico rostral bifurcou-se em ramo rostral, artéria medial do bulbo olfatório, presente em todas as peças em ambos os antímeros, e ramo caudal, artéria inter-hemisférica rostral. Esta originou-se conforme o padrão em 90% dos casos à direita e em 86,7% à esquerda, nas amostras restantes teve origem de um ramo emitido pelo antímero oposto. A artéria cerebral média, presente em todos os casos nos dois antímeros, apresentou os seguintes ramos colaterais: ramos ventrorostrais presentes em 66,7% preparações; ramo ventrocaudal medial para o lobo piriforme, presente em 100% das amostras à direita e em 83,3% à esquerda; ramo ventrocaudal lateral para o lobo piriforme, presente em 100% das peças nos dois antímeros; Iº ramo hemisférico dorsolateral presente em 93,3% dos encéfalos à direita, e em 90% à esquerda; IIº ramo hemisférico dorsolateral presente em 36,6% das peças bilateralmente; ramo dorso-rostral presente à direita em 56,7% das amostras e à esquerda em 80% das preparações. A artéria cerebral média continua-se como ramo terminal em todos os casos, bilateralmente. A artéria cerebral caudal, presente bilateralmente em todas as peças, originou: ramo hemisférico tentorial em 83,3% dos casos à direita e em 80% à esquerda, este ramo originou-se da artéria comunicante caudal em 10% das peças nos dois antímeros ou ainda originou-se do tronco cerebral médio-rostral em 6,7% à direita e em 10% das preparações à esquerda; artéria tectal rostral à direita ímpar em 96,7% e dupla em 3,3% dos casos e à esquerda ímpar em 100% das peças. O ramo terminal da artéria cerebral caudal originou a artéria inter-hemisférica caudal em 93,3% das amostras à direita e em 86,7% à esquerda. A artéria cerebelar rostral originou-se como continuação natural do ramo terminal da artéria basilar em 93,3% dos casos, bilateralmente, percorrendo as faces rostral e dorsal dos hemisférios cerebelares , formando por sobre estes uma alça. Em 6,7% das amostras não formou uma alça. A artéria tectal caudal apresentou-se à direita ímpar em 83,3% e dupla em 16,7% das peças, à esquerda ímpar em 86,7%, dupla em 10% e ausente em 3,3% das amostras. A artéria cerebelar média à direita mostrou-se dupla em 70% e ímpar em 30% das peças, à esquerda, dupla em 56,7% e ímpar em 43,3% dos achados. A artéria cerebelar caudal apresentou-se à direita ímpar em 93,3% e dupla em 6,7% das peças e à esquerda ímpar em 86,7% e dupla em 13,3% das amostras.

Análise integrada da carga física de trabalho para a prevenção da fadiga

A presente dissertação tem como objetivo a avaliação da relação entre as abordagens biomecânica, fisiológica e psicofisica da carga de trabalho. Para tanto, foi avaliado,dentro do enfoque macroergonômico, o trabalho do setor de pintura de uma empresa do ramo automotivo (JOHN DEERE BRASIL) , e o trabalho das centrais de armaçao, de carpintaria e de pré-moldados de uma empresa de construção civil pesada (ODEBRECHT). Os metodos utilizados no estudo foram a análise macroergonômica do trabalho- AMT (Guimarães,1999); a avaliação biomecânica das posturas assumidas pelo método WmOWAS@(Kivi&MattiI 1991); a avaliaçao fisiológica pela mensuração das variações da freqüência cardíaca, pulso de trabalho - PT (Grandjean, 1998) e a avaliação psicofisica da carga de trabalho (Corlett, 1995). Os dados foram tabulados, analisados e comparados com limites e critérios fisiológicos deteminados na literatura. A AMT, nos dois estudos de caso, pemitiu explicar os resultados na visão dos sujeitos e dos experts. As infomações coletadas nas entrevistas e questionários da fase de apreciação fornecem a base para entender os constrangimentos biomecânicos e as alterações fisiológicas identificadas. Além disso, pemitem apontar soluções e dar validade às propostas sugeridas, inclusive, quanto a fatores organizacionais, envolvendo as questões relativas ao leiaute, ritmo e rotinas de trabalho, visto que a avaliação biomecânica concentra-se nos constrangimentos biomecânicos geralmente relacionados com posturas de trabalho, a avaliação fisiológica em fatores que impactam no esforço, tais como o calor e falta de pausas no trabalho, e a avaliação psicofisica calca na percepção do sujeito quanto à carga de trabalho.

Estudo crítico das metodologias de cálculo para perfis dobrados a frio de vigas tipo canal sem enrijecedores de borda

Neste trabalho é realizada uma análise crítica sobre as formas de dimensionamento de estruturas de chapa dobrada, apresentando as metodologias utilizadas pelas principais normas modernas de cálculo, juntamente com metodologias não convencionais. Para realizar dita análise concentrou se a atenção num tipo de perfil simples (perfis não enrijecidos), e submetido a dois tipos de carregamento: carga axial concentrada no centro de gravidade; barra trabalhando como viga em balanço com carga aplicada no centro de corte. Foi realizado um estudo comparativo, entre as normas de cálculo NBR 14762, AISI e Eurocode 3, com outros métodos de dimensionamento, que são: o método direto desenvolvido por Schafer (1997); e duas formulações de análise por elementos finitos: uma considerando a carga de flambagem elástica; e outra empregando uma análise não linear física e geométrica, para perfis submetidos à compressão centrada e à flexão. Neste trabalho realizou-se um estudo paramétrico, registrando em gráficos e tabelas a influência de parâmetros tais como: comprimento do elemento analisado; da relação entre a alma e a mesa; da relação alma e a espessura. Através destes gráficos é possível observar as diversas formas de colapso (flambagem local, flambagem global, colapso plástico entre outras) se apresentando como críticas e interagindo entre si. Realizaram-se também dois experimentos de perfis submetidos à flexão, com o objetivo de verificar a grandeza dos valores obtidos através das diversas metodologias de cálculo. Para aplicação prática dos conceitos e métodos apresentados fez-se a verificação do dimensionamento de uma estrutura utilizada como prateleira em um almoxarifado de uma empresa metalúrgica situada na grande Porto Alegre, a qual apresentou flambagem lateral de seus elementos horizontais.

A sincronização de osciladores de Rössler acoplados

Neste trabalho utiliza-se como sistema dinâmico o circuito eletrônico que integra o sistema de equações acopladas de Rossler modificado. Este sistema possui uma nãolinearidade dada por uma função linear por partes e apresenta comportamento caótico para certos valores dos seus parâmetros. Isto e evidenciado pela rota de dobramento de período obtida variando-se um dos parâmetros do sistema. A caracterização experimental da dinâmica do sistema Rossler modificado e realizada através do diagrama de bifurcações. Apresenta-se uma fundamentação teórica de sistemas dinâmicos introduzindo conceitos importantes tais como atratores estranhos, variedades invariantes e tamb em uma análise da estabilidade de comportamentos assintóticos como pontos fixos e ciclos limites. Para uma caracterização métrica do caos, apresenta-se a definção dos expoentes de Lyapunov. São introduzidos também os expoentes de Lyapunov condicionais e transversais, que estão relacionados com a teoria de sincronizção de sistemas caóticos. A partir de uma montagem mestre-escravo, onde dois osciladores de Rossler estão acoplados unidirecionalmente, introduz-se a de nição de sincronização idêntica, sincronização de fase e variedade de sincronização. Demonstra-se a possibilidade de sincronização em uma rede de osciladores caóticos de Rossler, acoplados simetricamente via acoplamento de primeiros vizinhos. A rede composta por seis osciladores mostrou ser adequada pelo fato de apresentar uma rica estrutura espacial e, ao mesmo tempo, ser experimentalmente implementável. Além da sincronização global (osciladores identicamente sincronizados), obtém-se a sincronização parcial, onde parte dos osciladores sincronizam entre si e a outra parte não o faz. Esse tipo de sincronização abre a possibilidade da formação de padrões de sincronização e, portanto, exibe uma rica estrutura de comportamentos dinâmicos. A sincronização parcial e investigada em detalhes e apresentam-se vários resultados. A principal ferramenta utilizada na análise experimental e numérica e a inspeção visual do gráfico yi yj , fazendo todas as combinações entre elementos diferentes (i e j) da rede. Na análise numérica obtém-se como resultado complementar o máximo expoente de Lyapunov transversal, que descreve a estabilidade da variedade de sincronização global.

Análise do modelo t-J e sua aplicação aos compostos de óxidos de cobre

Neste trabalho estudamos modelos teóricos que descrevem sistemas eletrônicos fortemente correlacionados, em especial o modelo t-J, e suas aplicações a compostos de óxidos de cobre, notadamente os compostos que apresentam supercondutividade de alta temperatura crítica e o composto Sr2CuO2Cl2. No primeiro capítulo do trabalho, fazemos uma exposição de três modelos que envolvem o tratamento das interações elétron-elétron, que são os modelos de Hubbard de uma banda, o modelo de Heisenberg e o modelo t-J. Na dedução deste último fazemos uma expansão canônica do hamiltoniano de Hubbard, no limite de acoplamento forte, levando-nos a obter um novo hamiltoniano que pode ser utilizado para descrever um sistema antiferromagnético bidimensional na presen- ça de lacunas, que é exatamente o que caracteriza os compostos supercondutores de alta temperatura crítica na sua fase de baixa dopagem.Após termos obtido o hamiltoniano que descreve o modelo t-J, aplicamos à este uma descrição de polarons de spin, numa representação de holons, que são férmions sem spin, e spinons, que são bósons que carregam somente os graus de liberdade de spin. Utilizando uma função de Green para descrever a propagação do polaron pela rede, obtemos uma equação para a sua autoenergia somando uma série de diagramas de Feynman, sendo que para este cálculo utilizamos a aproxima ção de Born autoconsistente[1]. Do ponto de vista numérico demonstramos que a equação integral de Dyson resultante do tratamento anterior não requer um procedimento iterativo para sua solução, e com isto conseguimos trabalhar com sistemas com grande número de partículas. Os resultados mostram, como um aspecto novo, que o tempo de vida média do holon tem um valor bastante grande no ponto (π,0 ) da rede recíproca, perto da singularidade de Van Hove mencionada na literatura[2]. Este aspecto, e suas implicações, é amplamente discutido neste capítulo. No capítulo 3 estudamos o modelo estendido t-t'-J, com tunelamento à segundos vizinhos e a incorporação dos termos de três sítios[3]. Fazemos a mesma formulação do capítulo anterior, e discutimos as aplicações dos nossos resultados ao óxido mencionado anteriormente. Finalmente, no último capítulo apresentamos uma aplicação original do modelo t-J à uma rede retangular, levemente distorcida, e demonstramos que os resultados do capítulo 3 são reproduzidos sem necessidade de introduzir termos de tunelamento adicionais no hamiltoniano. Esta aplicação pode se tornar relevante para o estudo das fases de tiras encontradas recentemente nesses materiais de óxidos de cobre.

Controle de radiação sonora numa placa retangular através de atuadores piezelétricos discretos

Neste trabalho, é estudado o controle da transmissão do som numa placa retangular e fina. Para tanto, é encontrada a resposta dinâmica da placa, excitada por forças harmônicas pontuais e piezomomentos, obtida usando uma base não-clássica e uma análise modal. A radiação sonora emitida pela vibração da placa é encontrada. A potência sonora radiada pode ser calculada aplicando controle ativo diretamente na estrutura, na forma de uma entrada vibratória, uma vez conhecida a resposta na superfície da placa, obtendo-se uma signicativa redução analitica. Os piezocerâmicos, modelados como quatro momentos pontuais, são unidos a superfície da placa como atuadores. A potência sonora transmitida antes e depois do controle é comparada, usando diferentes número de atuadores. Uma estratégia clássica de controle linear quadrático (LQR) e empregada no contexto de um procedimento de otimização da posição dos atuadores do sistema.

Otimização em dois níveis da orientação e da topologia de cascas laminadas

Este trabalho é uma contribuição para o conhecimento de metodologias de projeto de estruturas de material composto, aplicando métodos de otimização estrutural a cascas laminadas e apresentando uma estratégia em dois níveis. No primeiro nível é realizada a minimização da flexibilidade da estrutura, tendo como variável de projeto a orientação de cada lâmina da estrutura. Utiliza-se Programação Linear Seqüencial (SLP) e direção de tensão principal para otimização da orientação. No segundo nível minimiza-se o volume de cada lâmina, usando a flexibilidade total da estrutura como restrição e a densidade relativa como variável de projeto, também através de SLP. Para evitar aparecimento de áreas com densidades intermediárias, utiliza-se um Método de Continuação, dividindo o nível de otimização topológica em duas ou mais etapas. As formulações desenvolvidas permitem a solução de problemas com múltiplos casos de carregamento. Para a solução da equação de equilíbrio de casca laminada, utiliza-se um elemento finito de casca degenerado de oito nós com integração explícita na direção da espessura. A implementação desse elemento é feita de modo a facilitar a obtenção das derivadas da matriz de rigidez, necessárias na linearização das funções objetivo e restrições. Evita-se assim o uso de derivadas numéricas. Resultados para vários tipos de estrutura são apresentados, incluindo comparações entre diferentes carregamentos, condições de contorno, número de lâminas, espessuras, etc. As soluções obtidas, formas de análise e possíveis aplicações são discutidas.