Aplicação da teoria de dano na análise do comportamento viscoelástico de estruturas de concreto

A presente dissertação de mestrado tem por assunto a representação do comportamento mecânico do concreto sob cargas de curta e longa duração, incluindo efeitos não-lineares. Para tal fim trabalha-se com equações baseadas na teoria do dano contínuo. São propostas equações para o caso triaxial e, baseado nelas, é implementado um programa computacional. Com diversos exemplos verifica-se que: a) A solução numérica aproxima bem os resultados teóricos. b) O comportamento do modelo representa bem as características qualitativas do concreto. c) O modelo permite aproximar bem alguns resultados experimentais, mas ainda deve ser aperfeiçoado, particularmente no que refere-se à identificação de parâmetros.

Solução da equação de transporte multidimensional em geometria cartesiana e meio infinito usando derivada fracionária

Neste trabalho, foi construída uma forma integral para a solução das equações de transporte em uma, duas e três dimensões, considerando o núcleo de espalhamento de Klein-Nishina, espalhamento isotrópico e o núcleo de espalhamento de Rutherford, respectivamente, seguindo a mesma idéia proposta em trabalhos recentes, nos quais foi construída uma solução para a equação de transporte de nêutrons em geometria cartesiana, usando derivada fracionária. A metodologia consiste em igualar a derivada fracionária do fluxo angular à equação integral, determinar a ordem da derivada fracionária comparando o núcleo da equação integral com o da definição de Riemann-Liouville. Essa formulação foi aplicada ao cálculo de dose absorvida. São apresentadas soluções geradas a partir do emprego do método da derivada fracionária e comparadas a resultados disponíveis na literatura.

Cálculo da complexidade exata de algoritmos do tipo divisão-e-conquista através das equações características

A equação de complexidade de um algoritmo pode ser expressa em termos de uma equação de recorrência. A partir destas equações obtém-se uma expressão assintótica para a complexidade, provada por indução. Neste trabalho, propõem-se um esquema de solução de equações de recorrência usando equações características que são resolvidas através de um "software" de computação simbólica, resultando em uma expressão algébrica exata para a complexidade. O objetivo é obter uma forma geral de calcular a complexidade de um algoritmo desenvolvido pelo método Divisão-e-Conquista.

Estudo e solução da equação de transporte de nêutrons bidimensionais pelo método LTSn para elevadas ordens de quadraturas angulares : LTSn2D - Diag e LTSn2D - DiagExp

O principal objetivo dessa tese consiste em determinar uma solução numéricada equação bidimensional do transporte de nêutrons para elevadas ordens de quadratura angular. Diagonalizando a matriz de transporte LTSN bidimensional , construímos dois algoritmos que se diferenciam pela forma de representar os termos de fuga transversal, que surgem nas equações LTSN integradas transversalmente. Esses termos no método LTSN2D − Diag são expressos como combinação linear dos autovetores multiplicados por exponenciais dos respectivos autovalores. No método LTSN2D − DiagExp os termos de fuga transversal são representados por uma função exponencial com constante de decaimento heuristicamente identificada com parâmetros materiais característicos do meio. A análise epectral desenvolvida permite realizar a diagonalização. Um estudo sobre o condicionamento é feito e também associamos um número de condicionamento ao termo de fuga transversal. Definimos os erros no fluxo aproximado e na fórmula da quadratura, e estabelecemos uma relação entre eles. A convergência ocorre com condições de fronteira e quadratura angular adequadas. Apresentamos os resultados numéricos gerados pelos novos métodos LTSN2D − Diag e LTSN2D − DiagExp para elevadas ordens de quadratura angular para um problema ilustrativo e comparamos com resultados disponíveis na literatura.

Uma metodologia unificada no domínio tempo para sistemas concentrados, discretos e distribuídos

A resposta impulso é utilizada como ferramenta padrão no estudo direto de sistemas concentrados, discretos e distribuídos de ordem arbitrária. Esta abordagem leva ao desenvolvimento de uma plataforma unificada para a obtenção de respostas dinâmicas. Em particular, as respostas forçadas dos sistemas são decompostas na soma de uma resposta permanente e de uma resposta livre induzida pelos valores iniciais da resposta permanente. A teoria desenvolve-se de maneira geral e direta para sistemas de n-ésima ordem, introduzindo-se a base dinâmica gerada pela resposta impulso na forma padrão e normalizada, sem utilizar-se a formulação de estado, através da qual reduz-se um sistema de ordem superior para um sistema de primeira ordem. Considerou-se sistemas de primeira ordem a fim de acompanhar-se os muitos resultados apresentados na literatura através da formulação de espaço de estado. Os métodos para o cálculo da resposta impulso foram classificados em espectrais, não espectrais e numéricos. A ênfase é dada aos métodos não espectrais, pois a resposta impulso admite uma fórmula fechada que requer o uso de três equações características do tipo algébrica, diferencial e em diferenças Realizou-se simulações numéricas onde foram apresentados modelos vibratórios clássicos e não clássicos. Os sistemas considerados foram sistemas do tipo concentrado, discreto e distribuído. Os resultados da decomposição da resposta dinâmica de sistemas concentrados diante de cargas harmônicas e não harmônicas foram apresentados em detalhe. A decomposição para o caso discreto foi desenvolvida utilizando-se os esquemas de integração numérica de Adams-Basforth, Strömer e Numerov. Para sistemas distribuídos, foi considerado o modelo de Euler-Bernoulli com força axial, sujeito a entradas oscilatórias com amplitude triangular, pulso e harmônica. As soluções permanentes foram calculadas com o uso da função de Green espacial. A resposta impulso foi aproximada com o uso do método espectral.

Desenvolvimento de um método de volumes finitos de alta ordem para a simulação de escoamentos de fluidos viscoelásticos

O uso da mecânica de fluidos computacional no estudo de processos envolvendo o escoamento de fluidos poliméricos está cada vez mais presente nas indústrias de transformação de polímeros. Um código computacional voltado a esta função, para que possa ser aplicado com sucesso, deve levar a predições mais próximas possível da realidade (modelagem), de uma forma relativamente rápida e eficiente (simulação). Em relação à etapa de modelagem, o ponto chave é a seleção de uma equação constitutiva que represente bem as características reológicas do fluido, dentre as diversas opções existentes. Para a etapa de simulação, ou seja, a resolução numérica das equações do modelo, existem diversas metodologias encontradas na literatura, cada qual com suas vantagens e desvantagens. Neste tópico se enquadra o trabalho em questão, que propõe uma nova metodologia para a resolução das equações governantes do escoamento de fluidos viscoelásticos. Esta se baseia no método dos volumes finitos, usando o arranjo co-localizado para as variáveis do problema, e na utilização de aproximações de alta ordem para os fluxos médios lineares e não-lineares e para outros termos não lineares que surgem da discretização das equações constitutivas. Nesta metodologia, trabalha-se com os valores médios das variáveis nos volumes durante todo o processo de resolução, sendo que os valores pontuais são obtidos ao final do procedimento via deconvolução. A solução do sistema de equações não lineares, resultante da discretização das equações, é feita de forma simultânea, usando o método de Newton São mostrados então, resultados da aplicação da metodologia proposta em problemas envolvendo escoamentos de fluidos newtonianos e fluidos viscoelásticos. Para descrever o comportamento reológico destes últimos, são usadas duas equações constitutivas, que são o modelo de Oldroyd-B e o modelo de Phan-Thien-Tanner Simplificado. Por estes resultados pode-se ver que a metodologia é muito promissora, apresentando algumas vantagens frente às metodologias convencionais em volumes finitos. A implementação atual da metodologia desenvolvida está restrita a malhas uniformes e, consequentemente, soluções para problemas com geometrias complexas, que necessitam de refinamento localizado da malha, foram obtidas somente para baixos números de Weissenberg, devido a limitação do custo computacional. Esta restrição pode ser contornada, tornando o seu uso competitivo.

Paralelizações de métodos numéricos em clusters empregando as bibliotecas MPICH, DECK e Pthread

Este trabalho tem como objetivo desenvolver e empregar técnicas e estruturas de dados agrupadas visando paralelizar os métodos do subespaço de Krylov, fazendo-se uso de diversas ferramentas e abordagens. A partir dos resultados é feita uma análise comparativa de desemvpenho destas ferramentas e abordagens. As paralelizações aqui desenvolvidas foram projetadas para serem executadas em um arquitetura formada por um agregado de máquinas indepentes e multiprocessadas (Cluster), ou seja , são considerados o paralelismo e intra-nodos. Para auxiliar a programação paralela em clusters foram, e estão sendo, desenvolvidas diferentes ferramentas (bibliotecas) que visam a exploração dos dois níveis de paralelismo existentes neste tipo de arquitetura. Neste trabalho emprega-se diferentes bibliotecas de troca de mensagens e de criação de threads para a exploração do paralelismo inter-nodos e intra-nodos. As bibliotecas adotadas são o DECK e o MPICH e a Pthread. Um dos itens a serem analisados nestes trabalho é acomparação do desempenho obtido com essas bibliotecas.O outro item é a análise da influência no desemepnho quando quando tulizadas múltiplas threads no paralelismo em clusters multiprocessados. Os métodos paralelizados nesse trabalho são o Gradiente Conjugação (GC) e o Resíduo Mínmo Generalizado (GMRES), quepodem ser adotados, respectivamente, para solução de sistemas de equações lineares sintéticos positivos e definidos e não simétricas. Tais sistemas surgem da discretização, por exemplo, dos modelos da hidrodinâmica e do transporte de massa que estão sendo desenvolvidos no GMCPAD. A utilização desses métodos é justificada pelo fato de serem métodos iterativos, o que os torna adequados à solução de sistemas de equações esparsas e de grande porte. Na solução desses sistemas através desses métodos iterativos paralelizados faz-se necessário o particionamento do domínio do problema, o qual deve ser feito visando um bom balanceamento de carga e minimização das fronteiras entre os sub-domínios. A estrutura de dados desenvolvida para os métodos paralelizados nesse trabalho permite que eles sejam adotados para solução de sistemas de equações gerados a partir de qualquer tipo de particionamento, pois o formato de armazenamento de dados adotado supre qualquer tipo de dependência de dados. Além disso, nesse trabalho são adotadas duas estratégias de ordenação para as comunicações, estratégias essas que podem ser importantes quando se considera a portabilidade das paralelizações para máquinas interligadas por redes de interconexão com buffer de tamanho insuficiente para evitar a ocorrência de dealock. Os resultados obtidos nessa dissertação contribuem nos trabalhos do GMCPAD, pois as paralelizações são adotadas em aplicações que estão sendo desenvolvidas no grupo.

Estimativas do erro nas aproximações de Galerkin para as equações de Navier-Stokes

Neste trabalho são provadas algumas estimativas de erro em espaços para as aproximações de Galerkin para a solução do sistema de equações de Navier-Stokes. Mostra-se que o erro decresce em proporção inversa aos autovalores do operador de Stokes.

Formulação semi-analítica para a equação transformada resultante da aplicação da GITT em problemas difusivos-advectivos

Neste trabalho se propõe um avanço para a Técnica Transformada Integral Generalizada, GITT. O problema transformado, usualmente resolvido por subrotinas numéricas, é aqui abordado analiticamente fazendo-se uso da Transformada de Laplace. Para exemplificar o uso associado destas duas transformadas integrais, resolvem-se dois problemas. Um de concentração de poluentes na atmosfera e outro de convecção forçada com escoamento laminar, entre placas planas paralelas, com desenvolvimento simultâneo dos perfis térmico e hidrodinâmico. O primeiro é difusivo, transiente e com coeficientes variáveis. Sua solução é obtida de forma totalmente analítica. Além de mostrar o uso da técnica, este exemplo apesar de ter coeficientes variáveis, é resolvido com o auxílio de um problema de autovalores associado com coeficientes constantes. No segundo, obtém-se a solução da Equação da Energia analiticamente. Já a Equação da Conservação do Momentum é linearizada e resolvida de forma iterativa. A solução de cada iteração é obtida analiticamente.

Respostas periódicas em sistemas lineares e fracamente não lineares não ressonantes e comportamento dinâmico de sistemas rotativos com o uso da base dinâmica

Este trabalho visa realizar o estudo do comportamento dinâmico de um eixo rotor flexível, modelado segundo a teoria de Euler-Bernoulli e caracterizar as respostas periódicas de sistemas LTI (sistemas lineares invariantes no tempo) e sistemas fracamente não lineares de ordem arbitrária. Para tanto, é utilizada a base dinâmica gerada pela resposta impulso ou solução fundamental. O comportamento dinâmico de um eixo rotor flexível foi discutido em termos da função de Green espacial e calculada de maneira não-modal. Foi realizado um estudo do problema de autovalor para o caso de um um eixo rotor biapoiado. As freqüências são obtidas e os modos escritos em termos da base dinâmica e da velocidade de rotação. As respostas periódicas de sistemas LTI, utilizadas nas aproximações com sistemas fracamente não lineares, são obtidas, independentemente da ordem do sistema, como um operador integral onde o núcleo é a função de Green T-periódica. Esta função é caracterizada em termos das propriedades de continuidade, periodicidade e salto da função de Green T-periódica, e da base dinâmica Simulações foram realizadas para sistemas concentrados, matriciais e escalares, com o objetivo de mostrar a validade da metodologia desenvolvida com as propriedades da função de Green T-periódica. Foi abordado um modelo não-linear para uma centrífuga utilizada na indústria textil [Starzinski, 1977].

Diferenciais de salários por raça e gênero: aplicação dos procedimentos de Oaxaca e Heckman em pesquisas amostrais complexas

Este artigo decompõe o diferencial de salários por cor e sexo dos trabalhadores brasileiros usando os microdados da Pesquisa Nacional por Amostra de Domicílio (PNAD). A metodologia consiste em estimar a equação de salários (Mincer, 1974) com a correção do viés de seleção das informações dos salários (Heckman, 1979). Em seguida, a decomposição do diferencial da média do logaritmo do salário/hora foi obtida pelo procedimento de Oaxaca (1973) apresentada em dois efeitos: características produtivas e discriminação. A análise empírica tem como foco o uso adequado de procedimentos de modelagem estatística em pesquisas, por amostragem complexa, conforme os trabalhos de Skinner e Smith (1989) e Pessoa e Silva (1998). Os resultados indicam a necessidade de se incorporar o plano amostral e a correção do viés de seleção da informação dos salários, visando melhorar a qualidade das estimativas das equações de salários e avaliar adequadamente as medidas de discriminação. Como exemplo, a estimativa do coeficiente de discriminação, D, entre homens e mulheres de cor branca é 0,37 sem a correção do viés e 0,30 com a correção do viés de seleção das informações dos salários.

Equações de advecção-difusão: propriedades e comportamento assintótico

Neste trabalho, examinamos em detalhe resultados recentes apresentados em [Zingano, 1999], [Zingano, 2004], [Zingano, 1996a] [T. Hagstrom, 2004] sobre o comportamento de soluções para equações (escalares) de ad vecção-difusão nãolineares, da forma Ut + div(f(u)) = div(A(u)V'u), x E ]Rn, t > O correspondentes a estados iniciais u(., O) E LI(]Rn) n DXJ(JRn).Aqui, A(u) E ]Rn é uniformemente positiva definida para todos os valores de u em questão, e f( u) = (f1(u),..., fn(u)) corresponde ao fluxo advectivo, com A, f suaves. Entre os vários resultados, tem-se em particular os limites assintóticos . !!. (I_l) Iml (47rÀ)~ 11mt2 p Ilu(" t)IILP(JRn) = (4 À)!!. - , t-++oo 7r 2 P para cada 1 :::;P :::;00, uniformemente em p, bem como lim t~(l-i) Ilu(" t) - u(',t)IILP(JRn) = O, t-++oo 1:::; p:::; 00 para duas soluçõesu(', t), u(', t) quaisquer correspondentesa estados iniciais u(', O),u(', O)E LI (]Rn) n Loo(]Rn) com a mesma massa, isto é, r u(x, O)dx = r u(x,O)dx JJRn JJRn Outra propriedade fundamental, válida em dimensão n ;:::2, é lim t%(l-~) Ilu(" t) - v(', t) IILP(JRn) = O t-++oo para cada 1 :::;p :::; 00, se v(', t) é solução da equação de advecção-difusão linear Vt + f (O) . V'v= div(A(O)V'v), x E ]Rn, t > O, com u(', O),v(', O) E U(]Rn) n Loo(JRn) tendo a mesma massa. Outros resultados de interesse são também discutidos.

Existência de soluções para equações elípticas semilineares com crescimento singular crítico

Neste trabalho estudamos uma equação diferencial parcial elíptica semilinear contendo uma singularidade e um termo de crescimento crítico. A existência de soluções depende da dimensão do espaço e do coeficiente da singularidade. Através da caracterização variacional e com o uso de seqüências de Palais-Smale provamos que o problema possui soluções não triviais.

Metodologia Samário/Neodímio : uma abordagem analítica simplificada e alguns exemplos aplicativos

O presente trabalho discorre sobre o desenvolvimento de uma nova técnica para análise isotópica de Sm-Nd com uso de um traçador combinado de Nd145/Sm149 e de somente a coluna primária, que envolve apenas razões que não requerem correção por interferência isobárica. O trabalho divide-se em duas partes: (i) na primeira parte está apresentado um breve histórico da aplicação da metodologia Sm/Nd e a sua prática rotineira nos laboratórios de geocronologia; (ii) na segunda parte, a nova metodologia Sm/Nd proposta envolvendo calibrações de um traçador combinado de Nd145 e Sm149 e os resultados obtidos nas análises de 3 amostras pré-selecionadas, incluindo uma amostra de referência (BHVO-2) e um teste da validade da técnica de se traçar uma amostra em alíquota. A razão Sm147/Sm149 e as razões de Nd (Nd143/Nd145, Nd143/Nd146 e Nd145/Nd146) nas amostras foram analisadas simultaneamente utilizando um arranjo triplo de filamento (Ta, Re, Ta). A presença de outros ETR não interferiu na eficiência de ionização, sendo as emissões do Sm e do Nd praticamente simultâneas. Esse fato implica que para corrigir o fracionamento observado na razão Sm147/Sm149 pode ser usada a razão Nd143/Nd146 = 0,7101232. A correção de fracionamento, como é empregado no método convencional, reduz os erros nos cálculos das concentrações de Sm e Nd. A razão Nd143/Nd144 normalizada, necessária para idades modelo ou isócrona, é derivada da multiplicação da razão Nd143/Nd145 obtida na amostra natural normalizada para Nd145/Nd146 = 0,482620 por um fator constante de 0,348404. Esse valor de normalização estabelecido é baseado em 45 análises independentes em duas amostras de referência laboratorial (SPEX e Johnson Mattey) obtidas no modo dinâmico de análise utilizando um espectrômetro de massa da VG Sector equipado com 9 coletores A razão Nd143/Nd144 calculada e a idade modelo TDM são plenamente concordantes quando comparadas aos valores obtidos pelo método convencional. Algumas restrições podem surgir quando os valores de Sm e Nd obtidos nas alíquotas forem comparados com os resultados obtidos da amostra total. Pequenas discrepâncias, provavelmente atribuíveis à não homogeneização da solução final, podem ser irrelevantes, uma vez que as razões Sm147/Nd144 e Nd143/N144d derivada estão dentro dos erros experimentais (2s). Para laboratórios de geocronologia com um espectrômetro de massa equipado com pelo menos 5 coletores e com o método Rb/Sr implementado, a nova metodologia proposta para o método Sm/Nd constitui-se em uma alternativa de aplicação quase imediata, desde que a solução de traçador combinado Nd145 e Sm149 esteja bem calibrada e a coluna de resina catiônica tenha resolução adequada para eficiente separação dos elementos terras raras.

Sobre convexidade de conjuntos de nível de soluções de certas equações elípticas

Neste trabalho é estudada a convexidade dos conjuntos de nível das soluções de dois problemas envolvendo equações elípticas. O primeiro desses problemas se refere a uma equação da forma 4u = °(u) em um anel convexo, com condições de fronteira u = 0 na fronteira externa e u = 1 na fronteira interna. Para provar a existência de solução do problema utiliza-se o método variacional. O problema de mostrar a convexidade dos conjuntos de nível é transformado em um problema de maximizar uma certa função. O segundo problema considerado é o de mostrar que é log-côncava a primeira autofunção do laplaciano, que tenha como peso uma função côncava.