Submitted by maria angelica Varella (angelica@sibi.ufrj.br) on 2017-12-15T16:47:23Z No. of bitstreams: 1 121849.pdf: 771845 bytes, checksum: e439d51023300ff6cfc6297b59c9fb17 (MD5) / Made available in DSpace on 2017-12-15T16:47:25Z (GMT). No. of bitstreams: 1 121849.pdf: 771845 bytes, checksum: e439d51023300ff6cfc6297b59c9fb17 (MD5) Previous issue date: 1969-03 / Introduzindo um jato de ar quente na secção de teste dum túnel de vento, no qual há um escoamento com turbulência gerada por uma grade, estudamos o efeito desta turbulência sobre o jato, medindo a distribuição de temperatura axial e transversalmente. Equacionando um balanço do fluxo de quantidade de movimento e de entalpia em secções normais do jato, e adotando a hipótese de similaridade dos perfis adimensionais e de axissimetria do jato, obtivemos uma expressão analítica - relacionando a sua espessura e quedas axiais de temperatura e velocidade. As medições diretas ainda nos levaram a aduzir grande influência da turbulência exterior na abertura do jato e no comprimento do seu núcleo, comparados com o caso de escoamento paralelo laminar. Usando a solução da equação da difusão da propriedade de uma fonte fixa e contínua, em ambiente móvel, calculamos a difusividade turbilhonar em vários pontos ao longo da secção de teste do túnel de vento. Observamos outrossim, uma interessante transição do valor da relação de difusividade turbilhonar térmica e de quantidade de movimento, de um valor próprio para jatos livres até um correspondente ao encontrado em esteiras, à medida que aumentávamos a velocidade do túnel de vento, fixando a do jato. Sugerimos assim, um estudo dessa transição que se obtém variando-se simplesmente a relação das velocidades do túnel de vento e do jato, sem necessidade de modificar a geometria.

Uma maneira de analisar a estabilidade de túneis escavados em maciços rochosos fraturados é a baseada na teoria da análise limite utilizando um critério de resistência macroscópico desenvolvido a partir da teoria de homogeneização. Este trabalho primeiramente descreve o comportamento mecânico dos constituintes do maciço rochoso fraturado e os métodos usuais de análise. Posteriormente é descrita a teoria da análise limite e da homogeneização em análise limite, utilizando estas teorias é apresentado o critério de resistência macroscópico para maciços rochosos fraturados em estado plano de deformações desenvolvido por Fréard (2000) e desenvolvido um critério de resistência macroscópico para maciços rochosos fraturados para o caso tridimensional. A partir destes critérios aplicaram-se mecanismos de ruptura para se encontrar limites superiores da função de estabilidade. Assim permitindo a análise da previsibilidade do colapso do túnel estação Pinheiros do metrô de São Paulo o qual os resultados apresentavam a possibilidade de prever o colapso. / One way of analyze the stability of tunnels excavated at fractured rock masses is based on the theory of yield design using a macroscopic strength criterion developed from the theory of homogenization. This dissertation first describes the mechanical behavior of the fractured rock mass constituents and the usual methods of analysis. Later is described the yield design theory and the theory of homogenization applied to the yield design using these theories is presented macroscopic strength criterion for fractured rock masses in a state plan developed by Fréard (2000) and developed a macroscopic strength criterion for rock masses fractured to the three-dimensional case. Based on these criteria were applied failure mechanisms to meet the upper bounding of the stability function. So it’s possible to analyze the predictability of the collapse of the tunnel-station Pinheiros of the Sao Paulo subway which the results showed the possibility of predicting the collapse.

O projeto estrutural de túneis envolve o conhecimento da variável convergência (fechamento radial) para definição da pressão atuante sobre o revestimento e do nível de descompressão do solo circundante (maciço). A adoção de modelos contínuos em elementos finitos consegue discretizar de forma mais adequada o modelo mecânico estrutural considerando a interação com o solo. Dessa forma, esta dissertação recai sobre a variação da razão de cobertura (profundidade) de um túnel, a fim de determinar a profundidade limite que o modelo numérico segue o padrão de tensões geostáticas crescentes com a profundidade versus um padrão geostático uniforme profundo. A validação do modelo numérico ocorre junto a soluções analíticas desenvolvidas por Park (2004) e Kirsch (1898) para túneis em maciços elásticos superficiais e profundos respectivamente. Nos exemplos analisados, quando a profundidade do eixo do túnel é dez vezes maior que seu diâmetro as convergências se aproximam daquelas dadas pelo GEOMEC 91 (Bernaud, 1991) e o túnel pode ser considerado profundo. A curva de assentamentos (recalques) foi considerada consequência importante a ser notada em túneis superficiais que mostraram um recalque máximo até duas vezes maior quando se reduz a profundidade Para esse estudo foram desenvolvidos scripts com profundidade variável adicionando a superfície de plasticidade de Drucker-Prager (utilizada em solos friccionais). O modelo numérico desenvolvido para túneis superficiais é validado junto a túneis reais instrumentados que tenham sido avaliados a partir de algum outro modelo contínuo, são eles: Heathrow Express, Green Park, Metrôs de Barcelona, Brasília e São Paulo com resultados variando entre 1,20% e 6,78%. Eventualmente, uma breve análise do desenvolvimento de tensões no revestimento para túneis circulares, elípticos e ferradura é realizada a fim de encontrar zonas de tração. Enquanto túneis em seção ferradura desenvolveram tração para qualquer profundidade devido ao formato retilíneo, túneis circulares precisaram se encontrar com baixa profundidade a fim de desenvolver zonas descomprimidas decorrentes da ovalização de túneis rasos. / The structural design of tunnels involves the knowledge of the variable convergence (radial closure) for the definition of the acting pressure upon the lining and for the decompression level of the involving soil (rockmass). The adoption of continuous models using finite element analysis (FEA) is able to assess the structural mechanical model accurately, concerning the soil interaction. This way, this work studies the covering ratio (depth) of tunnels, in a way to determine the threshold point which divides the geostatic growing stresses versus the deep uniform stresses pattern. The validation of the numerical model under an elastic rockmass is done by comparison with analytical solutions developed by Park (2004) and Kirsch (1898) for shallow and deep tunnels respectively. Regarding the analyzed examples, when the depth of the tunnel axis is ten times greater than its diameter; the convergences approach to the ones given by GEOMEC 91 (Bernaud, 1991) and the tunnel may be considered deep. The settlement through (ground movement) is considered an important consequence to be noticed upon shallow tunnels; these tunnels have led to a maximum settlement twice deeper when the covering ratio is shallow (H/D = 1). For this study, several scripts have been developed concerning the variable depth and adding the Drucker-Prager plasticity criterion for frictional soils. The numerical model developed for shallow tunnels is validated with other continuous models relied on real instrumentation of subway network tunnels, such as: Heathrow Express, Green Park, Barcelona, Brasília and São Paulo with the solutions ranging between 1,20% and 6,78%. Eventually, a brief analysis of the stress development in the lining for circular, ellipsoid and horseshoe tunnels is done in search of tension zones. In one hand, horseshoe tunnels develop tension for any depth due to its straight walls; in the other hand, circular tunnels need to have a small covering ratio in order to develop decompression zones caused by the shallow tunnels ovalization.

Este trabalho estuda os resultados obtidos da aplicação de um modelo tridimensional de previsão de comportamento de túneis profundos em rocha reforçados com tirantes passivos. Utiliza-se o método dos elementos finitos e simula-se a escavação pelo sistema de ativação e desativação de elementos. Os tirantes são tratados com uma lei de comportamento elastoplástica e o maciço reforçado com tirantes é considerado um meio homogêneo anisotrópico em escala macroscópica, através do método de homogeneização dos meios periódicos. Para a análise paramétrica foi verificada a influência dos parâmetros do sistema de reforço e do maciço. Resultados em termos de convergência do túnel e tensões são apresentados e indicam que, para o caso analisado, o comprimento dos tirantes é um fator de reduzida influência. A densidade de tirantes possui maior importância, porém, esta não é afetada por modificações na coesão e no módulo de elasticidade do maciço. A distância de aplicação dos tirantes à face de escavação se mostrou significativa somente para valores abaixo de 1/3 do raio do túnel, para valores acima, os deslocamentos já ocorreram e, portanto, não mobilizavam mais os tirantes. Verificou-se que a concentração de tensões σrr acontece na parte superior e inferior do túnel e, também, antes da frente de escavação. O ângulo de atrito demonstrou influência irrisória na convergência, diferentemente da coesão que diminui a convergência, tanto do teto como da parede do túnel, à medida que seu valor aumenta. O empuxo lateral exibiu ainda maior significância, revelando o quanto é importante sua correta determinação, especialmente na convergência lateral do túnel. Ainda, constatou-se a grande importância da combinação de tirantes com revestimento em concreto, que juntos, diminuíram em aproximadamente 25% os deslocamentos no teto do túnel. / This work studies the results obtained from the application of a three-dimensional model to predict the behavior of bolted reinforced rock mass deep tunnel. The finite element method is used to model and simulate the excavation through the activation and deactivation of elements. The bolts are treated with an elastoplastic law and the rock mass is considered as a homogenized anisotropic medium on the macroscopic scale through the homogenization method of periodic media. In the parametric analysis was studied the influence of the parameters of the reinforcement system and the rock mass. Results in terms of convergence of the tunnel and stresses are presented and indicate that, for the studied case, the length of the bolts is a factor of little influence. The density of bolts has a greater importance, however, its influence was not affected by changes in cohesion and elasticity modulus of the rock mass. The distance between the bolted zone and the facing of the tunnel was a significant parameter only for values below 1/3 of the tunnel’s radius, for distances greater than this, the displacements had already occurred and did not mobilize the bolts anymore. It was verified that stress concentration σrr occurs in the top and bottom of the tunnel, and also before the facing of the excavation. The friction angle demonstrated negligible influence on the convergence, unlike the cohesion, which decreases the convergence of the roof and tunnel’s wall as its value increases. The coefficient of lateral earth pressure at rest, k0, exhibited even greater significance, revealing how important its correct determination is, especially for the lateral convergence of the tunnel. Still, it was observed the great importance of the combination of bolts and concrete lining, which together decreased by approximately 25% the displacement of the tunnel’s roof.

Em túneis, a deformação do maciço e a pressão do solo no revestimento dependem das tensões e características do maciço, das tensões do maciço, da geometria, da rigidez, do instante da colocação do revestimento. A variação das tensões no revestimento e as deformações são causadas pelo avanço da escavação e as propriedades dependentes do tempo do maciço rochoso e revestimento. O objetivo desta dissertação é apresentar um modelo numérico para a análise tridimensional de túneis em maciços rochosos no programa ANSYS. Ele constitui uma ferramenta para calcular um túnel, revestido de concreto em diferentes situações de rigidez do revestimento, processo construtivo e características do maciço. O processo de escavação e colocação do revestimento é simulado pelo método da ativação/desativação de elementos em etapas sequenciais. O maciço apresenta comportamento elástico na primeira parte, e posteriormente uma lei viscoplástica. O concreto projetado e o concreto pré-moldado do revestimento são modelados como elásticos na primeira parte, e posteriormente, como materiais viscoelásticos utilizando o modelo da cadeia de Maxwell e de Kelvin, propostos por Bazant. Ambos os modelos foram calibrados por comparação com o Código Modelo fib 2010. A simulação do túnel é realizada através da incorporação de sub-rotinas para o concreto viscoelástico no programa ANSYS. Para validar o modelo são feitas comparações dos resultados obtidos no programa ANSYS com a solução analítica, com o programa GEOMEC91 e com o Novo Método Implícito, dependendo do caso. Um estudo paramétrico em elasticidade é realizado com o código ANSYS variando alguns parâmetros relevantes. Na segunda parte, o mesmo túnel revestido é modelado com dois diferentes tipos de revestimento e os dois modelos de cadeia citados. Finalmente, é comparado o comportamento com o túnel experimental Kielder, do qual se tem dados disponíveis da instrumentação. A comparação destes dados com os resultados do programa ANSYS mostra uma aproximação muito boa. / Stabilizing underground openings such as tunnels excavated in rock mass remain a major concern of geotechnical engineers dealing with this kind of structures. In tunnels, the rock mass strain and the ground pressure on lining depend on the stress and characteristics of the rock mass as well as of the geometry, the stiffness and the moment of the lining installation. Pressure variation on lining and strain are caused by the advance of excavation and the time-dependent properties of the rock mass and lining. This dissertation describes the numerical implementation of constitutive laws in the numerical code ANSYS. It is performed a numerical simulation with 3D finite elements of a tunnel. The rock mass presents elastic behavior in the first part and then a viscoplastic law. The shotcrete and precast concrete are modeled as elastic in the first part and then viscoelastic material with the constitutive laws of the Maxwell and the Kelvin chain model, proposed by Bazant. So the characteristics of the viscoelastic concrete changes with time until the stabilization. Both chain models had to be calibrated by comparing with de fib Model Code 2010. The finite elements simulation is performed by incorporating subroutines for the viscoelastic concrete model in the ANSYS code. The method to simulate the tunnel excavations is by activating and deactivating elements in sequential steps. In the first part of the dissertation two validations are performed. The analytical solution and the deformation achieved on the stabilization in the ANSYS code are compared for the validation of an unlined tunnel. For a lined tunnel, with an elastic and viscoplastic rock mass and elastic lining, validation is performed by comparing the results of the GEOMEC91 code with the ANSYS code. The results show a very good approximation in viscoplasticity and agree perfectly in elasticity. Also a parametric study is undertaken with the ANSYS code varying some relevant parameters. In the second part, it is compared the same tunnel with the two different concrete lining and the two chain models. Finally, it is compared with the Kielder experimental tunnel, which in situ measured data is available. The comparison of this data with the results from the ANSYS code shows a very good approximation.

No presente trabalho, é apresentado um modelo numérico para análise tridimensional de túneis escavados em maciços rochosos reforçados por tirantes aderidos ao longo de todo o seu comprimento. O método numérico utilizado é o método dos elementos finitos. A abordagem usada na modelagem da interação maciço/tirante é o método de homogeneização de meios periódicos. Essa abordagem substitui o meio compósito constituído de rocha e tirantes por um meio homogeneizado equivalente, dispensando a tarefa de discretizar cada tirante separadamente. O comportamento do maciço é considerado elastoplástico perfeito, caracterizado pelo critério de Drucker-Prager, com regra de fluxo associada. O revestimento é considerado elástico, enquanto os tirantes são tratados como inclusões lineares em tração/compressão com lei de comportamento elastoplástica. Com esse modelo foi implementado um código, tridimensional, que simula o processo de escavação pelo método da ativação/desativação de elementos. Esse código é capaz de levar em conta diferentes situações de profundidade, geometria e estado inicial de tensões. Com a finalidade de avaliar o desempenho da abordagem de homogeneização, foi implementado um segundo código, que modela a interação maciço/tirante por meio do método incorporado. Esse método, muito usado em análises de concreto armado, tem em comum com o método de homogeneização a capacidade de lidar com diferentes configurações de reforços utilizando a mesma malha base. Ao final, os dois códigos são utilizados em modelagens de um caso real, o túnel experimental Kielder (Inglaterra), para o qual se dispõe de dados medidos in situ. A comparação destes dados com os resultados das duas modelagens mostra uma aproximação muito boa. / A numerical model for three-dimensional analysis of tunnels excavated in reinforced rock mass by fully grouted bolts is presented in this work. The numeric method utilized is the finite element method and the approach used in the modeling of the interaction mass/bolt is the homogenization of periodic media method. This approach replaces the composite medium formed by rock and bolts by a homogenized equivalent medium eliminating the task of discretizing each bolt separately. The rockmass behavior is considered elastoplastic perfect caracterized by the Drucker-Prager criterion with associated flow rule. The lining is elastic, while the bolts are treated as linear inclusions in tension / compression with elastoplastic behavior law. A three-dimensional code was implemented with this model that simulates the excavation process by means of the activation-deactivation method. This code is able to consider different situations of depth, geometry and initial state of stress. In order to evaluate the performance of the homogenization approach, a second code was implemented to model the interaction rockmass/bolts by means of the embedded method. This method, widely used in analysis of reinforced concrete, has in common with the approach of the homogenization the capacity to deal with different configurations of reinforcements with the same base mesh. Finally, the two codes are used to model a real case, the Kielder experimental tunnel (England), for which in situ measured data are available. The comparison of these data with the results of the two modeling shows a very good approximation.

Novas tecnologias estão sendo desenvolvidas constantemente para coletar informações de superfícies ou de sólidos para diversas finalidades. Alguns métodos clássicos, como a Topografia e a Fotogrametria terrestre, com o passar dos anos, tiveram uma grande evolução. Na Fotogrametria terrestre todo o processo está sendo feito em meio digital. Na topografia, as estações totais automatizaram a medição de ângulos e distâncias. Essa evolução tecnológica fez com que os levantamentos se tornassem cada vez mais rápidos e precisos, aumentando a produtividade. O mais recente nessa evolução é o levantamento através do sistema de varredura a laser (Laser Scanner) 3D. São muitas as aplicações dessa tecnologia, dentre as quais pode-se citar: túneis, levantamento do como construído (as-built), mineração (principalmente subterrânea), arqueologia, levantamento de monumentos para restauração, refinarias e instalações industriais e outras, caracterizadas pela grande complexidade dos elementos envolvidos. A presente dissertação apresenta os conceitos envolvidos em todos os processos, desde a coleta de dados até o produto final. Desenvolve uma metodologia de uso que possa ser útil em diversas áreas, mostra uma aplicação completa na área de túneis, complementada por uma visão geral da área de plantas industriais e procura apresentar testes para quantificar a precisão que se obtém por essa tecnologia. / New technologies are constantly being developed in order to collect information of surfaces or solids for diverse purposes. Some classic methods such as topography and terrestrial photogrammetry have had a great evolution in the past. For example, all the processes of the terrestrial photogrammetry are made in digital way and the Total Stations have automated the measurements of angles and distances. This technical evolution made the surveying faster and accurate, increasing the productivity. However this evolution does not stop for there; in other words, the last technology in the area of topography is the surveying with the system known as Laser Scanner 3D. The Laser Scanner technology 3D has a lot of applications such as: tunnel, as-built, mining (mainly in the underground); archaeology (for restore monuments), refineries, industrial installations, etc., characterized by the great complexity of the involved elements. This work presents concepts involved in all the processes, since from data collection to the final product. It develops a methodology of use that can be applied in several areas, with emphasis in tunnels surveying area and presents some tests to quantization the accuracy.

A presente dissertação aborda o estudo dos efeitos de escavações realizadas acima de túneis já executados. Os trabalhos mais relevantes envolvendo a interação entre escavações e túneis já existentes são apresentados e discutidos. São realizadas simulações numéricas de casos hipotéticos, em estado plano de deformações utilizando o modelo de Mohr-Coulomb, permitindo a identificação dos parâmetros mais relevantes para esse tipo de problema e suas influências nos resultados. É discutido o comportamento mais rígido do maciço durante as trajetórias de descarregamento, bem como qual o módulo de elasticidade mais adequado para representar o fenômeno. Modelagens numéricas de um caso real em São Paulo (a escavação dos subsolos do edifício Stan Paulista, acima de túneis do Metrô) são realizadas, utilizando o modelo e os conceitos estudados. Os resultados obtidos das simulações são discutidos e comparados com os valores de instrumentação. Constata-se que a metodologia utilizada para a análise dos efeitos de escavações acima de túneis já existentes pode ser utilizada na avaliação de outros casos. / The presented research approaches the study of excavations effects above existing tunnels. Some of the most relevant papers related to interaction between excavation and tunnels are presented and discussed. Numerical simulations of hypothetic cases assuming plane strain condition and Mohr-Coulomb model for soil behavior are performed, allowing the study of influence of input parameters in the results. It is discussed behavior of the soil in stress paths that indicates reduction of the mean stress, wich is more rigid than the triaxial compression path. It is also discussed wich is an appropriate elastic modulus to be used in this situation. A real case (basement excavation of Stan Paulista building, above Metrô tunnels) is analyzed with numerical models, utilizing the concepts studied. The results of the simulations are compared with field instrumentation data. It is concluded that the methodology used for the analysis of the effect of excavation above existing tunnels can be utilized in other cases.

O projeto de submarinos é um tema que vem se desenvolvendo há muitos anos. No Brasil, vivemos um momento de pleno avanço na área, sendo encabeçado pela Marinha do Brasil com o desenvolvimento do primeiro submarino de propulsão nuclear brasileiro. No decorrer desse projeto, diversos estudos devem ser desenvolvidos, dentre eles, ensaios com modelos em escala, que serão realizados em um túnel de cavitação. Durante estes experimentos, espera-se que surjam diversas fontes de incerteza, que podem comprometer os resultados obtidos, destacando-se o efeito de blocagem hidrodinâmica. Nesse contexto, esta dissertação apresenta uma investigação do fenômeno de blocagem no modelo de submarino DARPA SUBOFF na seção de testes de um túnel de cavitação. Para isso, simulações tridimensionais de média de Reynolds (RANS) foram realizadas para o modelo de submarino posicionado no túnel (sob blocagem) e sob condição de águas profundas (referência), para diferentes ângulos de ataque e número de Reynolds (Re) constante. Coeficientes de força e momento hidrodinâmicos foram obtidos e os métodos de correção de blocagem clássico e WICS foram empregados para obter fatores que pudessem corrigir esses valores. A eficiência de cada método foi avaliada comparando os valores corrigidos com os de referência. Os resultados mostraram que ambos métodos de correção são eficientes na recuperação de coeficientes de força e momento hidrodinâmicos para razões de blocagem inferiores a 10%, mas ressalvas são feitas quanto à utilização de cada um desses métodos, sendo o método clássico indicado para o estudo de geometrias simples e o método WICS para modelos mais complexos. Paralelamente, a topologia do escoamento ao redor do submarino foi estudada através do método dos pontos críticos, para avaliar se a blocagem hidrodinâmica pode causar alterações nas estruturas presentes no escoamento. A análise topológica permitiu avaliar a evolução nas estruturas do escoamento para os ângulos de ataque e Re simulados. Foram capturadas as linhas de separação no corpo e avaliada sua relação com a tensão de cisalhamento ao longo do casco. A blocagem hidrodinâmica não causou alterações drásticas na topologia do escoamento. / The Project of submarines has been in development for years. In Brazil, we are facing a moment of full progress in this area, leaded by the Brazilian Navy with the development of the first submarine with nuclear propulsion in the country. During this project, many studies will be developed, among them, experiments with scaled models, that will be performed in a cavitation tunnel. Throughout these experiments, many sources of uncertainty are expected to arise that may compromise the results obtained, one effect that stands out is the hydrodinamic blockage. In this context, this master\'s thesis presents an investigation of the blockage phenomenon in the submarine model DARPA SUBOFF placed in the test section of a cavitation tunnel. To accomplish this, tridimensional simulations of Reynolds-Averaged Navier Stokes equations (RANS) were performed to the submarine model placed in the tunnel (under blockage) and in deep water condition (reference), to a variety of angles of attack and a constant Reynolds Number (Re). Cofficients of hydrodynamic forces and moment were measured and the classic and WICS correction methods were used to obtain factors that could correct these values. The efficiency of these methods were evaluated comparing the corrected values with the reference ones.The results showed that both correction methods are efficients in the recovery of hydrodynamic forces and moment to blockage ratios below 10%, but reservations are made about the use of each method. The classic method is suggested when simple geometries are tested and the WICS method when more complex models are used. At the same time, the flow topology around the submarine was studied, using the critical point theory, to evaluate if the hydrodynamic blockage can cause changes in the structures in the flow. The topology analysis permited to assess the evolution of structures in the flow to the conditions of angle of attack and Re simulated. Separation lines were captured in the body and their relationship with the shear stress in the hull was evaluated. The hydrodynamic blockage did not cause drastic changes in the flow topology.

Em túneis, a deformação do maciço e a pressão do solo no revestimento dependem das tensões e características do maciço, das tensões do maciço, da geometria, da rigidez, do instante da colocação do revestimento. A variação das tensões no revestimento e as deformações são causadas pelo avanço da escavação e as propriedades dependentes do tempo do maciço rochoso e revestimento. O objetivo desta dissertação é apresentar um modelo numérico para a análise tridimensional de túneis em maciços rochosos no programa ANSYS. Ele constitui uma ferramenta para calcular um túnel, revestido de concreto em diferentes situações de rigidez do revestimento, processo construtivo e características do maciço. O processo de escavação e colocação do revestimento é simulado pelo método da ativação/desativação de elementos em etapas sequenciais. O maciço apresenta comportamento elástico na primeira parte, e posteriormente uma lei viscoplástica. O concreto projetado e o concreto pré-moldado do revestimento são modelados como elásticos na primeira parte, e posteriormente, como materiais viscoelásticos utilizando o modelo da cadeia de Maxwell e de Kelvin, propostos por Bazant. Ambos os modelos foram calibrados por comparação com o Código Modelo fib 2010. A simulação do túnel é realizada através da incorporação de sub-rotinas para o concreto viscoelástico no programa ANSYS. Para validar o modelo são feitas comparações dos resultados obtidos no programa ANSYS com a solução analítica, com o programa GEOMEC91 e com o Novo Método Implícito, dependendo do caso. Um estudo paramétrico em elasticidade é realizado com o código ANSYS variando alguns parâmetros relevantes. Na segunda parte, o mesmo túnel revestido é modelado com dois diferentes tipos de revestimento e os dois modelos de cadeia citados. Finalmente, é comparado o comportamento com o túnel experimental Kielder, do qual se tem dados disponíveis da instrumentação. A comparação destes dados com os resultados do programa ANSYS mostra uma aproximação muito boa. / Stabilizing underground openings such as tunnels excavated in rock mass remain a major concern of geotechnical engineers dealing with this kind of structures. In tunnels, the rock mass strain and the ground pressure on lining depend on the stress and characteristics of the rock mass as well as of the geometry, the stiffness and the moment of the lining installation. Pressure variation on lining and strain are caused by the advance of excavation and the time-dependent properties of the rock mass and lining. This dissertation describes the numerical implementation of constitutive laws in the numerical code ANSYS. It is performed a numerical simulation with 3D finite elements of a tunnel. The rock mass presents elastic behavior in the first part and then a viscoplastic law. The shotcrete and precast concrete are modeled as elastic in the first part and then viscoelastic material with the constitutive laws of the Maxwell and the Kelvin chain model, proposed by Bazant. So the characteristics of the viscoelastic concrete changes with time until the stabilization. Both chain models had to be calibrated by comparing with de fib Model Code 2010. The finite elements simulation is performed by incorporating subroutines for the viscoelastic concrete model in the ANSYS code. The method to simulate the tunnel excavations is by activating and deactivating elements in sequential steps. In the first part of the dissertation two validations are performed. The analytical solution and the deformation achieved on the stabilization in the ANSYS code are compared for the validation of an unlined tunnel. For a lined tunnel, with an elastic and viscoplastic rock mass and elastic lining, validation is performed by comparing the results of the GEOMEC91 code with the ANSYS code. The results show a very good approximation in viscoplasticity and agree perfectly in elasticity. Also a parametric study is undertaken with the ANSYS code varying some relevant parameters. In the second part, it is compared the same tunnel with the two different concrete lining and the two chain models. Finally, it is compared with the Kielder experimental tunnel, which in situ measured data is available. The comparison of this data with the results from the ANSYS code shows a very good approximation.