Análise estática e dinâmica, linear e não-linear geométrica, através de elementos hexaédricos de oito nós com um ponto de integração

Para a análise estática e dinâmica, linear e não-linear de placas, cascas e vigas, implementa-se neste trabalho o elemento hexaédrico com integração reduzida, livre de travamento volumétrico e travamento de cisalhamento e que não apresenta modos espúrios. Na formulação do elemento, utiliza-se apenas um ponto de integração. Desta forma, a matriz de rigidez é dada de forma explícita e o tempo computacional é significativamente reduzido, especialmente em análise não-linear. Os modos espúrios são suprimidos através de um procedimento de estabilização que não exige parâmetros especificados pelo usuário. Para evitar o travamento de cisalhamento, desenvolve-se o vetor de deformações num sistema co-rotacional e remove-se certos termos não constantes nas componentes de deformações de cisalhamento. O travamento volumétrico é resolvido fazendo-se com que a parte dilatacional (esférica) da matriz gradiente seja avaliada apenas no ponto central do elemento. Como a eliminação do travamento de cisalhamento depende de uma abordagem no sistema local, emprega-se um procedimento co-rotacional para obter o incremento de deformação no sistema local e atualizar os vetores de tensões e forças internas na análise não-linear Para a solução das equações de equilíbrio na análise estática, utilizam-se métodos diretos baseados na eliminação de Gauss ou métodos iterativos de Gradientes Conjugados Precondicionado elemento-por-elemento (EBE). Para a análise dinâmica, as equações de equilíbrio são integradas através do método explícito de Taylor-Galerkin ou do método implícito de Newmark. Através de exemplos numéricos demonstra-se a eficiência e o potencial do elemento tridimensional na análise de casca, placas e vigas submetidas a grandes deslocamentos e grande rotações. Os resultados são comparados com trabalhos que utilizam elementos clássicos de placa e casca.

Defeitos pontuais nos compostos intermetálicos ZrNi e Zr/sub 2/Ni estudados por dinâmica molecular

Empregamos a técnica de Dinâmica Molecular para estudar propriedades de defeitos pontuais nos compostos intermetálicos ZrNi e Zr2Ni. Descrevemos as configurações estáveis de defeitos e mecanismos de migração, assim como as energias envolvidas. Os potenciais interatômicos foram derivados do Embedded Atom Model. No intuito de levar em conta a variação de estequiometria causada pela presença de alguns tipos de defeitos em intermetálicos, apresentamos um método numérico que fornece a energia efetiva de formação de defeitos e aplicamos o método ao ZrNi e Zr2Ni. Os resultados mostraram que vacâncias são mais estáveis na sub—rede do Ni, com energia de formação de 0,83 e-0,61 eV em ZrNi e Zr2Ni, respectivamente. Vacâncias de Zr são instáveis em ambos compostos; elas decaem espontaneamente em pares anti—sítio e vacância de Ni. Configurações e energias de formação de intersticiais também foram calculadas e mostraram comportamentos similares. Em ZrNi, a migração de vacâncias ocorre preferencialmente nas direções [025] e [100], com as respectivas energias de migração 0,67 e 0,73 eV, e é um processo essencialmente bidimensional no plano (001). Em Zr,Ni, a migração de vacâncias é unidimensional, ocorrendo na direção [001], com energia de migração de 0,67 eV. Em ambos compostos a presença de defeitos de anti—sítio de Ni diminui a energia de migração da vacância de Ni em até 3 vezes e facilita a movimentação em três dimensões. Mecanismos de anel não são energeticamente eficientes em comparação com saltos diretos. As configurações estáveis de intersticiais em ambos compostos consistem em um átomo de Ni sobre o plano (001) entre dois vizinhos de Zr fora do plano. Intersticiais de Zr são instáveis e tendem a deslocar um átomo de Ni, ocupando seu sítio. Energias de deslocamento foram estudadas através de simulações de irradiação de ambos compostos. Durante o processo de colisão binária, um potencial universal ZBL foi usado para colisões a curta distância. Para distâncias intermediárias usamos um potencial de união arbitrário. Zr mostrou—se mais difícil de ser arrancado de seu sítio do que Ni. Encontramos valores de energia de deslocamento no intervalo de aproximadamente 29 eV até 546 eV. Alguns resultados experimentais são mostrados e apresentam boa concordância com os cálculos.

Parâmetros bioecológicos de Gryon gallardoi (Brethes) (Hymenoptera: Scelionidae) e modelagem da dinâmica espaço-temporal da sua interação com Spartocera dentiventris (Berg) (Hemiptera: Coreidae) através da simulação de múltiplos agentes

Estudos referentes à dinâmica espaço-temporal da interação hospedeiro-parasitóide são de fundamental importância para o desenvolvimento de técnicas de manejo de insetos em agroecossistemas. Uma recente abordagem têm sido o estudo de populações através da simulação computacional dos indivíduos que a compõem. Este método tem possibilitado a incorporação de importantes características como individualidade e estrutura espacial nos modelos teóricos. Neste sentido, o presente estudo visou a elaboração de um modelo baseado no comportamento individual para simular a dinâmica espaço-temporal da interação Spartocera dentiventris (Berg) (Hemiptera: Coreidae), inseto associado à cultura do fumo, e de seu inimigo natural, o parasitóide de ovos Gryon gallardoi (Brethes) (Hymenoptera: Scelionidae). Para a elaboração do modelo, estudos sobre a bioecologia de G. gallardoi foram desenvolvidos, de modo a investigar (i) o efeito da temperatura no seu desenvolvimento e viabilidade, (ii) a mortalidade de imaturos em campo, (iii) os parâmetros reprodutivos e a longevidade dos adultos, (iv) os padrões de dispersão em campo e (v) a resposta funcional e a interferência mútua dos parasitóides em diferentes densidades de ovos de S. dentiventris. No intuito de avaliar a capacidade de predição do modelo elaborado, a dinâmica espaço-temporal de ovos sadios e parasitados de S. dentiventris foi acompanhada em um cultivo de fumo, estabelecido em área experimental do Departamento de Fitossanidade da Faculdade de Agronomia da UFRGS, Porto Alegre (30o01’S e 51o13’O), RS, Brasil A viabilidade do desenvolvimento dos parasitóides na faixa de 20 a 30oC não diferiu significativamente, alcançando 98,8%. O tempo de desenvolvimento ovo-adulto de machos e fêmeas foi inversamente proporcional ao aumento da temperatura. Os valores estimados para o limite térmico inferior de desenvolvimento e para a constante térmica foram 15,5oC e 185,19GD para machos e 15,6oC e 192,31GD para fêmeas. Diversos fatores afetam o sucesso de imaturos de G. gallardoi em campo, sendo que o malogro e a predação por sugadores são os principais responsáveis pela emergência de apenas 37,87% dos adultos. Foi observado um período médio de oviposição de 10,1 ± 1,74 dias, com o pico de oviposição no segundo dia, sendo depositados ao longo dos mesmos uma média de 67,5 ± 11,29 ovos por fêmea. Fêmeas de G. gallardoi foram significativamente mais longevas que os machos, vivendo, respectivamente, 13,7 ± 1,94 e 10,6 ± 1,78 dias. A razão sexual total observada foi de 0,79. A dispersão de G. gallardoi em cultivo de fumo não foi sensivelmente influenciada pelo vento, sendo estimada uma capacidade de locomoção diária média para fêmeas de no mínimo 7,62 m. O padrão de parasitismo de G. gallardoi ajustou-se perfeitamente à resposta funcional do tipo II (pseudo-r2= 0,99), sendo obtidos os valores de 0,0557 e 0,9989 h para os componentes a’ e Tm , respectivamente O aumento da densidade de parasitóides acarretou, de maneira geral, uma diminuição no número de ovos parasitados por parasitóide. Foi obtido para o índice de interferência mútua “m” o valor 0,626. A dinâmica temporal de S. dentiventris – G. gallardoi simulada ajustou-se muito bem aos dados obtidos em campo (r2= 0,82 para ovos sadios e r2= 0,72 para ovos parasitados). Da mesma forma, os arranjos espaciais dos ovos e da taxa de parasitismo observados em campo foram satisfatoriamente reproduzidos através das simulações computacionais. Os resultados indicam que G. gallardoi apresenta um bom potencial para controlar populações de S. dentiventris, seja pela sua viabilidade de desenvolvimento em uma larga faixa de temperatura ou pela sua capacidade reprodutiva e padrões de dispersão similares àqueles observados para outras espécies já utilizadas no controle biológico de insetos. Por outro lado, a sensibilidade às baixas temperaturas, a alta mortalidade de imaturos observada em campo e uma resposta funcional do tipo II podem ser fatores que se contrapõe ao sucesso do parasitóide no controle de populações de S. dentiventris. Considerando a independência entre os dados observados e os dados obtidos através de simulações, o modelo elaborado obteve pleno sucesso em simular a dinâmica da interação S. dentiventris – G. gallardoi a partir do comportamento individual dos diferentes agentes Duas importantes características associadas à estabilidade de sistemas hospedeiro-parasitóide foram observadas no estudo de campo e reproduzidas no modelo – um padrão de parasitismo inversamente dependente da densidade dos hospedeiros e uma resposta agregativa dos parasitóides em algumas regiões do sistema. Desta forma, futuros estudos podem avaliar, a partir de um enfoque baseado no indivíduo, quais mecanismos são responsáveis por tais comportamentos e qual é o seu efeito na dinâmica do sistema. Conclui-se, através deste trabalho, que o estudo da dinâmica de populações através de uma abordagem baseada no indivíduo, com base em modelos computacionais, pode ser uma alternativa viável para superar as limitações e a complexidade inerente ao uso de modelos populacionais analíticos. Desta forma, o maior realismo implícito nestes modelos os torna ferramentas muito úteis para diversas áreas da ecologia aplicada, particularmente para a entomologia agrícola.

Análise o comportamento hidrodinâmico e sedimentológico do estuario do Rio Piranji - CE (NE/Brasil)

Cerca de dois terços das maiores áreas metropolitanas mundiais estão nas proximidades dos estuários, portanto são regiões comumente sujeitas à situações de risco impostas pela pressão populacional. Instalações portuárias, efluentes diversos, uso indevido das suas margens, estão entre os riscos potenciais à grande diversidade de seus recursos naturais. A ação das marés, das ondas e o aporte fluvial interagem tornando bastante complexa a caracterização do seu funcionamento. A presente dissertação foca os seguintes aspectos: hidrodinâmica, transporte de sedimentos em suspensão, sedimentologia, morfologia e classificação do ambiente. A hidrodinâmica e o transporte de sedimentos foram estudados em duas etapas de campo em condições de maré de sizígia e de quadratura, cada qual com abrangência de dois ciclos de maré (25 hs). Estas etapas foram realizadas numa estação fixa com medições de: maré, velocidade da corrente, material em suspensão, sólidos totais dissolvidos, salinidade, pH, oxigênio dissolvido e temperatura. As características sedimentológicas foram avaliadas pela coleta de sedimentos de fundo em 24 pontos. Os aspectos morfológicos foram reconhecidos por uma batimetria detalhada e por aerofotografias. Os resultados revelaram a dominância das correntes de maré vazante em sizígia e de enchente em quadratura. O fluxo total de sedimentos por ciclo de maré na sizígia foi de 71 ton. (vazante), enquanto que na quadratura foi de 2,2 ton. (enchente). Os sedimentos de fundo variaram de areia muito grossa a silte. A profundidade máxima foi de 7,4 m (maré alta de sizígia). A comparação das aerofotografias revela a instabilidade morfológica da desembocadura. Relaciona-se este fato à influência variável da deriva litorânea e da descarga hidráulica. Este estuário foi classificado como: de frente de barreira arenosa; de circulação do tipo xvi parcialmente misturada e, por vezes, verticalmente homogênea; de Tipo 2b (quadratura); de tipo intermerdiário 2a e 2b (sizígia); de negativo; de mesomaré e de energia mista.

Estudo teórico-experimental do processo de microfiltração de partículas de sílica em suspensão aquosa

Cada vez mais muitos processos industriais convencionais de separação, tais como destilação, centrifugação e extração com solventes, entre outros, estão sendo substituídos pelos processos de separação por membranas em virtude desses apresentarem vantagens em relação aos processos convencionais, além de serem considerados uma tecnologia limpa. Entre estas vantagens pode-se citar a economia de energia, pois não envolvem mudança de fase; a possibilidade de processamento de substâncias termolábeis, pois podem operar à temperatura ambiente; a maior seletividade; a facilidade de scale – up, pois são modulares e não extensivos em mão de obra, entre outros. Uma característica limitante destes processos é a existência de uma camada de elevada concentração de partículas ou moléculas que se forma próxima à superfície da membrana e que ocasiona uma redução no fluxo permeado, podendo modificar as características de retenção da membrana. Este fenômeno é conhecido como polarização por concentração. O entendimento dos mecanismos de polarização é considerado de grande importância neste campo e o estudo e validação de modelos matemáticos preditivos torna-se necessário para se obter sucesso na aplicação dos processos com membranas no meio industrial O objetivo geral deste trabalho consiste na análise teórico-experimental do processo de microfiltração tangencial com suspensão aquosa de partículas de sílica. O objetivo específico é a utilização de dados experimentais para a estimação de um parâmetro do modelo matemático que descreve o perfil de concentração na superfície da membrana ( membrana). A configuração tangencial implica que a alimentação, contendo as partículas a serem removidas ou concentradas, flui paralelamente à superfície da membrana reduzindo, assim, o fenômeno de polarização por concentração. Os experimentos foram efetuados em uma unidade de bancada de microfiltração tangencial, também conhecido como fluxo cruzado, em um módulo para membrana plana com canal de escoamento retangular sob condições de escoamento laminar. Os seguintes parâmetros foram variados durante os experimentos: pressão transmembrana, velocidade tangencial e concentração de partículas de sílica Os dados experimentais foram utilizados em uma modelagem matemática, baseada nas equações fundamentais da mecânica dos fluidos e da transferência de massa, que considera as propriedades físicas da suspensão, viscosidade dinâmica e difusividade mássica, variáveis com a concentração de sílica para se estimar a concentração máxima de sílica junto à superfície da membrana. Esta modelagem engloba os parâmetros operacionais, bem como as características da membrana nesta única variável. Com a finalidade de se verificar como a concentração máxima sobre a membrana é afetada pelas condições operacionais e, conseqüentemente, como esta afeta o fluxo permeado, a concentração máxima sobre a superfície da membrana é apresentada em função das condições operacionais variadas durante os experimentos de microfiltração. Finalmente, após a estimação desta variável, pôde-se concluir que ela é influenciada pelas condições operacionais típicas do processo de microfiltração tangencial, tais como velocidade tangencial de escoamento, pressão através da membrana e concentração da alimentação.

Dinâmica estocástica de íons sujeitos a um conjunto quase-monocromático de ondas do tipo híbrida inferior

Neste trabalho, estudamos a interação de íons com um conjunto quase-monocromático de ondas eletrostáticas de frequência na faixa das frequências híbridas inferiores, propagando-se perpendicularmente a um campo magnético uniforme. Consideramos que as fases das ondas são aleatoriamente distribuídas (ondas incoerentes), tratando o caso de ondas de fases coerentes (ondas coerentes) como um caso particular. Derivamos o Hamiltoniano adequado a esse sistema, e deduzimos as equações de movimento, cujas soluções são analisadas numericamente, mostrando a ocorrência de difusão estocástica no espaçoo de fase ângulo-ação, para amplitudes de onda suficientemente grandes. Também fazemos estimativas sobre a amplitude mínima (threshold) para o aparecimento de ilhas de primeira ordem no espaço de fase. Estimamos, também, o limiar para as ilhas de segunda ordem e de ordens maiores, bem como o limiar de estocasticidade. A análise mostra que para o caso de várias ondas o comportamento estocástico ocorre antes do limiar de estocasticidade comparado com o caso de uma onda. No caso de ondas coerentes, observa-se que o limiar de estocasticidade diminui com o aumento do número de ondas que comp˜oem o conjunto de ondas, proporcionalmente ao inverso da raiz quadrada deste número, portanto, tendendo a ser nulo no limite em que o número de ondas no pacote tende a infinito. No caso de ondas incoerentes, observa-se também uma diminuição do limiar de estocasticidade com o aumento do número de ondas, mas nesse caso, saturando com valor até um terço do valor do limiar de estocasticidade para o caso de uma onda. Observa-se também que o limite superior da região de estocasticidade no espaço de fase aumenta com o aumento do número de ondas. No caso de ondas coerentes, esse aumento é proporcional à raiz cúbica do número de ondas que compõem o conjunto de ondas. No caso de ondas incoerentes o limite superior da região de estocasticidade têm um aumento de até o dobro em relação ao caso de uma onda. A análise também mostra que o mecanismo da estocasticidade para o caso de várias ondas é diferente do mecanismo atuante no caso de uma onda. No caso de uma onda, a estocasticidade ocorre por superposição de ilhas de ordens maiores do que um, com o aumento da intensidade da onda. No caso de várias ondas, a presençaa de ondas de frequências próximas à frequência de ressonância causa pequenas perturbações na trajetória principal das partículas, causada pela onda central, espalhando-a pelo espaço de fase de forma mais eficiente que o mecanismo de estocasticidade para o caso de uma onda.

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

Estimação e simulação dinâmica de um modelo econométrico para o mercado de alumínio no Brasil

O presente estudo pretende analisar o mercado brasileiro de alumínio primário, através da elaboração e simulação dinâmica de um modelo econométrico, com capacidade de explicar satisfatoriamente as relações típicas das variáveis econômicas consideradas neste mercado, conforme são postuladas pela teoria da empresa em competição imperfeita.

Dinâmica vegetacional e tipos funcionais em áreas excluídas e pastejadas sob diferentes condições iniciais de adubação

O entendimento da dinâmica vegetacional facilitará a tomada de decisão quanto ao tipo de manejo ou prática de melhoramento a ser adotada numa pastagem natural. Para descrição da vegetação são utilizadas espécies ou tipos funcionais. O uso de tipos funcionais promove melhor discernimento de padrões e processos que ocorrem no ecossistema. Este trabalho foi desenvolvido no Centro de Pesquisa de Forrageiras, da FEPAGRO, em São Gabriel, RS, numa área de pastagem natural, de janeiro de 1998 a dezembro de 2000, com os objetivos de estudar a dinâmica da vegetação e definir tipos funcionais para os ambientes avaliados. O delineamento experimental utilizado foi um fatorial completamente casualizado em parcela subdividida, com três repetições. Os tratamentos pastejo e exclusão foram alocados nas parcelas principais e os níveis de adubação ( com ou sem NPK) nas subparcelas. O pastejo foi realizado por bovinos, equinos e ovinos a partir de janeiro de 1998 e a adubação, anualmente, de fevereiro de 1994 a fevereiro de 1996. A exclusão iniciou em dezembro de 1996. Independente das condições iniciais, as trajetórias das áreas pastejadas, com ou sem adubação, foram convergentes, apresentando comunidades vegetais semelhantes. Nas áreas excluídas adubadas houve divergência, podendo estar relacionada à escala espacial utilizada. As espécies com maior cobertura nas áreas pastejadas são Paspalum notatum, Chevreulia acuminata, Hypoxis decumbens e Eryngium nudicaule e nas áreas excluídas do pastejo Paspalum plicatulum, Melica eremophila, Orthopappus angustifolius, Eryngium horridum e Coelorachis selloana. O pastejo é o fator determinante do tipo de comunidade vegetal resultante. Altura da planta, consistência e superfície ventral da lâmina foliar são atributos que melhor caracterizam os ambientes pastejados e excluídos.

Especificação e simulação dinâmica de um modelo macroeconômico de pequeno porte para a economia brasileira

None

Difusão sócio-geográfica do português em contato com o italiano no sul do Brasil

A Região Sul do Brasil, do ponto de vista da diversidade lingüística, caracteriza-se, entre diversos outros aspectos, pelo contato do português com as línguas dos imigrantes europeus que colonizaram a região desde o século XIX. Monolíngües no início, esses imigrantes tornaram-se bilíngües ao adquirir o português ao longo dos anos e, atualmente, a tendência é serem monolíngües em português. Em tal contexto, os italianos assumem posição de destaque, não só pelo número de falantes, mas também pelas áreas ocupadas e pela influência no contexto lingüístico, sociocultural e econômico. O português falado nas regiões em que ocorre o contato com o italiano assumiu traços específicos que refletem a constituição social e étnica dessas áreas, distinguindo-se, assim, do português falado em outras regiões e da variedadepadrão subjacente. Considerando esse cenário, o objetivo deste estudo é explicitar a dinâmica de difusão do português no espaço pluridimensional de contato com o italiano, mais especificamente em oito pontos (municípios) do Rio Grande do Sul e de Santa Catarina. A linha teórica da pesquisa segue a perspectiva da dialetologia pluridimensional e relacional, a qual busca constituir uma ciência da variação lingüística que corrija as deficiências da geolingüística tradicional e acrescente à sociolingüística uma importância maior ao valor do espaço no debate sobre a variação: “o conceito das interrelações no espaço”, segundo Harald Thun. Os dados foram coletados através de trinta e duas entrevistas, nos estilos conversa semidirigida, resposta a questionário e leitura, nas quais foram controladas dimensões sociais e geográficas, visando a verificar a pronúncia variável do ditongo nasal tônico [A)w)], do [r] forte, da vogal [a] seguida de consoante nasal, do alçamento das vogais átonas finais [e] e [o], da africação de [t] e [d] diante de [i], da realização das fricativas [S] e [Z]. Os resultados, demonstrados através de tabelas estatísticas e de mapas pluridimensionais, evidenciam que a difusão dos traços associados ao português varia no modo e na intensidade. No plano diatópico, ocorre difusão mais intensa em Orleans (SC) e Caxias do Sul (RS), ao passo que a maior resistência à inovação lingüística foi detectada em Rodeio (SC) e Sananduva (RS). Na perspectiva diassocial, o uso de variantes sem interferência do italiano é liderado, sucessivamente, pelos falantes urbanos, pelos mais jovens e mais escolarizados.

Dinâmica não-linear e sustentabilidade da dívida pública brasileira

This paper contributes to the debate on whether the Brazilian public debt is sustainable or not in the long run by considering threshold effects on the Brazilian Budget Deficit. Using data from 1947 to 1999 and a threshold autoregressive model, we find evidence of delays in fiscal stabilization. As suggested in Alesina (1991), delayed stabilizations reflect the existence of political constraints blocking deficit cuts, which are relaxed only when the budget deficit reaches a sufficiently high level, deemed to be unsustainable. In particular, our results suggest that, in the absence of seignorage, only when the increase in the budget deficit reaches 1.74% of the GDP will fiscal authorities intervene to reduce the deficit. If seignorage is allowed, the threshold increases to 2.2%, suggesting that seignorage makes government more tolerant to fiscal imbalances.