Η παρούσα διατριβή ασχολείται με την μελέτη της γραμμικής, ανιξωδικής ευστάθειας, συμπιεστού, υπερηχητικού οριακού στρώματος γύρω από κύλινδρο ή κώνο. Το υλικό που παρουσιάζεται μπορεί να χωριστεί σε δύο μέρη. i) Το αριθμητικό αέρος, όπου γίνεται επέκταση, βελτίωση και σύγκριση αριθμητικών τεχνικών που χρησιμοποιούνται για την μελέτη της γραμμικής, ανιξωδικής ευστάθειας, συμπιεστών, υπερηχητικών οριακών στρωμάτων. Σε αυτό το μέρος αναπτύσσονται δύο τεχνικές για την επίλυση των εξισώσεων ανιξωδικής ευστάθειας οι οποίες βασίζονται στους πίνακες παραγώγισης Chebyshev για την διακριτοποίηση: α) Επαναληπτική μέθοδος: Σύμφωνα με αυτή τη μέθοδο, οι εξισώσεις ευστάθειας διάκριτοποιούνται σε όλο το διάστημα ολοκλήρωσης με πίνακες παραγώγισης Chebyshev. Το διακριτοποιημένο, γραμμικό, αλγεβρικό σύστημα έχει τη μορφή AΧ - β, όπου ο πίνακας A περιέχει τους συντελεστές από τη διακριτοποίηση, το διάνυσμα Χ περιέχει τις (ζητούμενες) τιμές των ιδιοσυναρτήσεων στα σημεία του πλέγματος και το διάνυσμα b περιέχει τις συνοριακές συνθήκες. Δίνοντας μια αρχική εκτίμηση της ιδοτιμής το σύστημα επιλύεται επαναληπτικά μέχρι να ικανοποιηθεί η συνθήκη μη διαπερατότητας του στερεού ορίου. β) Ολική μέθοδος (QΖ αλγόριθμος): Σε αυτή τη μέθοδο, οι εξισώσεις ευστάθειας διάκριτοποιούνται σε όλο το διάστημα ολοκλήρωσης με ανάπτυγμα σε πολυώνυμα Chebyshe, ενώ οι ασυμπτωτικές τιμές των διαταραχών εκτός του οριακού στρώματος αντικαθίστανται από μηδενικές συνθήκες δίνοντας έτσι ένα γενικευμένο, αλγεβρικό πρόβλημα ιδιστιμώντης μορφής ΑΧ = cΒΧ. Το πρόβλημα επιλύεται από τον αλγόριθμο QΖ και δίνει μια προσέγγιση όλου του χάσματος ίδιο τιμών χωρίς την ανάγκη αρχικές εκτίμησης. Βασικό χαρακτηριστικό των εξισώσεων ανιξωδικής ευστάθειας και πηγή δυσκολιών στην αριθμητική τους επίλυση, είναι η ύπαρξη ιδιαζόντων (Singular) σημείων τα οποία κάνουν δύσκολο και σε κάποιες περιπτώσεις αδύνατο τον εντοπισμό και ακριβή υπολογισμό των ιδιοτιμών. Για την υπέρβαση αυτής της δυσκολίας, οι εξισώσεις ευστάθειας ολοκληρώνονται σε μια διαδρομή εντός του μιγαδικού επιπέδου, η οποία διέρχεται αρκετά μακριά από το ανώμαλο σημείο αυξάνοντας έτσι την ακρίβεια και την αποτελεσματικότητα των αριθμητικών σχημάτων. Το πρόβλημα που προκύπτει σε μια τέτοια περίπτωση είναι ο υπολογισμός των κατανομών της ταχύτητας και της θερμοκρασίας της βασικές ροής στο πλέγμα μιγαδικών σημείων όπου θα επιλυθούν οι εξισώσεις ευστάθειας. H συνήθης πρακτική για την αντιμετώπιση αυτού του προβλήματος είναι η επίλυση των εξισώσεων της βασικές ροής σε πλέγμα πραγματικών σημείων και στη συνέχεια το ανάπτυγμα κατά Taylor των λύσεων αυτών ώστε να υπολογιστούν οι τιμές σε ένα πλέγμα μιγαδικών σημείων. Στην εργασία αυτή 1) εξετάζεται η επίδραση του σφάλματος αποκοπής του αναπτύγματος Taylor στην ακρίβεια των υπολογισμών και 2) προτείνεται μια τεχνική όπου οι εξισώσεις του οριακού στρώματος επιλύονται απ’ ευθείας στο μιγαδικό επίπεδο. Για τις εξισώσεις ευστάθειας που χρησιμοποιούμε, υπάρχει αναλυτική έκφραση που δίνει την ασυμπτωτική μορφή των διαταραχών μακριά από το στερεό όριο. Στην έκφραση αυτή των διαταραχών εμπλέκονται μη-γραμμικά οι ιδιοτιμές. Υιοθετώντας λοιπόν αυτού του είδους τις συνθήκες θα πρέπει αναγκαστικά να χρησιμοποιήσουμε μια επαναληπτική μέθοδο όπως η (α) που προτείνουμε εδώ και να αποκλείσουμε τις ολικές τεχνικές όπως η (β). Στα πλαίσια της παρούσας εργασίας χρησιμοποιούμε άλλα δύο είδη συνοριακών συνθηκών, ομογενείς και γραμμικές μη-ομογενείς. Συγκρίνοντας τα αριθμητικά αποτελέσματα και από τα τρία είδη συνθηκών διερευνάται η επίδραση που έχει η επιλογή τους στην ακρίβεια και την αποτελεσματικότητα των υπολογισμών. Τέλος, τα αποτελέσματα όλων των παραπάνω τεχνικών συγκρίνονται με γνωστά από τη διεθνή βιβλιογραφία αποτελέσματα, τόσο για τη βασική ροή όσο και για τη χρονική Ευστάθεια, ώστε να διαπιστωθεί η ακρίβεια των υπολογισμών. ii) Στο δεύτερο αέρος επεκτείνουμε την μελέτη της φυσικής του προβλήματος με τη χρήση των αριθμητικών τεχνικών που αναπτύχθηκαν. H μελέτη επεκτάθηκε στην χωρική ευστάθεια, η οποία περιγράφει καλύτερα την φυσική του προβλήματος, και πιο συγκεκριμένα μελετήθηκε: α) Η επίδραση της επιλογής του νομού του ιξώδους και της τιμής του αριθμού Prandtl στην ακρίβεια των υπολογισμών. β) Η αποτελεσματικότητα διαφόρων τεχνικών για τον έλεγχο της ευστάθειας. Οι τεχνικές που μελετήσαμε είναι β1) η ομοιόμορφη θέρμανση/ψύξη των τοιχωμάτων, β2) η ομοιόμορφη έγχυση/αναρρόφηση και β3) κατανεμημένη έγχυση/αναρρόφηση ρευστού μέσω των τοιχωμάτων. / -

Σκοπός της παρούσας Διατριβής Ειδίκευσης είναι η παρακολούθηση του κατολισθητικού φαινομένου του οικισμού Καρυάς του Δήμου Πατρέων, με μετρήσεις αποκλισιομέτρων ώστε να προσδιοριστεί ποσοτικά ο ρυθμός μετακίνησης της περιοχής μελέτης και να εξαχθούν συμπεράσματα για την μελλοντική του εξέλιξη, όπως καθορίζεται μέσω αναλύσεων ευστάθειας πρανών. Η εν λόγω περιοχή παρουσιάζει αξιόλογο ενδιαφέρον, λόγω των ιδιαίτερων και σύνθετων γεωλογικών συνθηκών, των συχνών και σημαντικών κατολισθητικών φαινομένων των τελευταίων ετών που προκάλεσαν σημαντικές καταστροφές σε τεχνικά έργα και ανθρώπινες περιουσίες. Στα πλαίσια της Διατριβής και της παράλληλης Έρευνας του Εργαστηρίου Τεχνικής Γεωλογίας στην περιοχή ενδιαφέροντος πραγματοποιήθηκε λεπτομερής γεωτεχνική έρευνα που περιελάμβανε ανόρυξη δειγματοληπτικών γεωτρήσεων σε επιλεγμένες θέσεις με ταυτόχρονη τοποθέτηση αποκλισιομέτρων, συστηματική λήψη αποκλισιομετρικών μετρήσεων, σειρά εργαστηριακών και επί τόπου δοκιμών σε κατάλληλα διαμορφωμένα δείγματα, λεπτομερή καταγραφή των τεχνικογεωλογικών συνθηκών της περιοχής με βάση κυρίως τη σύσταση, τη δομή των γεωλογικών σχηματισμών. Έγινε διάκριση των σχηματισμών της περιοχής έρευνας σε επιμέρους Γεωτεχνικές Ενότητες και σχεδιάστηκαν αντιπροσωπευτικές γεωτεχνικές – τεχνικογεωλογικές διατομές με σκοπό τη σύνθεση και ερμηνεία του γεωτεχνικού περιβάλλοντος της περιοχής. Παράλληλα πραγματοποιήθηκαν αναλύσεις ευστάθειας πρανών με χρήση εξειδικευμένου λογισμικού σε επιλεγμένες διατομές κατά μήκος των γεωτρήσεων στην κατολισθαίνουσα ζώνη Καρυάς. Σκοπός των αναλύσεων ευστάθειας πρανών ήταν να γίνει εκτίμηση πιθανών μετακινήσεων σε επιλεγμένα σημεία της περιοχής έρευνας και να παρατηρηθεί εάν το είδος ολίσθησης εξακολουθεί να είναι το ίδιο με αυτό που είχε διαπιστωθεί κατά την τελευταία εκδήλωση του κατολισθητικού φαινομένου. / The Aim of this ΜSc. Project is the exploration of the landslide phenomeno in the region of Karya, Patras, Western Greece by inclinometer measurements in order to define quantitatively the movement rate of the study area and export conclusions on its future development, as it specified via slope stability analysis. The study area is an area with complex geological conditions and is well known especially for the important and large-scale landslide phenomena which have taken place during the last years near Karya village and caused important destructions in engineering works and human fortunes. This project was done along with survey of the Laboratory of Engineering Geology, Geology Department, University of Patras. The full survey program included geological mapping focused on the engineering properties of the soil and rock formations, fifteen (15) sampling boreholes, inclinometer installation and measurement and laboratory tests. Through laboratory testing and detailed description of soil and rock samples the formations were classified in several units and several geotechnical cross-sections were constructed to make a general geological-geotechnical model of the study area. The geological and geotechnical data were used to perform Limit Equilibrium Slope Stability analysis. The survey results were used to define probable movements of the study area and to establish the slipping type.

Σκοπός της παρούσας Διατριβής Ειδίκευσης είναι η ανάλυση, η παρουσίαση και η αξιολόγηση των τεχνικογεωλογικών – γεωτεχνικών συνθηκών των οικισμών Καρυάς και Ν. Σουλίου του Δήμου Πατρέων, σε σχέση με την ευστάθεια των φυσικών πρανών και τις συνθήκες γεωλογικής καταλληλότητας για τις προς δόμηση περιοχές. Η εν λόγω περιοχή παρουσιάζει αξιόλογο ενδιαφέρον, λόγω των ιδιαίτερων και σύνθετων γεωλογικών συνθηκών, των συχνών και σημαντικών κατολισθητικών φαινομένων των τελευταίων ετών που προκάλεσαν σημαντικές καταστροφές σε τεχνικά έργα και ανθρώπινες περιουσίες, αλλά και λόγω της ραγδαίας οικιστικής ανάπτυξης που παρατηρείται τα τελευταία χρόνια στα περίχωρα της πόλης της Πάτρας. Στα πλαίσια της Διατριβής και της παράλληλης Έρευνας του Εργαστηρίου Τεχνικής Γεωλογίας στην περιοχή ενδιαφέροντος πραγματοποιήθηκε κατά τα τη χρονική περίοδο 2003-2005 λεπτομερής γεωτεχνική έρευνα που περιελάμβανε ανόρυξη δειγματοληπτικών γεωτρήσεων σε επιλεγμένες θέσεις με ταυτόχρονη τοποθέτηση αποκλισιομέτρων, συστηματική λήψη αποκλισιομετρικών μετρήσεων, σειρά εργαστηριακών και επί τόπου δοκιμών σε κατάλληλα διαμορφωμένα δείγματα, λεπτομερή καταγραφή των τεχνικογεωλογικών συνθηκών της περιοχής με βάση κυρίως τη σύσταση, τη δομή των γεωλογικών σχηματισμών και τη γενικότερη φυσική τους κατάσταση και συμπεριφορά, αλλά και σε φαινόμενα ευστάθειας των φυσικών πρανών, ιδιαιτέρως λόγω του σημαντικού ιστορικού κατολισθητικών κινήσεων στην περιοχή της Καρυάς τα τελευταία χρόνια. Μέσω της κατασκευής αντιπροσωπευτικών εδαφικών προσομοιωμάτων έγινε διάκριση των σχηματισμών της περιοχής έρευνας σε επιμέρους Γεωτεχνικές Ενότητες και σχεδιάστηκαν αντιπροσωπευτικές γεωτεχνικές – τεχνικογεωλογικές διατομές με σκοπό τη σύνθεση και ερμηνεία του γεωτεχνικού περιβάλλοντος της περιοχής. Παράλληλα πραγματοποιήθηκαν αναλύσεις ευστάθειας πρανών με χρήση εξειδικευμένου λογισμικού σε επιλεγμένες διατομές κατά μήκος των γεωτρήσεων, τόσο στην κατολισθαίνουσα ζώνη της Καρυάς, όσο και στην οικιστική ζώνη. Η σύνθεση των αποτελεσμάτων της γεωλογικής και γεωτεχνικής έρευνας οδήγησε στον προκαταρκτικό διαχωρισμό της περιοχής σε ζώνες κατάλληλες, ακατάλληλες ή κατάλληλες για δόμηση υπό προϋποθέσεις με σκοπό τη διασφάλιση του δομημένου περιβάλλοντος από φυσικούς κινδύνους ή κινδύνους προερχόμενους από ανθρώπινες επεμβάσεις και δραστηριότητες. / The aim of this ΜSc. Project is the presentation, the analysis and assessment of the geological and geotechnical conditions of the sub-urban regions of N. Souli – Karya, Patras, Western Greece. The study area is an area with complex geological conditions and is well known especially for the important and large-scale landslide phenomena which have taken place during the last years near Karya village. This project was done along with survey of the Laboratory of Engineering Geology, Geology Department, University of Patras. The full survey program included geological mapping focused on the engineering properties of the soil and rock formations, fifteen (15) sampling boreholes, inclinometer installation and measurement and laboratory tests. Through laboratory testing and detailed description of soil and rock samples the formations were classified in several units and several geotechnical cross-sections were constructed to make a general geological-geotechnical model of the study area. The geological and geotechnical data were used to perform Limit Equilibrium Slope Stability analysis. The survey results were used to divide the regions of Karya and N. Souli in some zones according to soil suitability for construction purposes.

Αντικείμενο της παρούσας εργασίας είναι η πρόβλεψη τόσο της κίνησης, αλληλεπίδρασης και παραμόρφωσης δύο φυσαλίδων λόγω μεταβολής της πίεσης στο περιβάλλον ιξώδες υγρό, όσο και της ανοδικής κίνησης μιας φυσαλίδας λόγω άνωσης σε ένα Νευτωνικό ή ιξωδοπλαστικό ρευστό. Για τη μοντελοποίηση των αλληλεπιδρώντων φυσαλίδων, αναπτύχθηκε μιας νέα ελλειπτική μεθόδος κατασκευής του υπολογιστικού πλέγματος προκειμένου να αντιμετωπιστούν επιτυχώς τα ιδιάζοντα σημεία (πόλοι) των φυσαλίδων και οι μεγάλες παραμορφώσεις των διεπιφανειών τους. Με τη μέθοδο αυτή η πύκνωση του πλέγματος περιορίζεται μόνο στις περιοχές που είναι αναγκαίο, μειώνοντας έτσι το υπολογιστικό κόστος και αυξάνοντας την ακρίβεια των υπολογισμών. Για την επίλυση των παρακάτω προβλημάτων χρησιμοποιήθηκε η μέθοδος των μικτών πεπερασμένων στοιχείων κατά Galerkin. Στην περίπτωση των αλληλεπιδρώντων φυσαλίδων έχει εξετασθεί η επίδραση του σχετικού μεγέθους τους, της συχνότητας και του εύρους μεταβολής της επιβαλλόμενης πίεσης και πότε οδηγούν σε έλξη ή άπωση των φυσαλίδων. Στην περίπτωση ελκτικής δύναμης, ακολουθείται η κίνηση και η παραμόρφωσή τους μέχρι του σημείου που έρχονται σε επαφή, όπου αυτό είναι εφικτό. Για τη μελέτη του προβλήματος της φυσαλίδας που ανέρχεται λόγω άνωσης, υποθέτουμε αξονική συμμετρία και μόνιμη κατάσταση. Σύγκριση των προβλέψεών μας για το σχήμα των φυσαλίδων και το πεδίο ροής γύρω τους με προηγούμενα θεωρητικά και πειραματικά αποτελέσματα για Νευτωνικά ρευστά έδειξε άριστη συμφωνία. Στην περίπτωση του ιξωδοπλαστικού ρευστού εξετάστηκαν λεπτομερώς οι παραμορφώσεις των φυσαλίδων σαν συνάρτηση των αριθμών Bingham, Bond και Αρχιμήδη και υπολογίσθηκαν οι συνθήκες υπό τις οποίες είναι δυνατή η παγίδευση της φυσαλίδας μέσα σε αυτό. / The present study deals with the numerical simulation of the motion, interaction and deformation of two bubbles due to variation of the pressure of the ambient Newtonian fluid, and the buoyancy-driven rise of a bubble in a Newtonian or a viscoplastic fluid. A new elliptic mesh generation method is developed in order to deal with the singular points (poles) of the bubbles and the large deformations of their surface. This method permits us to increase the mesh resolution only in the regions that is necessary, decreasing thus the computational cost and increasing the precision of our calculations. The following problems are solved using the mixed finite element/Galerkin method. In the case of the interacting bubbles the effect of their relative size, the frequency and the width of the imposed pressure is examined as well as the conditions that lead in attraction or repulsion of the bubbles. In the case that attractive forces exist, the motion and the deformation of the bubbles followed up to the point that they come in contact, whenever this is possible. In order to study the problem of the bubble that rises due to buoyancy, axial symmetry and steady flow is assumed. Our results for the shape of the bubbles and the flow around them are in very good agreement with previous theoretical and experimental results for Newtonian fluids. The deformations of the bubbles rising in a viscoplastic material are also examined for various values of the Bingham, Bond and Archimedes numbers and the conditions under which entrapment of a bubble is possible are determined.

Στόχος της παρούσας διδακτορικής διατριβής είναι η μελέτη χρονικά μεταβαλλόμενων και ασταθών υπό συνθήκες ροών ιξωδοπλαστικών και ιξωδο-ελαστικών ρευστών λόγω της σημασίας που παρουσιάζουν τα υλικά αυτά σε διάφορες βιομηχανικές διεργασίες. Στα πλαίσια της εργασίας αυτής μελετήθηκε αρχικά το πρόβλημα της συμπίεσης ιξωδοπλαστικού ρευστού μεταξύ δύο δίσκων το οποίο είναι ένα βασικό πείραμα ρεολογικού χαρακτηρισμού μη Νευτωνικών υλικών. Η χρονικά μεταβαλλόμενη προσομοίωση επέτρεψε τον προσδιορισμό των σημαντικών διαφορών μεταξύ των δύο εκδοχών του συγκεκριμένου ρεολογικού πειράματος (δίσκοι κινούμενοι υπό σταθερή ταχύτητα ή υπό σταθερή δύναμη) που ήταν αδύνατος με τα προηγούμενα μοντέλα, ενώ παρουσιάζεται και πλήρης παραμετρική μελέτη της διεργασίας αυτής. Το υπόλοιπο μέρος της παρούσας εργασίας αφιερώθηκε στη μελέτη της διεργασίας εκβολής ιξωδοελαστικών υλικών η οποία συναντάται ευρέως στη βιομηχανία μορφοποίησης πολυμερών. Παρουσιάζεται μια πλήρης παραμετρική ανάλυση για την εκβολή ενός ιξωδοελαστικού ρευστού Phan-Tien & Tanner (PTT) από ένα δακτυλιοειδή αγωγό προκειμένου να εξετάσουμε την επίδραση της γεωμετρίας του αγωγού, και των ρεολογικών ιδιοτήτων του ρευστού. Επιπλέον ένα σημαντικό πρόβλημα που εμφανίζεται κατά την εκβολή από ένα κυλινδρικό, επίπεδο ή δακτυλιοειδή αγωγό είναι η εμφάνιση ασταθειών οι οποίες επηρεάζουν σημαντικά την ποιότητα του τελικού προϊόντος. Προκειμένου να διερευνηθεί ο μηχανισμός που προκαλεί την εμφάνιση αυτών των ασταθειών πραγματοποιήθηκε γραμμική ανάλυση ευστάθειας για τη αξονοσυμμετρική και καρτεσιανή δισδιάστατη ροή επικόλλησης-ολίσθησης (stick-slip flow) ενός ρευστού PTT. Γύρω από τη λύση μόνιμης κατάστασης πραγματοποιείται γραμμική ανάλυση ευστάθειας, και υπολογίζονται οι ιδιοτιμές του γενικευμένου προβλήματος με τη βοήθεια της μεθόδου Arnoldi. Έτσι προσδιορίζεται η επίδραση των ιδιοτήτων του υλικού στην ευστάθεια ή μη της ροής. / The purpose of this dissertation is to perform a study on time-dependent and conditionally unstable flows of viscoplastic and viscoelastic fluids because of their importance in various industrial processes. A typical experiment which is used widely for the rheological characterization of non-Newtonian fluids is squeeze flow which is the subject of the first chapter of this dissertation. The time-dependent simulation permitted the determination of the distinct differences between the two versions of this rheological experiment (disks moving with constant velocity or under constant force) which was impossible with the models used up to now, while a complete parametric analysis of this process is presented. The rest of this work focuses in the important problem of the extrusion process of a viscoelastic fluid, which is frequently encountered in the polymer industry. A complete parametric analysis is presented for the extrusion of a Phan-Thien Tanner (PTT) fluid from an annular die in order to examine the effect of the die geometry and the rheological properties of the fluid on this process. Moreover, it is widely known that during the extrusion from a cylindrical, planar or an annular die various flow instabilities may arise affecting significantly the quality of the final product. In order to investigate the mechanism, which causes the appearance of these instabilities, a linear stability analysis is performed for the cylindrical or planar stick-slip flow of a PTT fluid. The stability analysis is performed around the steady state solution and the eigenvalues of the generalized problem are calculated using the Arnoldi algorithm. With this method the effect of the various rheological properties of the fluid on the stability of the flow is determined.

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