Análise dinâmica de rotores de grupos turbo e hidrogeradores com o método das matrizes de transferência

Conhecer o comportamento vibratório de rotores de máquinas rotativas, já durante a fase de projeto, é uma necessidade cada vez maior nos dias atuais, quando modificações corretivas, após o início da operação, já não são mais aceitáveis, devido aos altos custos de material e execução empregados, se comparados com os lucros que são cada vez mais reduzidos, em função das regras impostas pelo mercado, e devido também ao tempo necessário para tais modificações, que implica em manter a máquina parada, sem gerar lucros nesse penodo, acarretando em mu1tas contratuais para o fabricante. o presente trabalho visa dar subsídiosà análise dinâmicadas linhas de eixos de turbo e hidrogeradores (conjuntos formados pela turbina mais o gerador, usados respectivamentenas centraistermelétricase hidrelétricas). Inicialmente é feita uma breve explanação sobre formas construtivas de hidrogeradores e tipos de turbinas hidráu1icas utilizados atualmente. O trabalho prossegue apresentando as propriedades dos mancais com filme de óleo, os fenômenos do empuxo magnético no rotor do gerador e do efeito giroscópico, bem como a norma ISO 1940, que define os níveis admissíveis de desbalanceamentos residuais para rotores de máquinas É apresentada uma formulação matemática pelo Método das Matrizes de Transferência para as análises das vibrações livres, bem como das vibrações forçadas com carregamentos cíclicos devido aos desbalanceamentos residuais das linhas de eixos. A análise das vibrações livres é feita para a obtenção das fteqüências naturais torcionais e das fteqüências tlexionais críticas em um plano. A análise das vibrações forçadas é elaborada para fornecer as respostas dinâmicas em dois planos ortogonais, devido ao desbalanceamento residual. Para aplicação da metodologia de cálculo apresentada, são elaboradas rotinas computacionais de cálculos, utilizando como plataforma o software MATLAB. Para testar o Método das Matrizes de Transferência, bem como os programas de computador elaborados, faz-se algumas análises com soluções analíticas exatas e comparam-se resultados obtidos. Como aplicação prática, faz-se um estudo aprofundado dos comportamentos dinâmicos algumas linhas de eixos de hidrogeradores, para variações de configurações de construção dessas, bem como para as diversas condições de operação.

Interação de espirais em 2D : redução da dinâmica à interação de defeitos e exploração de novas estruturas espaço-temporais

O presente trabalho apresenta um anova proposta de tratamento de estruturas espirais em meios contínuos oscilatórios na vizinhança de bifurcações de Hopf supercríticas. Tais estruturas são normalmente descritas pela Equação de Cinzburg-Landau Complexa a qual usa um campo complexo associado a essas oscilações. A proposta apresentada reduz a dinâmica de espirais à interação entre os centros das mesmas. Inicialmente, comparamos numericamente as duas descrições e com os ganhos computacionais decorrentes da abordagem reduzida caracterizamos finamente as estruturas espaço-temporais formadas nesses sistemas: em vez dos estados congelados mencionados anteriormente na literatura encontrou-se uma dinâmica espaço-temporal intermitente. Esse regime ocorre em duas fases distintas: Líquido de Vórtices e Vidros de Vórtices. Esta última evolui em escalas de tempo ultralentas como fenômenos semelhantes encontrados na Mecânica Estatística, apesar de sua origem puramente determinista.

A dinâmica de uma família de aplicações unidimensionais

Resumo não disponível.

Biologia e dinâmica populacional de Aegla leptodactyla Buckup & Rossi, 1977 (Crustacea, Anomura, Aeglidae)

O presente trabalho teve por objetivo avaliar alguns aspectos relacionados à zoogeografia e ecologia dos ostracodes da plataforma interna norte do Rio de Janeiro, próximo à cidade de Cabo Frio. O material de estudo resultou da análise de 43 amostras sedimentológicas coletadas durante a expedição GEOCOSTA RIO II. Este estudo permitiu o reconhecimento de 16 famílias, 49 gêneros, 66 espécies e dois gêneros que permaneceram em nomenclatura indeterminada. A descrição de uma nova espécie, Actinocythereis brasiliensis sp. n., também foi realizada no âmbito do presente trabalho. As famílias mais representativas quanto à abundância e o número de espécies foram Thaerocytheridae e Cytheruridae, respectivamente. A espécie mais constante na área foi Caudites ohmerti e, dentre as com maior dominância estão, Caudites ohmerti, Meridionalicythere? dubia, Callistocythere litoralensis, Paracytheridea bulbosa, Urocythereis dimorphica, Henryhowella inflata, Oculocytheropteron delicatum, Xestoleberis umbonata, Henryhowella macrocicatricosa e Brasilicythere reticulispinosa. As preferências sedimentológicas e batimétricas das espécies são também discutidas neste trabalho. É sugerido que se estenda a distribuição zoogeográfica, um pouco mais ao norte, para sete espécies: Cushmanidea variopunctata, Neocytherideis impudicus, Hemingwayella advena, Hemingwayella sp., Urocythereis dimorphica, Munseyella cornuta e Basslerites costata. A grande maioria das espécies são características da plataforma continental sul, sendo que destas o maior percentual pertence à Subprovíncia Sul-Brasileira. De acordo com a fauna registrada, predominantemente de águas frias, chegamos a conclusão que a sua presença se deve, em grande parte, à influência da ressurgência marinha na área de estudo A proporção sexual foi de 1,19 machos para 1 fêmea. Fêmeas ovígeras encontradas entre abril e setembro e o ingresso de juvenis na população (recrutamento) foi observado em novembro e dezembro de 2000.

Os ambientes de geometria dinâmica e o pensamento hipotético-dedutivo

O processo de demonstração é axial na construção do conhecimento matemático. Na geometria euclidiana, ele é um dos aspectos que apresenta grandes obstáculos aos alunos. Uma das dificuldades aparece na transição, necessária, entre o conhecimento de natureza empírica, já adquirido, e aquele a ser construído: a geometria euclidiana enquanto modelo teórico, organizado em axiomas, teoremas e demonstrações. Os recursos informáticos hoje disponíveis provocam a busca de estratégias pedagógicas favoráveis à construção deste conhecimento. Entender as suas potencialidades torna-se um objeto de investigação: o que acontece com os processos cognitivos quando ao sujeito em interação com a máquina é possibilitada a concretização de seus construtos e ações mentais, e quando, mediante realimentação imediata, ele é levado a novas reelaborações e construções mentais? E como tais processos concorrem para um novo conhecimento? Esta tese propõe uma engenharia didática, em ambiente de geometria dinâmica, que favorece a ascensão dos alunos em patamar de conhecimento — de empírico a hipotético-dedutivo. Toma-se como referencial a teoria piagetiana, bem como a teoria da situação didática em matemática desenvolvida pela escola francesa. A engenharia se desenrola em três níveis: no primeiro, o propósito é a compreensão do significado e da necessidade de demonstração por via de construções geométricas; no segundo nível, pretende-se o desenvolvimento das primeiras habilidades na produção de demonstrações; e, no terceiro, os problemas propostos ao alunos exigem mais de seus funcionamentos cognitivos no tratamento adequado de uma figura geométrica — trata-se das extensões de desenho e concomitantes apreensões operativas responsáveis pela identificação de subconfigurações geométricas que dão suporte à argumentação dedutiva. Análise a posteriori do desenrolar dos trabalhos dos alunos confirma as expectativas anunciadas na análise a priori apresentada na fase de concepção da situação didática cuja implementação é proposta: o progresso dos alunos na construção de conhecimento em geometria, como modelo matemático, foi expressivo.

Novos modelos para cálculo de energia livre de solvatação em simulações de dinâmica molecular

Combinando métodos computacionais da eletrostática e dinâmica molecular, este trabalho teve como objetivo descrever processos de solvatação em sistemas quimicamente importantes. Foram obtidas propriedades termodinâmicas necessárias para o entendimento do processo de solvatação. Estão descritos e testados modelos para descrever a interação soluto-solvente, possibilitando, assim, aprimorar a descrição físico-química dos processos de solvatação. Utilizaram-se programas desenvolvidos em nosso grupo e programas comerciais que permitem os cálculos de dinâmica molecular e química quântica. Uma nova abordagem para o cálculo de energia livre de solvatação foi desenvolvida proporcionando a obtenção acurada e eficiente dessa propriedade, dentro do enfoque da dinâmica molecular. Nessa nova abordagem, novas metodologias para a geração de cavidades moleculares foram propostas e avaliadas. As energias livres de solvatação obtidas estão em boa concordância com os valores experimentais.

Derivação e solução de equações modelo da dinâmica de gases rarefeitos

Neste trabalho, duas equações modedlo na área da dinâmica de gases rarefeitos, são derivadas a partir de algumas soluções exatas da equação linearizada de Boltzmann homogênea e não homogênea. Em adição, uma versão analítca do método de ordenadas discretas é usado para resolver problemas clássicos nesta área, descritos pelo "Modelo S". Resultados numéricos são apresentados para os problemas de fluxo de Couette, fluxo de Poiseuille, "Creep" Térmico, Deslizamento Térmico e problema de Kramers.

Soluções unificadas para modelos com freqüência de colisão variável da dinâmica de gases rarefeitos

Neste trabalho, uma versão analítica do método de ordenadas discretas é usada para desenvolver soluções para alguns problemas da dinÂmica de gases rarefeitos, baseado em um modelo com freqüência de colisão variável (modelo CLF) da equação de Boltzmann linearizada. Em particular, resultados numéricos obtidos para os problemas de salto de temperatura, fluxo de Poiseuille, fluxo de Couette, Kramers, creep-térmico e deslizamento térmico são apresentados e discutidos.

Dinâmica de gases rarefeitos e transferência radiativa : aplicações em geometria cilíndrica

Neste trabalho são investigados problemas formulados em geometria cilíndrica na área da dinâmica de gases rarefeitos bem como na área de transferência radiativa. Com relação á dinâmica de gases rarefeitos, primeiramente são abordadas duas formas diferenciadas de se avaliar numericamente as funções de Chapmann-Enskog e de Burnett, necessárias na composição de soluções gerais nessa geometria. Em seguida é apresentada a derivação de uma equação integral baseada no modelo BGK para descrever o fluxo de um gás rarefeito em um tubo cilíndrico. Problemas relacionados á transferência radiativa, incluindo o caso não-linear acoplado radiação-condução, são solucionados com a aplicação de uma versão reformulada do método de ordenadas discretas, sendo que resultados numéricos relevantes a estes problemas são também apresentados.

Dinâmica de paredes de domínios magnéticos : um estudo através da impedanciometria

Neste trabalho é apresentado um estudo sobre a dinâmica de paredes de domínios através de medidas de impedanciometria. É proposto um método que permite a obtenção dos seguintes parâmetros das amostras estudadas: mobilidade e velocidade crítica das paredes, largura dos domínios e de suas paredes, densidade de energia de parede e constante efetiva de troca. Todas essas informações são obtidas a partir do espectro em freqüência da permeabilidade e de relações apropriadas entre a permeabilidade e a impedância complexa. O elo de ligação entre essas quantidades é feita através do efeito da profundidade de penetração, cuja definição inclui a freqüência, a resistividade e a permeabilidade do material. O método foi aplicado ao estudo de dois tipos diferentes de materiais, (i) (110)[001]FeSi3%, policristalino e altamente texturizado com tamanho de grãos bastante grande (~ 5 mm) e (ii) amostras nanocristalinas obtidas através do recozimento de fitas amorfas de Fe73.5Cu1Nb3Si16.5B6. Enquanto o primeiro sistema foi utilizado para se fazer uma comparação entre os parâmetros aqui obtidos com aqueles de outros autores e técnicas, o segundo foi estudado em termos das modificações da anisotropia magnética associadas ao alívio das tensões internas com a temperatura de recozimento.

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.

Adaptação dinâmica do timeout de detectores de defeitos através do uso de séries temporais

Uma aplicação distribuída freqüentemente tem que ser especificada e implementada para executar sobre uma rede de longa distância (wide-área network-WAN), tipicamente a Internet. Neste ambiente, tais aplicações são sujeitas a defeitos do tipo colapso(falha geral num dado nó), teporização (flutuações na latência de comunicação) e omissão (perdas de mensagens). Para evitar que este defeitos gerem comseqüências indesejáveis e irreparáveis na aplicação, explora-se técnicas para tolerá-los. A abstração de detectores de defeitos não confiáveis auxilia a especificação e trato de algoritmos distribuídos utilizados em sistemas tolerantes a falhas, pois permite uma modelagem baseada na noção de estado (suspeito ou não suspeito) dos componentes (objetos, processo ou processadores) da aplicação. Para garantir terminação, os algoritmos de detecção de defeitos costumam utilizar a noção de limites de tempo de espera (timeout). Adicionalmente, para minimizar seu erro (falasas suspeitas) e não comprometer seu desempenho (tempo para detecção de um defeito), alguns detectores de defeitos ajustam dinamicamente o timeout com base em previsões do atraso de comunicação. Esta tese explora o ajuste dinâmico do timeout realizado de acordo com métodos de previsão baseados na teoria de séries temporais. Tais métodos supõem uma amostragem periódica e fornececm estimativas relativamente confiáveis do comportamento futuro da variável aleatória. Neste trabalho é especificado uma interface para transformar uma amostragem aperiódica do atraso de ida e volta de uma mensagem (rtt) numa amostragem periódica, é analisado comportamento de séries reais do rtt e a precisão dee sete preditores distintos (três baseados em séries temporais e quatrro não), e é avaliado a influência destes preditores na qualidade de serviço de um detector de defeitos do estilopull. Uma arquitetura orientada a objetos que possibilita a escolha/troca de algoritmos de previsão e de margem de segurança é também proposta. Como resultado, esta tese mostra: (i) que embora a amostragem do rtt seja aperiódica, pode-se modelá-la como sendo uma série temporal (uma amostragem periódica) aplciando uma interface de transformação; (ii) que a série temporal rtt é não estacionária na maioria dos casos de teste, contradizendo a maioria das hipóteses comumente consideradas em detectores de defeitos; (iii) que dentre sete modelos de predição, o modelo ARIMA (autoregressive integrated moving-average model) é o que oferece a melhor precisão na predição de atrasos de comunicação, em termos do erro quadrático médio: (iv) que o impacto de preditores baseados em séries temporais na qualidade de serviço do detector de defeitos não significativo em relação a modelos bem mais simples, mas varia dependendo da margem de segurança adotada; e (v) que um serviço de detecção de defeitos pode possibilitar a fácil escolha de algoritmos de previsão e de margens de segurança, pois o preditor pode ser modelado como sendo um módulo dissociado do detector.

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.

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.