O método dos elementos de contorno (MEC) surgiu como uma poderosa alternativa ao método dos elementos finitos (MEF) principalmente em casos como problemas de concentração de tensões ou onde o domínio se estende para o infinito. Em virtude das potencialidades já identificadas do MEC para a solução de problemas da geotecnia, em especial para problemas de túneis, este trabalho tem como objetivo desenvolver um programa que seja capaz de analisar as variáveis envolvidas na construção de túneis profundos através de um modelo numérico bidimensional baseado no MEC, implementando técnicas numéricas tais como: subelementação, técnica da sub-região e modelagem de inclusão e enrijecedores. O modelo numérico bidimensional foi calibrado para considerar o efeito tridimensional do problema de túneis no que se refere ao avanço da frente de escavação, para dois casos a saber: i) túneis sem suporte e ii) túneis com suporte. Os resultados mostraram grande precisão quando comparados com os resultados analíticos mesmo utilizando um número pequeno de elementos, provocando uma redução significativa no tempo de processamento se comparado com outros métodos. A técnica da subelementação produziu uma suavização nos resultados dos pontos internos localizados muito próximos do contorno. A técnica da sub-região, bem como a modelagem de inclusão e enrijecedores apresentaram resultados consistentes dando ao programa uma versatilidade maior. Na calibração dos parâmetros para a consideração do efeito tridimensional na escavação de túneis sem suporte, foi proposto o método da redução do carregamento com a construção do perfil de deformação longitudinal do túnel - LDP (Longitudinal Deformation Profile). Para a escavação de túneis com suporte foram propostos quatro métodos de análise: i) Método da redução do carregamento sobre o túnel, ii) Método da redução de rigidez do suporte, iii) Método do acréscimo do carregamento sobre o túnel e iv) Método do alívio de carga sobre o suporte. Todos esses métodos foram desenvolvidos a partir do modelo Kappa ($\kappa$), elaborado neste trabalho a partir dos resultados encontrados na simulação numérica tridimensional realizado nos programas BEFE e BEFE++, e comparado com o modelo de Schwartz e Einstein. Por fim, o método para a construção do gráfico de deslocamento radial para túneis circulares suportados, considerando o atraso na instalação do suporte, utilizando um método numérico ou resultado análitico do estado plano de deformação se mostra como uma alternativa simples para análise do efeito tridimensional contido no problema de túneis.

The boundary element method (BEM) has appeared as a powerful alternative to the finite element method (FEM) mainly in the cases where a good accuracy is required, as for problems with strain or stress concentration and problems with domain extending to infinite. The objective of this work is to develop a formulation and the corresponding computational code to analyse the variables in a design of deep tunnels, using a improved BEM two-dimensional numerical model, in which the following techniques were implemented: sub-elementation, sub-region technique, reinforcements introduced by modifying locally the domain rigidity. The two-dimensional model was calibrated to take into account the three-dimensional effects appearing around the tunnel face advance for two cases: i) tunnel without support and ii) tunnel with support. The results showed good accuracy when compared with analytical results even when obtained by using coarse discretizations and therefore requiring less computer time in comparison with other numerical procedures. The sub-elementation technique has smoothed the results for internal points near the boundary. The sub-region technique and the reinforcement inclusions lead to accurate making the computer code reliable. For the parameter calibration to take into account the three-dimensional effects applied to non lined tunnels the method of loading reductions was proposed obtaining a tunnel longitudinal deformation profile - LDP. For the excavation of lined tunnels four methods of analysis were proposed: i) load reduction model, ii) reduction support stiffness model, iii) additional load model, and iv) decrease of lining load model. All these methods were developed from the kappa ($\kappa$) model, developed in this work using three-dimensional results obtained by using the computational systems BEFE and BEFE++ and compared with the Schwartz-Einstein method. Finally the method used to build the radial displacement graphic for lined circular tunnels, taking into account the support insertion delay, using either a numerical method or plane strain analytical solutions, was developed.