Neste trabalho foi desenvolvida uma formulação alternativa para o acoplamento entre o método dos elementos de contorno (MEC) e o método dos elementos finitos (MEF) para análise não linear geométrica de estruturas reticuladas ligadas a meios contínuos bidimensionais heterogêneos, aplicado a problemas de interação solo-estrutura. O solo foi considerado com comportamento elástico linear e modelado via MEC por meio de uma formulação alternativa à clássica técnica de sub-região permitindo a consideração de múltiplas inclusões mais ou menos rígidas do que o material padrão e de linhas de carga internas aos domínios. Este código foi então acoplado ao programa AcadFrame, baseado no MEF posicional para análise não linear geométrica de pórticos com consideração de cinemática exata. O acoplamento numérico foi realizado por meio de uma formulação algébrica onde a matriz de rigidez do solo e a força de contato são condensadas e somadas à matriz e ao vetor de forças internas da estrutura a cada iteração no processo de Newton-Raphson. Em ambos os programas foi utilizada uma generalização do grau de aproximação dos elementos através dos polinômios de Lagrange, o que permite a utilização de elementos curvos de alta ordem. Foi utilizada ainda a técnica dos mínimos quadrados para reduzir as oscilações de forças de superfície no contato. Os resultados obtidos de forma geral são bastante satisfatórios e comprovam a eficiência da formulação. O trabalho permite a análise de problemas de edificações apoiadas sobre solos estratificados com múltiplas inclusões e linhas de carga. Permite tanto a análise de elementos apoiados diretamente sobre o solo (sapatas, radies) quanto de elementos internos e em qualquer direção, como no caso de estacas verticais ou inclinadas. Pode-se inclusive considerar as estacas passando por diferentes camadas de solo. A aplicação pode ser estendida ainda a outros problemas elásticos, acoplamento entre peças mecânicas e análise de materiais compostos. / This work presents an alternative coupling of the boundary element method (BEM) and the finite element method (FEM) to create a computer program for non linear geometric analysis of frames coupled to continuous domains, applied to soil-structure interaction. A linear elastic behavior is considered for the soil, modeled by BEM. An alternative formulation is adopted for the classic sub-region technique, allowing the consideration of multiple inclusions and load lines inside the soil domain. The BEM computational code is coupled to the AcadFrame software, based on positional FEM for non linear geometric analysis of frames, considering exact kinematics. The numerical coupling is made by an algebraic formulation where the soil stiffness matrix and contact forces are condensed and added to the structure matrix and internal forces for each iteration on Newton-Raphson process. On both programs it is adopted a generalization of the element degree assuming the Lagrange polynomials, which allows the use of curved high order elements. It was also implemented the least square method in order to obtains better and smoother results of surface forces in the contact interface. The obtained results are satisfactory and prove the formulation efficiency. The program allows the analysis of buildings supported by layered soils with multiples inclusions and load lines. It considers directly supported elements over the soil (footing foundations, radies) and internal elements in any direction, like vertical and diagonal piles. It can also consider piles going through different layers of the soil. This formulation can be applied to other elastic problems like coupling between mechanic pieces and composite material analysis.

Materiais enrijecidos são encontrados cotidianamente na engenharia, como na construção de aeronaves, veículos, navios e manufaturas diversas. Na Engenharia Civil, são exemplos de materiais enrijecidos o concreto armado e o solo reforçado. A solução aproximada para esse tipo de problema pode ser encontrada pelo uso de métodos numéricos como o Método dos Elementos Finitos (MEF) e o Método dos Elementos de Contorno (MEC). Particularmente, o MEC é muito vantajoso em casos onde se tem elevados gradientes de tensão, como em problemas da mecânica da fratura, além de modelar de forma natural meios infinitos, sendo muito utilizado na análise de interação solo-estrutura ou de túneis. No entanto, como no MEC somente o contorno é discretizado, torna-se necessária a utilização de um acoplamento numérico para que o efeito mecânico dos enrijecedores seja considerado. Para esse acoplamento, podem ser utilizados diversos métodos numéricos, sendo neste trabalho realizado com o MEF e o MEC em sua forma unidimensional (MEC 1D). Dentro desse contexto, o acoplamento com o MEC 1D destaca-se por vantagens, como a compatibilidade dos métodos e a redução de aproximações. Com isso, este trabalho tem como principal objetivo o desenvolvimento e comparação de soluções numéricas para o problema de domínios enrijecidos, no âmbito do MEC aplicado a problemas bidimensionais. Os enrijecedores foram considerados como elementos de treliça, sendo realizado inicialmente na sua forma mais usual, modelando-os por meio do MEF. A seguir foi implementada nova formulação para o acoplamento, na qual os enrijecedores são modelados por meio do MEC 1D. A implementação do MEC 1D foi validada pela comparação de resultados com soluções analítica e do programa computacional FTOOL. No caso do acoplamento, os resultados de ambas formulações foram comparados com resposta do programa computacional ANSYS. Foram avaliados quatro exemplos, sendo dois isotrópicos e dois anisotrópicos. Foram aplicadas diferentes condições de carregamentos, apoios e materiais, sendo utilizadas diversas discretizações e graus de aproximação nos enrijecedores. Verificou-se os efeitos da mudança do grau de aproximação, mantendo-se constante o número de graus de liberdade. Os resultados obtidos foram similares às respostas do ANSYS, foram mecanicamente equivalentes, não havendo diferença relevante em custo computacional. O acoplamento com o MEC 1D levou a resultados estáveis e, em geral, melhores do que com o MEF. No caso de aproximação quadrática, o acoplamento com o MEF frequentemente levou a resultados instáveis. Com o aumento do grau de aproximação, os resultados de ambos métodos se aproximam, se tornando muito próximos com aproximação do quarto grau. Verificou-se que próximo às pontas dos enrijecedores ocorrem concentrações de tensão, havendo maior influência da discretização e grau de aproximação adotado nestas áreas. / The reinforced materials are commonly encountered in engineering applications, as in the construction of airplanes, vehicles, ships and general manufactures. In Civil Engineering, the reinforced concrete and the reinforced soil are examples of reinforced materials. The solution of this kind of problem can be found by the use of numerical methods as the Finite Element Method (FEM) and the Boundary Element Method (BEM). Particularly, the BEM use in cases with high stress gradients, such as in fracture mechanics problems, is very advantageous. The BEM models infinite domains in a natural way, being largely used in analysis as soil-structure interaction and tunnel modelling. However, as in the BEM only the boundary is discretized, became necessary the use of a numerical coupling to consider the mechanical effect of the stiffeners. For the coupling many numerical methods can be used, in this study it was realized with the FEM and with the BEM in its unidimensional way (BEM 1D). In this context, the coupling BEM/BEM 1D stands out for its advantages, like the compatibility between the methods and the reduction of approximations. Therefore, the main objective of this study is the development and the comparison of numerical solutions for reinforced media problems, considering the framework of the bidimensional BEM problems. The coupling was implemented considering the stiffeners as truss elements, initially it was implemented in the most usual way, by modelling the stiffeners through the FEM. In the following, a new coupling formulation was proposed, in which the stiffeners are modelled through the BEM 1D method. The BEM 1D implementation was validated by comparing the results with analytical and numerical solutions, the last one obtained with the software FTOOL. Both coupling strategies were compared with the solution obtained by the software ANSYS. Four examples were evaluated, with two isotropic and two anisotropic domains. Different boundary conditions, supports and materials were applied, as many discretization and the approximation degree of the stiffeners were tested. The effects of changing the approximation degree in both coupling strategies keeping the degrees of freedom constant were analysed. The results obtained with the implemented algorithms were mechanically coherent, and they were similar to ANSYS results. For all examples, there is no significant computational costs differences between the two coupling strategies. However, the coupling with the BEM 1D leaded to more stable results and better approximations. It was observed that the MEF results were instable for many results, mainly in the quadratic approximations. When the approximation degree rises, the methods tend to converge to equivalent results, becoming very close in fourth degree approximation. Lastly, it was observed stress concentration in the stiffeners ends. In these regions, the discretization and the approximation degree have large influence to the numerical response.

Neste trabalho foi desenvolvida uma formulação alternativa para o acoplamento entre o método dos elementos de contorno (MEC) e o método dos elementos finitos (MEF) para análise não linear geométrica de estruturas reticuladas ligadas a meios contínuos bidimensionais heterogêneos, aplicado a problemas de interação solo-estrutura. O solo foi considerado com comportamento elástico linear e modelado via MEC por meio de uma formulação alternativa à clássica técnica de sub-região permitindo a consideração de múltiplas inclusões mais ou menos rígidas do que o material padrão e de linhas de carga internas aos domínios. Este código foi então acoplado ao programa AcadFrame, baseado no MEF posicional para análise não linear geométrica de pórticos com consideração de cinemática exata. O acoplamento numérico foi realizado por meio de uma formulação algébrica onde a matriz de rigidez do solo e a força de contato são condensadas e somadas à matriz e ao vetor de forças internas da estrutura a cada iteração no processo de Newton-Raphson. Em ambos os programas foi utilizada uma generalização do grau de aproximação dos elementos através dos polinômios de Lagrange, o que permite a utilização de elementos curvos de alta ordem. Foi utilizada ainda a técnica dos mínimos quadrados para reduzir as oscilações de forças de superfície no contato. Os resultados obtidos de forma geral são bastante satisfatórios e comprovam a eficiência da formulação. O trabalho permite a análise de problemas de edificações apoiadas sobre solos estratificados com múltiplas inclusões e linhas de carga. Permite tanto a análise de elementos apoiados diretamente sobre o solo (sapatas, radies) quanto de elementos internos e em qualquer direção, como no caso de estacas verticais ou inclinadas. Pode-se inclusive considerar as estacas passando por diferentes camadas de solo. A aplicação pode ser estendida ainda a outros problemas elásticos, acoplamento entre peças mecânicas e análise de materiais compostos. / This work presents an alternative coupling of the boundary element method (BEM) and the finite element method (FEM) to create a computer program for non linear geometric analysis of frames coupled to continuous domains, applied to soil-structure interaction. A linear elastic behavior is considered for the soil, modeled by BEM. An alternative formulation is adopted for the classic sub-region technique, allowing the consideration of multiple inclusions and load lines inside the soil domain. The BEM computational code is coupled to the AcadFrame software, based on positional FEM for non linear geometric analysis of frames, considering exact kinematics. The numerical coupling is made by an algebraic formulation where the soil stiffness matrix and contact forces are condensed and added to the structure matrix and internal forces for each iteration on Newton-Raphson process. On both programs it is adopted a generalization of the element degree assuming the Lagrange polynomials, which allows the use of curved high order elements. It was also implemented the least square method in order to obtains better and smoother results of surface forces in the contact interface. The obtained results are satisfactory and prove the formulation efficiency. The program allows the analysis of buildings supported by layered soils with multiples inclusions and load lines. It considers directly supported elements over the soil (footing foundations, radies) and internal elements in any direction, like vertical and diagonal piles. It can also consider piles going through different layers of the soil. This formulation can be applied to other elastic problems like coupling between mechanic pieces and composite material analysis.

Este trabalho busca aplicar técnicas de confiabilidade ao problema de grupo de estacas utilizadas como fundação de estruturas correntes. Para isso, lança-se mão de um modelo tridimensional de interação estaca-solo onde estão presentes o Método dos Elementos de Contorno (MEC) e o método dos Elementos Finitos (MEF) que atuam de forma acoplada. O MEC, com as soluções fundamentais de Mindlin (meio semi-infinito, homogêneo, isotrópico e elástico-linear é utiliza), é utilizado para modelar o solo. Já o MEF é utilizado para modelar as estacas. Definido o modelo de funcionamento estrutural das estacas, parte-se para a aplicação de métodos trazidos da confiabilidade estrutural para avaliação da adequabilidade em relação aos estados limite de serviço e estados limites últimos. Os métodos de confiabilidade utilizados foram o Método de Monte Carlo, o método FOSM (First-Order Second-Moment) e o método FORM (First-Order Reliability Method). / This work seeks to apply reliability techniques to the problem of piles groups used as current structures foundation. For this, makes use of a three-dimensional model of pile-soil interaction with the boundary element method (BEM) and the finite element method (FEM) working coupled. The BEM, with Mindlin fundamental solutions (semi-infinite medium, homogeneous, isotropic and linear elastic) is used to model the soil. The MEF is used to model the piles. Defined the model of structural functioning of the piles, the aim goes to the application of structural reliability for assessing the adequacy of the serviceability limit states and ultimate limit states. Reliability methods used were the Monte Carlo method, the FOSM (First-Order Second-Moment) method and the FORM method (First-Order Reliability Method).

Materiais enrijecidos são encontrados cotidianamente na engenharia, como na construção de aeronaves, veículos, navios e manufaturas diversas. Na Engenharia Civil, são exemplos de materiais enrijecidos o concreto armado e o solo reforçado. A solução aproximada para esse tipo de problema pode ser encontrada pelo uso de métodos numéricos como o Método dos Elementos Finitos (MEF) e o Método dos Elementos de Contorno (MEC). Particularmente, o MEC é muito vantajoso em casos onde se tem elevados gradientes de tensão, como em problemas da mecânica da fratura, além de modelar de forma natural meios infinitos, sendo muito utilizado na análise de interação solo-estrutura ou de túneis. No entanto, como no MEC somente o contorno é discretizado, torna-se necessária a utilização de um acoplamento numérico para que o efeito mecânico dos enrijecedores seja considerado. Para esse acoplamento, podem ser utilizados diversos métodos numéricos, sendo neste trabalho realizado com o MEF e o MEC em sua forma unidimensional (MEC 1D). Dentro desse contexto, o acoplamento com o MEC 1D destaca-se por vantagens, como a compatibilidade dos métodos e a redução de aproximações. Com isso, este trabalho tem como principal objetivo o desenvolvimento e comparação de soluções numéricas para o problema de domínios enrijecidos, no âmbito do MEC aplicado a problemas bidimensionais. Os enrijecedores foram considerados como elementos de treliça, sendo realizado inicialmente na sua forma mais usual, modelando-os por meio do MEF. A seguir foi implementada nova formulação para o acoplamento, na qual os enrijecedores são modelados por meio do MEC 1D. A implementação do MEC 1D foi validada pela comparação de resultados com soluções analítica e do programa computacional FTOOL. No caso do acoplamento, os resultados de ambas formulações foram comparados com resposta do programa computacional ANSYS. Foram avaliados quatro exemplos, sendo dois isotrópicos e dois anisotrópicos. Foram aplicadas diferentes condições de carregamentos, apoios e materiais, sendo utilizadas diversas discretizações e graus de aproximação nos enrijecedores. Verificou-se os efeitos da mudança do grau de aproximação, mantendo-se constante o número de graus de liberdade. Os resultados obtidos foram similares às respostas do ANSYS, foram mecanicamente equivalentes, não havendo diferença relevante em custo computacional. O acoplamento com o MEC 1D levou a resultados estáveis e, em geral, melhores do que com o MEF. No caso de aproximação quadrática, o acoplamento com o MEF frequentemente levou a resultados instáveis. Com o aumento do grau de aproximação, os resultados de ambos métodos se aproximam, se tornando muito próximos com aproximação do quarto grau. Verificou-se que próximo às pontas dos enrijecedores ocorrem concentrações de tensão, havendo maior influência da discretização e grau de aproximação adotado nestas áreas. / The reinforced materials are commonly encountered in engineering applications, as in the construction of airplanes, vehicles, ships and general manufactures. In Civil Engineering, the reinforced concrete and the reinforced soil are examples of reinforced materials. The solution of this kind of problem can be found by the use of numerical methods as the Finite Element Method (FEM) and the Boundary Element Method (BEM). Particularly, the BEM use in cases with high stress gradients, such as in fracture mechanics problems, is very advantageous. The BEM models infinite domains in a natural way, being largely used in analysis as soil-structure interaction and tunnel modelling. However, as in the BEM only the boundary is discretized, became necessary the use of a numerical coupling to consider the mechanical effect of the stiffeners. For the coupling many numerical methods can be used, in this study it was realized with the FEM and with the BEM in its unidimensional way (BEM 1D). In this context, the coupling BEM/BEM 1D stands out for its advantages, like the compatibility between the methods and the reduction of approximations. Therefore, the main objective of this study is the development and the comparison of numerical solutions for reinforced media problems, considering the framework of the bidimensional BEM problems. The coupling was implemented considering the stiffeners as truss elements, initially it was implemented in the most usual way, by modelling the stiffeners through the FEM. In the following, a new coupling formulation was proposed, in which the stiffeners are modelled through the BEM 1D method. The BEM 1D implementation was validated by comparing the results with analytical and numerical solutions, the last one obtained with the software FTOOL. Both coupling strategies were compared with the solution obtained by the software ANSYS. Four examples were evaluated, with two isotropic and two anisotropic domains. Different boundary conditions, supports and materials were applied, as many discretization and the approximation degree of the stiffeners were tested. The effects of changing the approximation degree in both coupling strategies keeping the degrees of freedom constant were analysed. The results obtained with the implemented algorithms were mechanically coherent, and they were similar to ANSYS results. For all examples, there is no significant computational costs differences between the two coupling strategies. However, the coupling with the BEM 1D leaded to more stable results and better approximations. It was observed that the MEF results were instable for many results, mainly in the quadratic approximations. When the approximation degree rises, the methods tend to converge to equivalent results, becoming very close in fourth degree approximation. Lastly, it was observed stress concentration in the stiffeners ends. In these regions, the discretization and the approximation degree have large influence to the numerical response.

Neste trabalho, é desenvolvido um programa de computador para a análise estática e tridimensional de problemas de interação solo-estrutura. O programa permite considerar várias camadas de solo, cada qual com características físicas diferentes. Sobre este solo, o qual pode conter estacas, podem ser apoiados diversos tipos de estruturas, tais como placas e até um edifício. Todos os materiais considerados são homogêneos, isotrópicos, elásticos e lineares. O solo tridimensional é modelado com o método dos elementos de contorno (MEC), empregando as soluções fundamentais de Kelvin e uma técnica alternativa na consideração do maciço não-homogêneo. Esta técnica, que é uma contribuição original deste trabalho, é baseada no relacionamento das soluções fundamentais de deslocamento dos diferentes domínios, permitindo que sejam analisados como um único sólido sem a necessidade de equações de equilíbrio e compatibilidade. Isso reduz o sistema de equações final e melhora a precisão dos resultados, conforme comprovado nos exemplos apresentados. Para reduzir o custo computacional sem prejudicar a precisão dos resultados, é utilizada uma malha de elementos de contorno infinitos (ECI) nas bordas da malha de ECs para modelar o comportamento das variáveis de campo em longas distâncias. A formulação do ECI mapeado utilizado é outra contribuição original deste trabalho, sendo baseado em um EC triangular. É demonstrado por meio de exemplos que tal formulação é eficiente para a redução de malha, contribuindo de forma significativa na redução do custo computacional. Todas as estruturas que interagem com o solo, incluindo as de fundação, são simuladas empregando o método dos elementos finitos (MEF). Cada estaca é modelada como uma linha de carga empregando um único elemento finito com 14 parâmetros nodais, o qual utiliza funções de forma do quarto grau para aproximar os deslocamentos horizontais, do terceiro grau para as forças horizontais e deslocamentos verticais, do segundo grau para as forças cisalhantes verticais e constantes para as reações da base. Este elemento é empregado em outros trabalhos, no entanto os autores utilizam as soluções fundamentais de Mindlin na consideração da presença da estaca no solo. Desta forma, a formulação desenvolvida neste trabalho com as soluções fundamentais de Kelvin pode ser considerada mais uma contribuição original. No edifício, que pode incluir um radier como estrutura de fundação, são utilizados dois tipos de EFs. Os pilares e vigas são simulados com elementos de barra, os quais possuem dois nós e seis graus de liberdade por nó. As lajes e o radier são modelados empregando elementos planos, triangulares e com três nós. Nestes EFs triangulares são superpostos efeitos de membrana e flexão, totalizando também seis graus de liberdade por nó. O acoplamento MEC/MEF é feito transformando as cargas de superfície do MEC em carregamentos nodais reativos no MEF. Além de exemplos específicos nos Capítulos teóricos, um Capítulo inteiro é dedicado a demonstrar a abrangência e precisão da formulação desenvolvida, comparando-a com resultados de outros autores. / In this work, a computer code is developed for the static analysis of three-dimensional soil-structure interaction problems. The program allows considering a layered soil, which may contain piles. This soil may support several structures, such as shells or even an entire building. All materials are considered homogeneous, isotropic, elastic and linear. The three-dimensional soil is modeled with the boundary element method (BEM), employing Kelvin fundamental solutions and an alternative multi-region technique. This technique, which is an original contribution of this work, is based on relating the displacement fundamental solution of the different domains, allowing evaluating them as an unique solid and not requiring compatibility or equilibrium equations. In such a way, the final system of equations is reduced and more accurate results are obtained, as demonstrated in the presented examples. In order to reduce the computational cost maintaining the accuracy, an infinite boundary element (IBE) mesh is employed at the BE mesh limits to model the far field behavior. The mapped IBE utilized, based on a triangular EC, is another original contribution of this work. In the presented examples it is demonstrated that this IBE formulation is efficient for mesh reduction, implying on a significant computational cost reduction. All structures that interact with the soil, including the foundations, are simulated with de finite element method (FEM). The piles are modeled using a one-dimensional 14 parameter finite element, with forth degree shape functions for horizontal displacement approximation, third degree shape functions for horizontal forces and vertical displacement, second degree shape functions for vertical share force, and constant for the base reaction. This element is employed in other works, however the authors utilize Mindlin fundamental solutions for the pile presence consideration in the soil. In such a way, the formulation developed in this work with Kelvin fundamental solutions may be considered one more original contribution. The building, which may include a radier as a foundation structure, is modeled using two types os FEs. Piles and beams are simulated using bar FEs with two nodes and six degrees of freedom per node. The radier and pavements are modeled employing plane triangular three-node FEs. In these FEs plate and membrane effects are superposed, totalizing six degrees of freedom per node. FEM/BEM coupling is made by transforming the BEM tractions in nodal reactions in the FEM. Even though specific examples are presented in the theoretical Chapters, a role Chapter is dedicated for demonstrating the formulation accuracy and coverage. In most examples, the results are compared with the ones obtained by other authors.

Apresenta-se neste trabalho, uma formulação de cascas laminadas anisotrópicas enrijecidas ou não, considerando-se não-linearidade física com lei de fluxo não-associativa e acoplamento com meio contínuo tridimensional viscoelástico. Para tanto, são desenvolvidos elementos finitos triangulares planos com aproximação cúbica de variáveis para modelagem das cascas e elementos de barra de mesma aproximação para os elementos de barra geral (enrijecedores). A cinemática de laminados, ou Reissner geral, é utilizada para ambos possibilitando a representação de estruturas enrijecidas excentricamente e consideração de elementos compostos de camadas com diferentes propriedades físicas e espessuras, tornando-se assim a formulação aplicável a um grande número de problemas. Com relação à plasticidade na casca, adota-se o critério de Tsai-Wu para materiais anisotrópicos gerais, obtendo-se expressões fechadas para o multiplicador plástico com fluxo não-associativo. Nas barras, critérios uniaxiais são considerados, desprezando-se a contribuição do cisalhamento na plastificação. Para estes elementos, permite-se a utilização de diagrama multilinear para a relação tensão x deformação. A modelagem do meio contínuo viscoelástico é realizada utilizando-se elementos de contorno triangulares com aproximação linear de variáveis. As soluções fundamentais de Kelvin e de Mindlin são apresentadas e implementadas. O acoplamento foi realizado utilizando-se técnica de matriz de rigidez equivalente, proporcionando uma contribuição direta das matrizes do MEC na matriz de rigidez do MEF. Exemplos gerais são resolvidos para a verificação e validação da formulação proposta e implementada / This work presents an anisotropic laminated stiffened shell formulation, considering physical non-linearity with non-associative law, coupled to viscoelastic three-dimensional continuum medium. Plane triangular finite elements with cubic approximation for nodal variables are developed to model the shell. Bar elements with the same approximation are derived for the general bar element. Laminated kinematics is used for both elements, making possible the representation of eccentrically stiffened structures and the consideration of composed elements with different properties and thickness for each layer. Therefore, the formulation is applicable for a large number of problems. In order to model plasticity in shell, the Tsai-Wu criterion for general anisotropic materials is adopted. Closed expression for the plastic multiplier using non-associative law is founded. For bars, uniaxial criterion is considered, and shear contribution for plasticity is neglected. For these elements, the use of multilinear stress x strain relation is developed. The viscoelastic continuum is modeled by triangular boundary elements with linear approximation of variables. The fundamental solutions of Kelvin and Mindlin are presented and implemented. The coupling is made by the equivalent stiffness matrix method, making possible a direct contribution of the BEM matrix on the FEM stiffness matrix. General examples are presented to verify and validate the proposed formulation

O objetivo central deste trabalho é o estudo da interação do solo com a estrutura. Para tanto, são introduzidos mais recursos na ferramenta numérica desenvolvida no trabalho de Almeida (2003a). O solo é modelado pelo método dos elementos de contorno (MEC) tridimensional, aplicando a solução fundamental de Kelvin. É possível analisar problemas nos quais o solo é composto por camadas de diferentes características físicas, apoiadas em uma superfície de deslocamento nulo e enrijecidas por elementos de fundação, também modelados pelo MEC tridimensional. A superestrutura tridimensional, diferentemente do modelo utilizado em Almeida (2003a), é simulada pelo método dos elementos finitos (MEF), sendo composta por elementos planos e reticulares com seis graus de liberdade por nó. Também é introduzido no programa o recurso de simular um número qualquer de blocos, modelados pelo MEC tridimensional, apoiados sobre o solo. Estes blocos podem ser utilizados como elementos de fundação para o edifício, permitindo estudar a interação do solo em conjunto com os blocos e o edifício. São analisados alguns exemplos, nos quais é validada a formulação empregada e é demonstrada a necessidade de se considerar a interação do solo com a estrutura em problemas práticos de engenharia / The main objective of this work is to study the soil structure interaction problem. For such, more resources in the numerical tool developed in Almeida (2003a) are introduced. The soil is simulated by the three-dimensional boundary element method (BEM), applying Kelvin’s fundamental solution. It is possible to analyze problems in which the soil is composed by layers of different physical characteristics, supported by a rigid and adhesive interface and reinforced by foundation elements, also simulated by the three-dimensional BEM. The three-dimensional superstructure is simulated using the finite element method (FEM), with shell and frame elements with six degrees of freedom by node. This model is different of the one used in Almeida (2003a). It is also introduced in the program the resource to consider blocks, simulated by the three-dimensional BEM and supported by the soil. These blocks can be used as foundation elements for the building, coupling the non-homogeneous soil-foundation-blocks-superstructure system as a whole. Some examples are analyzed, in order to validate the theory employed and demonstrate the necessity of considering the soil structure interaction in practical problems of engineering

Neste trabalho, uma combinação do Método dos Elementos de Contorno (MEC) com o Método dos Elementos Finitos (MEF) é apresentada para análise bidimensional de sólidos elastostáticos reforçados, sendo considerados modelos de aderência no acoplamento. O elemento de contorno é adotado para modelar o comportamento do domínio, enquanto que o modelo por elementos finitos é utilizado para modelar o enrijecedor. Devido às singularidades nas equações integrais do MEC, estudou-se o erro ocasionado pelos integrandos de ordem e e como conseqüência sugerem-se, neste trabalho, equações mais simples para representar o erro das integrações. Para a formulação do acoplamento, um polinômio do terceiro grau é adotado para aproximar tanto o campo de deslocamento quanto a rotação do enrijecedor, enquanto aproximações lineares são usadas para representar a força de contato entre o domínio e o enrijecedor. Modelos de escorregamento, apresentados, são lineares e governados em função do carregamento escrito em termos das forças de contato e o deslocamento relativo. A partir da combinação entre o MEC e o MEF obtém-se uma matriz retangular contendo duas equações para o MEC e uma para o MEF. O resultado das equações algébricas redundantes é eliminado pela aplicação do procedimento dos mínimos quadrados. Exemplos ilustram o bom ajuste e os melhores resultados proporcionados pelo controle do erro das equações integrais, mostrando ainda através de exemplos, a potencialidade e as limitações no acoplamento entre os dois materiais, considerando modelos de aderência ou não. / In this work it is presented a coupling between the Boundary Element Method and the Finite Element Method for two-dimensional elastostatic analysis of reinforced bodies considering adherence. The Boundary Element is used to model the matrix while the reinforcement is modeled by the Finite Element. Due to the inherent singularities present in Boundary Element formulations the quadrature rules, used to develop the necessary integrals may present undesired errors. In this sense the behavior of this integration error is studied and a simple way to control it is proposed along the work. Regarding the coupling formulation a third degree polynomial is adopted to describe the displacements and rotations of the reinforcement, while a linear polynomial is used to describe the contact forces among the continuum and the reinforcement. Adherence (or sliding) models are presented and implemented in the computer code. A linear relation between relative displacement and transmitted force is adopted. From difference of approximation regarding contact forces and displacements a rectangular matrix arrises from the BEM/FEM coupling. The additional equations are eliminated by the use of a least square method based on the multiplication of transpose matrices. Examples are shown to demonstrate the good behavior of error control applied on gaussian quadratures regarding Boundary Element simulations for coupled or not situations, considering or not adherence models.

O estudo do comportamento mecânico do complexo sistema advindo da interação entre solo/subestrutura/superestrutura é o tema do trabalho. Neste contexto, a representação do maciço é feita usando-se o método dos elementos de contorno (MEC) em abordagem 3D, de maneira que se possa simular o maciço com características mecânicas não-homogêneas, além de se considerar uma camada de apoio indeslocável a distâncias prescritas a priori e condição de aderência perfeita. A subestrutura também é representada via MEC tridimensional, a qual está imersa dentro deste meio heterogêneo. A infra e a superestrutura são modeladas empregando o método dos elementos finitos (MEF), com o uso de elementos estruturais reticulares e elementos laminares. São apresentados alguns exemplos em que se valida a formulação e outros que demonstram a potencialidade e a necessidade de se empregar a formulação para a melhor análise do complexo fenômeno em estudo. Por fim, demonstra-se a obrigatoriedade de se otimizar a formulação, empregando-se duas grandes ferramentas numéricas: o paralelismo e o emprego de um adequado método de resolução de sistemas esparsos. / The analysis of the soil-structure system interaction is a vast field of interest in the area of civil engineering. A realistic representation of its behaviour. Thus, in the present research, the soil is considered a non-homogeneous continuum supported by a rigid and adhesive interface and modelled by boundary element method via Kelvin solution in 3D space. The foundation is also modelled by this above-mentioned modelling technique. The raft foundation and the superstructure are represented by finite shell and 3D frame elements. In order to estimate the accuracy and the potentiality of the proposed numerical formulation, some examples are validated when compared to similar approaches, and others simulations are presented to stress the necessity of coupling the non-homogeneous soil-foundation-radier-superstructure system as a whole. Finally, to acquire numerical time efficiency, it is shown that it is imperative to apply parallel processing and sparse techniques for the solution of the final system.