Estuda-se uma combinação do método dos elementos finitos (MEF) com o método dos elementos de contorno (MEC) no acoplamento de uma estrutura reticulada em um domínio bidimensional. Para o caso em análise, os elementos uniaxiais são tratados através do MEF, enquanto que o MEC é utilizado na modelagem do meio contínuo que pode ser homogêneo ou não-homogêneo. Em problemas geomecânicos, é possível simular a sequência de escavação com as modificações estruturais necessárias. As equações do domínio bidimensional, com a influência das estruturas reticuladas, são agrupadas em um sistema que é resolvido pelo algoritmo proposto por Crotty. Utilizando-se o método dos elementos de contorno, a plasticidade do meio contínuo é analisada com um procedimento incremental e iterativo baseado no processo das tensões iniciais. / A combination of finite element method boundary element method (MEC) is studied (MEF) with the for the frame structure and two-dimensional domain. In the analyzed case, the uniaxial elements are treated by MEF, while the MEC be used to model the continuum media which can is homogeneous or not. In geomechanical problems, it is possible to simulate a sequence excavation with the necessary structural modifications. The equations of the two-dimensional domain with the influence of the frame structure, are assembled in the system which is solved by algorithm proposed by Crotty. Using the boundary element method, it is analyzed the plastic behavior of a continuum media by an incremental and iterative procedure, which is based in the initial stress process.

A busca por novos modelos matemáticos e técnicas inovadoras para análises numéricas tem sido tema de muitas pesquisas. Em análises de modelos que possuem domínios infinitos e semi-infinitos, o Método dos Elementos de Contorno (MEC) sobressai-se como uma das mais eficientes ferramentas numéricas. Por outro lado, em análises não-lineares o MEC requer a avaliação de integrais de domínio, diminuindo as vantagens de uma discretização apenas do contorno do modelo analisado. Neste trabalho apresenta-se uma técnica inovadora que trata as integrais de domínio, não adequadas para uma representação pura do contorno, em análises de modelos com materiais viscoplásticos. Na abordagem proposta, utiliza-se um novo algoritmo de visualização proposto por Noronha & Pereira para detectar as regiões de plastificação automaticamente. Este procedimento de detecção é realizado de forma incremental por meio de predições (gradiente como direção de busca) e iterações (Newton-Raphson). Uma vez que as regiões sejam obtidas, torna-se possível transformar as integrais de domínio em integrais de contorno de forma direta. Obtém-se assim uma abordagem baseada apenas na discretização do contorno dos modelos, mantendo uma das principais vantagens da utilização do MEC. Foram realizados neste trabalho alguns exemplos numéricos que apresentaram excelentes resultados em comparação com o Método dos Elementos Finitos (MEF) e com resultados encontrados na literatura. / The search for new mathematical models and innovative techniques for numerical analyses has been subject of many research studies. In analysis of models with infinite and semi-infinite domains, the Boundary Element Method (BEM) has been proved to be one of the most efficient numerical tools. On the other hand, in nonlinear analyses the BEM requires the evaluation of domain integrals, diminishing the advantages of a discretization only of the boundary of the model. This work presents an innovative technique that treats the domain integrals, not suitable for pure boundary representations, in analyses of models with viscoplastic materials. The proposed approach is based on a new post-processing algorithm developed by Noronha & Pereira to detect the plastic regions automatically. The detection procedure herein proposed is an incremental technique that uses prediction (along the gradient direction) and iteration (Newton-Raphson) loops. Once the plastic regions are found, it becomes possible to transform the domain integrals in boundary integrals in a straightforward manner. The proposed approach results in a pure boundary discretization, preserving the main advantage of the BEM. The numerical examples presented in this work are in good agreement with the Finite Element Method (FEM) and with results found in the literature.

Neste trabalho, desenvolve-se uma formulação linear de placas através do Método dos Elementos de Contorno, baseada na teoria clássica de Kirchhoff, onde a integração numérica, sobre os elementos do contorno, é feita considerando-se a técnica de subelementos. Estende-se essa formulação à análise não-linear de placas de concreto armado, através da inclusão de um campo de momentos iniciais, onde as integrais de domínio são calculadas aproximando-se o campo de momentos iniciais em células internas. Consideram-se dois modelos constitutivos para o concreto: um elasto-plástico, onde o critério utilizado é o de Von Mises, sem considerar resistência à tração, enquanto que o outro é o modelo de dano de Mazars. A distribuição das tensões é aproximada, em uma seção qualquer da placa, por pontos discretos, que seguindo um esquema gaussiano, permite a integração numérica para o cálculo dos esforços. Em cada ponto, considerado ao longo da espessura, verifica-se o modelo constitutivo adotado. Numa primeira aproximação, considera-se que a linha neutra é definida pela superfície média da placa e, numa aproximação seguinte, a posição da mesma é calculada de tal forma que a força normal resultante seja nula. / In this work a linear boundary element formulation for Kirchhoff plate in bending is presented, including a precise numerical scheme based on sub-element integration to compute the boundary matrices. Then, the formulation is extended to perform non-linear analysis of reinforced concrete slabs by incorporating initial moment fields. The domain integral required to evaluate the initial moment influences are performed by using the well known cell sub-division. The non-linear behaviour of concrete slabs is included by considering two particular models: the Von Mises criterion modified to consider no strength in tension and the damage model proposed by Mazars. Those criteria are verified at points along the plate thickness, appropriately placed to allow performing numerical integration to approach moments and normal forces using Gauss point schemes. A first model was developed assuming that the plate middle surface is coincident with the plate neutral axes. Then, this model is modified to find the proper neutral axis positions by enforcing the normal force resultants to be zero.

Neste trabalho, a formulação linear do método dos elementos de contorno - MEC, baseada nas hipóteses de Kirchhoff, é adaptada à análise de estruturas de pavimentos de edifícios, considerando-se as interações entre elementos lineares e de superfície. Leva-se em conta, além da flexão, o comportamento dos elementos como membrana. A representação integral deduzida contempla todos os elementos estruturais envolvidos, portanto garantido a monoliticidade do conjunto sem a necessidade de impor compatibilizações de deslocamentos e equilíbrio das forças generalizadas de superfícies ao longo das interfaces. A formulação integral é deduzida a partir da primeira identidade de Betti, onde a placa é considerada com variação de espessura, quer seja contínua ou abrupta. Porém, nesse trabalho apenas o caso de placas e vigas com rigidez constante são tratados. A partir dessa formulação, a fim de reduzir o número de graus de liberdade do problema, apresenta-se um modelo alternativo, onde as vigas são representadas por seus eixos médios. Estende-se essa formulação à análise não-linear, através da inclusão de campos de esforços iniciais, onde as integrais de domínio são calculadas aproximando-se o campo de esforços iniciais em células internas. A solução não-linear é obtida a partir da formulação implícita, na qual as correções que devem ser dadas aos estados de curvatura e das deformações de chapa em uma determinada iteração, são obtidas através do operador tangente consistente, que é atualizado a cada iteração, e da correção dos esforços nos pontos da placa. O critério elasto-plástico utilizado é o de Von Mises e a distribuição das tensões é aproximada, em uma seção qualquer da placa, por pontos discretos, que seguindo um esquema gaussiano, permite a integração numérica para o cálculo dos esforços. / In this work, the plate bending linear formulation of the boundary element method - BEM, based on the Kirchhoff\'s hypothesis, is extended to incorporate beam elements. The final objective of the work is to obtain a numerical model to analyse building floor structures, in which stiffness is further increased by the presence of membrane effects. From the boundary integral representations of the bending and the stretching problems a particular integral equation to represent the equilibrium of the whole body is obtained. Using this integral equation, no approximation of the generalized forces along the interface is required. Moreover, compatibility and equilibrium conditions along the interface are automatically imposed by the integral equation. An alternative formulation where the number of degrees of freedom is further reduced is also investigated. In this case, the kinematics Navier-Bernoulli hypothesis is assumed to simplify the strain field for the thin sub-regions (beams). Then, the formulation is extended to perform non-linear analysis by incorporating initial effort fields. Then non-linear solution is obtained using the concept of the local consistent tangent operator. The domain integral required, to evaluate the initial effort influences, are performed by using the well-known cell sub-division. The non-linear behaviour is evaluated by the Von Mises criterion, that is verified at points along the plate thickness, appropriately placed to allow performing numerical integration to approach moments and normal forces using Gauss point schemes.

A capacidade de prever os movimentos de uma plataforma de petróleo sujeita a ondas é bastante importante no contexto da engenharia naval e oceânica, já que esses movimentos terão diversas implicações no projeto deste sistema, com impactos diretos nos custos de produção e tempo de retorno do investimento. Esse trabalho apresenta os fundamentos teóricos sobre o problema de comportamento no mar de corpos flutuantes sujeitos a ondas de gravidades e um método numérico para solução do problema 2D no domínio do tempo. A hipótese básica adotada é a de escoamento potencial, que permitiu a utilização do método de elementos de contorno para descrever a região fluida. Optou-se pela utilização de fontes de Rankine como função de Green no desenvolvimento do método, o qual será abordado somente no contexto linear do problema matemático, delimitado através de um procedimento combinado entre expansão de Stokes e série de Taylor. As simulações são realizadas no domínio do tempo sendo, portanto, resolvido o problema de valor inicial com relação às equações do movimento e equações que descrevem a superfície-livre combinadas com dois problemas de valor de contorno, um para o potencial de velocidades e outro para o potencial de aceleração do escoamento. As equações integrais de contorno permitem transformar o sistema de equações diferenciais parciais da superfície livre num sistema de equações diferenciais ordinárias, as quais são resolvidas através do método de Runge-Kutta de 4a. ordem. As equações integrais são tratadas de forma singularizada e o método utilizado para discretizar as mesmas é de ordem baixa tanto para a função potencial quanto para a aproximação geométrica, sendo as integrações necessárias realizadas numericamente através de quadratura Gauss-Legendre. O algoritmo numérico é testado e validado através de comparações com soluções analíticas, numéricas e experimentais presentes na literatura, considerando os problemas de geração de ondas, cálculo de massa adicional e amortecimento potencial através de ensaios de oscilação forçada, testes de decaimento e, por último, resposta em ondas. Os resultados obtiveram boa concordância com aqueles adotados como paradigma. / The ability to predict the seakeeping characteristics of an offshore structure (such as an oil platform) is very important in offshore engineering since these motions have important consequences regarding its design and therefore its cost and payback period. This work presents the theoretical and numerical aspects concerning the evaluation of the 2D seakeeping problem under the potential flow hypothesis, which allows the use a Boundary Elements Method to describe the fluid region with Rankine sources as Green function. The linearized version of the mathematical problem is built by a combined Stokes expansion and Taylor series procedure and solved in time domain. The initial value problem concerning the motion and free surface equations are solved combined to the boundary value problems considering the velocity and acceleration flow potentials, which transform the partial differential equations of the free surface into ordinary differential equations, that are solved using the 4th order Runge-Kutta method. The integral equations are solved in its singularized version using a low order method both for the potential function and the geometrical approximation, with the terms of the linear system evaluated using Gauss Legendre quadrature. The numerical scheme is tested and validated considering analytical, numerical and experimental results obtained in the literature, concerning wave generation, added mass and potential damping evaluation, decay tests and response to waves. The results achieved good agreement with respect to those used as paradigm.

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.

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.

Made available in DSpace on 2014-06-11T19:25:21Z (GMT). No. of bitstreams: 0 Previous issue date: 2008-05-15Bitstream added on 2014-06-13T18:53:17Z : No. of bitstreams: 1 konda_dh_me_ilha.pdf: 2260135 bytes, checksum: 7b50e08a2458da851841fed8236b0b31 (MD5) / Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES) / Esse trabalho trata da análise linear de flexão de pavimentos de edifícios utilizando o Método dos Elementos de Contorno. A formulação de placas baseada nas hipóteses de Reissner, que possibilita a análise de placas finas e moderadamente espessas, é estendida a fim de se modelar pavimentos de edifícios. Nesse modelo, o pavimento é representado por uma placa composta de sub-regiões, sendo cada subregião a representação de uma placa ou viga. Assim, a viga é modelada por uma subregião de pequena largura e maior rigidez. As equações integrais dos deslocamentos são deduzidas utilizando-se o Método dos Resíduos Ponderados. Inicialmente as variáveis são definidas ao longo das interfaces e contorno externo da placa. Então, a fim de reduzir o número de graus de liberdade do problema, apresenta-se um modelo alternativo onde as variáveis ao longo das interfaces são escritas em função dos valores nos eixos das vigas. A validação da formulação proposta é feita através de exemplos, cujos resultados numéricos são comparados com valores previstos analiticamente ou com valores obtidos numericamente a partir de um software comercial baseado no Método dos Elementos Finitos. / This work deals with linear bending analysis of building floor structures by the Boundary Element Method. The plate formulation based on Reissner’s hypothesis, which can be used either for thin or thick plates, is extended to model building floor structures which is represented by a zoned plate, where each sub-region defines a slab or a beam. The beams are modelled as narrow sub-regions with larger thickness. The integral equations are derived by applying the weighted residual method to each subregion and summing them to obtain the equation for the whole body. Initially the values are defined on the interfaces and external boundary. Then, in order to reduce further the degrees of freedom, the interfaces values are translated to the the beam axis. The accuracy of the proposed model is confirmed by several numerical examples whose results are compared with the analytical solution if it exists or with a wellknown finite element code.

Made available in DSpace on 2014-06-11T19:25:22Z (GMT). No. of bitstreams: 0 Previous issue date: 2007-12-07Bitstream added on 2014-06-13T19:12:02Z : No. of bitstreams: 1 denipotti_gj_me_ilha.pdf: 1071673 bytes, checksum: 6d761992a1a3009f2f3628701a9ccbe8 (MD5) / Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP) / Este trabalho trata da inclusão de pilares numa formulação do Método dos Elementos de Contorno para análise linear de pavimentos de edifícios. Inicialmente é apresentada a formulação para análise do pavimento sujeito à flexão simples, onde o pavimento é modelado por uma placa composta por sub-regiões de diferentes rigidezes, sendo cada sub-região a representação de uma viga ou laje. Nesse modelo as forças de superfície são eliminadas das interfaces e todas as sub-regiões são representadas por suas superfícies médias. Então, a fim de reduzir o número de graus de liberdade do problema, aproximam-se os deslocamentos ao longo da largura da viga, obtendo-se um novo modelo em que as variáveis ficam definidas ao longo dos eixos das vigas ao invés de seus contornos. Em seguida, introduz-se o efeito de membrana nessa formulação para que se possa fazer a análise do pavimento sujeito à flexão composta, isto é, representando todas as sub-regiões por uma mesma superfície de referência. Nesse modelo considera-se aproximação linear para as forças de membrana ao longo da largura das vigas. Finalmente, estendem-se essas formulações para análise de pavimentos de edifícios com pilares, através da inclusão de pontos do domínio com carregamentos prescritos, que são escritos em termos das reações do pilar sobre a placa. Os esforços de flexão do pilar sobre a placa são escritos em termos da tensão uniformemente distribuída sobre a seção transversal do pilar e os esforços de membrana são escritos em termos de linhas de carga atuando na seção transversal do pilar. / This work deals with the inclusion of columns in a formulation of the Boundary Element Method (BEM) for linear analysis of building floors structures. Initially it is presented the BEM formulation for analyzing building floors subjected to simple bending, where the floor is modeled by a plate divided into sub-regions of different rigidities, being each sub-region the representation of either a beam or slab. In this model the tractions are eliminated on the interfaces and all the sub-regions are represented by their middle surfaces. Then, in order to reduce the number of degrees of freedom, some approximations for the displacements are adopted along the beam width, leading to a model where the values are defined on the beam axis, instead of its boundaries. Then the membrane effects are introduced into this formulation in order to represent all sub-regions by a same reference surface. In this model are assumed linear approximation for the in-plane tractions along the beam cross section. Finally, those formulations are extended in order to consider columns in the plate domain. This is done by considering domain points with prescribed loadings which are written in terms of the column reactions on the plate. The columns bending reactions are written in terms of a stress uniformly distributed over the column cross section and the membrane forces are written in terms of loading lines acting in the column cross section.

Resumo: Esse trabalho trata da análise linear de flexão de pavimentos de edifícios utilizando o Método dos Elementos de Contorno. A formulação de placas baseada nas hipóteses de Reissner, que possibilita a análise de placas finas e moderadamente espessas, é estendida a fim de se modelar pavimentos de edifícios. Nesse modelo, o pavimento é representado por uma placa composta de sub-regiões, sendo cada subregião a representação de uma placa ou viga. Assim, a viga é modelada por uma subregião de pequena largura e maior rigidez. As equações integrais dos deslocamentos são deduzidas utilizando-se o Método dos Resíduos Ponderados. Inicialmente as variáveis são definidas ao longo das interfaces e contorno externo da placa. Então, a fim de reduzir o número de graus de liberdade do problema, apresenta-se um modelo alternativo onde as variáveis ao longo das interfaces são escritas em função dos valores nos eixos das vigas. A validação da formulação proposta é feita através de exemplos, cujos resultados numéricos são comparados com valores previstos analiticamente ou com valores obtidos numericamente a partir de um software comercial baseado no Método dos Elementos Finitos. / Abstract: This work deals with linear bending analysis of building floor structures by the Boundary Element Method. The plate formulation based on Reissner's hypothesis, which can be used either for thin or thick plates, is extended to model building floor structures which is represented by a zoned plate, where each sub-region defines a slab or a beam. The beams are modelled as narrow sub-regions with larger thickness. The integral equations are derived by applying the weighted residual method to each subregion and summing them to obtain the equation for the whole body. Initially the values are defined on the interfaces and external boundary. Then, in order to reduce further the degrees of freedom, the interfaces values are translated to the the beam axis. The accuracy of the proposed model is confirmed by several numerical examples whose results are compared with the analytical solution if it exists or with a wellknown finite element code. / Orientador: Gabriela Rezende Fernandes / Coorientador: Luiz Carlos Facundo Sanches / Banca: Renato Bertolino Júnior / Banca: Leandro Palermo Júnior / Mestre