Simulação numérica da cavitação em turbomáquinas usando uma formulação Euler-Lagrange.

Turbomáquinas são máquinas operacionais que transferem energia mecânica entre um rotor e um fluido. Estas máquinas têm muitas aplicações industriais. Um dos componentes de uma turbomáquina responsável pela transferência da energia, ou receber a rotação do eixo e transformar em energia de fluido em caso de bomba ou transferir a energia do fluido para o eixo em caso de uma turbina, é o impelidor ou rotor. O fenómeno da cavitação envolve escoamento bifásico: o líquido a ser bombeado e as bolhas de vapor que são formadas durante o processo de bombeamento. O processo de formação dessas bolhas é complexo, mas ocorre principalmente devido a presença de regiões de pressões muito baixas. O colapso dessas bolhas pode muitas vezes levar a deterioração do material, dependendo da intensidade ou da velocidade de colapso das bolhas. O principal objetivo deste trabalho foi estudar o comportamento hidrodinâmico do escoamento nos canais do impelidor de uma turbomáquina do tipo radial usando recursos de fluidodinâmica computacional (CFD). Uma abordagem Euler-Lagrange acoplada com o modelo da equação de Langevin foi empregada para estimar a trajetória das bolhas. Resultados das simulações mostram as particularidades de um escoamento líquido-bolha de vapor passando em um canal de geometria curva, fornecendo assim informações que podem nos ajudar na prevenção da cavitação nessas máquinas.

Identificação da fase sigma em uma junta soldada pelo processo TIG manual autógeno de um aço hiperduplex SAF 2707HD (UNS S32707) por difração de raios-x.

O aço inoxidável hiperduplex possui alta resistência a corrosão por pite em ambientes contendo cloretos, quando comparado a outros aços inoxidáveis comercialmente conhecidos. Possui boas propriedades mecânicas, com limite de escoamento superior a 700MPa e limite de resistência a tração em torno de 1000MPa. Essas propriedades o tornam muito atrativos para aplicações em ambientes contendo cloretos, e por isso tem tido destaque na indústria de óleo e gás, refinarias, plataformas offshore, etc. A liga hiperduplex é composta por uma estrutura bifásica, contendo proporções aproximadamente iguais de ferrita e austenita. Esse material possui boa soldabilidade, mas por ser termodinamicamente metaestável, em altas temperaturas pode ocorrer a precipitação de fases intermetálicas não desejáveis, o que resulta em perda de propriedades mecânicas e diminuição da resistência a corrosão. A fase sigma tem sido fortemente estudada, pois é comum sua precipitação nos aços inoxidáveis da família duplex durante o procedimento de soldagem se este não for muito bem controlado. A fase sigma precipita preferencialmente na fase ferrítica, devido a maior concentração de Cr e Mo, que são os elementos formadores da fase. A resistência a corrosão é reduzida e as propriedades mecânicas do material são alteradas o tornando frágil devido a presença da fase sigma. É formada entre 600C e 1000C e possui uma estrutura tetragonal complexa. O objetivo do trabalho foi identificar a possível presença da fase sigma na junta soldada do aço inoxidável hiperduplex SAF 2707 HD (UNS S32707) pelo processo TIG autógeno manual através da difração de raios-x. Nessa pesquisa, foram analisadas uma junta soldada do material pelo processo TIG autógeno manual com arco pulsado. Complementando o estudo foram analisadas seis amostras do aço inoxidável superduplex, sendo que cinco amostras sofreram tratamento térmico para a proposital formação da fase sigma. O refinamento do resultado da difração das amostras foi feito utilizando o método de Rietveld no software Topas Academic versão 4.1. O resultado da amostra soldada de hiperduplex apresentou as fases austenita, ferrita e alguns prováveis óxidos. Os resultados das amostras de superduplex tratadas termicamente apresentaram a fase sigma, conforme esperado na pesquisa, e as fases austenita e ferrita.

Implementação paralela de um código de elementos finitos em 2D para as Equações de Navier-Stokes para fluidos incompressíveis com transporte de escalares.

O estudo do fluxo de água e do transporte escalar em reservatórios hidrelétricos é importante para a determinação da qualidade da água durante as fases iniciais do enchimento e durante a vida útil do reservatório. Neste contexto, um código de elementos finitos paralelo 2D foi implementado para resolver as equações de Navier-Stokes para fluido incompressível acopladas a transporte escalar, utilizando o modelo de programação de troca de mensagens, a fim de realizar simulações em um ambiente de cluster de computadores. A discretização espacial é baseada no elemento MINI, que satisfaz as condições de Babuska-Brezzi (BB), que permite uma formulação mista estável. Todas as estruturas de dados distribuídos necessárias nas diferentes fases do código, como pré-processamento, solução e pós-processamento, foram implementadas usando a biblioteca PETSc. Os sistemas lineares resultantes foram resolvidos usando o método da projeção discreto com fatoração LU por blocos. Para aumentar o desempenho paralelo na solução dos sistemas lineares, foi empregado o método de condensação estática para resolver a velocidade intermediária nos vértices e no centróide do elemento MINI separadamente. Os resultados de desempenho do método de condensação estática com a abordagem da solução do sistema completo foram comparados. Os testes mostraram que o método de condensação estática apresenta melhor desempenho para grandes problemas, às custas de maior uso de memória. O desempenho de outras partes do código também são apresentados.

Análise vibratória de fundações de máquinas sobre estacas.

A análise de fundações sob solicitações dinâmicas é algo sempre presente em projetos na área industrial. É um campo pouco explorado na área de engenharia geotécnica, onde existem relativamente poucas informações no Brasil, de maneira geral. O método mais comum de realizar essas análises é a simplificação de modelos estruturais a partir do uso de molas. Sabe-se que esses coeficientes de reação têm uma variação relativamente grande e que esse enfoque de projeto pode, em alguns casos, mostrar-se contra a segurança ou levar a superdimensionamentos desnecessários. Verifica-se, então, a necessidade de uma avaliação mais criteriosa, utilizando a interação solo x estrutura, onde as molas comumente utilizadas nas análises vibratórias convencionais são substituídas pela rigidez real do solo quando concebido como um meio contínuo, através de sua discretização pelo método dos elementos finitos. A presente dissertação analisa o problema através do módulo de dinâmica do programa Plaxis 2D. Neste tipo de análise, além da modelagem do solo como um meio contínuo, torna-se possível introduzir condições de contorno específicas ao problema em estudo, múltiplas camadas de solo, sejam horizontais ou inclinadas, além da introdução de amortecedores capazes de evitar a reflexão espúria das ondas incidentes nos limites da malha de elementos finitos e assim modelar mais adequadamente a perda de energia por radiação. A presente dissertação compara medições experimentais e soluções eficientes de métodos vibratórios clássicos com a resposta obtida pelo MEF, mostrando resultados bastante satisfatórios tanto pelos métodos clássicos quanto pelo MEF.

Aplicação do modelo hidrológico SWMM na gestão das águas pluviais urbanas: estudo de caso da bacia hidrográfica do Rio Morto, Rio de Janeiro.

O crescimento da população e dos núcleos urbanos durante o século XX, sobretudo nos países em desenvolvimento, contribuiu para o aumento das áreas impermeáveis das bacias hidrográficas, com impactos importantes nos sistemas de drenagem urbana e na ocorrência de enchentes associadas. As enchentes trazem prejuízos materiais, na saúde e sociais. Recentemente, têm sido propostas práticas conservacionistas e medidas compensatórias, que buscam contribuir para o controle das enchentes urbanas, através do retardo do pico e amortecimento dos hidrogramas. Modelos matemáticos hidrológicos-hidráulicos permitem a simulação da adoção destas medidas de controle, demonstrando e otimizando sua localização. Esta dissertação apresenta os resultados da aplicação do modelo hidrológico Storm Water Management Model (SWMM) à bacia hidrográfica de estudo e representativa do rio Morto localizada em área peri-urbana em Jacarepaguá na cidade do Rio de Janeiro, com área de 9,41 km. O processamento do modelo SWMM foi realizado com o apoio da interface Storm and Sanitary Analysis (SSA), integrada ao sistema AutoCAD Civil 3D. Além da verificação da adequabilidade do modelo à representação dos sistemas hidrológico e hidráulico na bacia, foram desenvolvidos estudos para dois cenários como medidas de controle de enchentes: cenário 1, envolvendo implantação de um reservatório de detenção e, cenário 2, considerando a implantação de reservatórios de águas pluviais nos lotes. Os hidrogramas resultantes foram comparados ao hidrograma resultante da simulação nas condições atuais. Além disso, foram avaliados os custos associados a cada um dos cenários usando o sistema de orçamento da Empresa Rio Águas da PCRJ. Nas simulações foram adotadas a base cartográfica, e os dados climatológicos e hidrológicos previamente observados no contexto do projeto HIDROCIDADES, Rede de Pesquisa BRUM/FINEP, na qual este estudo se insere. Foram representados os processos de geração e propagação do escoamento superficial e de base. Durante o processo de calibração, realizou-se a análise de sensibilidade dos parâmetros, resultando como parâmetros mais sensíveis os relativos às áreas impermeáveis, especialmente o percentual de área impermeável da bacia (Ai). A calibração foi realizada através do ajuste manual de sete parâmetros do escoamento superficial e cinco do escoamento de base para três eventos. Foram obtidos coeficientes de determinação entre 0,52 e 0,64, e a diferença entre os volumes escoados e observados entre 0,60% e 4,96%. Para a validação do modelo foi adotado um evento pluviométrico excepcional observado na cidade em abril de 2010, que à época causou enchentes e grandes transtornos na cidade. Neste caso, o coeficiente de determinação foi igual a 0,78 e a diferença entre volumes foi de 15%. As principais distorções entre hidrogramas observados e simulados foram verificados para as vazões máximas. Em ambos os cenários as enchentes foram controladas. A partir destes estudos, pôde-se concluir que o melhor custo-benefício foi o cenário 2. Para este cenário, foi observado maiores amortecimento e retardo da vazão de pico do hidrograma, igual a 21,51% da vazão simulada para as condições atuais da bacia. Os custos de implantação orçados para os reservatórios de lote ficaram 52% a menos do que o do reservatório de detenção.

Jato transversal de gotas: simulações por ALE/FEM e efeitos interfaciais.

Um código computacional para escoamentos bifásicos incorporando metodologia híbrida entre oMétodo dos Elementos Finitos e a descrição Lagrangeana-Euleriana Arbitrária do movimento é usado para simular a dinâmica de um jato transversal de gotas na zona primária de quebra. Os corpos dispersos são descritos por meio de um método do tipo front-tracking que produz interfaces de espessura zero através de malhas formadas pela união de elementos adjacentes em ambas as fases e de técnicas de refinamento adaptativo. Condições de contorno periódicas são implementadas de modo variacionalmente consistente para todos os campos envolvidos nas simulações apresentadas e uma versão modificada do campo de pressão é adicionada à formulação do tipo um-fluido usada na equação da quantidade de movimento linear. Simulações numéricas diretas em três dimensões são executadas para diferentes configurações de líquidos imiscí veis compatíveis com resultados experimentais encontrados na literatura. Análises da hidrodinâmica do jato transversal de gotas nessas configurações considerando trajetórias, variação de formato de gota, espectro de pequenas perturbações, além de aspectos complementares relativos à qualidade de malha são apresentados e discutidos.

Uma nova descrição para a transferência de massa em meios porosos com transição saturado-insaturado.

Esse texto trata do problema de um fluido contaminado escoando por um meio poroso, tratando os componentes na mistura como meios contínuos. Na primeira parte, desenvolvemos a teoria de misturas de meios contínuos e discutimos equações da continuidade, momento linear e momento angular. A seguir, descrevemos o problema em detalhe e fazemos hipóteses para simplificar o escoamento. Aplicamos as equações encontradas anteriormente para encontrarmos um sistema de equações diferenciais parciais. Desse ponto em diante, o problema se torna quase puramente matemático. Discutimos o caso insaturado, e depois a saturação do meio poroso. Finalmente, adicionamos um contaminante à mistura e, em seguida, N contaminantes.

Determinação por ultrassom de constantes elásticas de pastilhas de alumina utilizadas em varetas combustíveis de reatores nucleares.

Pastilhas de dióxido de urânio (UO2) são empregadas como combustível nos reatores nucleares de potência, onde as condições operacionais introduzem elevados gradientes térmicos nas pastilhas. Potências elevadas propiciam a fusão da parte central das pastilhas. O inchamento das pastilhas de dióxido de urânio (UO2), decorrente dos produtos de fissão, pode causar o trincamento do material, em função disso, o estudo do seu comportamento mecânico é importante. Esse trabalho avaliou a aplicação de método ultrassônico na obtenção de constantes elásticas. Em função das dificuldades no manuseio do material nuclear, optou-se por um estudo comparativo em pastilhas de alumina (Al2O3). Foram fabricados e usados dois conjuntos de 10 pastilhas de Al2O3 com densidades de 92% e 96%. Foi desenvolvida técnica ultrassônica por transmissão obtendo medidas do tempo de percurso de ondas ultrassônicas, longitudinais e transversais, usadas para a determinação das constantes elásticas do material. Equações relacionando a velocidade da onda ultrassônica ao módulo de elasticidade, módulo de cisalhamento e coeficiente de Poisson, permitiram obter essas constantes elásticas, que apresentaram excelente concordância com a literatura disponível para o Al2O3.

Desenvolvimento e implementação de algoritmo para análise de dutos fabricados ou reparados com material compósito carregados com pressão interna utilizando o método dos elementos finitos.

A integridade de dutos ganha importância à medida em que o desenvolvimento da indústria expande a malha para transportar líquidos e gases num contexto de aumento das exigências ambientais. Um importante aliado para manutenção da integridade de dutos são reparos de materiais compósitos. Estes materiais apresentam baixa densidade, capacidade de direcionar resistência de acordo com as direções das solicitações, execução de reparo a frio sem necessidade de interromper produção ou grande maquinário. Este trabalho tem como objetivo desenvolver e implementar um algoritmo de elementos finitos que permita avaliar os esforços e a resistência das paredes de um tubos fabricados ou reparados com laminados de material compósito carregados com pressão interna. Entre as vantagens de desenvolver um programa tem-se: agilidade de avaliação, menor custo com licença, menores exigências computacionais, possibilidade de desenvolver o programa e o melhor entendimento da modelagem dos fenômenos. Utiliza-se como entrada do programa o diâmetro do duto, pressão interna e parâmetros do laminado. A modelagem em elementos finitos é realizada a partir da teoria clássica de laminados. Aplicando o carregamento resultante da pressão interna, determina-se os deslocamentos e são calculadas as tensões e aplicado o critério de falha de Tsai-Hill em cada camada. Estudos experimentais e numéricos encontrados na literatura foram simulados com o programa gerado e os resultados para propriedades do laminado, tensões nos dutos e pressão de ruptura apresentam concordância com os resultados da literatura.O programa ainda tem sua estrutura modificada para encontrar a pressão de falha a partir dos dados do laminado. O programa implementado permite uma avaliação rápida de resistência do reparo e possibilita avaliar rapidamente a resposta a mudanças nos parâmetros de projeto do laminado.

Análise de vibrações e estudo de conforto humano sobre pisos mistos (aço-concreto) submetidos a ações humanas rítmicas.

Este trabalho tem por objetivo avaliar o comportamento dinâmico de pisos mistos (aço-concreto) sob a ação de cargas provenientes das atividades humanas rítmicas, especificamente a prática de ginástica aeróbica, sob o ponto de vista do conforto humano. Tal avaliação torna-se necessária por crescentes problemas estruturais associados às vibrações excessivas, decorrentes da concepção de sistemas estruturais com baixos níveis de amortecimento e com frequências naturais cada vez mais baixas e bastante próximas das faixas de frequência das excitações associadas às atividades humanas rítmicas. O modelo estrutural investigado baseiase em um piso misto (aço-concreto) submetido a aulas de ginástica aeróbica. A modelagem numérica do piso misto investigado foi realizada com base no emprego do programa ANSYS e foram utilizadas técnicas de discretização por meio do método dos elementos finitos (MEF). As cargas aplicadas sobre o piso, oriundas das atividades aeróbicas, são simuladas através de dois modelos de carregamentos dinâmicos distintos. Uma extensa análise paramétrica foi desenvolvida sobre o modelo estrutural investigado e a resposta dinâmica do sistema foi obtida, em termos dos deslocamentos e das acelerações, e comparada com os limites recomendados por normas e critérios de projeto. A resposta dinâmica do piso estudado viola os critérios de projeto relativos ao conforto humano e indica níveis de vibrações excessivas nos casos de carregamento dinâmicos analisados nesta dissertação.

Simulação de problemas de transferência de calor em regime permanente com uma relação entre condutividade térmica e temperatura constante por partes.

Este trabalho estuda a transferência de calor por condução considerando a condutividade térmica como uma função constante por partes da temperatura. Esta relação, embora fisicamente mais realista que supor a condutividade térmica constante, permite obter uma forma explícita bem simples para a inversa da Transformada de Kirchhoff (empregada para tratar a não linearidade do problema). Como exemplo, apresenta-se uma solução exata para um problema com simetria esférica. Em seguida, propôe-se uma formulação variacional (com unicidade demonstrada) que introduz um funcional cuja minimização é equivalente à solução do problema na forma forte. Finalmente compara-se uma solução exata obtida pela inversa da Transformada de Kirchhoff com a solução obtida via formulação variacional.

Comportamento mecânico de cerâmica produzida a partir de reciclagem de para-brisa de automóvel

É crescente o interesse nos materiais cerâmicos, devido as suas características como baixa massa específica e maior resistência a ambientes agressivos do que a maioria das ligas metálicas. Este trabalho tem o objetivo de produzir a partir de material descartado, cerâmicas em diferentes temperaturas de sinterização e avaliar a sua tensão de ruptura em ensaio de flexão de três pontos e a confiança desta medida. Devido ao custo de produzir pó de vidro, tanto pelo alto gasto de energia para fundir a matéria-prima como pelo consumo de minerais industriais, foi proposto utilizar vidro de para-brisa obsoleto reduzindo despesas operacionais e definindo um destino econômico e ambiental viável para estes rejeitos. A metodologia consistiu-se na obtenção do pó de vidro com características adequadas para ser conformado e sinterizado. Foram usadas duas composições e quatro tratamentos térmicos para obter oito materiais. Uma composição com apenas o pó oriundo da moagem de para-brisa e outra com este pó mais 4% de óxido de nióbio. A resistência à flexão dos produtos obtidos foi avaliada. Utilizou-se a estatística de Weibull para caracterizar estes resultados. Os resultados obtidos indicam que o material de composição pó de vidro e temperatura final de sinterização de 650C obteve a maior resistência mecânica entre os materiais sintetizados. A adição do óxido de nióbio provoca um decréscimo na resistência mecânica se comparada com o material sem a adição deste óxido. Entretanto, comparando as duas composições na mesma temperatura final de sinterização, a adição de óxido de nióbio provocou um aumento no módulo de Weibull, excetuando-se dois de oito materiais obtidos. As diferentes composições e temperaturas de sinterização afetaram as propriedades mecânicas dos materiais obtidos.

Estudo paramétrico para um muro de gravidade submetido a carregamentos sísmicos.

A presente dissertação apresenta um estudo do comportamento de um muro de gravidade flexível submetido a carregamentos sísmicos. A influência do carregamento sísmico, e mais precisamente da variação da aceleração horizontal de pico é avaliada num estudo paramétrico, através da análise dos resultados obtidos para os deslocamentos e acelerações horizontais que ocorreram no paramento do muro de gravidade. Optou-se pela análise ao longo do tempo, introduzindo o carregamento sísmico ao modelo numérico do qual fazia parte o muro de gravidade através de 5 históricos temporais de acelerogramas horizontais normalizados em 0,05g, 0,10g, 0,15g, 0,20g e 0,25g. O evento sísmico de origem dos acelerogramas normalizados foi o terremoto ocorrido no Chile em 27 de fevereiro de 2010 e medido na estação em San Jose de Maipo, a 332,7km do epicentro do terremoto. Utiliza-se um software computacional aplicado à Geotecnia, o Plaxis, para a calibração de um modelo numérico em situação estática e posterior desenvolvimento das análises paramétricas em condições sísmicas. As análises realizadas para avaliação do comportamento do muro de gravidade em condição sísmica foram a verificação da influência da variação do acelerograma horizontal normalizada introduzido na base do modelo nos valores das acelerações e deslocamentos horizontais atuantes no muro e a verificação da influência do carregamento sísmico em comparação com a situação estática, comparando-se os valores dos deslocamentos horizontais obtidos na análise estática e nas análises dinâmicas. Os fatores que influenciam potencialmente nos resultados numéricos em condições sísmicas foram apresentados, destacando-se os aspectos relacionados à modelagem numérica em condição de sismo. Os modelos constitutivos oferecidos, a possibilidade de incorporação do carregamento sísmico na base do modelo e a possibilidade de consideração do amortecimento visco-elástico de Rayleigh nos materiais são as principais qualidades do Plaxis, utilizado na modelagem numérica. Avaliou-se positivamente a potencialidade do Plaxis, visto que o programa mostrou-se uma ferramenta capaz de simular o comportamento de muros de gravidade sujeitos a carregamentos sísmicos.

Uma Metodologia de Estudo de Simulação Tridimensional de Escoamento Turbulento Estratificado no Reservatório de Plantas Hidrelétricas.

Uma simulação numérica que leva em conta os efeitos de estratificação e mistura escalar (como a temperatura, salinidade ou substância solúvel em água) é necessária para estudar e prever os impactos ambientais que um reservatório de usina hidrelétrica pode produzir. Este trabalho sugere uma metodologia para o estudo de escoamentos ambientais, principalmente aqueles em que o conhecimento da interação entre a estratificação e mistura pode dar noções importantes dos fenômenos que ocorrem. Por esta razão, ferramentas de simulação numérica 3D de escoamento ambiental são desenvolvidas. Um gerador de malha de tetraedros do reservatório e o modelo de turbulência algébrico baseado no número de Richardson são as principais ferramentas desenvolvidas. A principal dificuldade na geração de uma malha de tetraedros de um reservatório é a distribuição não uniforme dos pontos relacionada com a relação desproporcional entre as escalas horizontais e verticais do reservatório. Neste tipo de distribuição de pontos, o algoritmo convencional de geração de malha de tetraedros pode tornar-se instável. Por esta razão, um gerador de malha não estruturada de tetraedros é desenvolvido e a metodologia utilizada para obter elementos conformes é descrita. A geração de malha superficial de triângulos utilizando a triangulação Delaunay e a construção do tetraedros a partir da malha triangular são os principais passos para o gerador de malha. A simulação hidrodinâmica com o modelo de turbulência fornece uma ferramenta útil e computacionalmente viável para fins de engenharia. Além disso, o modelo de turbulência baseado no número de Richardson leva em conta os efeitos da interação entre turbulência e estratificação. O modelo algébrico é o mais simples entre os diversos modelos de turbulência. Mas, fornece resultados realistas com o ajuste de uma pequena quantidade de parâmetros. São incorporados os modelos de viscosidade/difusividade turbulenta para escoamento estratificado. Na aproximação das equações médias de Reynolds e transporte de escalar é utilizando o Método dos Elementos Finitos. Os termos convectivos são aproximados utilizando o método semi-Lagrangeano, e a aproximação espacial é baseada no método de Galerkin. Os resultados computacionais são comparados com os resultados disponíveis na literatura. E, finalmente, a simulação de escoamento em um braço de reservatório é apresentada.

Métodos sem malha: aplicações do Método de Galerkin sem elementos e do Método de Interpolação de Ponto em casos estruturais.

Apesar de serem intensamente estudados em muitos países que caminham na vanguarda do conhecimento, os métodos sem malha ainda são pouco explorados pelas universidades brasileiras. De modo a gerar uma maior difusão ou, para a maioria, fazer sua introdução, esta dissertação objetiva efetuar o entendimento dos métodos sem malha baseando-se em aplicações atinentes à mecânica dos sólidos. Para tanto, são apresentados os conceitos primários dos métodos sem malha e o seu desenvolvimento histórico desde sua origem no método smooth particle hydrodynamic até o método da partição da unidade, sua forma mais abrangente. Dentro deste contexto, foi investigada detalhadamente a forma mais tradicional dos métodos sem malha: o método de Galerkin sem elementos, e também um método diferenciado: o método de interpolação de ponto. Assim, por meio de aplicações em análises de barras e chapas em estado plano de tensão, são apresentadas as características, virtudes e deficiências desses métodos em comparação aos métodos tradicionais, como o método dos elementos finitos. É realizado ainda um estudo em uma importante área de aplicação dos métodos sem malha, a mecânica da fratura, buscando compreender como é efetuada a representação computacional da trinca, com especialidade, por meio dos critérios de visibilidade e de difração. Utilizando-se esses critérios e os conceitos da mecânica da fratura, é calculado o fator de intensidade de tensão através do conceito da integral J.