67 resultados para Modelos de elementos finitos

em Lume - Repositório Digital da Universidade Federal do Rio Grande do Sul


Relevância:

100.00% 100.00%

Publicador:

Resumo:

o objetivo desta tese apresentar dois modelos distintos para representar as fissuras em peas de concreto armado: um do tipo distribudo e outro do tipo incorporado. Os modelos de fissura incorporada se baseiam no conceito de descontinuidades incorporadas dentro do campo de deslocamento do elemento finito padro. J nos modelos de fissura distribuda a descontinuidade do campo de deslocamentos causada pela fissura espalhada ao longo do elemento. O modelo distribudo proposto apresenta um diagrama tenso-deformao multilinear e sua calibragem feita, segundo o Cdigo Modelo CEB-FIP 1990, atravs de ensaios de tirantes de concreto armado. O modelo incorporado implementado baseado no modelo de Dvorlcin, CuiMo e Gioia, o qual no inclui a contribuio da armadura no equilbrio interno de foras do elemento. A incluso da parcela da armadura feita atravs de um modelo de transferncia de tenso por aderncia, conforme Russo, Zingone e Romano, Russo e Romano e FIE - Bulletin 10. Para representar o comportamento do concreto intacto, utiliza-se o modelo constitutivo de Ottosen. Trata-se de um modelo elstico no-linear, tridimensional, que utiliza valores secantes dos parmetros do material. J para simular o comportamento das barras de ao da armadura, emprega-se o modelo incorporado desenvolvido por Elwi e Hrudey. Neste modelo, permite-se uma disposio arbitrria das barras de ao no interior dos elementos de concreto. O modelo constitutivo adotado para a armadura do tipo elasto-plstico com endurecimento. Por fim, alguns exemplos numricos so analisados com o objetivo de comprovar a eficcia dos dois modelos propostos.

Relevância:

100.00% 100.00%

Publicador:

Resumo:

O objetivo deste trabalho desenvolver um modelo computacional, baseado no mtodo dos elementos finitos, para o estudo de peas de concreto armado e protendido submetidas a estados planos de tenso. O estudo abrange situaes de carga de curta e longa durao, onde consideram-se fluncia e retrao do concreto e relaxao do ao. So utilizados modelos constitutivos elasto-viscoplsticos para descrever o comportamento dos materiais. Implementou-se um modelo de camadas superpostas para melhor representar o comportamento do concreto, onde o material composto de diversas camadas que sofrem a mesma deformao. Cada camada possui diferentes caractersticas materiais e a tenso total obtida pela soma das diferentes contribuies de cada camada. Para a fissurao da concreto, utilizou-se um modelo de fissuras distribudas, que leva em conta a contribuio do concreto entre fissuras. Tanto a amadura passiva como a de pratenso so introduzidas no modelo como uma linha de material mais rgido dentro do elemento de concreto. Os deslocamentos ao longo da armadura so referenciados aos deslocamentos nodais do elemento de concreto. Deste modo, obtm-se uma matriz de rigidez para a armadura com as mesmas dimenses que a matriz de rigidez do elemento de concreto, A matriz de rigidez do elemento concreto-ao a soma das duas matrizes. Considera-se aderncia perfeita entre o concreto e o ao. Os resultados obtidos com esse programa computacionai so comparados com valores experimentais disponveis.

Relevância:

100.00% 100.00%

Publicador:

Resumo:

A corroso da annadura um dos principais problemas para um grande nmero de es1ruturas de concreto annado por causa da reduo substancial da rea de ao, da perda de aderncia, da fonnao de fissuras e lascamento do concreto. O modelo tridimensional proposto por este trabalho baseado no processo de corroso iniciado por ataque de cloretos e propagado a partir de micropilhas de corroso, utilizando o mtodo dos elementos finitos. A influncia da temperatura, da umidade relativa do concreto, da ligao dos cloretos com compostos do concreto, da fonnao de fissuras provenientes do carregamento e da resistncia hmica do concreto so levados em considerao. O modelo dividido em 4 mdulos que interagem entre si, sendo eles: anlise estrutural, transferncia de calor, difuso de cloretos e corroso eletroqumica. O modelo computacional para a anlise estrutural utiliza dois procedimentos. No primeiro busca-se a resposta da estrutura para um carregamento instantneo, considerando-se o material com comportamento elasto-viscoplstico, associado com modelos de annadura incorporada e fissuras distribudas.No segundo, detennina-se a resposta da estrutura ao longo do tempo com a teoria da solidificao proposta por Bazant e Prasannan (1989), utilizando as cadeias Kelvin, apresentando timos resultados com uma grande economia computacional. O processo de transferncia de calor baseado na fonnulao tridimensional da equao de Helmholtz, com condies de contorno de fluxos provenientes da radiao solar, conveco e radiao, apresentando bons resultados em relao a fonnulaes analticas. Para anlise de vigas com a variao da temperatura mdia diria, constatou-se que as mesmas se encontram em equilbrio tnnico com o ambiente A anlise de difuso de cloretos implementada apresentou um comportamento idntico ao da segunda lei de Fick para casos unidimensionais. Para anlises de vigas carregadas se constatou que as fissuras provenientes de cargas de servio no influenciam o processo de difuso e que o coeficiente de difuso referencial o principal parmetro na anlise. O modelo de corroso eletroqumica baseado na fonnulao proposta por Liu (1996), obtendo-se como resultado a taxa de corroso em funo do tempo.

Relevância:

100.00% 100.00%

Publicador:

Resumo:

Este trabalho apresenta um modelo computacional para anlise tridimensional de peas de concreto armado, utilizando o mtodo dos elementos finitos. So utilizados modelos constitutivos elasto-viscoplsticos para representar o comportamento dos materiais, quando submetidos a cargas de curta ou longa durao. Se a estrutura est submetida a cargas de curta durao, o modelo desenvolvido fornece, quando alcanado o estado estvel, a soluo do problema elastoplstico. Em contrapartida, se a estrutura est submetida a cargas de longa durao, o modelo pode representar o comportamento viscoelstico com envelhecimento do concreto. Foi utilizado um modelo de camadas superpostas para assegurar uma representao adequada do comportamento real do concreto ao longo do tempo. Deste modo, admite-se que o concreto constitudo por um nmero conveniente de camadas, que apresentam a mesma deformao. Para a fissurao do concreto, utilizou-se um modelo de fissuras distribudas, que leva em considerao a contribuio do concreto entre fissuras O objetivo deste estudo a modelagem da armadura e do reforo nas peas de concreto armado, atravs de um modelo de aderncia entre o concreto e o material adjacente. A armadura introduzida no modelo como uma linha de material mais rgido dentro do elemento de concreto. E, a considerao da degradao da aderncia realizada atravs de um modelo onde os efeitos da aderncia imperfeita so includos pela introduo de graus de liberdade para os deslocamentos relativos entre os materiais. Os resultados obtidos com este programa computacional so comparados com valores experimentais disponveis na literatura, visando validar o modelo matemtico e a metodologia numrica.

Relevância:

100.00% 100.00%

Publicador:

Resumo:

O objetivo deste trabalho o desenvolvimento de um programa computacional, elaborado atravs do mtodo dos elementos finitos, para anlise tridimensional de situaes de puno em lajes de concreto armado com e sem armadura de cisalhamento. So utilizados modelos constitutivos elasto-viscoplsticos para representar o comportamento do concreto. Para a simulao do concreto submetido a um carregamento de longa durao, empregou-se o modelo de camadas superpostas. Este modelo admite que o concreto constitudo por um nmero conveniente de camadas com a mesma deformao, porm caractersticas materiais diferentes. A tenso total obtida pela soma das contribuies de cada camada. A armadura introduzida no modelo como uma linha de material mais rgido no interior do elemento de concreto. Considera-se aderncia perfeita entre o concreto e o ao. Assim, os deslocamentos ao longo das barras de ao so determinados a partir dos deslocamentos nodais dos elementos de concreto. Desta forma, a matriz de rigidez da armadura tem a mesma dimenso da matriz de rigidez do elemento de concreto, as quais so somadas para a obteno da matriz de rigidez total Para representar o comportamento aps fissurao do concreto, utilizado o modelo de fissuras distribudas, onde apenas as propriedades dos materiais so modificadas e a malha de elementos finitos no necessita de atualizao. Sugerem-se duas formulaes para a considerao da contribuio do concreto entre fissuras e da reduo na capacidade de transferncia de corte devido ao aumento da abertura da fissura. Este ltimo includo atravs de uma simplificao, que consiste em atribuir ao mdulo de elasticidade transversal, correspondente ao plano fissurado, um valor reduzido. Para comprovar a validade do programa computacional desenvolvido, comparam-se os resultados numricos com valores experimentais de ensaios realizados por Coelho (1999), Trautwein (2001) e Fusco (1988). Obtiveram-se boas aproximaes para lajes com e sem armadura de cisalhamento. Desta forma, o programa computacional desenvolvido pode ser empregado para generalizar resultados experimentais em lajes cogumelo com variao na resistncia do concreto, na bitola das barras de armadura, no espaamento entre as barras e no nmero de camadas de armadura de cisalhamento.

Relevância:

100.00% 100.00%

Publicador:

Resumo:

O objetivo deste trabalho desenvolver um programa computacional, baseado no mtodo dos elementos finitos, para simular situaes de reforo e recuperao de peas de concreto armado. A linguagem Matlab a ferramenta utilizada para a elaborao do programa. uma linguagem de alta performance para computao tcnica. Integra computao, visualizao e programao em um fcil ambiente para uso, onde problemas e solues so expressos em familiar notao matemtica. A caracterstica principal deste programa de permitir alteraes na modelagem numrica durante a anlise do problema, podendo-se retirar ou acrescentar elementos da estrutura, aumentar ou diminuir rigidezes, alterar materiais, etc, viabilizando-se, assim, uma avaliao das diversas etapas ou hipteses dos procedimentos de recuperao ou reforo de estruturas. O programa permite a mudana de vinculaes do elemento estrutural, durante a anlise. Alm disto, permitido no s incrementos de foras como incrementos de deslocamentos impostos estrutura. O programa computacional utiliza duas etapas distintas, baseadas em um modelo elasto-viscoplstico. Na primeira etapa, faz-se a determinao da resposta da estrutura ao longo do tempo. Considera-se, nesta etapa, que o material possui comportamento viscoelstico. Na segunda, busca-se determinar a resposta da estrutura para um carregamento instantneo, considerando-se o material com comportamento elastoplstico Para melhor representar o comportamento do concreto, quando submetido a carregamento de longa durao, utilizou-se um modelo de camadas superpostas. A armadura introduzida no modelo como uma linha de material mais rgido dentro do elemento de concreto. Considera-se aderncia perfeita entre o concreto e o ao. Assim, os deslocamentos ao longo das barras de ao so determinados a partir dos deslocamentos nodais dos elementos de concreto. Para verificar a preciso do programa computacional desenvolvido, comparam-se os resultados numricos com resultados experimentais, cujos ensaios foram realizados no Laboratrio de Ensaios e Modelos Estruturais do Programa de Ps-Graduao em Engenharia Civil da Escola de Engenharia da Universidade Federal do Rio Grande do Sul, e no Laboratrio de Estruturas da Escola de Engenharia de So Carlos- USP. Nessas comparaes, obtiveram-se timas aproximaes entre os resultados numricos e experimentais.

Relevância:

100.00% 100.00%

Publicador:

Resumo:

Este trabalho apresenta um estudo de fluxo de gua em barragens de terra, em regimes permanente e transiente, com a utilizao do Mtodo de Elementos Finitos. No estudo de fluxo em regime permanente duas formas de abordar o problema so apresentadas e comparadas. A primeira considera, para a discretizao da malha de elementos finitos, somente a regio saturada, de maneira que a linha fretica obtida atravs de ajustes desta malha de elementos finitos. A segunda considera toda a regio saturada-insaturada, sendo discretizado todo o domnio fsico da barragem. A malha de elementos finitos no modificada ao longo das iteraes e a linha fretica obtida por interpolao dentro dos elementos, em funo dos valores nodais do potencial de presses. O desenvolvimento terico das equaes utilizadas para as duas formas de abardagem apresentado, mostrando onde elas diferem entre si. No estudo de fluxo em regime transiente utilizado apenas o esquema de malha fixa de elementos finitos.

Relevância:

100.00% 100.00%

Publicador:

Resumo:

o presente trabalho tem por objetivo desenvolver um cdigo computacional destinado simulao numrica tridimensional e de escoamentos compressveis hipersnicos em especial para analisar o processo de reentrada de veculos espaciais na atmosfera terrestre. Devido s temperaturas elevadas encontradas neste tipo de escoamento, efeitos de altas temperaturas so includos, como excitao vibracional das molculas e reaes qumcas. Utiliza-se o modelo de 2 temperaturas (translacional e vibracional);5 espces (N2, O2, NO, N, O) e 17 reaes qumicas de Park [1986 para o modelamento do fenmeno termo-qumico. Introduz-se uma formulao baseada no mtodo dos elementos finitos, atravs do esquema de Taylor-Calerkin, para a discretizao do tempo e do espao e para a soluo do sistema de equaes. Alguns exemplos e aplicaes de problemas hipersnicos bi e tridimensionais, difusivos e no difusivos, so apresentados, incluindo comparaes entre os resultados obtidos no presente trabalho e na literatura.

Relevância:

100.00% 100.00%

Publicador:

Resumo:

o confinamento uma tcnica que tem sido cada vez mais aplicada no reforo e reabilitao de peas comprimidas de concreto, com o objetivo de aumentar sua capacidade resistente e recuper-las de eventuais deterioraes. Este trabalho apresenta um modelo no-linear fisico e geomtrico, baseado no Mtodo dos Elementos Finitos, para a anlise de pilares esbeltos de concreto armado, sob estado triplo de tenses. Neste estudo, so utilizados elementos finitos isoparamtricos tridimensionais para o concreto e um modelo incorporado de armadura. O comportamento do concreto descrito por um modelo elasto-viscoplstico, que permite a anlise da estrutura sob cargas de curta e longa durao. Para representao do concreto ao longo do tempo, utiliza-se um modelo de camadas superpostas. A fissurao do concreto representada atravs de um modelo de fissuras distribudas, onde se considera a contribuio do concreto entre fissuras. admitida aderncia perfeita entre os materiais. O ao modelado como um material elastoplstico perfeito ou com endurecimento linear, e os materiais utilizados para confinamento, como os compsitos de fibras de carbono, so modelados como um caso particular do modelo elastoplstico definido para o ao. O modelo no-linear geomtrico foi desenvolvido com base na formulao Lagrangeana Total. Consideram-se grandes deslocamentos e pequenas deformaes. Os resultados obtidos atravs do modelo computacional so cOnITontados com resultados experimentais disponveis na literatura, de modo a validar o modelo matemtico e a metodologia numrica. Boa correlao obtida entre os resultados numricos e experimentais, que confirmam estudos prvios realizados por diversos autores, na medida em que so verificados significativos ganhos de resistncia e/ou ductilidade em pilares comprimidos de concreto, devido ao confinamento.

Relevância:

100.00% 100.00%

Publicador:

Resumo:

A anlise da iterao solo-estrutura em fundaes um importante campo de pesquisa que ainda tem um grande progresso a ser feito. No presente trabalho foi desenvolvido um programa computacional para a anlise da interao solo-estrutura de fundaes de concreto armado. Este tema abrange duas reas da engenharia civil: estruturas e geotecnia. O mtodo dos elementos finitos foi usado no trabalho na seqncia para resolver o problema considerando estado plano de defonnao e comportamento elastoplsti.co dos materiais estudados (solo, concreto e ao). A linguagem de programao MATLAB foi usada em toda esta pesquisa como alternativa ao FORTRAN. O MATLAB foi escolhido uma vez que uma linguagem de programao que permite facilmente construir uma interfce de pr e psprocessamento amigvel. Os passos para a soluo completa do problema foram os seguintes: Primeiramente um programa foi desenvolvido considerando o comportamento elastoplstico com critrios de plastificao e ruptura especficos para o concreto e solo. Solues analticas fechadas foram usadas para checar a preciso do programa. O segundo passo foi a introduo do reforo de ao no concreto por meio de um modelo para armaduras. Logo aps, um modelo de fissuras para o concreto 1racionado foi in1roduzido no programa. Na seqncia o programa de pr e ps-processamento foi desenvolvido para gerar a malha de elementos finitos (pr-processamento), distribuio tenses e deformaes, mapa de fissuras, etc (psprocessamento). Finalmente, os parme1ros constitutivos do solo, concreto e ao foram calibrados e vrias situaes reais de interao do solo-concreto de fundaes de concreto armado foram simuladas. Nesta dissertao so encontrados resultados para as presses de contato sapata-solo. Diferentes diagramas de tenses de interfce foram obtidos em funo rigidez relativa do elemento estrutural de concreto armado-solo. Na anlise mnnrica, rigidez relativa desempenhou uma relevante funo no comportamento mecnico do elemento estrutural de concreto armado (sapata) e da base de assentamento (solo), uma vez ruptura em ambos os casos esteve diretamente relacionada a esta grandeza. So encon1rados, em funo da rigidez relativa, resultados indicativos dos modos de falha da fundao, excessiva plastificao do solo em fundaes com rigidez relativa alta, e a plastificao armaduras, esmagamento do concreto, formao de fissuras, bielas e confinamento concreto para fundaes de rigidez relativa baixa. Na anlise numrica, obteve-se resultados importantes com relao ao projeto de fundaes. Estes resultados foram cOnITontadoscom normas, destacando-se as discordncias com relao s recomendaes da nonna brasileira Projeto e Execuo de Fundaes" NBR-6122 (1996) para os diagramas de tenses interface sapata-solo usados no dimensionamento de fundaes de concreto armado.

Relevância:

100.00% 100.00%

Publicador:

Resumo:

O objetivo deste trabalho apresentar o desenvolvimento de um modelo computacional, baseado no mtodo dos elementos finitos, para o estudo de estruturas de concreto armado e protendido, sob estados planos de tenso. O modelo computacional utiliza dois procedimentos distintos, baseados em um modelo elasto-viscoplstico. No primeiro, determina-se a resposta da estrutura ao longo do tempo. Considera-se, nesta etapa, que o material possui comportamento viscoelstico. No segundo, busca-se a resposta da estrutura para um carregamento instantneo, considerando-se o material com um comportamento elastoplstico. Um modelo de camadas superpostas adotado para melhor representar o comportamento viscoelstico do concreto e da armadura protendida, onde estes materiais so compostos de diversas camadas que sofrem a mesma deformao. Cada camada possui diferentes propriedades mecnicas e a tenso total obtida pela soma das diferentes contribuies de cada camada. Para representar o concreto fissurado, utiliza-se um modelo de fissuras distribudas, que leva em conta a contribuio do concreto entre fissuras. Tanto a armadura passiva quanto a armadura de protenso so introduzidas no modelo como uma linha de material mais rgido dentro do elemento de concreto. Considera-se aderncia perfeita entre o concreto e o ao. Os deslocamentos ao longo das barras de ao so determinados a partir dos deslocamentos nodais dos elementos de concreto. Desta forma, as matrizes de rigidez das armaduras tm as mesmas dimenses das matrizes de rigidez dos elementos de concreto e ambas so somadas para a obteno das matrizes de rigidez totais. Para verificar a preciso do modelo computacional, compararam-se resultados numricos com valores experimentais disponveis.

Relevância:

100.00% 100.00%

Publicador:

Resumo:

O presente trabalho apresenta o estudo e implementao de um algoritmo numrico para anlise de escoamentos turbulentos, tridimensionais, transientes, incompressveis e isotrmicos, atravs da Simulao de Grande Escalas, empregando o Mtodo de Elementos Finitos. A modelagem matemtica do problema baseia-se nas equaes de conservao de massa e quantidade de movimento de um fluido quase-incompressvel. Adota-se um esquema de Taylor-Galerkin, com integrao reduzida e frmulas analticas das funes de interpolao, para o elemento hexadrico de oito ns, com funes lineares para as componentes de velocidade e constante no elemento para a presso. Para abordar o problema da turbulncia, emprega-se a Simulao de Grandes Escalas, com modelo para escalas inferiores resoluo da malha. Foram implementados o modelo clssico de Smagorinsky e o modelo dinmico de viscosidade turbulenta, inicialmente proposto por Germano et al, 1991. Uma nova metodologia, denominada filtragem por elementos finitos independentes, proposta e empregada, para o processo de segunda filtragem do modelo dinmico. O esquema, que utiliza elementos finitos independentes envolvendo cada n da malha original, apresentou bons resultados com um baixo custo computacional adicional. So apresentados resultados para problemas clssicos, que demonstram a validade do sistema desenvolvido. A aplicabilidade do esquema utilizado, para anlise de escoamentos caracterizados por elevados nmeros de Reynolds, discutida no captulo final. So apresentadas sugestes para aprimorar o esquema, visando superar as dificuldades encontradas com respeito ao tempo total de processamento, para anlise de escoamentos tridimensionais, turbulentos e transientes .

Relevância:

100.00% 100.00%

Publicador:

Resumo:

Neste trabalho apresenta-se um algoritmo para a simulao de problemas tridimensionais de interao fluido-estrutura utilizando a tcnica de elementos finitos. Um esquema de Taylor-Galerkin de dois passos e elementos tetradricos lineares so empregados para o fluido, que pode ser compressvel ou incompressvel. adotada uma formulao lagrangeana-euleriana arbitrria (ALE), compatvel com o movimento da interface fluidoestrutura. Um mtodo ftacionado de correo de velocidade utilizado para os fluidos incompressveis. A estrutura analisada usando elementos triangulares com trs ns e seis graus de liberdade por n (trs componentes de deslocamentos e trs componentes de rotao). Os efeitos da no-linearidade geomtrica so includos. O mtodo de Newmark empregado para integrar no tempo as equaes dinmicas de equilbrio, usando-se uma descrio lagrangeana atualizada. O sistema de equaes alge'bricas solucionado atravs do mtodo dos gradientes conjugados e o sistema no-linear, resultante de deslocamentos e rotaces finitas da estrutura, solucionado com um esquema incremental-iterativo. O cdigo otimizado para aproveitar as vantagens do processamento vetorial.

Relevância:

100.00% 100.00%

Publicador:

Resumo:

Aduelas protendidas externamente utilizadas em pontes contrastam com as clssicas construes monolticas, consistindo de pequenos segmentos pr-moldados protendidos por cordoalhas externas. H muitas vantagens na utilizao deste tipo de composio, como rapidez e versatilidade na construo, alto controle de qualidade e economia. No presente trabalho, discute-se uma formulao, baseada no mtodo dos elementos finitos, para simular o comportamento de estruturas compostas por aduelas externamente protendidas. Esta formulao permite tanto a anlise de situaes de servio, onde as sees permanecem completamente comprimidas, como tambm de situaes ltimas, onde se considera a abertura das juntas e a transferncia de esforos atravs das mesmas. O modelo computacional utiliza dois procedimentos distintos, baseados em um modelo elasto-viscoplstico. No primeiro, determina-se a resposta da estrutura ao longo do tempo, considerando-se, nesta etapa, que o material possui comportamento viscoelstico. No segundo, busca-se a resposta da estrutura para um carregamento instantneo, considerando-se o material com um comportamento elastoplstico. As juntas entre as aduelas so representadas atravs de elementos de contato unidimensionais. Os ns de cada elemento de contato coincidem com os ns do elemento de concreto que a ele foi associado.Um modelo de camadas superpostas utilizado para representar o comportamento viscoelstico do concreto e da armadura protendida Para representar o concreto fissurado, utiliza-se um modelo de fissuras distribudas, que leva em conta a contribuio do concreto entre fissuras. A armadura passiva introduzida na anlise atravs de um modelo incorporado e a armadura de protenso no aderente disposta em elementos fictcios, paralelos aos elementos reais da estrutura. Para verificar a adequao do modelo computacional, compararam-se resultados numricos com valores experimentais disponveis na literatura.

Relevância:

100.00% 100.00%

Publicador:

Resumo:

Neste trabalho, apresenta-se um estudo numrico de um modelo convectivo-difusivo-reativo em combusto baseado no Mtodo de Elementos Finitos. Primeiramente, apresenta-se o desenvolvimento das equaes de balano (quantidade de movimento, massa, espcie e energia) que modelam um processo de mistura molecular e reao qumica, irreversvel, de passo nico e exotrmica entre duas espcies qumicas F (Combustvel) e O (Oxidante). Tais espcies reagem e formam um produto P, conforme vFF +vOO ! vPP + calor, onde vF , vO e vP so os coeficientes estequiomtricos molares. No modelo, considera-se que a reao de primeira ordem com respeito a cada um dos reagentes e que a taxa de reao especfica segue a cintica de Arrhenius. Em seguida, o modelo estudado numericamente considerando-se um domnio retangular e condies de contorno do tipo Neumann. Tanto a Tcnica das Diferenas Finitas como a Tcnica de Elementos Finitos so utilizadas na discretizao espacial das equaes do modelo. Para a integrao no tempo, utiliza-se a mtodo de Runge-Kutta simplificado de trs estgios. Os diferentes cdigos computacionais obtidos, tanto pela Tcnica de Diferenas Finitas como de Elementos Finitos, so comparados frente ao problema de interesse. Observa-se que ambas as tcnicas apresentam resultados equivalentes. Alm disso, os cdigos desenvolvidos so robustos (capazes de lidar com vrios conjuntos de parmetros), de baixo custo e precisos. Por fim, apresenta-se uma reviso do trabalho de Zavaleta [48], no qual obtem-se uma estimativa local do erro na aproximao do problema estudado pela Tcnica de Elementos Finitos.