Non linear geometric analysis of soil-structure interaction via BEM/FEM coupling

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