A Costa Peruana apresenta alta atividade sísmica, tornando imprescindível a execução de análises de estabilidade que considerem os eventos sísmicos. Com o desenvolvimento de novas ferramentas numéricas, a análise dinâmica está se tornando cada vez mais importante e usual na fase de projeto, não se justificando mais a execução de análises estáticas em locais tão vulneráveis a sismos. A presente dissertação tem como objetivo realizar a análise da estabilidade dos Penhascos de Lima, Peru, considerando quatro taludes distintos, em uma região bastante afetada por abalos sísmicos. As análises foram executadas pelos métodos de equilíbrio limite, elementos finitos e pseudo-estático, buscando-se uma comparação entre os diferentes métodos. O trabalho apresenta uma descrição completa dos taludes em estudo, uma avaliação da condição sismológica da região e finalmente a estabilidade dos Penhascos de Lima, fazendo uso dos programas computacionais Slide (método do equilíbrio limite) e Plaxis (método dos elementos finitos) em 2D. Os resultados mostraram que os três métodos adotados forneceram fatores de segurança compatíveis, principalmente quando se considera perfis menos estratificados. Para perfis homogêneos, as diferenças obtidas foram da ordem de 0,5 a 1,0 %. As análises ressaltaram a importância de considerar a condição dinâmica, e mostraram-se bastante sensíveis aos valores de carga sísmica adotado. / The Peruvian coast has high seismic activity, making essential the implementation of stability analyzes that consider the seismic events. With the development of new numerical tools, dynamic analysis is becoming increasingly important and usual in the design phase is not justified over the execution of static analysis in places as vulnerable to earthquakes. This thesis aims to conduct the analysis of the stability of the cliffs of Lima, Peru, considering for different slopes, in a region highly affected by earthquakes. The analyzes were performed by equilibrium methods limit finite elements and pseudo-static, seeking a comparison between the different methods. The dissertation presents a complete description of the slopes in the study, an evaluation of the seismic condition of the region and finally the stability of Lima cliffs, making use of computer programs Slide (limit equilibrium method) and Plaxis (finite element method). The results showed that all three methods provided adopted gave safety factors, especially when considering less laminated profiles. For homogeneous profiles, obtained differences were of the order of 1.0 at 0.5%. The analysis highlighted the importance of considering the dynamic condition, and were quite sensitive to seismic load values adopted.

Main goal of thesis entitled "Landslide remediation design" is - based on acquired knowledge and submitted documents about the problematic of landslides - to determine the causes of the landslides on II. class road between villages Koryčany and Jestřabice and to propose optimal redevelopment.

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior / O perfeito entendimento de deslizamentos de terra tem sido objeto de grande interesse na engenharia geotécnica há muitos anos. Encostas naturais representam um desafio constante para os engenheiros, mais do que os taludes artificiais. Não é incomum que os escorregamentos causam perdas de vidas humanas e também danos materiais. As informações obtidas a partir de retro-análise de casos históricos, têm servido como uma ferramenta valiosa para a compreensão da estabilidade de taludes. Vale a pena notar que o primeiro método de cálculo de estabilidade de taludes, usando divisão de fatia, foi desenvolvido após um acidente na Noruega. Métodos de estabilidade por Equilíbrio Limite, desenvolvidos no século passado são atualmente utilizados na prática da engenharia. Apesar de suas deficiências, eles provaram ser confiáveis na determinação do Fator de Segurança. A geometria 2D e a abordagem completa da condição saturada são provavelmente as questões mais relevantes dos métodos atuais. A condição do solo não saturado está geralmente presente em taludes naturais e os métodos ignoram a influência da sucção do solo sobre a estabilidade do talude. O desenvolvimento de computadores no meio do século passado, tornaram os cálculos Equilíbrio Limite fáceis de serem executados. Por outro lado, eles não podem explicar a natureza do problema ou prever o comportamento do solo antes da ruptura. A década de 1960 marcou o início do desenvolvimento de programas de computador baseados em métodos numéricos que, entre outros, destaca-se, em geotecnia, a aplicação dos métodos de Elementos Finitos (FEM). Este método é ideal para muitos fins, uma vez que permite a simulação da sequência de construção e também a incorporação de diferentes modelos constitutivos. Este trabalho mostra a influência da sucção sobre a estabilidade de encostas naturais. Uma retroanálise de um escorregamento no Morro dos Cabritos, no Rio de Janeiro, foi feita através de um programa de análise de tensão (Plaxis 2D). Os resultados foram comparados com o FS obtido por programa de Equilíbrio Limite (Slide v5). O efeito de sucção do solo foi incorporada nas análises, através de sua influência sobre os parâmetros de resistência do solo e por diferentes possibilidades de infiltração de água na encosta. Este caso histórico foi estudado anteriormente por Gerscovich (1994), que na época, fez analise de fluxo e estabilidade utilizando um abordagem de Equilíbrio Limite 3D. Os resultados mostraram boa concordância entre a análise de Elementos Finitos e a simulação de Equilíbrio Limite, sem a presença de nível de água. Os resultados confirmaram que foi necessário haver pressões positivas para explicar a ruptura do talude.

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior / O perfeito entendimento de deslizamentos de terra tem sido objeto de grande interesse na engenharia geotécnica há muitos anos. Encostas naturais representam um desafio constante para os engenheiros, mais do que os taludes artificiais. Não é incomum que os escorregamentos causam perdas de vidas humanas e também danos materiais. As informações obtidas a partir de retro-análise de casos históricos, têm servido como uma ferramenta valiosa para a compreensão da estabilidade de taludes. Vale a pena notar que o primeiro método de cálculo de estabilidade de taludes, usando divisão de fatia, foi desenvolvido após um acidente na Noruega. Métodos de estabilidade por Equilíbrio Limite, desenvolvidos no século passado são atualmente utilizados na prática da engenharia. Apesar de suas deficiências, eles provaram ser confiáveis na determinação do Fator de Segurança. A geometria 2D e a abordagem completa da condição saturada são provavelmente as questões mais relevantes dos métodos atuais. A condição do solo não saturado está geralmente presente em taludes naturais e os métodos ignoram a influência da sucção do solo sobre a estabilidade do talude. O desenvolvimento de computadores no meio do século passado, tornaram os cálculos Equilíbrio Limite fáceis de serem executados. Por outro lado, eles não podem explicar a natureza do problema ou prever o comportamento do solo antes da ruptura. A década de 1960 marcou o início do desenvolvimento de programas de computador baseados em métodos numéricos que, entre outros, destaca-se, em geotecnia, a aplicação dos métodos de Elementos Finitos (FEM). Este método é ideal para muitos fins, uma vez que permite a simulação da sequência de construção e também a incorporação de diferentes modelos constitutivos. Este trabalho mostra a influência da sucção sobre a estabilidade de encostas naturais. Uma retroanálise de um escorregamento no Morro dos Cabritos, no Rio de Janeiro, foi feita através de um programa de análise de tensão (Plaxis 2D). Os resultados foram comparados com o FS obtido por programa de Equilíbrio Limite (Slide v5). O efeito de sucção do solo foi incorporada nas análises, através de sua influência sobre os parâmetros de resistência do solo e por diferentes possibilidades de infiltração de água na encosta. Este caso histórico foi estudado anteriormente por Gerscovich (1994), que na época, fez analise de fluxo e estabilidade utilizando um abordagem de Equilíbrio Limite 3D. Os resultados mostraram boa concordância entre a análise de Elementos Finitos e a simulação de Equilíbrio Limite, sem a presença de nível de água. Os resultados confirmaram que foi necessário haver pressões positivas para explicar a ruptura do talude.

Submitted by Franciele Moreira (francielemoreyra@gmail.com) on 2018-09-21T14:50:47Z No. of bitstreams: 2 Dissertação - Nícolas Rodrigues Moura - 2018.pdf: 4165797 bytes, checksum: 86460ffa29b4d346005f977c2f75693e (MD5) license_rdf: 0 bytes, checksum: d41d8cd98f00b204e9800998ecf8427e (MD5) / Approved for entry into archive by Luciana Ferreira (lucgeral@gmail.com) on 2018-09-24T13:21:40Z (GMT) No. of bitstreams: 2 Dissertação - Nícolas Rodrigues Moura - 2018.pdf: 4165797 bytes, checksum: 86460ffa29b4d346005f977c2f75693e (MD5) license_rdf: 0 bytes, checksum: d41d8cd98f00b204e9800998ecf8427e (MD5) / Made available in DSpace on 2018-09-24T13:21:40Z (GMT). No. of bitstreams: 2 Dissertação - Nícolas Rodrigues Moura - 2018.pdf: 4165797 bytes, checksum: 86460ffa29b4d346005f977c2f75693e (MD5) license_rdf: 0 bytes, checksum: d41d8cd98f00b204e9800998ecf8427e (MD5) Previous issue date: 2018-08-27 / Coordenação de Aperfeiçoamento de Pessoal de Nível Superior - CAPES / Análises de estabilidade de taludes são predominantemente realizadas em duas-dimensões (2D) por métodos de equilíbrio limite. Os métodos convencionais de equilíbrio limite são alvos de muitas críticas pelo o uso do equilíbrio estático, pressuposições de forças na análise e por possuir o mesmo fator de segurança ao longo de toda a superfície de deslizamento. Para superar essas deficiências, o método proposto por Kulhawy (1969) incorpora o comportamento tensão-deformação do solo com tensões oriundas do Método de Elementos Finitos (FEM) e tem a mesma estrutura dos métodos de equilíbrio dos métodos de equilíbrio de equilíbrio limite clássicos. É também chamado de método aperfeiçoado. Devido à existência de taludes com geometrias complexas, a abordagem tridimensional (3D) é a única capaz de incorporar todos os aspectos geométricos do problema. Esta dissertação é dividida em dois artigos com o tema de estabilidade de taludes pelo o método de Kulhawy em abordagens 2D e 3D. O objetivo primário dessa dissertação é avaliar os principais fatores que influenciam a estabilidade de taludes. O primeiro artigo propôs geometrias variando de convexo a côncavo em um conjunto de propriedades típicas de solos. Duas trajetórias de geração de tensões in- situ foram utilizadas, i.e., construção e escavação. O método GLE foi calculado para fins de comparação. A superfície de deslizamento é muito dependente do histórico de tensão e da forma da face do talude. A diferença relativa no estado de tensões inicial típico alcançou 15% nas análises 2D e 34% nas análises 3D. A análise de sensibilidade do efeito do Poisson no fator de segurança mostrou-se muito relevante. A relação entre os fatores de segurança 2D e 3D são quase lineares independente do estado de tensões, das propriedades de resistência, formato da face do talude, na ordem de 30 a 40%. O segundo artigo propõe análises em barragens construídas em vales estreitos. Três abordagens de análises foram utilizadas: redução de resistência ao cisalhamento (SSR), métodos de Kulhawy e GLE. Condições típicas de fundações são muito variáveis como topografia e deformabilidade. A análise de sensibilidade do módulo de fundação de Young no fator de segurança foi insignificante em 2D e relevantes em 3D. A relação entre o comprimento do aterro (W) e a altura () é uma variável representativa do grau de estreitamento da barragem e na estabilidade de taludes. A análise de deformações planas tende a altas relações W/H e o oposto mostra elevados efeitos 3D. A inclinação da ombreira mostrou-se relevante no fator de segurança apenas em aterros em vales estreitos. Fundações mais rígidas são mais beneficias para o fator de segurança devido à produção de altos efeitos de arqueamento neste cenário. O método de Kulhawy provou-se ser vantajoso pela a incorporação do comportamento mecânico do solo, o uso de tensões mais realísticas nas análises e avaliação da estabilidade local na superfície de deslizamento. Isso aumenta a quantidade de informações úteis e confiáveis para o profissional de engenharia no projeto de taludes. Todas as análises foram realizadas usando o software SoilVision Systems, Saskatchewan, Canadá. Os módulos utilizados foram o SVLOPE, SVSOLID e o SVFLUX (SoilVision Systems Ltd., 2018). A análise de tensão para a simulação de trajetórias de escavação e construção utilizou do solucionador de Equações Diferenciais Parcias (PDE) o FlexPDE 6 (PDE Solutions Inc., 2014). / Slope stability analyses are predominantly performed using the two-dimensional (2D) conventional limit equilibrium methods. Conventional limit equilibrium methods are targets of many criticisms for using static equilibrium, assumptions of forces in the analysis and the same factor of safety along the slip surface. In order to overcome these shortcomings, the method proposed by Kulhawy (1969) incorporates the stress-strain behavior of the soil with stresses from the Finite Element Method (FEM) with the same structure as limit equilibrium methods. It is also called enhanced method. Due to the existence of slopes with complex geometries, the three-dimensional (3D) approach is only able to incorporate all the geometric aspects of the problem. This thesis is divided into two papers with the subject of slope stability analyses by the Kulhawy’s method in 2D and 3D approaches. The primary objective of this thesis is to evaluate major factors that influence slope stability. The first paper proposed geometries ranging from convex to concave in a set of soil properties. Two stress paths were considered, i.e., construction and excavation. The GLE method was computed for comparison purposes. The slip surface is very dependent on the stress history and the shape of slope face. The relative difference in typical initial stress state reached 15% in 2D analyses and 34% in 3D analyses. Sensitivity analysis of effect of Poisson’s ratio on factor of safety was very relevant. The relationship between 2D and 3D factors of safety are nearly linear independent of stress state on order of 30 to 40%. The second paper propose analyses on embankment dams built in narrow valleys. Three approaches of slope stability were used: shear strength reduction (SSR), Kulhawy’s and GLE methods. The typical conditions of foundation sites are very variable as topography and deformability. Sensitivity analysis of Young’s modulus of foundation on factor of safety was insignificant in 2D and relevant in 3D analyses. The ratio between embankment length (W) and height (H) was very important on slope stability. The plane strain analysis tends to high W/H ratio and the opposite shows high 3D effects. The inclination of abutment showed to be relevant on factor of safety only in embankments in narrow valleys. Stiffer foundations are beneficial to the factor of safety due producing high arching effects in this scenario. The Kulhawy’s method proved to be advantageous by the incorporation of mechanical behavior of soil, more realistic stress and local stability evaluation on slip surface. It increases the amount of useful and reliable information to the engineering professional in design of slopes. All analyses were performed using SoilVision Systems suite software, Saskatchewan, Canada. The modules used were SVSLOPE, SVSOLID and SVFLUX (SoilVision Systems Ltd., 2018). The stress analysis for the simulation of excavation and construction trajectories used partial differential equation (PDE) solver, FlexPDE 6 (PDE Solutions Inc., 2014).

La stabilité des exploitations à ciel ouvert (mines ou carrières) excavées dans des massifs rocheux dépend de leur géométrie, des caractéristiques géométriques de la fracturation (orientation et espacement) du massif et des caractéristiques mécaniques du massif et des discontinuités. L'évaluation du risque d'instabilité rocheuse dépend de la quantité et la qualité des informations disponibles sur le massif rocheux et des méthodes d'analyse mises en œuvre pour évaluer le comportement mécanique de l'assemblage de blocs rocheux. Différents sites de talus naturels, d'excavation routière et de carrières et mines à ciel ouvert présentés dans le mémoire montrent que les informations connues peuvent beaucoup varier d'un site à l'autre. Les étapes de collecte et de synthèse des données, puis de modélisation afin d'analyser la stabilité d'une pente rocheuse conduisent à des incertitudes. Le mémoire traite plus particulièrement de l'influence des incertitudes géométriques et différentes méthodes de regroupement des fractures en familles sont explorées. La combinaison du travail proposé par différents auteurs conduit à une nouvelle méthode nommée PSMY. Cette méthode ainsi que la méthode spectrale ont été programmées sous l'environnement Mathematica et les résultats sont comparés à des regroupements "manuels" réalisé avec le logiciel DIPS. Les méthodes de regroupement sont présentées selon le pourcentage d'orientations classé par la méthode. L’orientation et l’espacement des familles de fractures sont ajustées à des lois statistiques. Les paramètres statistiques sont comparés en fonction de la méthode regroupement. Ces paramètres influencent la construction d'un modèle géométrique du massif rocheux appelé DFN (Discrete Fracture Network). L'influence des regroupements sur l'analyse de stabilité est étudiée à partir de modélisation stochastiques à l'équilibre limite utilisant les logiciels SWEDGE et RESOBLOK. Ces logiciels ne prennent pas en compte les mêmes incertitudes et en sortie ils fournissent des indicateurs de stabilité différents. Une analyse de sensibilité des indicateurs (nombre de blocs instables, volume moyen du bloc instable, volume total instable) aux méthodes de regroupement et aux orientations de talus est réalisée. Une analyse de variance permet de préciser l'influence de ces deux facteurs. Une évaluation de l’état de stabilité global du massif, en fonction de la cohésion et de l'angle de frottement est proposée. L'influence de la méthode de modélisation est évoquée en comparant des calculs tridimensionnels à l'équilibre limite et des modélisations bidimensionnel et tridimensionnels en éléments discrets rigides ou déformables. Un couplage entre le logiciel RESOBLOK (équilibre limite) et LMGC90 (éléments discrets) permet de comparer les analyses pour une même géométrie. Pour le cas de la déviation d'Ax-les-Thermes et certaines configurations géométriques instables l’influence des paramètres de modèles sur stabilité du talus est testée. Plusieurs variantes sont comparées. La simulation d’une excavation en 3D multi phases est réalisée et l’indice de mobilisation est étudié pour comparer les divers types de contacts dans LMGC90 par rapport à l’éventuel glissement des blocs dans les étapes successives / The stability of open-cast operations (mines or quarries), excavated in rock mass depends on their geometry, the geometrical characteristics of the rock mass fractures (orientation and spacing) and the mechanical characteristics of the rock mass and the discontinuities. The assessment of the rock instability risk depends on the quantity and quality of the available information on the rock mass and the analysis methods used for the evaluation of the mechanical behavior of the rock block collection. Different sites of natural rock slopes, of road cuts, and of open-cast mines and quarries, described in the present document, show that the known information can vary a lot from one site to the other. The steps of data gathering and analysis, then the modelling step used to analyze the stability of the rock slope, lead to uncertainties. The present work deals in particular with the influence of the geometrical uncertainties, and the different clustering methods, to define families of fractures, are examined. The combination of works suggested by several authors leads to a new method called PSMY. This method, together with the spectral method, was coded in the Mathematica platform, and the obtained results are compared with "hand-made" clusters, done with the DIPS software. The clustering methods are presented according to the ratio of classified orientations. The orientation and spacing of fracture families are fitted by statistical law. The statistical parameters are compared according to the clustering method used. These parameters have an influence on the construction of the rock mass geometrical model, called DFN (Discrete Fracture Network). The influence of the clusters on the stability analysis is studied from the stochastic models based on the limit equilibrium analysis, in the SWEDGE and RESOBLOK software. These programs do not take into account the same uncertainties, and as a result, they give different stability indicators. A sensitivity analysis of these indicators (number of unstable blocks, average volume of unstable blocks, and total volume of unstable blocks) versus the clustering methods used, and the orientation of the slopes, is carried out. A variance analysis allows an evaluation of the influence of these factors. The assessment of the global stability condition of the rock mass, depending on the cohesion and the friction angle, is proposed. The influence of the modeling method is analyzed by comparing 3D calculations using the limit equilibrium calculations, and 2D and 3D models using discrete rigid and deformable blocks. A coupling between RESOBLOK (limit equilibrium) and LMGC90 (discrete elements) allows the comparison of results on the same original geometry. For the case of Ax-les-Thermes road-cut, and for various unstable geometries, the influence of the model parameters is tested. Several cases are compared. The 3D simulation of an excavation, at different steps, is performed, and the mobilization index is studied, in order to compare several types of contacts, within the LMGC model, in relation to the possible sliding of blocks, at those different steps

This document is focused on the elucidation of embankment deformation nearby žst. Valašská Polanka, with subsequent proposal of remediation measures, for which the own software and commercial software Plaxis 2D were used.

[pt] O principal objetivo deste trabalho é a aplicação do método dos elementos finitos (MEF) na análise da estabilidade de taludes em solo, para identificação da superfície crítica de deslizamento e cálculo do correspondente fator de segurança. Ainda que o método dos elementos finitos tenha se tornado ferramenta popular para análise de muitos problemas da engenharia geotécnica, existem relativamente poucas aplicações específicas publicadas na literatura sobre estabilidade de taludes. Dois métodos podem ser empregados para a utilização do MEF na área da estabilidade de taludes: o método direto, onde os resultados numéricos são exclusivamente usados para localização da potencial superfície de deslizamento e cálculo do fator de segurança, e o método indireto, onde seus resultados são acoplados a uma técnica para auxiliar a pesquisa da superfície crítica. Comparações entre resultados obtidos pelo método dos elementos finitos e de equilíbrio limite são feitas em vários exemplos numéricos investigados nesta dissertação, comentando-se a aplicabilidade, vantagens e limitações de ambas as metodologias em estudos da estabilidade de taludes em solos. / [en] The main objective of this research is the application of the finite element method (FEM) to analyze the stability of soil slopes, searching the position of the failure slip surface in the soil mass and its corresponding safety factor. The finite element method has become a popular tool for the study of many geotechnical engineering problems but, in the specific area of slope stability analysis, there are relatively few applications reported in the literature. Basically, two methods have evolved for application of the FEM results in slope stability problems: the direct method, that uses the numerical results only for establishing the slip surface and computation of the related safety factor, and the indirect method, where the FEM results are combined with another numerical technique to improve the search process of slip surfaces. Several numerical examples have been studied in this dissertation, making comparisons among the results obtained using the finite element and limit equilibrium methods. The applicability, advantages and shortcomings of both approaches are also discussed. / [es] 1276 El objetivo principal de este trabajo es la aplicación del método de los elementos finitos (MEF) al análisis de estabilidad de taludes de suelo, con el próposito de identificar la superfície crítica de deslizamento y cálculo del correspondiente factor de seguridad. Aunque el método de los elementos finitos se haya convertido en una herramienta popular en el análisis de muchos problemas de la ingeniería geotécnica, existen relativamente pocas aplicaciones específicas publicadas en la literatura sobre la estabilidad de taludes. Dos métodos pueden ser empleados para la utilización del MEF en el área de estabilidad de taludes: el método directo, donde los resultados numéricos son exclusivamente usados para localización de la potencial superficie de deslizamento y el cálculo del factor de seguridad; y el método indirecto, donde los resultados se acoplan a una técnica para auxiliar la investigación de la superficie crítica. Se utilizan varios ejemplos con el fin de comparar los resultados obtenidos por el método de los elementos finitos y por el método de equilibrio límite investigados en esta disertación. Se incluyen comentarios sobre la aplicabilidad, ventajas y limitaciones de ambas metodologías en estudios de la estabilidad de taludes en suelos.

The objective of this thesis was to identify the mechanisms associated with the recurrence of rock-slope instability along the R518 and R523 roads in Limpopo. Advanced geotechnical methods and ASTER imagery were used for the purpose while a predictive rockfall hazard rating matrix chart and rock slope stability charts for unsaturated sensitive clay soil and rock slopes were to be developed. The influence of extreme rainfall on the slope stability of the sensitive clay soil was also evaluated. To achieve the above, field observations, geological mapping, kinematic analysis, and limit equilibrium were performed. The latter involved toppling, transitional and rotational analyses. Numerical simulation was finally resorted to. The following software packages were employed: SWEDGE, SLIDE, RocData, RocFall, DIPS, RocPlane, and Phase 2. The simulation outputs were analyzed in conjunction with ASTER images. The advanced remote sensing data paved the way for landslide susceptibility analysis. From all the above, rockfall hazard prediction charts and slope stability prediction charts were developed. Several factors were also shown by numerical simulation to influence slope instability in the area of study, i.e. sites along the R518 and R523 roads in the Thulamela Municipality. The most important factors are extreme rainfall, steep slopes, geological features and water streams in the region, and improper road construction. Owing to the complexity of the failure mechanisms in the study area, it was concluded that both slope stability prediction charts and rock hazard matrix charts are very useful. They indeed enable one to characterize slope instability in sensitive clay soils as well as rockfall hazards in the study area. It is however recommended that future work is undertaken to explore the use of sophisticated and scientific methods. This is instrumental in the development of predictive tools for rock deformation and displacement in landslide events. / Electrical and Mining Engineering / D. Phil. (Mining Engineering)

Η παρούσα διατριβή επανεξετάζει την οριακή ανάλυση ευστάθειας των λίθινων αψίδων. Η οριακή ανάλυση ευστάθειας χρησιμοποιείται σήμερα ως το βασικό εργαλείο αποτίμησης της ευστάθειας τόξων και θολωτών κατασκευών από τοιχοποιία, όπως ακριβώς συνέβαινε και τους τελευταίους τέσσερις αιώνες. Παρά την τόσο μακρόχρονη ιστορία της μεθόδου, δεν έχουν πλήρως διασαφηνιστεί στην επιστημονική κοινότητα θεμελιώδης έννοιες και δεν έχουν σαφώς απαντηθεί ερωτήματα όπως: Ποιες είναι οι φυσικά πραγματοποιήσιμες γραμμές ώθησης και ποιες όχι; Ποια είναι η επίδραση της στερεοτομίας ενός τόξου στην οριακή του ευστάθεια; Ποιος είναι ο ρόλος της αλυσοειδούς καμπύλης και κατά πόσο αυτή είναι μία φυσικά αποδεκτή γραμμή ώθησης; Τι σχέση υπάρχει ανάμεσα στην κλίση της συνισταμένης θλιπτικής δύναμης και στην κλίση της γραμμής ώθησης στο σημείο εφαρμογής της; Η παρούσα διατριβή αναζητά απαντήσεις στα ερωτήματα αυτά, και έχει ως στόχο τη βαθύτερη κατανόηση της οριακής ανάλυσης ευστάθειας των τόξων, με παράλληλη ανάδειξη νέων υπολογιστικών διαδικασιών. Η δομή της παρουσιάζεται συνοπτικά παρακάτω. Στο πρώτο κεφάλαιο γίνεται ιστορική ανάλυση της μεθόδου μέσα από παρουσίαση και σχολιασμό των εργασιών με τη σημαντικότερη συμβολή, από τα μέσα του 17ου αιώνα μέχρι σήμερα. Στο δεύτερο κεφάλαιο επανεξετάζεται ένα από τα πιο κλασικά προβλήματα της μηχανικής: ποιο είναι το ελάχιστο επιτρεπτό πάχος ενός ημικυκλικού τόξου υπό τη δράση του ιδίου βάρους του για να είναι ευσταθές. Παράλληλα απαντώνται τα ερωτήματα που τέθηκαν παραπάνω αναπτύσσοντας νέες κλειστές μαθηματικές εκφράσεις των γραμμών ώθησης μέσω γεωμετρικής προσέγγισης, αλλά και μέσω του λογισμού των μεταβολών. Στο τρίτο κεφάλαιο χρησιμοποιείται παρόμοια διαδικασία για την ανάλυση της γενικής περίπτωσης των ελλειπτικών τόξων, οποιουδήποτε γεωμετρικού λόγου ύψος προς βάση, καθώς δεν είναι διαθέσιμα αναλυτικά αποτελέσματα στη διεθνή βιβλιογραφία, όπως συμβαίνει για τα κυκλικά τόξα. Στο τέταρτο κεφάλαιο εξετάζεται η οριακή ευστάθεια κυκλικών τόξων οποιασδήποτε γωνίας εναγκαλισμού, υπό την ταυτόχρονη δράση του ιδίου βάρους τους και σταθερής οριζόντιας εδαφικής επιτάχυνσης, ενώ υπολογίζεται με ακρίβεια η μορφή που θα έχει ο επικείμενος μηχανισμός κατάρρευσης μαζί με το οριακό πάχος, συναρτήσει της σεισμικής φόρτισης. Τα αποτελέσματα της μαθηματικής ανάλυσης (Κεφ. 2-4) επιβεβαιώνουν την ακρίβεια του λογισμικού που αναπτύχθηκε για τις ανάγκες της διατριβής, καθώς και τα αποτελέσματα που προκύπτουν από εμπορικό λογισμικό της μεθόδου των διακριτών στοιχείων. Στο πέμπτο κεφάλαιο γίνεται εφαρμογή και σύγκριση των πιο αντιπροσωπευτικών υπολογιστικών μεθόδων που απαντώνται σήμερα στη βιβλιογραφία για την αποτίμηση της ευστάθειας και φέρουσας ικανότητας της υπόγειας Θολωτής Διόδου του Σταδίου της Αρχαίας Νεμέας, ενώ η οριακή ανάλυση ευστάθειας αναδεικνύεται ως ένα μοναδικό εργαλείο για την κατανόηση της αλληλεπίδρασης της κατασκευής με το περιβάλλον έδαφος. Επιπλέων των συμπερασμάτων στο τέλος κάθε κεφαλαίου (Κεφ. 2 έως 5), στο έκτο κεφάλαιο παρουσιάζονται τα πιο σημαντικά συμπεράσματα και η συνεισφορά της παρούσας διατριβής. / This doctoral thesis revisits the limit equilibrium analysis of masonry arches. Limit equilibrium analysis is used today as the main analysis method for the assessment of the stability of masonry arches and vaulted structures, and is the outcome of important contributions that happened during the last four centuries. Although this method has a long history and a rich literature, there are still fundamental concepts that have not been thoroughly clarified, such as: What are the physically admissible thrust lines of an arch? How the stereotomy of an arch affects its limit stability? What is the role of the catenary curve (the alysoid)? Is the catenary curve a physically admissible thrust line? What is the relation between the direction of the thrust force and the slope of the thrust line at the point of application of the force? This thesis investigates these questions and aims to a better understanding of the limit equilibrium analysis of masonry arches, and at the same time, to present innovative methodologies and new analysis tools. Chapter 1 presents the work of other authors that have contributed the most to the stability analysis of masonry arches and vaulted structures over the last centuries. Chapter 2 revisits one of the most classical problems of Mechanics—what is the minimum thickness of a semicircular masonry arch subjected to its own weight. At the same time, the analysis presented in this chapter answers to the aforementioned questions through the development of closed-form expressions of the thrust line and the application of calculus of variation. Chapter 3 is focused on the limit equilibrium state of elliptical masonry arches, using the same approaches that were used in Chapter 2. This analysis was motivated from the fact that numerical results have been available in literature only for circular and not for elliptical masonry arches. Chapter 4 computes the location of the imminent hinges and the minimum thickness of circular masonry arches, for every given embrace angle, which can just sustain their own weight, together with a given level of horizontal ground acceleration. The numerical results presented in Chapters 2 to 4 confirm the accuracy of the in-house software that was developed for the needs of this thesis and the results obtained with a representative, commercially available software of the distinct element method. Chapter 5 present a comprehensive structural analysis of the Tunnel-Entrance to the Stadium of Ancient Nemea which ranges from the thrust line limit analysis and the discrete element method, to a 3-dimensional finite-element analysis. Limit equilibrium analysis emerges as a unique analysis method for the assessment of the stability of the structure and its interaction with the surrounding soil. While at the end of every chapter (Chapters 2 to 5) are presented detailed comments and conclusions, Chapter 6 is focused on outlining the most important conclusions and the main contribution of this thesis.