Με την τομογραφία επιχειρείται η δημιουργία μιας εικόνας-τομής των ιδιοτήτων ενός αντικειμένου χρησιμοποιώντας μετρήσεις της ενέργειας που διαδόθηκε μέσω του αντικειμένου, σε δέκτες στην περιφέρειά του. Έχει αποδειχθεί ότι ένα αντικείμενο είναι δυνατόν να απεικονισθεί με ακρίβεια από ένα πλήρες σύνολο περιμετρικών προβολών του. Η τομογραφική ανάλυση είναι δυνατό να γίνει αξιοποιώντας κυματική ενέργεια που έχει διαδοθεί μέσω του υπό μελέτη αντικειμένου (ακτινική τομογραφία) ή προέρχεται από ανάκλαση, περίθλαση κ.ο.κ. σε ανομοιογενή τμήματα του αντικειμένου (περιθλαστική τομογραφία). Στην παρούσα εργασία γίνεται πειραματική εφαρμογή της μεθόδου της σεισμικής τομογραφικής απεικόνισης σε αρχαιολογικά προβλήματα, όπως στην ανίχνευση αρχαιολογικών αντικειμένων στο υπέδαφος ή σε περιβάλλον τεχνητών λόφων (τύμβων -ταφικών μνημείων). Η εφαρμογή αυτή συνιστά ενδιαφέρουσα πρόκληση λόγω του ιδιαίτερα ανομοιογενούς περιβάλλοντος στο οποίο επιχειρείται η ανίχνευση των αντικειμένων καθώς και των μικρών γενικά, αλλά και συχνά απροσδιόριστων, διαστάσεων των αντικειμένων αυτών. Ο μη καταστρεπτικός χαρακτήρας της εφαρμογής είναι ιδιαίτερα ευνοϊκός και αποδεκτός στην αρχαιολογική έρευνα, τόσο σε αναγνωριστικό όσο και σε ανασκαφικό/ σωστικό επίπεδο. Παρουσιάζονται αναλυτικά οι αρχές και το θεωρητικό μέρος της τομογραφικής μεθοδολογίας, με ιδιαίτερη αναλυτική αναφορά στις δύο κύριες μεθόδους ανασύνθεσης των χαρακτηριστικών του εσωτερικού του ερευνώμενου μέσου (χώρου), την ακτινική τομογραφία και περίθλασης, καθώς και τις σχετικές τεχνικές και αλγόριθμους αντιστροφής. Τέλος γίνεται εφαρμογή της ακτινικής τομογραφίας σε 3 διαστάσεις, σε σεισμικές καταγραφές από αρχαιολογικούς-μνημειακούς χώρους, με τη χρήση του αλγόριθμου SIRT (τεχνική ταυτόχρονης επαναληπτικής ανασύνθεσης). Τα δεδομένα ήταν χρόνοι πρώτης άφιξης καθώς και πλάτη σε μια περίπτωση. Kατά την εφαρμογή της ακτινικής τομογραφίας, γίνονται δοκιμές ρύθμισης των παραμέτρων και μεταβλητών κατά τη διαδικασία της αντιστροφής, ώστε να διερευνηθεί η επίδρασή τους στο αποτέλεσμα. Επίσης, διερευνάται η αποτελεσματικότητα της μεθόδου με κατασκευή μοντέλων με βάση πληροφορίες για τους χώρους έρευνας, υπολογισμό συνθετικών χρόνων και αντιστροφή των συνθετικών χρόνων για σύγκριση με τα αποτελέσματα των πραγματικών χρόνων και εξαγωγή συμπερασμάτων. Λαμβάνοντας υπόψη τα στοιχεία από τις ανωτέρω δοκιμές, επιλέγονται τελικές τιμές παραμέτρων για την αντιστροφή των πραγματικών δεδομένων και γίνεται αξιολόγηση και ερμηνεία των αποτελεσμάτων. Γίνεται επίσης δοκιμαστική εφαρμογή τεχνικών επεξεργασίας / ενίσχυσης εικόνας στα αποτελέσματα της αντιστροφής, για την ενίσχυση κυρίως των μικρών, τοπικών, μεταβολών της σεισμικής ταχύτητας και την πληρέστερη απεικόνιση των χώρων έρευνας. Παρουσιάζονται τα αποτελέσματα από την έρευνα του προϊστορικού οικισμού Παλαμαρίου Σκύρου, του χώρου θεμελίωσης του Καθολικού (Ναού) της Ι.Μ. Δαφνίου Αττικής, του εσωτερικού του τεχνητού λόφου – πιθανού τύμβου Καστά Μεσολακκιάς Σερρών και των νεκροταφείων της περιοχής Chatby της Αλεξάνδρειας Αιγύπτου. / With tomography is attempted the creation of an image – slice of the internal properties of a medium / object with use of measurements of the energy which was transmitted through the medium, at receivers on the outside of it. It has been proved that a medium can be imaged with precision by a set of perimetric projections. Tomographic analysis is possible by employing wave energy that has been transmitted through the medium under investigation (ray tomography) or has been reflected or diffracted at various inhomegeneities within the medium (diffraction tomography). In the present study, an experimental application of the seismic tomographic imaging method is made in the field of archaeological exploration, such as detection of archaeological remains in the subsurface or in artificial hills (tumuli). The application poses an interesting challenge, because of the highly inhomogeneous environment where the detection is tried and the small and unknown, in general, dimensions of the target objects. The non-destructive character of the application makes it quite attractive and acceptable in archaeological exploration, both at reconnaissance and later stages of excavation. The principles and theoretical foundations of the tomographic methodology are presented, with particular reference to the two main methods of reconstruction of the internal character of the investigated medium and included objects, that is ray and diffraction tomography, as well as the relevant inversion techniques and algorithms. Finally, an application of ray tomography on seismic recordings from archaeological / monumental sites is made, with use of the simultaneous iterative reconstruction technique (SIRT). In the application of ray tomography, adjustment tests of inversion parameters and variables are made, in order to study their effect on the inversion result. The effectiveness of the method and of the algorithm is investigated with the aid of synthetic models based on possible information from the test sites and the bibliography. Using the models as starting models, synthetic travel times are calculated and are then inverted, while the resulting image is compared with the synthetic starting model and conclusions are drawn. Based on the results of the above tests, final parameter values are set for the inversion of real field data and an assessment and an interpretation of the results is made. Image processing and enhancement is applied to the resulting tomograms, in order to enhance and clarify the small (hidden) variations of seismic velocity, which may reflect important subsurface information, mainly related to particular archaeological targets. There are presented the results from the investigation of the subsurface of a part of the prehistoric settlement of Palamari (Skyros isl., Greece) of the foundations of the church of Dafni Monastery (Attica, Greece), of the interior of the artificial hill / possible tumulus at Kasta Messolakia (Serrae, Greece) and of the subsurface of the cemeteries of Chatby area (Alexandria, Egypt).

Στην εργασία αυτή διερευνάται η δυναμική απόκριση κρυογενικών δεξαμενών υγροποιημένου φυσικού αερίου υποκείμενων σε δράση σεισμού. Επειδή οι εν λόγω δεξαμενές αποτελούν κρίσιμα στοιχεία για τη λειτουργία ενός συστήματος διανομής φυσικού αερίου και επειδή ενδεχόμενη καταστροφική αστοχία τους θα είχε δυσβάστακτες κοινωνικές επιπτώσεις, ιδιαίτερα αυστηρές απαιτήσεις ισχύουν σχετικά με τον αντισεισμικό σχεδιασμό τους. Για το σχεδιασμό έναντι σεισμού προδιαγράφονται συνήθως ένα σεισμικό γεγονός με μέση περίοδο επαναφοράς της τάξης των 500 ετών, κατά το οποίο οι δεξαμενές απαιτείται να παραμένουν πλήρως λειτουργικές, και ένα σεισμικό γεγονός με μέση περίοδο επαναφοράς της τάξης των 5000 έως 10000 ετών, κατά το οποίο απαιτείται η εξασφάλιση της ασφαλούς διακοπής της λειτουργίας των δεξαμενών, ενώ ελάχιστη μόνο βλάβη είναι αποδεκτή στα δομικά στοιχεία τους. Οι εν λόγω δεξαμενές αποτελούνται από δύο κελύφη: Ένα εξωτερικό κυλινδρικό κέλυφος από προεντεταμένο σκυρόδεμα με θολωτή στέγη από οπλισμένο σκυρόδεμα και ένα εσωτερικό κέλυφος από κρυογενικό χάλυβα, εντός του οποίου αποθηκεύεται το προϊόν. Τα δύο κελύφη εδράζονται σε κοινή πλάκα βάσης. Στην εργασία αυτή εξετάζονται δύο εναλλακτικοί τρόποι έδρασης της πλάκας βάσης: α) απλή έδραση επί του διαμορφωμένου εδάφους και β) εφαρμογή κάποιας διάταξης σεισμικής μόνωσης επί της οποίας εδράζεται η πλάκα βάσης. Ειδικότερα, εξετάζονται τα ακόλουθα συστήματα σεισμικής μόνωσης: α) μόνωση με ελαστομεταλλικά εφέδρανα με ελαστομερές υψηλής απόσβεσης, β) μόνωση με ελαστομεταλλικά εφέδρανα με πυρήνα μόλυβδου, και γ) μόνωση με συνδυασμό γραμμικών ελαστομεταλλικών εφεδράνων και μη γραμμικών αποσβεστήρων. Στις ιδιαιτερότητες του δυναμικού προβλήματος περιλαμβάνονται, εκτός από την επίδραση της μόνωσης, η δυναμική αλληλεπίδραση των δύο κελυφών, η δυναμική αλληλεπίδραση υγρού-δεξαμενής, ο κυματισμός της ελεύθερης επιφάνειας του υγρού και η δυναμική αλληλεπίδραση εδάφους-δεξαμενής. Τα φαινόμενα αυτά λαμβάνονται υπόψη στην ανάλυση. Για τη μελέτη του κυματισμού της ελεύθερης επιφάνειας χρησιμοποιείται η γραμμική θεωρία κυμάτων. Για τη μελέτη της αλληλεπίδρασης υγρού-δεξαμενής εφαρμόζεται μια προσέγγιση τύπου Lagrange με χρήση πεπερασμένων στοιχείων υγρού. Για τη μελέτη της αλληλεπίδρασης εδάφους-δεξαμενής χρησιμοποιούνται συγκεντρωμένα στοιχεία δυσπαραμορφωσιμότητας και απόσβεσης. Επίσης, συγκεντρωμένα στοιχεία δυσπαραμορφωσιμότητας και απόσβεσης χρησιμοποιούνται για την προσομοίωση των συστημάτων μόνωσης. Τα χαρακτηριστικά των παραπάνω συγκεντρωμένων στοιχείων για την προσομοίωση της αλληλεπίδρασης εδάφους-δεξαμενής και των συστημάτων μόνωσης επιλέγονται με βάση ρεαλιστικές υποθέσεις για τις δυναμικές ιδιότητες του εδάφους θεμελίωσης και τα χαρακτηριστικά των συστημάτων μόνωσης. Δύο περιπτώσεις υφιστάμενων δεξαμενών υιοθετούνται προς ανάλυση. Οι δύο δεξαμενές διαφέρουν κυρίως ως προς τη χωρητικότητα και το λόγο ύψους προς ακτίνα. Και οι δύο δεξαμενές έχουν μελετηθεί και κατασκευαστεί χωρίς σεισμική μόνωση. Ένα συνθετικό επιταχυνσιογράφημα συμβατό με το φάσμα απόκρισης σχεδιασμού του Ευροκώδικα 8 για λόγο απόσβεσης 5%, για την αντίστοιχη κατηγορία εδάφους και για μέγιστη εδαφική επιτάχυνση λίγο μεγαλύτερη από αυτή για την οποία είναι σχεδιασμένες οι δεξαμενές περιγράφει τη διέγερση της κατασκευής στην οριζόντια διεύθυνση. Για την ανάλυση αναπτύσσεται ένα συνολικό παραμετρικό προσομοίωμα πεπερασμένων στοιχείων το οποίο περιλαμβάνει τα δύο κελύφη και το αποθηκευμένο υγρό. Στο προσομοίωμα συμπεριλαμβάνονται, ανάλογα με την περίπτωση που κάθε φορά αναλύεται, τα κατάλληλα στοιχεία για την προσομοίωση της αλληλεπίδρασης εδάφους-δεξαμενής ή των διαφόρων συστημάτων μόνωσης. Το προσομοίωμα χρησιμοποιείται για τον υπολογισμό των ιδιομορφών και την πραγματοποίηση δυναμικών γραμμικών και μη γραμμικών αναλύσεων με ολοκλήρωση στο πεδίο του χρόνου. Οι μη γραμμικές αναλύσεις είναι αναγκαίες για να ληφθεί υπόψη με ρεαλιστικό τρόπο η έντονα μη γραμμική συμπεριφορά των συστημάτων μόνωσης. Στην περίπτωση αυτή η μη γραμμικότητα είναι συγκεντρωμένη στα στοιχεία που προσομοιώνουν τη μόνωση ενώ για την υπόλοιπη κατασκευή γίνεται σε κάθε περίπτωση η παραδοχή γραμμικής συμπεριφοράς. Η παραδοχή αυτή προκύπτει από το γεγονός ότι για τις δεξαμενές αυτού του τύπου η αποδεκτή βλάβη, και επομένως μη γραμμικότητα στη συμπεριφορά, των δομικών στοιχείων είναι ελάχιστη. Παρουσιάζονται επιλεγμένα αποτελέσματα των αναλύσεων που πραγματοποιήθηκαν με χρήση του προσομοιώματος πεπερασμένων στοιχείων. Τα αποτελέσματα αυτά αναφέρονται κυρίως στην τέμνουσα βάσης, τη ροπή ανατροπής, τις οριζόντιες μετακινήσεις της κατασκευής σε διάφορες στάθμες και το ύψος του κυματισμού. Έμφαση δίνεται στην ποσοτικοποίηση της επίδρασης της αλληλεπίδρασης εδάφους-δεξαμενής και της σεισμική μόνωσης στα ανωτέρω μεγέθη. Η εργασία τελειώνει με γενικά συμπεράσματα. / In the present work the dynamic response of liquefied natural gas cryogenic tanks subjected to earthquake action is investigated. Since these tanks are critical elements for the function of a natural gas distribution system and because a potential failure of them could lead to a major disaster, very severe requirements concerning the seismic design of these tanks are imposed. For the design against earthquake action two seismic events are generally considered: an event with a mean return period of 475 years during which the tanks are expected to remain fully operational and an event with a mean return period of the order of 5000 to 1000 years during which the safe shut down of the tanks is to be ensured, while minimum damage is accepted at the structural parts of the tanks. The type of tanks under consideration is constituted of two shells: an outer one made of prestressed concrete and an inner one made of cryogenic steel in which the product is stored. Both shells rest on a common base slab. In these work two alternative ways for the support of this slab are examined: a) the slab lays directly on the ground and b) the slab lays on a number of devices which provide seismic isolation. The following isolation systems are examined. a) high damping rubber bearings, b) lead core rubber bearings and c) low damping rubber bearing in conjunction with non-linear viscous dampers. The examined structural dynamic problem is characterised by certain particularities, besides the one of the implementation of an isolation system. The most important of them are sloshing of the free surface, dynamic fluid-structure interaction, dynamic soil-structure interaction and interaction between the inner and outer shell. These phenomena are taken into consideration in the preformed analyses. For the modeling of the free surface effect the realistic assumption of small wave height is made and the linear wave theory is employed. For the modeling of the fluid-structure interaction a Lagrangian approach is applied using finite elements for the modeling of the liquid and the solid domain. The soil-structure interaction is simulated by using concentrated stiffness and damping elements. Concentrated stiffness and damping elements are used for the modeling of the isolation systems as well. Two cases of existing tanks are adopted for analysis. The two tanks differ mainly in capacity and in the height to radius ratio. Both tanks are designed and constructed without seismic isolation. A global parametric finite element model is developed for the analyses. In this model the appropriate elements are incorporated for the modeling of the soil-structure interaction and the seismic isolation system, depending on the analysed case. The model is used for the calculation of eigenmodes and eigenfrequencies and for performing linear and non-linear transient analyses in time domain. Non-linear analyses are necessary in order for the highly non-linear behavior of the isolation devices to be properly simulated. In these cases the non-linearity is concentrated at the elements which simulate the seismic isolation system while the rest of the structure is considered elastic in all cases. This consideration is justified by the fact that for these tanks the acceptable damage, and consequently non-linearity in the behavior of the structural part, is minimum. Selective results of the performed analyses are presented. These results refer mainly to base shear forces, overturning moments, horizontal displacements at different levels of the tank and sloshing heights. Emphasis is on the quantification of the impact of the soil-structure interaction and the implementation of the examined seismic isolation systems.

Η παρούσα διατριβή ανήκει στο γενικότερο θεματικό πεδίο της σεισμικής αποτίμησης, σχεδιασμού ή ανασχεδιασμού κατασκευών οπλισμένου σκυροδέματος με βάση τις μετακινήσεις. Οι σύγχρονες μέθοδοι αυτού του τύπου, στηρίζονται σε έλεγχο και σύγκριση της σεισμικής απαίτησης με την ικανότητα των μελών της κατασκευής σε όρους μετακινήσεων παρά σε όρους δυνάμεων. Δημιουργείται επομένως η ανάγκη για απλό και αξιόπιστο υπολογισμό της συμπεριφοράς μελών οπλισμένου σκυροδέματος σε κάμψη και διάτμηση, σε όρους μετακινήσεων. Το αντικείμενο της παρούσης διατριβής είναι η ανάπτυξη προσομοιωμάτων για τον υπολογισμό των βασικών χαρακτηριστικών της συμπεριφοράς καμπτόμενων μελών οπλισμένου σκυροδέματος και συγκεκριμένα: της ροπής διαρροής, της παραμόρφωσης στη διαρροή, της ενεργού δυσκαμψίας, της παραμόρφωσης στην αστοχία, της διατμητικής αντοχής σε ανακυκλιζόμενη φόρτιση, της αντοχής μελών με χαμηλό λόγο διάτμησης και της συμπεριφοράς υπό διαξονική καταπόνηση. Εξετάζονται μέλη διαφόρων τύπων και διαφορετικής διατομής, μέλη με ενίσχυση μανδύα οπλισμένου σκυροδέματος ή μανδύα σύνθετων υλικών, καθώς επίσης και μέλη με μάτιση του διαμήκους οπλισμού στην περιοχή πλαστικής άρθρωσης. Για την ανάπτυξη των προσομοιωμάτων, καθώς και για τον έλεγχο άλλων παλαιότερων, αναπτύχθηκε και αξιοποιήθηκε βάση πειραματικών δεδομένων μελών οπλισμένου σκυροδέματος με περισσότερα από 2800 πειράματα από τη διεθνή βιβλιογραφία. Για τον υπολογισμό της ροπής και της καμπυλότητας στη διαρροή, αναπτύσσονται απλές σχέσεις υπολογισμού, βασιζόμενες σε ανάλυση σε επίπεδο διατομής και καθορίζονται τα κατάλληλα κριτήρια διαρροής. Αναπτύσσονται ακολούθως σχέσεις υπολογισμού της παραμόρφωσης στη διαρροή, και συγκεκριμένα της γωνίας στροφής χορδής του μέλους στη διαρροή, θy, ως άθροισμα τριών όρων: καμπτικής παραμόρφωσης, διατμητικής παραμόρφωσης και παραμόρφωσης λόγω ολίσθησης των ράβδων διαμήκους οπλισμού από την περιοχή αγκύρωσης. Προτείνονται δε δύο εναλλακτικοί τρόποι υπολογισμού της ενεργού δυσκαμψίας, ένας θεωρητικός και ένας καθαρά εμπειρικός. Στη συνέχεια εξετάζεται η παραμόρφωση στην αστοχία και προτείνονται δύο εναλλακτικοί μέθοδοι υπολογισμού της γωνίας στροφής χορδής στην αστοχία, θu. Η 1η βασίζεται στον υπολογισμό της καμπυλότητας στην αστοχία, φu, με εφαρμογή του κατάλληλου προσομοιώματος περίσφιγξης του σκυροδέματος, και στην εφαρμογή της φu σε μήκος πλαστικής άρθρωσης ίσο με Lpl, ενώ η 2η σε καθαρά εμπειρικές εξισώσεις. Εξετάζεται ακολούθως η διατμητική αντοχή σε ανακυκλιζόμενη φόρτιση και προτείνονται προσομοιώματα για αστοχία σε διαγώνιο εφελκυσμό ή αστοχία σε λοξή θλίψη, μετά την καμπτική διαρροή. Στη συνέχεια εξετάζεται η συμπεριφορά μελών οπλισμένου σκυροδέματος υπό διαξονική καταπόνηση. Εξετάζονται επίσης μέλη με χαμηλό λόγο διάτμησης και προτείνονται νέα αντιπροσωπευτικότερα κριτήρια για τον χαρακτηρισμό ενός μέλους ως “κοντό μέλος”, καθώς και νέα μεθοδολογία υπολογισμού της αντοχής των μελών αυτών, με κατάλληλο συνδυασμό του προσομοιώματος των Shohara and Kato, 1981 και των Φαρδής και συνεργάτες 1998. Ακολούθως εξετάζονται μέλη ενισχυμένα με μανδύα σύνθετων υλικών και προτείνονται προσομοιώματα υπολογισμού της γωνίας στροφής χορδής στη διαρροή και την καμπτική αστοχία, καθώς και προσομοίωμα υπολογισμού της διατμητικής αντοχής. Στη συνέχεια εξετάζεται η συμπεριφορά μελών με μάτιση του διαμήκους οπλισμού στην περιοχή πλαστικής άρθρωσης, καθώς και η εφαρμογή μανδύα σύνθετων υλικών για την ενίσχυση της περιοχής αυτής. Τέλος εξετάζεται η συμπεριφορά στη διαρροή και στην αστοχία, μελών ενισχυμένων με μανδύα οπλισμένου σκυροδέματος. Η ανάπτυξη όλων των προτεινόμενων προσομοιωμάτων της διατριβής βασίζεται στην καλύτερη δυνατή συμφωνία με τα πειραματικά αποτελέσματα της βάσης δεδομένων, χωρίς όμως να θυσιάζεται η απλότητα και η ευχρηστία αυτών. / The present Thesis belongs in the general field of seismic assessment, design and redesign of concrete structures with displacement based procedures. Modern methods of this kind are based in controlling and comparing seismic demand with structural elements capacity in terms of displacements rather than forces. This leads in the need of estimating reinforced concrete elements performance under bending and shear, in terms of displacements. The object of the Thesis is development of models for calculating the basic performance characteristics of reinforced concrete elements under bending, in particular: yield moment, deformation at yielding, effective stiffness, deformation at ultimate, shear strength under cyclic loading, maximum strength of members with low shear ratio and behavior under biaxial loading. Members with various types of section and various characteristics are included, as also members retrofitted with FRP jacket or concrete jacket and members with lap-splice of longitudinal reinforcement in plastic hinge region. In order to develop new models and check older ones, a database of more than 2800 experiments from international literature on reinforced concrete elements was created and used here. Simple equations and procedures are suggested for calculating yield moment and corresponding curvature, based on section analysis, by specifying the appropriate yield criteria. Equations for calculating deformation at yielding, in particular chord rotation at yielding, θy as the sum of deformations due to bending, due to shear and due to slippage of longitudinal reinforcement from anchorage zone, are also developed. Calculation of effective stiffness is based on two alternative models, one theoretical and one purely empirical. Deformation at ultimate is then examined where two methods for calculating chord rotation at ultimate are suggested. 1st one is based on ultimate curvature, φu, where an appropriate concrete confinement model is used, and plastic hinge length Lpl, while 2nd one is based on purely empirical equations. Shear strength under cyclic loading is also examined and new models for calculating shear strength for shear tension and shear compression failure after flexural yield are developed. Behavior of reinforced concrete elements under biaxial loading is then examined. Elements with low shear ratio are also covered and new, more representative, criteria to characterize an element as a “short element” are suggested. A procedure based on an appropriate combination of Shohara and Kato 1981 model and Fardis et al. 1998 model is then suggested for calculating maximum strength of such “short elements”. Retrofitted members with FRP jacket are then examined and models for chord rotation at yielding and ultimate, as well as for shear strength are suggested. Behavior of members with lap-splice of longitudinal reinforcement inside plastic hinge region is then examined, including also retrofitting of this region with FRP jacket. Performance at yielding and ultimate of retrofitted members with concrete jacket is also examined. Development of all the suggested models of the Thesis is based on best fit with experimental results of the database, without sacrificing simplicity and applicability of the models.

Η παρούσα διατριβή επανεξετάζει την οριακή ανάλυση ευστάθειας των λίθινων αψίδων. Η οριακή ανάλυση ευστάθειας χρησιμοποιείται σήμερα ως το βασικό εργαλείο αποτίμησης της ευστάθειας τόξων και θολωτών κατασκευών από τοιχοποιία, όπως ακριβώς συνέβαινε και τους τελευταίους τέσσερις αιώνες. Παρά την τόσο μακρόχρονη ιστορία της μεθόδου, δεν έχουν πλήρως διασαφηνιστεί στην επιστημονική κοινότητα θεμελιώδης έννοιες και δεν έχουν σαφώς απαντηθεί ερωτήματα όπως: Ποιες είναι οι φυσικά πραγματοποιήσιμες γραμμές ώθησης και ποιες όχι; Ποια είναι η επίδραση της στερεοτομίας ενός τόξου στην οριακή του ευστάθεια; Ποιος είναι ο ρόλος της αλυσοειδούς καμπύλης και κατά πόσο αυτή είναι μία φυσικά αποδεκτή γραμμή ώθησης; Τι σχέση υπάρχει ανάμεσα στην κλίση της συνισταμένης θλιπτικής δύναμης και στην κλίση της γραμμής ώθησης στο σημείο εφαρμογής της; Η παρούσα διατριβή αναζητά απαντήσεις στα ερωτήματα αυτά, και έχει ως στόχο τη βαθύτερη κατανόηση της οριακής ανάλυσης ευστάθειας των τόξων, με παράλληλη ανάδειξη νέων υπολογιστικών διαδικασιών. Η δομή της παρουσιάζεται συνοπτικά παρακάτω. Στο πρώτο κεφάλαιο γίνεται ιστορική ανάλυση της μεθόδου μέσα από παρουσίαση και σχολιασμό των εργασιών με τη σημαντικότερη συμβολή, από τα μέσα του 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.