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 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.

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.

Características do sinal mecanomiográfico em atletas velocistas, fundistas e indivíduos sedentários

A mecanomiografia (MMG) foi utilizada para estudar o comportamento mecânico e fisiológico do músculo vasto lateral de atletas velocistas, de atletas fundistas e de indivíduos sedentários. Partindo-se do pressuposto que o músculo acima apresenta diferentes composições de fibras musculares nos atletas velocistas e nos atletas fundistas, esperava-se que o sinal MMG produzido durante a contração fosse diferente entre esses dois grupos. A amostra foi constituída por 30 sujeitos (10 atletas velocistas, 10 atletas fundistas e 10 indivíduos sedentários), do sexo masculino (18 a 30 anos de idade) sem história de lesão ou doença neuromuscular. Os sujeitos foram submetidos a um teste de esforço voluntário e a um teste de contrações produzidas por meio de estimulação elétrica artificial. Paralelamente aos sinais MMG, foram também coletados os sinais eletromiográficos (EMG), através dos quais se verificou um crescimento na ativação elétrica do músculo vasto lateral dos três grupos da amostra, à medida que aumentou o nível de esforço voluntário. Os sinais MMG obtidos não apresentaram diferenças significativas entre os três grupos da amostra, tanto na sua magnitude quanto em seu conteúdo de freqüência, na maior parte dos testes realizados. Esses resultados sugerem, ao contrário das idéias de alguns autores, que a MMG não é uma técnica que possibilite a fácil detecção do comportamento mecânico e fisiológico de músculos com diferentes percentuais de unidades motoras (UMs).Entretanto, verificou-se uma tendência no sentido de que as respostas dos sinais MMG do músculo vasto lateral dos atletas fundistas são menores que aquelas dos velocistas e indivíduos sedentários.

Analgesia em pós-operatório de toratomias : estudo prospectivo e randomizado

Introdução: A dor é um importante fator de incremento da morbidade e mortalidade em pacientes submetidos a procedimentos cirúrgicos que incluem toracotomias. Diversos fatores contribuem para que esses pacientes apresentem um alto grau de dor no pós-operatório, entre os quais a secção da pele, músculos e pleura, retração dos músculos e ligamentos pelo afastador de Finochietto, irritação da pleura e nervos intercostais pelos drenos tubulares torácicos e fraturas ocasionais dos arcos costais. O aumento das taxas de morbidade e mortalidade é dado principalmente à respiração superficial decorrente da pouca mobilidade da parede torácica e conseqüente à dor e pela perda da efetividade do principal mecanismo de eliminação de secreções da árvore traqueobrônquica (tosse), resultando em atelectasias, inadequado gradiente ventilação / perfusão, hipoxemia e pneumonia. Uma vez caracterizada a necessidade de atenuação da dor como fator primordial na melhora dos índices de morbidade e mortalidade no período pós-operatório de cirurgia torácica, torna-se imperiosa uma análise das terapêuticas disponíveis na atualidade para tanto. Objetivos: Avaliar a utilização de três diferentes métodos de analgesia: 1. bloqueio peridural com morfina (BPM); 2. morfina parenteral (MP); e 3. bloqueio intercostal extrapleural contínuo com lidocaína” (BIC), em pacientes submetidos a procedimentos que incluíram toracotomias em sua execução, além de analisar o custo financeiro desses métodos. Materiais e métodos: Trata-se de um estudo prospectivo, randomizado, no qual foram analisados 79 pacientes, submetidos a toracotomias, subdivididos de forma aleatória em três grupos, de acordo com a modalidade terapêutica instituída: 25 pacientes no grupo BIC, 29 pacientes no grupo BPM e 25 pacientes no grupo MP. Cada paciente foi observado e analisado por profissionais de enfermagem previamente treinados. As variáveis analisadas foram a dor e a sedação. (quantificadas através de escores e analisadas através do método de Kruskal-Wallis com correção pelo teste de Dunn), além do custo financeiro de cada método e da necessidade de administração de opióides adicionais. Resultados: As variáveis dor e sedação foram obtidas através das seguintes medianas, respectivamente: grupos BIC (2,5 e 0); BPM (4 e 0) e MP (3,5 e 0). O custo financeiro foi de US$ 78,69 para o grupo BIC; US$ 28,61 para o grupo BPM e US$ 11,98 para o grupo MP. A necessidade adicional de opióide foi de 4,2 mg/dia para o grupo BIC; 5,7 mg/dia para o grupo BPM e 10,7 mg/dia para o grupo MP. Conclusões: A intensidade da dor foi significativamente menor no grupo BIC, quando comparado ao grupo MP. Não foram identificadas diferenças significativas de intensidade da dor quando comparados os grupos BIC versus BPM e BPM versus MP. A intensidade de sedação foi significativamente maior no grupo MP quando comparado aos grupos BIC e BPM. Não foram evidenciadas diferenças significativas quanto à sedação entre os grupos BIC e BPM. O custo financeiro do grupo MP foi sensivelmente menor quando comparado aos grupos BIC e BPM. A necessidade adicional de morfina foi significativamente maior no grupo MP, quando comparados aos grupos BIC e BPM.

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.

Um estudo da fatoração incompleta LU e Cholesky como pré-condicionadores nos métodos iterativos

Neste trabalho procuramos analisar alguns métodos iterativos e os processos de aceleração na solução lineares grandes e esparsos, associando o uso de alguns pré-condicionadores, tais como os métodos de fatoração incompleta. De forma mais específica, nos detivemos no estudo deos métodos de fatoração incompleta LU, ou ILU, e o método de Cholesky incompleto. Para isso procuramos antes definir algumas especificidades sobre esses métodos, tais como, crtérios de existência, limitação. Alguns fatores analisam tais problemas e sugerem algumas técnicas de conserto, ou seja, algumas maneiras de eliminar tais falhas para que os métodos de iteração possam ser utlizados para determinar soluções mais próximas da solução real. Procedemos a uma revisão teórica de alguns dos métodos iterativos, dos pré-condicionadores. Jacobi, fatoração incompleta LU e fatoração incompleta de Cholesky e a sua associação com os métodos iterativos GMRES e Gradiente Conjugado. Utilizando os pré-condionadores associados aos métodos iterativos citados e fixando alguns parâmetros de parada, aplicamos algusn testes. Os resultados e a análise dos mesmos encontram-se neste trabalho.

Algoritmos adaptativos para o método GMRES(m)

Nesse trabalho apresentamos algoritmos adaptativos do M´etodo do Res´ıduo M´ınimo Generalizado (GMRES) [Saad e Schultz, 1986], um m´etodo iterativo para resolver sistemas de equa¸c˜oes lineares com matrizes n˜ao sim´etricas e esparsas, o qual baseia-se nos m´etodos de proje¸c˜ao ortogonal sobre um subespa¸co de Krylov. O GMRES apresenta uma vers˜ao reinicializada, denotada por GMRES(m), tamb´em proposta por [Saad e Schultz, 1986], com o intuito de permitir a utiliza¸c˜ao do m´etodo para resolver grandes sistemas de n equa¸c˜oes, sendo n a dimens˜ao da matriz dos coeficientes do sistema, j´a que a vers˜ao n˜ao-reinicializada (“Full-GMRES”) apresenta um gasto de mem´oria proporcional a n2 e de n´umero de opera¸c˜oes de ponto-flutuante proporcional a n3, no pior caso. No entanto, escolher um valor apropriado para m ´e dif´ıcil, sendo m a dimens˜ao da base do subespa¸co de Krylov, visto que dependendo do valor do m podemos obter a estagna¸c˜ao ou uma r´apida convergˆencia. Dessa forma, nesse trabalho, acrescentamos ao GMRES(m) e algumas de suas variantes um crit´erio que tem por objetivo escolher, adequadamente, a dimens˜ao, m da base do subespa¸co de Krylov para o problema o qual deseja-se resolver, visando assim uma mais r´apida, e poss´ıvel, convergˆencia. Aproximadamente duas centenas de experimentos foram realizados utilizando as matrizes da Cole¸c˜ao Harwell-Boeing [MCSD/ITL/NIST, 2003], que foram utilizados para mostrar o comportamento dos algoritmos adaptativos. Foram obtidos resultados muito bons; isso poder´a ser constatado atrav´es da an´alise das tabelas e tamb´em da observa ¸c˜ao dos gr´aficos expostos ao longo desse trabalho.