Made available in DSpace on 2016-12-23T14:08:13Z (GMT). No. of bitstreams: 1 Dissertacao de Felipe Patricio das Neves.pdf: 512132 bytes, checksum: d42e56e29214989211674c822342bbea (MD5) Previous issue date: 2009-12-04 / In this work is implemented a numerical model to simulate computationally the distribution of pressures, velocities, temperatures and heat flows in two-dimensional stationary control volumes. The relation between temperatures and velocities is established by the advective-diffusive Equation, using the Dual Reciprocity Boundary Element Method formulation... / Neste trabalho é desenvolvido um modelo numérico para simular computacionalmente a distribuição de pressões, velocidades, temperaturas e fluxos de calor estacionários em volumes de controle bidimensionais. A relação do campo de temperaturas e velocidades é governada pela equação da Difusão-Advecção, resolvida através da formulação com Dupla Reciprocidade do Método dos Elementos de Contorno. Admite-se a lei de Darcy para associar pressão e velocidade, resultando num modelo matemático dado pela Equação de Laplace, no caso linear. Na análise não-linear insere-se a dependência entre do campo de velocidades e as temperaturas, resultando num campo matematicamente representado pela Equação de Poisson. Os resultados da solução desse problema são então implementados no modelo difusivo-advectivo, gerando temperaturas e fluxos de calor

Made available in DSpace on 2016-12-23T14:08:15Z (GMT). No. of bitstreams: 1 Atila Lupim Cruz.pdf: 1394501 bytes, checksum: 0954b2c5b1fdcb864ee81cef7d14e9e5 (MD5) Previous issue date: 2012-10-19 / A pesquisa envolvida na presente dissertação se baseou no uso de funções de base radial para gerar uma nova formulação integral, que interpola diretamente o termo não homogêneo da equação diferencial de governo, no contexto do Método dos Elementos de Contorno (MEC). Emprega-se o uso de funções primitivas das funções de interpolação originais no núcleo da integral de domínio, permitindo a transformação desta última numa integral de contorno, evitando assim a discretização do domínio por meio de células, semelhante ao realizado na Dupla Reciprocidade. Para melhor avaliação das potencialidades da formulação, os testes numéricos apresentados abordaram apenas a solução de problemas governados pela Equação de Poisson. Os problemas escolhidos dentro desta categoria possuem solução analítica, o que permitiu aferir com mais rigor a precisão dos resultados. Para melhor balizamento da eficiência da formulação proposta, todos os problemas abordados também foram resolvidos pela formulação com Dupla Reciprocidade. O custo computacional dispendido para cada uma dessas formulações também foi comparado. Para ambas as formulações também foram testados esquemas de ajuste da interpolação realizada, visando avaliar seus efeitos na precisão dos resultados e também propositando obter economia computacional em futuras aplicações em simulações na área de propagações de ondas / This research was based on the use of radial basis functions to generate a new integral formulation that interpolates directly the domain action, related to the inhomogeneous term of the governing differential equation, using the Boundary Element Method (BEM). The use of primitive functions of the original interpolation functions in the kernel of the inhomogeneous integral is proposed, allowing its transformation into a boundary integral, thus avoiding the domain discretization through cells, similar to that conducted in the Dual Reciprocity. To better evaluation of the capability of the proposed formulation, the numerical tests presented only solved problems governed by the Poisson Equation. Test problems chosen have known analytical solution, which allowed a better evaluation of the numerical accuracy. To better check the efficiency of the proposed formulation, all the problems were also solved by the Dual Reciprocity Boundary Element Formulation. The computational cost expended for each of these formulations was also compared. Fitting interpolation schemes for both formulations were also tested in order to evaluate their effects on the accuracy of the results and also looking for economy in future computational applications related to wave propagation problems

This dissertation examines the nature of jet-boundary interactions and the role of barotropic instabilities in the Brazil Current system while still attached to the western boundary. The motivation was due to the frequent observation of bipolar features associated to the Brazil Current flow south of Cape Frio (RJ) and off Santos Bight (22º-28ºS). Such observations were mainly obtained from sea surface temperature images, and also from few \"in situ\" hydrographic and direct velocity measurements data. Therefore, our main focus is on the formation of bipolar features associated with the flow, and the main hypothesis is that barotropic instability is responsible for the generation of such dipoles along the western boundary current jet, while it has to deal with topographic variations along its path poleward. We address the system dynamics from a semi-theoretical perspective, and also through the application of numerical modeling on a process study approach. So it includes semi-theoretical studies of jet-lateral boundary interactions in idealized configurations relevant to the Brazil Current system off Cape Frio, considering the quasi-geostrophic theory as an appropriate approximation of the system we want to investigate. Additionally, numerical modeling is used through the construction of idealized scenarios where we simplify the physics in order to isolate the dynamical process of interest. All the dynamical analyzes were initialized from synoptic hydrographic data set which comprised the Brazil Current system off Cape Frio region. We validated the quasi-geostrophic theory we assumed and obtained the kinematics characteristics of the jet. In the dynamical analyses, we started the investigation from the simplest framework applied here, where we evaluate a piecewise constant potential vorticity field in a quasi-geostrophic contour dynamics model. It dealt with a meridionally-oriented jet flowing southward along a straight western boundary. Next, in a second model, we added more complexity in the system, idealizing western boundary coastline scenarios considering a quasi-continuous potential vorticity field in a quasi-geostrophic numerical model. Finally, in a third model, we constructed a more complex scenario for the Brazil Current jet that incorporated real topography and stratification of the water column in a primitive equation numerical model. Among our main finds, we verified that purely barotropic instability is able to generate vortex dipoles. Moreover, variations in the western boundary can indeed trigger perturbations in the jet and dipoles form. Therefore, sites with abrupt change in bathymetry and coastline orientation are preferred to the formation of the BC vortical dipoles. Vortex streets and instability trains can also develop downstream of such locations as consequence of perturbing a potentially batropically unstable jet. Hence, the horizontal shear is key to the generation of vortex dipoles. A weak shear does not allow these features to form, instead, frontal eddies are possibly generated, with the predominance of anticyclones. Although the three different model strategies applied here differ in dynamical configurations and approximations, they still kept nearly the same regime for the Brazil Current vortex-dipole formation. Topographic variations such as those associated with the change of coastline orientation near Cape Frio (23ºS) account for those vortex dipoles and streets be dominantly observed within the Santos Bight. / A presente tese examina a natureza das interações entre jato e contorno e o papel de instabilidades barotrópicas no sistema Corrente do Brasil (CB) quando este ainda se encontra fluindo junto à margem continental oeste. A motivação se deu através da frequente observação de feições bipolares associadas ao escoamento da CB ao sul de Cabo Frio (RJ) e ao largo da Bacia de Santos (entre 22º e 28ºS). Tais observações se devem principalmente a imagens termais de temperatura da superfície do mar, a algumas raras observações \"in situ\" através de dados hidrográficos e medições diretas de velocidade. Assim, o principal foco é na formação de feições bipolares associadas ao fluxo médio, sendo a principal hipótese a de que instabilidade barotrópica é responsável pela formação destes dipolos ao longo da corrente quando esta tem que lidar com variações da topografia. Este estudo tenta abordar os problemas em uma perspectiva semi-teórica, e também através de modelagem numérica em uma abordagem de estudos de processo. Assim, inclui estudos semi-teóricos em configurações idealizadas relevantes para o sistema Corrente do Brasil ao largo de Cabo Frio, considerando a teoria quase-geostrófica como a aproximação apropriada para a dinâmica do sistema que queremos avaliar. Adicionalmente, a modelagem numérica é usada através da construção de cenários idealizados onde simplificamos a física a fim de isolar os processos que queremos investigar. Toda a análise dinâmica partiu de um conjunto de dados que compreendeu o sistema Corrente do Brasil ao largo de Cabo Frio. Validamos a teoria quase-geostrófica que estamos considerando e obtivemos as características cinemáticas do jato. Nas análises dinâmicas, começamos a investigação do problema partindo do cenário mais simples utilizado aqui, onde consideramos um campo de vorticidade potencial discretizado em camadas horizontais em um modelo quase-geostrófico de dinâmica de contornos. O modelo incorporou a presença de uma linha de costa retilínea orientada meridionalmente no contorno oeste. Posteriormente, incluímos mais complexidade no sistema, idealizando diferentes cenários de linha de costa e considerando um campo de vorticidade potencial quase-contínuo em um modelo numérico quase-geostrófico. Finalmente, construímos um cenário ainda mais complexo para a Corrente do Brasil, o qual incorporou a topografia real da região e a estratificação da coluna de água em um modelo numérico de equações primitivas. Dentre as principais conclusões, pudemos comprovar que instabilidade barotrópica pode promover a formação de dipolos. Além disso, variações no contorno podem consequentemente ser gatilhos para gerar perturbações no jato e dipolos se formam. Assim, locais de mudança abrupta de batimetria e orientação de linha de costa são preferidos para a formação de dipolos vorticais. Rua de vórtices e trens de instabilidade também podem se desenvolver à juzante de tais locais como consequência de se perturbar um jato potencialmente barotropicamente instável. Consequentemente, um cisalhamento horizontal é chave para a geração de dipolos vorticais. Um cisalhamento relativamente fraco não permite a formação de tais feições, e em vez disso, vórtices frontais são possivelmente gerados, com a predominância de anticiclones. Apesar de as três diferentes estratégias aplicadas aqui diferirem em termos de configurações dinâmicas e aproximações, estas ainda mantiveram aproximadamente o mesmo regime para a formação de dipolos. Variações na topografia tais como aquelas associadas com mudanças de orientação de costa próximo a Cabo Frio (23ºS) contam com o fato de dipolos vorticais e rua de vórtices serem frequentemente observados dentro da Bacia de Santos.

Submitted by Erika Demachki (erikademachki@gmail.com) on 2014-09-24T20:35:00Z No. of bitstreams: 2 Santos, Claudia Gomes de Oliveira - Dissertação - 2013.pdf: 1950939 bytes, checksum: 050c57553672656134c6b1264cb562a6 (MD5) license_rdf: 23148 bytes, checksum: 9da0b6dfac957114c6a7714714b86306 (MD5) / Approved for entry into archive by Jaqueline Silva (jtas29@gmail.com) on 2014-09-24T20:42:50Z (GMT) No. of bitstreams: 2 Santos, Claudia Gomes de Oliveira - Dissertação - 2013.pdf: 1950939 bytes, checksum: 050c57553672656134c6b1264cb562a6 (MD5) license_rdf: 23148 bytes, checksum: 9da0b6dfac957114c6a7714714b86306 (MD5) / Made available in DSpace on 2014-09-24T20:42:50Z (GMT). No. of bitstreams: 2 Santos, Claudia Gomes de Oliveira - Dissertação - 2013.pdf: 1950939 bytes, checksum: 050c57553672656134c6b1264cb562a6 (MD5) license_rdf: 23148 bytes, checksum: 9da0b6dfac957114c6a7714714b86306 (MD5) Previous issue date: 2013-07-31 / Coordenação de Aperfeiçoamento de Pessoal de Nível Superior - CAPES / The Boundary Element Method (BEM) has been successfully employed in the analysis of various engineering problems. The BEM consists in a mathematical modeling, for a numerical solution of a system of integral equations, and in their cores may appear singularities. This paper presents the Classical and Hypersingular formulation of the Boundary Element Method for dimensional elastostatic problems with smooth boundary geometry. The improper integrals arising from the singularities of the core in the hypersingular formulation are treated by Hadamard finite parts. In the discretization process two types of interpolation are used, one traditional and the other special. Traditional interpolation is used in all bondary elements that have no point , special interpolation ensures the continuity of the tangential derivative of displacements on the element that contains the point . To accomplish this, a theoretical mathematics study of related topics was performed. The hypersingular formulation developed in this work was implemented through the Intel Visual Fortran compiler. Some problems were analyzed and the obtained results were compared with those of analytical solution or through the Finite Element Method. The results achieved were satisfactory validating the proposed formulation / O Método dos Elementos de Contorno (MEC) vem sendo empregado com sucesso na análise de diversos problemas de engenharia. O MEC consisti em uma modelagem matemática, para resolução numérica de um sistema de equações integrais, e que em seus núcleos podem aparecer singularidades. Nesse trabalho apresenta a formulação Clássica e Hipersingular do Método dos Elementos de Contorno para problemas de elastostática bidimensional com geometria de contornos não suaves. As integrais impróprias que surgem da singularidade do núcleo na formulação hipersingular são tratados por partes finitas de Hadamard. No processo de discretização utiliza-se de dois tipos de interpolação, uma tradicional e outra especial. A interpolação tradicional é utilizada em todos os elementos de contorno que não tem o ponto , a interpolação especial garante a continuidade da derivada tangencial dos deslocamentos no elemento que contém o ponto . Para a realização deste, foi realizado um estudo teórico-matemático dos tópicos afins. Implementou-se a formulação hipersingular desenvolvidas no trabalho através do compilador Intel Visual FORTRAN. Foram analisados alguns problemas e os resultados obtidos comparados àqueles de solução analítica ou através do Método dos Elementos Finitos. Os resultados alcançados mostraram-se satisfatórios validando a formulação proposta.

A Corrente de Contorno Oeste Profunda (CCP) é um dos principais componentes da Célula de Revolvimento Meridional do Atlântico e, portanto, é determinante para a manutenção do clima global. Ela é gerada nas regiões subpolares do Atlântico Norte e transporta as águas formadas neste sítio para o Hemisfério Sul. A sua trajetória bem como a sua variabilidade têm sido estudada ao longo dos tempos, mas informações acerva desta corrente no Atlântico Sul ainda são escassas. Estudos recentes mostram que ao sul da latitude de 8ºS a CCP se quebra em vórtices anticiclônicos e se propagam para sul. Reunimos um amplo conjunto de dados observacionais de cruzeiros oceanográficos e os utilizamos para observar estas feições ao largo da costa leste brasileira. Utilizamos o cálculo geostrófico referenciado em 4000 dbar a partir dos dados hidrográficos destes. Os resultados do cálculo geostrófico confirmam os achados na literatura, identificando estruturas vorticais com diâmetros entre 162 e 220 m e velocidades azimutais máximas ultrapassando 25 cm s-1. Um fundeio correntográfico instalado nas imediações da latitude 11ºS pela Universidade de Kiel no inicio do século foi utilizado para avaliar a regularidade e sazonalidade na formação dos anticiclones. Foi conferido que a formação de vórtices possui uma modulação sazonal, ocorrendo predominantemente em períodos que se observa uma maior intensidade da CCP, ou seja, de abril a setembro. Os resultados de simulação global do modelo HYCOM nos permitiram analisar a região de estudo de forma mais ampla, tanto espacial como temporalmente. A análise dos resultados de modelo sugere que a ocorrência de vórtices anticiclônicos da CCP entre 11ºS e 20ºS é bem regular e se propaga com velocidade velocidade de translação média de 5 cm s-1, porém é severamente interrompido ao se deparar com a complexa topografia ao sul de 17ºS. Entre os obstáculos topográficos, podemos destacar a Cadeia de Abrolhos (principalmente), o Monte Submarino Hot Spur e a Cadeia Vitoria Trindade, sendo que este último deflete para leste o escoamento associado provavelmente a uma CCP restabelecida. Embora os anticiclones da CCP tenham sido descritos há mais de uma década, acreditamos ser este o trabalho que primeiro mostra o destino e a destruição dessas feições ao largo da margem continental leste brasileira. / The Deep Western Boundary Current (DWBC) is one of the main components of the Atlantic Meridional Overturning Circulation and is therefor determinant in the maintainance of the global climate. This current has its origins in the subpolar region of the North Atlantic and transports the water mass formed at deep levels to the Southern Hemisphere. Its path as well as its variability have been studied in the last twenty years but information about the DWBC in the South Atlantic are still sparse. Recent investigations showed that the DWBC breaks into large anticyclonic eddies at around 8ºS. In this wrok, we gathered an ample observation data set of oceanographic cruises and used them to describe the vortical rings off the Eastern Brazilian coast. We emplyed the classic dynamics method referenced to 4,000 dbar to infer geostrophic velocity patterns from the data set. We indentify vortical structures of typically 162-220 km of diameter and azimuthal velocities higher than 25 cm s-1. A currentmeter mooring deployed by the University of Kiel in the begining of the century was used to investigate regularity and sazonality of the eddy formation and their passing of the 11ºS paralell. We found by inspecting the 5-year time series that there are about 3,7 events/year and a seasonal modulation dictated by the strength of the DWBC. More rings are shed during the April-September season when the DWBC is more intense and transports more. We also analyzed the output of a assimilative global simulation by the Hybrid Ocean Circulation Model (HYCOM) Consortium after validating it with our observational data set. Our results indicate that the anticyclones propagate at a characteristic speed of about 5cm s-1 until interacts with complex seamount topography at 17ºS. The topography is formed by th Abrolhos Ridge (17ºS), the Seamount Hot Sput (18ºS) and the Vitória-Trindade Ridge (20ºS). Moreover, it seems that the vast majority of the vortical features are destroted on the northernmost obstacle - the Abrolhos Ridge. Only a few survive and are destructed south. It seems from the analysis that a recoalesced DWBC existis nearby 20ºS, where it veers east and flows toward the interior of the Atlantic Ocean basin. We should emphasize that despite the DWBC anticyclones original description is more than a decade old, this is the first work to describe the fate of these large rings off the Eastern Brazil Continental Margin.

This work is devoted to the numerical analysis of saturated porous media, taking into account the damage phenomenon on the solid skeleton. The porous media is taken into poroelastic framework, in full-saturated condition, based on the Biot\'s Theory. A scalar damage model is assumed for this analysis. An implicit Boundary Element Method (BEM) formulation, based on time-independent fundamental solutions, is developed and implemented to couple the fluid flow and the elasto-damage problems. The integration over boundary elements is evaluated by using a numerical Gauss procedure. A semi-analytical scheme for the case of triangular domain cells is followed to carry out the relevant domain integrals. The non-linear system is solved by a Newton-Raphson procedure. Numerical examples are presented, in order to validate the implemented formulation and to illustrate its efficiency. / Este trabalho trata da análise numérica de meios porosos saturados, considerando danificação na matriz sólida. O meio poroso é admitido em regime poroelástico, em condição saturada, com base na teoria de Biot. Um modelo de dano escalar é empregado nesta análise. Uma formulação implícita do Método dos Elementos de Contorno (MEC), baseada em soluções fundamentais independentes do tempo, é desenvolvida e implementada de forma a acoplar os problemas de difusão de fluido e de elasto-dano. A integração sobre os elementos de contorno é feita através da quadratura de Gauss. Um esquema semi-analítico é aplicado sobre células triangulares para avaliar as integrais de domínio do problema. A solução do sistema não linear é obtida através de um procedimento do tipo Newton-Raphson. Apresentam-se exemplos numéricos a fim de validar a formulação implementada e demonstrar sua eficiência.

Neste trabalho apresenta-se um modelo numérico para a análise de problemas tridimensionais envolvendo a interação mecânica estaca-solo, acoplando-se o Método dos Elementos de Contorno (MEC) ao Método dos Elementos Finitos (MEF). O solo é modelado com o MEC utilizando-se as soluções fundamentais de Mindlin, assumindo um meio semi-infinito, homogêneo, isotrópico e elástico-linear. As estacas, modeladas com o MEF, consistem em um elemento único, com quatro nós e 14 parâmetros nodais (três deslocamentos em cada nó e mais duas rotações no topo da estaca). Cada uma das estacas é levada em consideração no MEC como uma linha de carga. Considera-se o escorregamento das estacas em relação ao maciço, empregando modelos de aderência para a definição da evolução das tensões tangenciais ao longo do comprimento das estacas. São empregados, como funções de forma, polinômios do quarto grau para os deslocamentos horizontais, cúbicos para os deslocamentos verticais e tensões horizontais ao longo do fuste e quadráticos para as tensões verticais do fuste e escorregamento. A reação da ponta da estaca é calculada assumindo tensão constante na base. / This work presents a method for tri-dimensional pile-soil interaction problems, by coupling the Boundary Element Method (BEM) to the Finite Element Method (FEM). The soil is modeled with BEM, using the Mindlins fundamental solutions, supposing a semi-infinite, homogeneous, isotropic, elastic and linear space. Piles are modeled with FEM and are represented by one element with four nodes and 14 nodal parameters (three displacements in each node and two rotations at the top node). Each pile is represented in BEM as a line load. The pile slip is considered using adherence models to evaluate the evolution of shaft tractions. There are employed fourth grade polynomial shape functions for horizontal displacements, cubic polynomial functions for vertical displacements and horizontal tractions along shaft and quadratic polynomial functions for vertical tractions and slip. Tip reaction is calculated supposing constant traction at the base.

O trabalho tem como objetivo trazer contribuições à simulação numérica pelo método dos elementos de contorno (MEC) de formação e propagação de descontinuidades no campo de deslocamentos (descontinuidades fortes) em sólidos bidimensionais. A formação de descontinuidades fortes caracteriza o processo de falha material, que pode estar associado ao fraturamento em materiais quase frágeis ou a superfícies de deslizamentos de materiais dúcteis. Apresenta-se uma formulação do MEC baseada na incorporação de interfaces de descontinuidade no interior de células internas, que possibilita propagação arbitrária de descontinuidades usando uma malha de células internas fixa, definida antes da análise. Comparam-se diferentes alternativas provenientes do relaxamento dos requisitos de consistência estática e analisa-se a influência do alinhamento da malha. Apresenta-se também um possível esquema de construção adaptativa de células internas com interface incorporada para capturar a trajetória arbitrária da descontinuidade que se propaga durante o processo de carregamento. Este esquema visa aumentar a robustez e reduzir o esforço computacional. As características geométricas das células internas geradas são estabelecidas em função da orientação da descontinuidade fornecida pelo critério de falha, de maneira a proporcionar melhor eficiência numérica. Os estudos são levados a cabo através da simulação numérica de testes experimentais colhidos da literatura. / This work has the objective of bringing contributions to the numeric simulation using the boundary elements method (BEM) to model the initiation and propagation of strong discontinuities in the displacement field in bidimensional solids. The initiation process of strong discontinuities characterizes the failure process of material, which can be associated with the fracture of quasi-brittle materials and slip lines in ductile materials such as metals. The effect of the displacement jump of a discontinuity interface embedded in an internal cell is provided by an equivalent strain field over the cell. This model allows the study of arbitrary crack growth using a fixed mesh defined before the analysis. The dissipative process in the cell interface is described by an isotropic damage model in the continuum approach of strong discontinues. Alternatives that come from relaxing the static consistencies and the influence of the mesh alignment are analyzed. An adaptative algorithm for internal cells creation is also presented to capture the path of the crack growth during the loading process. This algorithm intends to overcome some convergence problems found in models with predefined meshes and also to reduce the computational efforts. The geometric characteristics of the generated internal cells are defined using the crack orientation, given by the failure criterion, to provide a better numerical efficiency. The results obtained with the proposed formulation are compared with the ones obtained with other numerical methods and also from experiments.

Orientador: German Jesus Lozada Cruz / Banca: Luiz Augusto Fernandes de Oliveira / Banca: José Marcio Machado / Resumo: O objetivo deste trabalho é estudar problemas de valor de contorno do tipo {ÿ + f(t) =0 y(0)=0˙ y(1)= ky(η), (1) onde η ∈ (0, 1), k ∈ R e f ∈C([0, 1],R). Para antingirmos nosso objetivo usamosas funções de Green G(t,s)que nos permitem escrever a solução do problema(1)na seguinte forma: w(t)= ∫ 1 0 G(t,s)f(s)ds. Usando esta solução, investigamos através do ponto fixo de Schauder a solvabilidade do problema não linear { y + f(t,y)=0 y(0)=0˙ y(1)= ky(η). / Abstract: The main goal of this work is study the following boundary value problems {ÿ + f(t) = 0 =0 y(0)=0˙ y(1)= ky(η), (1), where η ∈ (0, 1), k ∈ R e f ∈C([0, 1],R). To achieve our goal we use the Green's function G(t,s) which allow us to write the solution of the problem (2) in the form: w(t)= ∫ 1 0 G(t,s)f(s)ds. Using this solution and the Schauder point theory, also we study the solvability of a nonlinear problem { y + f(t,y)=0 y(0)=0˙ y(1)= ky(η). / Mestre

Tese (Doutorado) - Universidade Federal de Santa Catarina, Centro Tecnológico. / Made available in DSpace on 2012-10-17T13:30:56Z (GMT). No. of bitstreams: 0Bitstream added on 2014-09-25T19:15:18Z : No. of bitstreams: 1 161846.pdf: 91791177 bytes, checksum: f51a3f5dd4b7584432ec2bf16b4b082c (MD5) / Apresenta-se uma metodologia de modelagem global para caixas de engrenagem, tendo por base os métodos de elementos finitos e elementos de contorno, cujo enfoque principal é para as vibrações e ruídos transportados pela estrutura e que têm como origem a excitação provocada pelo engrenamento, de modo a permitir a predição de parâmetros vibro-acústicos. Assim, são apresentadas as formulações matemáticas que permitem a determinação de matrizes de rigidez generalizadas, para a representação de mancais de elementos rolantes e do engrenamento; o desenvolvimento de um modelo numérico de sistema de eixos engrenados, que permite estimar as reações nos mancais devido ao engrenamento e a influência da modificação do perfil dos dentes das engrenagens; a concepção e criação de uma malha de elementos finitos para representação do alojamento da caixa de engrenagem, bem como a malha acústica correspondente usada para a estimativa de parâmetros acústicos pelo método de elementos de contorno. A metodologia é aplicada a uma caixa de engrenagem de uso veicular real e o comportamento vibro-acústico é predito para condições operacionais distintas.