Fundações do tipo radier estaqueado são aquelas formadas pelos elementos estruturais de placa e estacas (elementos de barras) e o solo . Ao contrário de outras tipos de fundações, onde a carga da superestrutura é transferida ao solo pelo radier ou pelas estacas apenas, no radier estaqueado a contribuição das estacas, bem como a do radier são consideradas. As estacas transferem as cargas da superestrutura ao solo e, assim, permitem a redução dos recalques de uma forma muito econômica. O objetivo do presente trabalho é a análise da interação solo-estrutura através do acoplamento MEC-MEF. O solo é considerado um semi-espaço homogêneo, elástico e linear governado pela equação de Navier e modelado pelo Método dos Elementos de Contorno (MEC), admitindo a solução fundamental de Mindlin. As estacas são modeladas pelo Método dos Elementos Finitos (MEF) e cada elemento possui quatro nós. Além disso, as estacas podem receber forças horizontais, verticais e momentos. A tensão de cisalhamento ao longo da estaca é aproximada por um polinômio do segundo grau e as forças na direção horizontal são aproximadas por um polinômio do quarto grau. O elemento de fundação que faz a ligação do pilar com a estaca é representado por uma placa de grande rigidez, que apresenta o comportamento de um bloco. A interação entre o radier estaqueado e o solo é feita através da reação resultante da interação estaca-solo, nos nós com estaca. A interface radier-solo é dividida em elementos triangulares e para a reação do solo considera-se a variação linear ao longo de cada elemento. A superestrutura é modelada pelo MEF. Vários exemplos de interação solo-estrutura são estudados nesta tese, e mostram que as soluções obtidas a partir do programa computacional desenvolvido no presente trabalho denominado SSI estão de acordo com outros autores. / Piled raft foundations are structures consisting of piles, the raft and the soil. Unlike classical foundation design where the building load is either transferred by the raft or the piles alone, in a piled raft foundation the contribution of the piles as well as the raft is taken into account. The piles transfer a part of the building loads into the soil and thereby allow the reduction of settlement in a very economic way. The objective of the present work is the analysis of soil-structure interaction using BEM-FEM coupling. The soil, assumed to be an elastic linear homogeneous half space is governed by Navier\'s equation and it is modeled by the Boundary Elements Method (BEM) using Mindlin\'s fundamental solution. The piles are modeled by the Finite Element Method (FEM) with four nodes each. In addition, the piles can received horizontal and vertical forces and bending moments. The shear traction along the pile is approximated by a second-degree polynomial and the tractions in the horizontal direction are approximated by a fourth degree polynomial. The cap of the pile group is assumed to be rigid. The interaction between the raft and soil is made through the subgrade reaction. The soil-cap interface is divided into triangular elements and the subgrade reaction is assumed to vary linearly across each element. The building\'s structure is modeled by FEM. Several soil structure interaction examples are studied in this thesis, and they show that the solutions obtained from program SSI are in good agreement with others authors.

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

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

Fundações do tipo radier estaqueado são aquelas formadas pelos elementos estruturais de placa e estacas (elementos de barras) e o solo . Ao contrário de outras tipos de fundações, onde a carga da superestrutura é transferida ao solo pelo radier ou pelas estacas apenas, no radier estaqueado a contribuição das estacas, bem como a do radier são consideradas. As estacas transferem as cargas da superestrutura ao solo e, assim, permitem a redução dos recalques de uma forma muito econômica. O objetivo do presente trabalho é a análise da interação solo-estrutura através do acoplamento MEC-MEF. O solo é considerado um semi-espaço homogêneo, elástico e linear governado pela equação de Navier e modelado pelo Método dos Elementos de Contorno (MEC), admitindo a solução fundamental de Mindlin. As estacas são modeladas pelo Método dos Elementos Finitos (MEF) e cada elemento possui quatro nós. Além disso, as estacas podem receber forças horizontais, verticais e momentos. A tensão de cisalhamento ao longo da estaca é aproximada por um polinômio do segundo grau e as forças na direção horizontal são aproximadas por um polinômio do quarto grau. O elemento de fundação que faz a ligação do pilar com a estaca é representado por uma placa de grande rigidez, que apresenta o comportamento de um bloco. A interação entre o radier estaqueado e o solo é feita através da reação resultante da interação estaca-solo, nos nós com estaca. A interface radier-solo é dividida em elementos triangulares e para a reação do solo considera-se a variação linear ao longo de cada elemento. A superestrutura é modelada pelo MEF. Vários exemplos de interação solo-estrutura são estudados nesta tese, e mostram que as soluções obtidas a partir do programa computacional desenvolvido no presente trabalho denominado SSI estão de acordo com outros autores. / Piled raft foundations are structures consisting of piles, the raft and the soil. Unlike classical foundation design where the building load is either transferred by the raft or the piles alone, in a piled raft foundation the contribution of the piles as well as the raft is taken into account. The piles transfer a part of the building loads into the soil and thereby allow the reduction of settlement in a very economic way. The objective of the present work is the analysis of soil-structure interaction using BEM-FEM coupling. The soil, assumed to be an elastic linear homogeneous half space is governed by Navier\'s equation and it is modeled by the Boundary Elements Method (BEM) using Mindlin\'s fundamental solution. The piles are modeled by the Finite Element Method (FEM) with four nodes each. In addition, the piles can received horizontal and vertical forces and bending moments. The shear traction along the pile is approximated by a second-degree polynomial and the tractions in the horizontal direction are approximated by a fourth degree polynomial. The cap of the pile group is assumed to be rigid. The interaction between the raft and soil is made through the subgrade reaction. The soil-cap interface is divided into triangular elements and the subgrade reaction is assumed to vary linearly across each element. The building\'s structure is modeled by FEM. Several soil structure interaction examples are studied in this thesis, and they show that the solutions obtained from program SSI are in good agreement with others 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

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.

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.

O tema desta pesquisa refere-se ao tratamento numérico de estruturas de pavimentos, dando ênfase à utilização do método dos elementos de contorno - MEC para o tratamento de elementos planos, enquanto os elementos lineares serão tratados através do método dos elementos finitos - MEF. Busca-se contribuir com uma série de trabalhos realizados nesta área, principalmente com a consideração do efeito da excentricidade do eixo neutro das barras em relação à superfície neutra da placa somando-o, portanto, ao fenômeno de flexão desta última. A técnica de acoplamento dos elementos estruturais utilizada foi a técnica das sub-regiões, que facilita a visualização do problema da combinação e viabiliza o uso da técnica de condensação estática na resolução do sistema de equações. / This research refers to the numeric analysis of structures, emphasizing the use of the boundary element method - BEM - in the discretization of the plates, while the linear elements wil be formulated through the finite element method - FEM. The aim is to contribute with a number of works related to this subject, specially in the consideration of the effect of the stiffeners eccentricity in relation to the plates neutral surface, in addition to its bending state. The tool adopted to fulfill the coupling among these kinds of structural elements is the well-known sub region or multi-domain technique, which works well with the static condensation technique for solving the sparse system of equation generated.

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.