Este trabalho trata da formulação do Método dos Elementos de Contorno para o problema de flexão de placas através da teoria clássica de Kirchhoff. Ênfase especial é dada a placas com variação de espessura uma vez que este tema é muito pouco abordado. O estudo de flexão de placas com variação de espessura é exercida em vários campos, tais como engenharia civil, engenharia aeroespacial e projeto de máquinas. Adotou-se uma variação linear da rigidez, resultando nas equações integrais de deslocamentos com termos de domínios, que serão tratados por discretização do domínio em células. / This work deals with the formulation of the Boundary Element Method applied to plate bending problem using the Kirchhoff\'s theories. Special emphasis is given to plates with varying thickness since this subject is not much tackled. The study of the bending of plates of variable thickness is pursued in various fields, such as civil engineering, aerospace engineering, and the design of machines. It was adopted linearly varying rigidity, giving integral equations of displacements with domain terms, that will be treated by domain discretization into cells.

Este trabalho trata da formulação do método dos elementos de contorno para a análise não linear de pavimentos de concreto armado. A teoria utilizada é a de Reissner, que mostrou-se eficiente tanto para placas esbeltas quanto para as moderadamente espessas. Considera-se a ocorrência de cargas concentradas, distribuídas em sub-regiões da placa e em linha. Admite-se também a possibilidade de um campo de momentos iniciais, que viabiliza o estudo da não linearidade física nos problemas. Foram utilizados campos de momentos iniciais aplicados apenas em pontos internos ao domínio. As integrais que envolvem as células de domínio foram modificadas, eliminando-se os núcleos complexos e as aproximações através de séries. Foi desenvolvida uma formulação para a análise de placas vinculadas a estruturas quaisquer em seu domínio, com o uso de cargas aplicadas incógnitas atuando como enrijecedores. O acoplamento MEC/MEF foi empregado utilizando-se modelos simples, porém robustos. O sistema de equações algébricas foi otimizado com a utilização da técnica dos mínimos quadrados. O concreto foi modelado adotando-se o modelo de dano de Mazars; para as armaduras um modelo elastoplástico uniaxial com endurecimento isótropo. A análise não linear do problema é efetuada utilizando-se procedimento incremental-iterativo. São apresentados alguns exemplos simples que mostram a precisão da técnica usada. / This work refers to the formulation of the boundary element method for non-linear analysis of building floor structures. The plate bending theory adopted to develop the work wad due to Reissner, which has demonstrated to be efficient for thick, moderated thick and thin plates. The kinds of load applied on the plate medium surface have been taken into account: concentrated loading, distributed over sub-domains; distributed along internal lines. The presence of initial moment fields convenient to model temperature effects and to be used to build up non-linear solutions has also been considered in the formulation. The domain integrals containing complex kernels to take into account the initial moment field influences were modified by introducing their primitive functions, avoiding therefore using series expansions. To integrate the initial moments fields only approximations based on internal nodal points were used. The resulting cell integrals have been transformed to the cell boundary which results into regular integral only. A boundary element formulation to treat structural system defined by combining plates with other structural element was developed, using interface force as unknowns. The BEM/FEM coupling developed to treat this case is simple but robust; only displacements have been coupled avoiding important singularities that may happen when coupling rotations. The resulting system of algebraic equations has been regularized by using the least square method. The concrete material was modeled by using the Mazar\'s damage model, while the steel reinforcement was assumed to behave as elastoplastic material with isotropic hardening. Finally, some examples are shown to illustrate the accuracy of the presented formulation and the numerical schemes proposed in this work.

Um dos principais problemas que o Método de Elementos de Contorno (MEC) apresenta encontra-se na avaliação de integrais singulares e quase-singulares que envolvem as soluções fundamentais de Kelvin em deslocamento e força. O processo de integração numérica em MEC tem sido o objetivo de inúmeras pesquisas, pois dele depende a qualidade das respostas quando se deseja obter uma excelente precisão numérica em uma análise. Este trabalho apresenta uma nova proposta de integração numérica para análise tridimensional com MEC. Esta técnica possui três características importantes. A primeira é a determinação da parcela efetiva de singularidade que ocorre na função raio, distância entre o ponto fonte e o elemento de contorno bidimensional. A correta obtenção desta parcela permite representar sem aproximações o comportamento da singularidade da função raio, que é a verdadeira fonte de singularidade e quase-singularidade nas soluções fundamentais. A segunda característica da técnica proposta é que ela baseia-se em um Método Semi-Analítico de avaliação de integrais, onde, para cada parcela efetiva de singularidade, utiliza-se uma quadratura numérica cujos pesos específicos são calculados analiticamente. A terceira característica da técnica proposta é apresentar um tratamento unificado para todos os tipos de integrais singulares, quasesingulares e regulares. Esta técnica foi implementada na plataforma computacional desenvolvida pelo grupo GoBEM, utilizando o conceito de Programação Orientada a Objetos através da Linguagem de programação Java. Com a implementação da nova técnica de integração na plataforma computacional torna-se possível realizar o desenvolvimento de vários tipos de pesquisa sobre análise tridimensional com o MEC como, por exemplo: visualização de isosuperfícies em análise tridimensional sem discretização do domínio, automatização do cálculo elasto-plástico, modelagem de problemas geotécnicos de forma precisa, etc. Para a validação da técnica proposta três procedimentos foram considerados: análise da eficiência da parcela efetiva de singularidade, testes de convergência da integração numérica específica e exemplos numéricos utilizando o MEC em problemas de engenharia. / One of the main problems with the Boundary Elements Method (BEM) is the evaluation of singular and quasi-singular integrals due to the Kelvin\'s Fundamental Solutions in displacement and traction. Today there is an increasing body of research work that focuses on numerical evaluation of BEM integrals, for this is a crucial issue in order to achieve highly accurate results. This work presents an innovative numerical integration procedure for threedimensional analysis with BEMs. The proposed technique encompasses three important features. First, it corresponds to an accurate representation of the effective term of singularity in the radius function, which measures the distance between the source point and a twodimensional boundary element. The correct evaluation of this term leads without approximation to the actual singularity behavior of the radius function, which is the true source of singularity and quasi-singularity in the fundamental solutions. Second, the proposed technique is based on a semi-analytical procedure, i.e, for each singularity effective term it uses a quadrature scheme with specific weights analytically evaluated. Last, the proposed technique represents a unified procedure to singular, quasi-singular and regular integrals. This technique was implemented in the computational platform developed by the group GoBEM, using Object Oriented Programming and the Java programming language. The implementation of the proposed technique into this computational plataform opens new possibilities for future researches on three-dimensional BEM, e.g.: visualization of isosurfaces in three-dimensional analysis without any domain discretization, automatic elasto-plastic analysis, accurate modeling of geomechanical problems, etc. Validation of the proposed technique was carried out using three procedures: efficiency analysis of the singularity effective term, convergence tests for the specific numerical integration and numerical examples using BEM in engineering problems.

Neste trabalho a formulação direta do Método dos Elementos de Contorno é utilizada para o estudo da interação de placas com o meio contínuo. O solo, considerado como um meio contínuo, tem a sua reação representada pelo acréscimo de uma integral de domínio nas equações integrais usuais de placas. Essa integral de domínio é tratada utilizando-se células internas, o processo da reciprocidade dual e uma formulação alternativa. Inicialmente a reação do solo é aproximada utilizando-se a teoria de Winkler. Em seguida, o solo é considerado como sendo um sólido tridimensional elástico de domínio semi-infinito, e analisado pelo Método dos Elementos de Contorno utilizando-se as soluções fundamentais de Boussinesq-Cerruti e Mindlin, neste caso a interação entre a placa e o solo, é feita impondo-se o equilíbrio dos esforços e a compatibilidade de deslocamentos transversais, em todos os pontos da interface. Adotando-se um critério de plastificação simples e bilinear, é considerada a não-Iinearidade da reação do solo. Em cada caso, são apresentadas aplicações numéricas utilizando-se as três formulações, cujos resultados são comparados entre si e com valores teóricos, mostrando a eficiência das mesmas. / In this work the direct formulation of Boundary Element Method is adopted to study interaction of plates in bending with the continuum medium. The soil material, assumed as a continuum medium, applies on the plate surface what is represented in the usual integral equations by a domain integral. That domain integral is treated by approaching the subgrade reaction using cells, the dual reciprocity method or an alternative procedure. Initially, the reaction of soil is given by assuming the Winkler\'s theory. Then, the soil is assumed as a semi-infinite three-dimensional elastic solid for which the Boundary Element Method is applied using the Boussineq-Cerruti and Mindlin\'s fundamental solutions. In this case, the interaction of the plate with the soil is made by enforcing displacement compatibility and equilibrium at all interface points define by the plate surface discretization. Non-linear behaviour is also assumed to govern the interaction reaction between plate and the soil medium. For that, a simple non symmetric stress-strain curve is taken to represent the elastoplastic responses. In each case, numerical examples, using the three subgrade reaction approximations discussed here, are shown to illustrate the accuracy and efficiency of the proposed models.

O sistema de correntes de contorno oeste que flui ao largo da costa leste brasileira entre 15°S e 22°S, é composto pela Corrente do Brasil (CB) fluindo para sul, a Sub-Corrente Norte do Brasil (SNB) fluindo para norte e a Corrente de Contorno Oeste Profunda (CCP) também fluindo para sul. Vigorosos meandros e vórtices são observados à jusante dos Bancos de Abrolhos (BA) e Royal Charlotte (BRC) e da Cadeia Vitória-Trindade. O objetivo central desta dissertação é o estudo da estabilidade deste sistema de correntes, aqui denomidado Sistema CB-SNB-CCP, utilizando o modelo oceânico da Universidade de Princeton (POM) em um cenário idealizado. Buscamos responder o quão instável é este sistema e quais seriam as características das ondas instáveis geradas a fim de contribuir para o conhecimento acerca da rica dinâmica de meso-escala observada nesta região. Objetivamos também, elucidar o papel dos BA e BRC nesta dinâmica. Embasados na alta baroclinicidade do sistema de correntes de contorno ao largo do sudeste brasileiro, optamos por representar o escoamento CB-SNB-CCP através de um modelo paramétrico do campo de massa, calibrado com os dados hidrográficos oriundos dos Cruzeiros Abrolhos [Silveira et al., 2006]. Mantendo o caráter idealizado do estudo, também empregamos topografia de fundo analítica, onde representamos o talude da região por uma função tangente hiporbólica. O BA e o BRC foram aproximados através de funções gaussianas devidademente ajustadas aos contornos da isóbata de 80 m extraídas do conjunto ETOPO 2. Para identificar os mecanismos de crescimento das possíveis ondas instáveis, calculou-se o balanço de energia das simulações realizadas de acordo com o método de [Xue & Bane, 1997]. Os resultados de três experimentos numéricos realizados sugerem que o sistema CB-SNB-CCP é instável. Ciclones quase-estacionários do lado costeiro da CB surgem como o principal modo de variabilidade desta corrente. De acordo com a análise do balanço energético, o crescimento dessas feições resulta primariamente de instabilidade baroclínica do escoamento. A escala horizontal típica das ondas e vórtices instáveis modelados é dada pelo raio de deformação interno, como esperado pela teoria de instabilidade baroclínica de escoamentos realisticamente estratificados. Em particular, os resultados do experimento com o BA e o BRC idealizados comprovam que estes funcionam como gatilhos para o desenvolvimento de ondas instáveis, favorecendo amplamente o crescimento das estruturas verticais. Os trens de onda instáveis quase-estacionários aqui obtidos sugerem que provavelmente o meandramento da CB observado em latitudes que se estendem até 28°S pode ser parte de um único sistema que se origina na região dos BA e BRC. / The western boundary currents system that flows off the eastern brazilian coast between 15°S e 22°S is composed by the southward-flowing Brazil Current (BC), the northward-flowing North Brazil Under Current (NBUC) and the Deep Western Boundary Current (DWBC) that flows south. Vigorous meanders and eddies are observed downstream of the Abrolhos (AB) and Royal Charlotte Banks (RCB) and the Vitória-Trindade Ridge. The main goal of this dissertation is to study the stability of this currents system, hereby named BC-NBUC-DWBC system, using the Princeton University Ocean Model (POM) in an idealized scenerio. We seek to answer how unstable is this system and what are the characteristics of the unstable waves in order to contribute to the understanding of the rich mesoescale dynamics observed in this region. We also aim to elucidate the role of the AB and of the RCB on this dynamics. Based on the high degree of baroclinicity of the western boundary currents system off the southeastern brazilian coast, we opted to represent the BC-NBUC-DWBC system through a parametric model of the mass field, calibrated with hydrographic data from the Abrolhos Cruises [Silveira et al., 2006]. Maintaining the idealized character of the study, we have also employed an analytical bottom topography, in which the region\'s continental slope is approximated by hyperbolic tangent function. The AB and RCB were approximated by gaussian functions properly adjusted to the 80 m isobath extracted from the ETOPO 2 database. To identify the growth mecanism of the unstable waves, the energy budget of the simulations was calculated according to [Xue & Bane, 1997]. The results from the three experiments conducted here suggest that the BC-NBUC-DWBC system is indeed unstable. Quasi-stationary cyclones in the coastal side of BC arise as the main mode of variability of this current. According to the energy budget analisys, the growth of these features results primarily from baroclinic instability of BC-NBUC-DWBC flow. The typical horizontal scale of the modeled unstable waves and eddies is given by the internal radius of deformation, as expected by baroclinic instability theory of realistically stratified flows. The quasi-stationary unstable wave trains modeled in the present study also suggest that the meandering of the BC observed down to 28°S are probably part of a single system that originates at the AB and RCB region.

Variabilidades de longo período têm sido detectadas no padrão de ventos em larga escala no hemisfério sul, com base em dados observacionais, reanálises e modelos. Estudos recentes têm sugerido algumas respostas oceânicas a mudanças no padrão atmosférico, desde a migração da Frente Subtropical, ao aumento do Vazamento das Agulhas e à intensificação dos giros subtropicais (GS). Apesar da importância do Atlântico Sul (AS) na distribuição global de calor, ainda não está claro a variabilidade de sua circulação e suas respostas às mudanças no campo do vento em larga escala. O presente estudo tem como objetivo investigar, usando resultados de uma simulação com o modelo numérico HYCOM (Hybrid Coordinate Ocean Model), a variabilidade do transporte de volume da Corrente do Brasil (CB) em sua região central (∼30ºS) frente ao campo de vento da Reanálise I do NCEP/NCAR (National Centers for Environmental Prediction/National Centers for Atmospheric Research). Os produtos do NCEP/NCAR apresentam intensificação e deslocamento em direção ao polo do cinturão de ventos de oeste no hemisfério sul (HS), o que ocasiona alteração do rotacional da tensão de cisalhamento do vento, ao longo de toda a bacia do AS, e a migração para sul das linhas de rotacional nulo e máximo em latitudes subtropicais. Mesmo com a intensificação do rotacional médio sobre o AS, a série computada do fluxo de retorno do transporte integrado de Sverdrup em 30ºS, através dos dados da reanálise, demonstram uma diminuição entre os anos de 1960 e 2010. Este comportamento é coerente com a tendência à redução do rotacional mediado zonalmente ao longo da mesma latitude. O transporte da CB computado com resultados do HYCOM apresenta uma tendência próxima à encontrada para o transporte de Sverdrup integrado, de cerca de 0,1 Sv por década, e estas apresentam correlação máxima de 0,6 com defasagem de 2 anos. Ao longo de 30°S, a leste da CB modelada, é encontrado um fluxo para norte identificado como o retorno de uma célula de recirculação, com transporte médio de 4,25±2,87Sv. Este fluxo também apresenta tendência à redução, a qual poderia estar relacionada com uma migração para sul de toda a estrutura da recirculação, acompanhando o deslocamento da confluência Brasil-Malvinas e da Frente Subtropical. O transporte residual na borda oeste (CB+Recirculação) apresenta uma intensificação do fluxo para sul, a qual poderia indicar um fortalecimento do GS, o que estaria em concordância com a intensificação e maior abrangência da Alta Subtropical observada no campo de Reanálise I. No entanto não foi encontrada correlação significativa entre o campo de vento e a variabilidade do transporte residual ao longo da fronteira oeste. / Long term variability has been detected in the South Hemisphere large-scale wind pattern, using models, reanalysis and observations. Recent studies have suggested an oceanic response to changes in the atmospheric pattern, such as the Subtropical Front migration, increase in the Agulhas leakage and intensification of the subtropical gyre (SG). Despite the importance of the South Atlantic in the global heat distribution, the variability of its circulation and response to changes in the large-scale wind field remains unclear. The current study aims to investigate, from a numerical simulation using HYCOM (Hybrid Coordinate Ocean Model), the variability of the Brazilian Current (BC) volume transport at its central region (∼30ºS) according to the wind field from Reanalysis I - NCEP/NCAR (National Center for Environmental Prediction/National Center for Atmospheric Research). The NCEP/NCAR results show an intensification and a poleward shift of the Southern Hemisphere westerly winds, leading changes in the wind curl along the SA, and a southern migration of the zero and maximum curl lines in subtropical latitudes. Even with the intensification of the mean wind curl over the SA, the time series of integrated Sverdrup return flow transport at 30ºS show a decrease between 1960 and 2010. This pattern is in agreement with the reduction trend in the zonal mean curl along the same latitude. The computed BC transport by HYCOM results present a similar trend to the integrated Sverdup transport (approximately 0.1 Sv per decade), with a 0.6 maximum correlation with a 2 year delay. Along the 30ºS, close to the western boundary, a north flow was identified as the return of a recirculation cell, with a 4,25±2,87 Sv mean transport. This flow also presents a reduction trend, which could be related to southern migration of the entire recirculation structure, following the Brazil-Malvinas confluence and Subtropical Front displacement. The residual transport at the western boundary (BC+Recirculation) shows an increase in the south flow that could indicates an increase of the SG, which is in agreement with the intensification and expansion of the Subtropical Ridge observed in the Reanalysis I. However, no significant correlation between the wind field and residual transport variability along the western boundary was found.

Desenvolvem-se a formulação do Método dos Elementos de Contorno e correspondente algoritmo (para implementação em microcomputador) para a análise de propagação de fraturas em domínios bidimensionais. São utilizados elementos lineares isoparamétricos, tanto para discretizar o contorno quanto para simular a fratura. Os elementos de fratura são descontínuos. A formulação é baseada em equações integrais de tensões e de deslocamentos, onde o termo que considera tensões iniciais concentradas na linha de fratura é formulado a partir da definição de dipolos. O critério adotado é o modelo de fratura coesiva. Os termos singulares e hiper-singulares da formulação são tratados analiticamente e os termos quase-singulares são calculados através de um esquema numérico baseado na utilização de sub-elementos. Os valores dos dipolos são estimados ponto a ponto. Ao longo das fraturas, o valor máximo da tensão normal de tração permite definir novos elementos. As tensões de cisalhamento são removidas para manter a direção principal durante o processo. / The Boundary Element Method formulation and corresponding algorithrn (for microcomputer implementation) are developed for crack growth analysis in two-dimensional domains. Linear isoparametric elements are used to discretize both boundary and crack path. Fracture elements are assumed to be discontinuous. The formulation is based on stress and displacement integral equations, where the term that takes into account initial stresses concentrated along fracture line is formulated from dipoles definition. The coesive fracture rnodel is the criterium adopted. Singular and hipersingular formulation terms are anallitically treated and quasi-singular terms are computed by a numerical scheme based on element subdivision. Dipole values are estirnated point by point. Along fractures, the maximum normal tensile strenght is used to define new elements. Shear stresses are also removed to maintain the principal direction during the process.

Submitted by Guilherme Avino Ohland null (guilherme01261@aluno.feis.unesp.br) on 2017-03-04T18:34:57Z No. of bitstreams: 1 DISSERTAÇÃO - GUILHERME AVINO OHLAND.pdf: 5200722 bytes, checksum: 2b9af38c712c37f2ad2bc2bd13a092c5 (MD5) / Approved for entry into archive by Juliano Benedito Ferreira (julianoferreira@reitoria.unesp.br) on 2017-03-09T13:21:06Z (GMT) No. of bitstreams: 1 ohland_ga_me_ilha.pdf: 5200722 bytes, checksum: 2b9af38c712c37f2ad2bc2bd13a092c5 (MD5) / Made available in DSpace on 2017-03-09T13:21:06Z (GMT). No. of bitstreams: 1 ohland_ga_me_ilha.pdf: 5200722 bytes, checksum: 2b9af38c712c37f2ad2bc2bd13a092c5 (MD5) Previous issue date: 2017-02-07 / Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq) / Este trabalho tem como objetivo apresentar uma formulação do Método dos Elementos de Contorno (MEC) para análises de microestruturas heterogêneas, onde dentro da matriz podem ser definidos vazios ou inclusões com diferentes propriedades elásticas. A microestrutura é modelada por uma chapa em sub-regiões, onde diferentes valores de coeficientes de Poisson e módulo de Young podem ser definidos para cada sub-região. Para resolver as integrais de domínio escritas em termos de deslocamentos no plano, o domínio da matriz e das inclusões são discretizadas em células, onde os deslocamentos têm que ser aproximados. Assim, neste modelo, além de valores no contorno para deslocamentos e forças no plano, valores nodais de deslocamentos no plano são definidos também no domínio. Então, adotando-se técnicas de homogeneização, os valores homogeneizados para o tensor constitutivo e o tensor das tensões são calculados. A formulação é proposta dentro do contexto de análise em multi-escala de estruturas, onde a microestrutura do material é denominada de EVR (Elemento de Volume Representativo), sendo seu problema de equilíbrio definido em termos de flutuação dos deslocamentos. Neste trabalho será adotado comportamento elástico linear para os diferentes materiais (ou fases) do EVR, porém tal formulação pode ser facilmente estendida no futuro a fim de considerar deformações residuais. Nos exemplos numéricos, os resultados são comparados com uma formulação desenvolvida com o Método dos Elementos Finitos (MEF), a fim de validar o modelo proposto. / A formulation of the boundary element method (BEM) to perform elastic analysis of heterogeneous microstructures is presented. The microstructure is modelled as a zoned plate where voids or inclusions can be considered inside a matrix. Thus, each sub-region represents either the matrix or an inclusion, where different Poisson’s ratio and Young’s modulus can be defined. In the proposed model domain integrals in terms of in-plane displacements arise in the formulation, which are solved by discretizing the sub-regions into cells where the displacements are approximated. Thus, besides the boundary values for in-plane displacements and tractions, nodal values of in-plane displacements are defined in the domain. Although the proposed model can be used to analyse the stretching problem of plates composed by different materials, in this paper the formulation is proposed in the context of a multi-scale analysis, i. e. it will be used to model the RVE (Representative Volume Element), which equilibrium equation is solved in terms of displacement fluctuations. In this paper, only elastic behaviour will be considered for matrix and inclusions, although the proposed model can be extended to consider dissipative phenomena in the microstructure. To make the micro-to-macro transition necessary in a multi-scale analysis, the homogenized values for stress and constitutive tensor have to be computed adopting homogenization techniques. In a future work this formulation will be coupled to a BEM formulation to model the macro-continuum in order to perform the multi-scale using only the Boundary Element Method. Some numerical examples of heterogeneous microstructures are presented and compared to a formulation of the Finite Element Method to show the accuracy of the proposed model. / CNPq: 130642/2015-0

Resumo: No projeto estrutural de qualquer obra de engenharia busca-se o menor custo, a alta resistência e confiabilidade. Portanto, com o avanço da técnica e do desenvolvimento dos materiais empregados, são projetadas estruturas cada vez mais esbeltas, que conduzem à necessidade de conhecimentos nas áreas de instabilidade, dinâmica e reologia. Neste trabalho são realizados estudos de problemas de segunda ordem em placas elásticas delgadas através do Método dos Elementos de Contorno (MEC), considerando-se para a montagem do sistema de equações algébricas apenas a equação integral do deslocamento transversal da placa e as variáveis de canto como incógnitas. Apresenta-se uma nova Metodologia para tratar dos efeitos de segunda ordem de placas pelo MEC. Primeiramente é feito um estudo da teoria de primeira ordem em placas e também em chapas, pois estes estão relacionados com o estudo de segunda ordem. Este trabalho trata do problema desde o início mostrando as diferenças entre a teoria de primeira ordem e a de segunda ordem, após explicitar as hipóteses e as deduções das equações que conduzem a uma teoria de segunda ordem simplificada, da qual uma das aplicações é a equação da flambagem de Von Kármán. Outra contribuição deste trabalho foi dar um tratamento novo para as integrais de domínio, cujo integrando é função da curvatura "w,ij" que passou a ser função apenas do deslocamento transversal "w", permitindo assim escrever na formulação do problema pelo MEC apenas equações integrais de "w", tanto no contorno quanto no domínio / Abstract: In structural design of any engineering project seeks to lower cost, high strength and reliability. Therefore, with the advance of technique and strength of materials used, more slender structures are designed, which lead to the need for knowledge in the areas of instability, dynamics and rheology. This work studies the problems of second order and buckling in thin elastic plates through the Boundary Element Method (BEM), considering, for mounting the system of algebraic equations, only the integral equation of the transverse displacement of the plates and unknown corner variables. It presents a new methodology to adress the effects of second order plates by the BEM. First it is studied the theory of first order in bending plates, and also in plane stress, as they are related to the study of second order. This paper addresses the problem from the beginning showing the differences between the theory of first and second order, after explaining the hypotheses and deductions of equations that lead to a simplified theory of second order, which is one of the applications of the equation Von Kármán buckling. Another contribution of this paper was to give a new treatment for the domain integrals, whose integrand is function of curvature "w,ij" which is now a function only of transversal displacement, "w", enabling to write in the problem's formulation by BEM, only the integral equations of "w", both in boundary and in the domain / Orientador: Edson Luiz França Senne / Coorientador: Webe João Mansur / Banca: Celso Pinto Morais Pereira / Banca: Fernando de Azevedo Silva / Banca: Lutgardes de Oliveira Neto / Banca: Valdir Schalgh / Doutor

Neste trabalho, estudam-se problemas bidimensionais de propagação de ondas acústicas e elásticas, no domínio da freqüência, formulados através do Método dos Elementos de Contorno. A formulação é baseada nas representações integrais das equações diferenciais que governam os fenômenos de propagação de ondas acústicas num meio fluido e de ondas elásticas numa estrutura elástica. Analisa-se também a interação entre o fluido e a estrutura com o uso de sistemas de equações acoplados. As soluções fundamentais utilizadas são expressões exatas e não há necessidade de subdivisão dos domínios em células de integração. São aplicadas técnicas de integração alternativas na escolha das equações algébricas no domínio do fluido, visando a melhora das respostas globais do conjunto. Apresentam-se ainda exemplos numéricos, com o objetivo de possibilitar a modelagem numérica de problemas de acoplamento fluido-estrutura e de radiação e difusão acústica. / In this work, acoustic and elastic wave propagation problems in 2D, in frequency domain, are studied and formulated with the Boundary Element Methods. The formulation is based on the integral representations derived from the differential equations that govern the phenomena of acoustic wave propagation in a fluid medium and elastic wave propagation inside an elastic domain. The fluid-structure interaction is also formulated by coupling appropriately the corresponding systems of equations. The fundamental solutions adopted in this work are conveniently chosen to avoid the mass integral terms in the elastic wave integral representation and the equivalent terms in the acoustic integral equation. Thus, the algebraic representations of both problems are written only in terms of boundary values. Subdivisions of the domain to perform integrals over cells are not required. In an attempt to improve the global answers of the fluid problem, several integration techniques have been experimented to build alternative algebraic matrix equations. Numerical examples are presented in order to shown the accuracy of the studied acoustic radiation and scattering problems and also to verify the proposed fluid-structure coupling.