Αποτίμηση ανελαστικής συμπεριφοράς 16- EC8 EC2 Evaluation...

Αποτίμηση της ανελαστικής συμπεριφοράς στρεπτικά ευαίσθητου 16-ώροφου κτιρίου

σχεδιασμένο με EC8 και EC2 Evaluation of the inelastic behaviour of a 16-storey torsionally sensitive building

designed according to EC8 and EC2

Γρηγόρης Γ. ΠΕΝΕΛΗΣ1, Βασίλης Κ. ΠΑΠΑΝΙΚΟΛΑΟΥ2

Λέξεις κλειδιά : Οπλισμένο σκυρόδεμα, ανελαστική ανάλυση, πεπερασμένα στοιχεία ΠΕΡΙΛΗΨΗ : Ένα δεκαεξααώροφο κτίριο που βρίσκεται υπό κατασκευή στο Βουκουρέστι σχεδιάστηκε σύμφωνα με τις διατάξεις των ΕC8 και ΕC2, χρησιμοποιώντας ελαστική δυναμική φασματική ανάλυση. Δεδομένου ότι το κτίριο είναι στρεπτικά ευαίσθητο, ελέγχθηκε περαιτέρω με τη χρήση ανελαστικής δυναμικής και στατικής ανάλυσης, βάσει ενός λεπτομερούς ραβδωτού προσομοιώματος στο χώρο. Συγκεκριμένα, εκτελέστηκαν ανελαστικές δυναμικές αναλύσεις για τρεις διαφορετικές διεγέρσεις, καθώς και αντίστοιχες στατικές ανελαστικές αναλύσεις, λαμβάνοντας υπόψη και τα στρεπτικά φαινόμενα. Τα αναλυτικά αποτελέσματα αξιολογήθηκαν τόσο σε καθολικό επίπεδο (φορέα) όσο και σε τοπικό επίπεδο (στοιχείου) και εξετάστηκε μια σειρά θεμάτων που αφορούν την επάρκεια του αρχικού ελαστικού σχεδιασμού και την εφαρμοσιμότητα των παραπάνω εξελιγμένων μεθόδων ανάλυσης. ABSTRACT : A sixteen storey building under construction in Bucharest has been designed according to the provisions of EC2 and EC8, using elastic spectral modal analysis. Considering that the building is torsionally sensitive, it was further checked and verified using nonlinear dynamic and static procedures, on a detailed space frame model. Specifically, time history analysis for three different excitations, as well as respective inelastic static analysis, taking into account torsional effects, were performed. The results were evaluated regarding both structural (global) and member (local) response and various issues concerning the adequacy of the original elastic design and the applicability of the above advanced analysis methods are discussed. 1 Δρ. Πολιτικός Μηχανικός ΑΠΘ, MSc DIC, Πενέλης Σύμβουλοι Μηχανικοί Α.Ε.,

email : [email protected] 2 Δρ. Πολιτικός Μηχανικός ΑΠΘ, MSc DIC, Πενέλης Σύμβουλοι Μηχανικοί Α.Ε.,

email : [email protected]

16ο Συνέδριο Σκυροδέματος, ΤΕΕ, ΕΤΕΚ, 21-23/10/ 2009, Πάφος, Κύπρος 1

ΕΙΣΑΓΩΓΗ

Ο στόχος της παρούσας εργασίας είναι η εφαρμογή των πλέον εξελιγμένων μη-γραμμικών αναλυτικών μεθόδων σε ένα σύνθετο πολυώροφο κτίριο από οπλισμένο σκυρόδεμα (Ο/Σ), το οποίο θεωρείται υψηλό (high-rise) κατά τα πρότυπα της νοτίου Ευρώπης και έχει σχεδιαστεί (βρίσκεται υπό κατασκευή επί του παρόντος) σύμφωνα με τους σύγχρονους ευρωπαϊκούς Κανονισμούς ΕC2 και ΕC8 (CEN, 2004α, 2004β). Το κτίριο είναι ένα δεκαεξαώροφο μεικτό σύστημα 16 ορόφων (πλαίσια και τοιχώματα) και έχει σχεδιαστεί για συντελεστές : PGA = 0.24·g, q = 4.0, β = 2.75, ψ2 = 0.4, S = 1.0, TB = 0.16 sec, TC = 1.60 sec και TD = 2.00 sec, σύμφωνα με τις διατάξεις του ΕC8 (CEN, 2004β). Η κάτοψη ξυλοτύπου και η τομή του κτιρίου φαίνονται στο Σχήμα 1 και είναι φανερό ότι η έκκεντρη θέση του ισχυρού πυρήνα στο άνω κεντρικό τμήμα της κάτοψης καθιστά την κατασκευή στρεπτικά ευαίσθητη, γεγονός που αποτυπώνεται στη δεύτερη και τρίτη ιδιομορφή, όπως θα παρουσιαστεί στην επόμενη ενότητα.

Σχήμα 1. Κάτοψη ξυλοτύπου (αριστερά) και τομή (δεξιά) του 16-ώροφου κτιρίου. Ο αρχικός ελαστικός σχεδιασμός έγινε με τη μέθοδο της δυναμικής φασματικής ανάλυσης στο πρόγραμμα πεπερασμένων στοιχείων ETABS (Computers and Structures Inc., 2008), ενώ ο ικανοτικός σχεδιασμός και η διαστασιολόγηση των δομικών στοιχείων έγινε με το λογισμικό eTools (Πενέλης Λογισμικά ΕΠΕ, 2009).


Για τη στατική και δυναμική ανελαστική ανάλυση του κτιρίου χρησιμοποιήθηκε το πρόγραμμα πεπερασμένων στοιχείων Zeus-NL (Elnashai et al., 2009), το οποίο περιλαμβάνει ραβδωτά πεπερασμένα στοιχεία κατανεμημένης πλαστικότητας με προσομοίωση ινών και κατάλληλους καταστατικούς νόμους για σκυρόδεμα και χάλυβα. Περισσότερες λεπτομέρειες όσον αφορά τη διαδικασία προσομοίωσης θα δοθούν στην επόμενη ενότητα, ωστόσο αξίζει να σημειωθεί ότι οι υπολογιστικές απαιτήσεις της παρούσας ανάλυσης ξεπέρασαν, σύμφωνα με την άποψη των συγγραφέων, προηγούμενες αντίστοιχες εφαρμογές του ιδίου ή παρόμοιων προγραμμάτων, γεγονός που απαίτησε κατάλληλες προσαρμογές στον πηγαίο κώδικα του εν λόγω λογισμικού. Από τη διαδικασία εφαρμογής των ανελαστικών μεθόδων, εξετάστηκαν διάφορα θέματα σχετικά με την επάρκεια των διατάξεων του EC8 σχετικά με τα στρεπτικά ευαίσθητα κτίρια, όπως η διαθέσιμη υπεραντοχή, η απαιτούμενη πλαστιμότητα (καθολική και τοπική) και το αναμενόμενο σχήμα παραμόρφωσης του φορέα (σχετικές παραμορφώσεις στο κέντρο μάζας και στα όρια της κάτοψης). Επιπλέον, αντιμετωπίστηκαν αρκετά θέματα που αφορούν την ακριβή προσομοίωση ενός πραγματικού κτιριακού φορέα στο πλαίσιο της μη γραμμικής ανάλυσης και συγκεκριμένα η προσομοίωση των διαφραγμάτων, των σύνθετων πυρήνων Ο/Σ και της κατανομής μαζών. Τα παραπάνω θέματα αναδεικνύουν την ανάγκη χρήσης εμπορικών πακέτων λογισμικού με απλοποιημένες προσομοιώσεις (όπως για παράδειγμα στοιχεία με σημειακές πλαστικές αρθρώσεις) για τις ανάγκες ελέγχου πραγματικών κτιρίων σε περιβάλλον γραφείου, σε αντίθεση με την ανάλυση πρωτότυπων προσομοιωμάτων με εξελιγμένα προγράμματα πεπερασμένων στοιχείων για ερευνητικούς σκοπούς. ΠΡΟΣΟΜΟΙΩΣΗ ΠΕΠΕΡΑΣΜΕΝΩΝ ΣΤΟΙΧΕΙΩΝ ΓΙΑ ΜΗ

ΓΡΑΜΜΙΚΗ ΑΝΑΛΥΣΗ Στην ενότητα αυτήν, περιγράφεται συνοπτικά η διαδικασία προσομοίωσης του κτιριακού φορέα με ραβδωτά πεπερασμένα στοιχεία. Ο φορέας προσομοιώθηκε εξαιρώντας τα τέσσερα υπόγεια και θεωρήθηκε πακτωμένος στο επίπεδο του δαπέδου ισογείου. Αποτελείται από 16 ορόφους, έναν ημιώροφο και το θάλαμο απόληξης του ανελκυστήρα. Οι διαστάσεις της κάτοψης είναι 32.2 × 18.8 m και το συνολικό ύψος ανέρχεται σε 58.25 m. Ο φορέας προσομοιώθηκε με το πρόγραμμα πεπερασμένων στοιχείων Zeus-NL (Elnashai et al., 2009) (Σχήμα 2), το οποίο περιλαμβάνει ραβδωτά στοιχεία κατανεμημένης πλαστικότητας καθ’ ύψος της διατομής και κατά μήκος του στοιχείου, βασισμένα στη προσομοίωση ινών (fiber model). Επιπλέον, έχει τη δυνατότητα προσδιορισμού της ανελαστικής συμπεριφοράς της κατασκευής υπό στατικά ή δυναμικά φορτία, λαμβάνοντας υπόψη και τα φαινόμενα δευτέρας τάξης (γεωμετρική μη-γραμμικότητα). Για την προσομοίωση του σκυροδέματος, θεωρήθηκε ένας μονοαξονικός καταστατικός νόμος σταθερής περίσφιξης, με θλιπτική αντοχή fc = 32.0 MPa, εφελκυστική αντοχή ft = 3.0 MPa, παραμόρφωση στη μέγιστη αντοχή ίση με


εc = 0.0022 και σταθερό συντελεστή περίσφιξης Κ = 1.2. Για την προσομοίωση του χάλυβα, θεωρήθηκε ένας διγραμμικός ελαστοπλαστικός νόμος με κινηματική κράτυνση, μέτρο ελαστικότητας Es = 200 GPa, τάση διαρροής fy = 500 MPa και συντελεστή κράτυνσης 5 ‰. Οι παραπάνω ιδιότητες των υλικών βρίσκονται σε συμφωνία με τα αντίστοιχα μεγέθη σχεδιασμού, ήτοι σκυρόδεμα C32/40 και χάλυβας B500C.

Ραβδωτά στοιχεία

Όψη στερεού

Σχήμα 2. Προσομοίωση του φορέα με ραβδωτά πεπερασμένα στοιχεία.


Οι διατομές οπλισμένου σκυροδέματος (γεωμετρία και τοπολογία ράβδων οπλισμού) ορίστηκαν σύμφωνα με τα σχέδια ξυλοτύπων και ανατέθηκαν σε κυβικά ελαστο-πλαστικά ραβδωτά στοιχεία, των οποίων το εφαπτομενικό μητρώο δυσκαμψίας προκύπτει από ολοκλήρωση 2 σημείων Gauss. Κατά συνέπεια, το μήκος των στοιχείων αυτών είναι καθοριστικό για την αξιόπιστη αποτύπωση της μη γραμμικής συμπεριφοράς στις κρίσιμες ζώνες. Για το λόγο αυτό, όλα τα φυσικά δομικά στοιχεία του φορέα (δοκοί, υποστυλώματα και τοιχώματα) διακριτοποιήθηκαν σε 4 ραβδωτά στοιχεία με διαίρεση 15 – 35 – 35 – 15 % κατά μήκος, γεγονός που οδηγεί σε πυκνότερη διακριτοποίηση στις περιοχές σύνδεσης δοκών-υποστυλωμάτων/τοιχωμάτων, όπου οι δράσεις και οι μετακινήσεις αναμένονται να είναι αυξημένες. Τα μεμονωμένα τοιχώματα και οι τοιχωματικοί πυρήνες προσομοιώθηκαν ανά σκέλος με κάθετα ραβδωτά στοιχεία (στο κέντρο μάζας κάθε σκέλους) και η οριζόντια κινηματική δέσμευση μεταξύ σκελών και μεταξύ σκέλους-δοκού προσομοιώθηκε με άκαμπτα στοιχεία (ελαστικό υλικό με ERigid = 10·Es και διατομή 1.0×1.0 m) (Σχήμα 3, αριστερά). Θεωρήθηκε επίσης δράση άκαμπτων διαφραγμάτων σε όλα τα επίπεδα ορόφων, ωστόσο η έμμεση προσομοίωσή τους με κινηματικές δεσμεύσεις τύπου master-slave δεν ήταν διαθέσιμη στο λογισμικό, οπότε προσομοιώθηκε με άμεσο τρόπο χρησιμοποιώντας χιαστί άκαμπτους αρθρωτούς συνδέσμους (Σχήμα 3, δεξιά).

Σχήμα 3. Προσομοίωση σκέλους τοιχώματος (αριστερά) και διαφραγμάτων (δεξιά). Η συνολική φόρτιση βαρύτητας σε κάθε όροφο υπολογίστηκε από το συνδυασμό G+ψ2·Q του αρχικού σχεδιασμού και κατανεμήθηκε στους εσώτερους 3 κόμβους κάθε δοκού και στους ακραίους κόμβους κάθε σκέλους τοιχώματος, σύμφωνα με τις αντίστοιχες επιφάνειες επιρροής (Σχήμα 4, αριστερά). Ειδικότερα για τη δυναμική ανάλυση, η συνολική μάζα κάθε ορόφου m = (G+ψ2·Q)/g κατανεμήθηκε στους ακραίους κόμβους όλων των κατακόρυφων στοιχείων (υποστυλώματα και τοιχώματα), σύμφωνα με τις αντίστοιχες επιφάνειες επιρροής (Σχήμα 4, δεξιά). Η οριζόντια φόρτιση για την ανελαστική στατική ανάλυση


εφαρμόστηκε στο κέντρο μάζας κάθε ορόφου (C.M., Σχήμα 2), ενώ για τη δυναμική ανάλυση τοποθετήθηκαν τα αντίστοιχα επιταχυνσιογραφήματα ταυτόχρονα και στις δύο διευθύνσεις (x και y) σε όλους τους κόμβους της βάσης.

Σχήμα 4. Κατανομή φορτίων βαρύτητας (αριστερά) και μαζών (δεξιά). Πρέπει να σημειωθεί στο σημείο αυτό ότι το παρόν προσομοίωμα, κατά την άποψη των Συγγραφέων, είναι το μεγαλύτερο που έχει υλοποιηθεί στο πλαίσιο της μη γραμμικής ανάλυσης ραβδωτών φορέων με πεπερασμένα στοιχεία κατανεμημένης πλαστικότητας. Για το λόγο αυτό απαιτήθηκαν κατάλληλες προσαρμογές στον πηγαίο κώδικα του προγράμματος Zeus-NL για την αντιμετώπιση των ιδιαίτερα αυξημένων υπολογιστικών απαιτήσεων. Επιπλέον, αναπτύχθηκε ένας ειδικός γραφικός προ/μετα-επεξεργαστής για την απεικόνιση του στερεού προσομοιώματος στο χώρο (βλ. Σχήμα 2), για λόγους εποπτείας και ελέγχου της ορθότητας της γεωμετρίας του προσομοιώματος, αλλά και της αξιοπιστίας των αναλυτικών αποτελεσμάτων. Μερικά στατιστικά στοιχεία του μεγέθους της παρούσας ανάλυσης φαίνονται στον Πίνακα 1. Για τον έλεγχο της αξιοπιστίας του προσομοιώματος, αρχικά εκτελέστηκε μια ιδιομορφική ανάλυση και τα αποτελέσματα αυτής συγκρίθηκαν ποιοτικά και ποσοτικά με τα αντίστοιχα αποτελέσματα του προγράμματος ΕΤΑBS, με το οποίο έγινε ο αρχικός ελαστικός σχεδιασμός. Παρατηρείται ικανοποιητική σύγκλιση τόσο ως προς τις ιδιοπεριόδους (Πίνακας 2) όσο και ως προς τις ιδιομορφές. Η πρώτη ιδιομορφή προέκυψε μεταφορική στην ασθενή διεύθυνση (y) του κτιρίου, η δεύτερη προέκυψε μεικτή (μεταφορική/στρεπτική) και η τρίτη μόνο στρεπτική (Σχήμα 5). Λαμβάνοντας υπόψη ότι τα προγράμματα Zeus-NL και ETABS είναι τελείως διαφορετικά ως προς τους καταστατικούς νόμους υλικών και τις αρχές προσομοίωσης (π.χ. στο ETABS τα τοιχώματα και οι πλάκες προσομοιώθηκαν με επιφανειακά στοιχεία κελύφους, τα διαφράγματα με κινηματικές δεσμεύσεις τύπου master-slave και τα φορτία βαρύτητας εφαρμόζονται ως επιφανειακές φορτίσεις στις πλάκες), θεωρείται ότι η παρούσα προσομοίωση είναι επιτυχής και αξιόπιστη, ώστε στη συνέχεια να εφαρμοστεί για τη στατική και δυναμική ανελαστική ανάλυση του κτιρίου, όπως θα παρουσιαστεί στις επόμενες ενότητες.


Πίνακας 1. Στοιχεία προσομοιώματος

Τύπος ανάλυσης Στοιχείο Ιδιομορφική Στατική Δυναμική Κόμβοι 4897

Διατομές Ο/Σ 339 Ραβδωτά στοιχεία 5424 Στοιχεία αρθρώσεων 968

Στοιχεία μάζας 508 − 508 Επικόμβιες φορτίσεις βαρύτητας − 1102 1102 Επικόμβιες οριζόντιες φορτίσεις − 51 −

Επικόμβιες φορτίσεις επιταχύνσεων − − 60 Απαιτούμενη μνήμη RAM ≈ 160,000,000 32-bit words ≈ 640 MB

Πίνακας 2. Σύγκριση ιδιοπεριόδων μεταξύ των προγραμμάτων Zeus-NL και ETABS

Ιδιομορφή Zeus-NL (sec) ETABS (sec) 1 1.03 0.98 2 0.81 0.73 3 0.53 0.56 4 0.25 0.22 5 0.21 0.21 6 0.16 0.16 7 0.14 0.12 8 0.13 0.09

1η ιδιομορφή (1.03 sec)



Σχήμα 5. Ιδιομορφές κτιρίου όπως προέκυψαν από το πρόγραμμα Zeus-NL


ΑΝΕΛΑΣΤΙΚΗ ΔΥΝΑΜΙΚΗ ΑΝΑΛΥΣΗ Γενικά Για τη σύγκριση μεταξύ δυναμικών και στατικών ανελαστικών αναλύσεων, η αρχική ιδέα ήταν η κατασκευή της ‘δυναμικής καμπύλης αντίστασης’, με την εκτέλεση πολλαπλών δυναμικών αναλύσεων (Papanikolaou & Elnashai, 2005). Ωστόσο, το απαιτούμενο υπολογιστικό κόστος για το παρόν προσομοίωμα αποδείχτηκε απαγορευτικό και τελικά επιλέχθηκε η εφαρμογή τριών ζευγών διεγέρσεων (στις διευθύνσεις x και y), ώστε να προσδιοριστεί η μέγιστη μετακίνηση κορυφής, το σχήμα παραμόρφωσης του φορέα και η απαιτούμενη πλαστιμότητα. Σεισμικές διεγέρσεις Οι σεισμικές διεγέρσεις που επιλέχθηκαν ήταν ο σεισμός Vrancea του 1977 (συνιστώσες ΝS & EW) και ένα συνθετικό επιταχυνσιογράφημα με συμβατό ως προς τον EC8 συχνοτικό περιεχόμενο (Σχήμα 6). Ο σεισμός του 1977 επιλέχθηκε αυτούσιος γιατί αντιστοιχεί στο πραγματικό συμβάν που έπληξε την πόλη του Βουκουρεστίου και θεωρείται σημαντικός για την αξιολόγηση της συμπεριφοράς του εν λόγω υψηλού κτιρίου επειδή το ελαστικό του φάσμα παρουσιάζει σημαντική ενίσχυση στην περιοχή του 1.0 sec, η οποία αντιστοιχεί στη θεμελιώδη ιδιομορφή της κατασκευής (Σχήματα 7 & 8). Απεναντίας, το συνθετικό επιταχυνσιογράφημα εισήχθη ανηγμένο σε PGA = 0.24·g, σύμφωνα με το ρουμανικό εθνικό προσάρτημα του EC8 για την περιοχή του Βουκουρεστίου. Διαδικασία ανάλυσης Διενεργήθηκαν τρεις δυναμικές αναλύσεις, δύο με τη χρήση του συμβατού προς τον EC8 συνθετικού επιταχυνσιογραφήματος και μια με τη διέγερση Vrancea-1977 ως εξής : (α) 100 % του EC8 κατά τη διεύθυνση x και 30 % του ΕC8 κατά τη διεύθυνση y (EC8 x + 0.3·y) (β) 100 % του ΕC8 κατά τη διεύθυνση y και 30 % του EC8 κατά τη διεύθυνση x (EC8 0.3·x + y) (γ) Συνιστώσα Vrancea NS κατά τη διεύθυνση x και EW κατά τη διεύθυνση y

(Vrancea x:NS + y:EW) Τα αναλυτικά αποτελέσματα που θα παρουσιαστούν στη σχετική ενότητα (συγκρινόμενα με τα αντίστοιχα της στατικής ανάλυσης) επιλέχθηκαν ως εξής : (α) Τέμνουσα βάσης – μετακίνηση κορυφής (καμπύλη αντίστασης P-δ). (β) Προφίλ μέγιστων μετακινήσεων ορίων κάτοψης για όλες τις διεγέρσεις. (γ) Διαγράμματα ροπών-στροφών χορδής για το υποστύλωμα C13 (βλ. Σχήμα 2).


-0.3

-0.2

-0.1

0

0.1

0.2

0.3

0 2 4 6 8 10

EC8 Συνθετικό (PGA = 0.24 sec)

Χρόνος (sec)

a (g)

0

0.1

0.2

0.3

0.4

0.5

0.6

0.7

0.8

0 0.5 1 1.5 2 2.5 3 3.5 4

EC8 Συνθετικό (PGA = 0.24 sec)

T (sec)

Sa (g)

Ελαστικό φάσμα 5 % απόσβεση

Σχήμα 6. Συνθετικό επιταχυνσιογράφημα EC8, ανηγμένο σε 0.24·g (αριστερά) καιελαστικό φάσμα με απόσβεση 5 % (δεξιά).

-0.3

-0.2

-0.1

0

0.1

0.2

0.3

0 5 10 15 20

Vrancea N-S (PGA = 0.20 sec)

a (g)

Χρόνος (sec) 0

0.1

0.2

0.3

0.4

0.5

0.6

0.7

0.8

0 0.5 1 1.5 2 2.5 3 3.5 4

Vrancea N-S (PGA = 0.20 sec)

T (sec)

Sa (g)


Σχήμα 7. Επιταχυνσιογράφημα Vrancea-NS (αριστερά) και ελαστικό φάσμα με απόσβεση 5 % (δεξιά).

-0.3

-0.2

-0.1

0

0.1

0.2

0.3

0 5 10 15 20

Vrancea E-W (PGA = 0.165 sec)

a (g)

Χρόνος (sec) 0

0.1

0.2

0.3

0.4

0.5

0.6

0.7

0.8

0 0.5 1 1.5 2 2.5 3 3.5 4

Vrancea E-W (PGA = 0.165 sec)

T (sec)

Sa (g)


Σχήμα 8. Επιταχυνσιογράφημα Vrancea-EW (αριστερά) και ελαστικό φάσμα με απόσβεση 5 % (δεξιά).

ΑΝΕΛΑΣΤΙΚΗ ΣΤΑΤΙΚΗ ΑΝΑΛΥΣΗ Γενικά Για τη στατική ανελαστική ανάλυση, εφαρμόστηκε η μεθοδολογία των Penelis & Kappos (2002, 2005) για τον υπολογισμό των φασματικών φορτίσεων καθώς και του ισοδύναμου μονοβάθμιου ταλαντωτή, συμπεριλαμβάνοντας τόσο τις μεταφορικές όσο και τις στρεπτικές συνιστώσες απόκρισης του κτιρίου.


Ειδικότερα, τα στατικά φορτία (P) σε κάθε όροφο προσδιορίζονται από την παρακάτω σχέση (σε μητρωική μορφή) : i i aiP Γ φ S= ⋅ ⋅∑ (1) όπου : Γi συντελεστής συμμετοχής της μάζας για την ιδιομορφή (i) φi ιδιοδιάνυσμα της ιδιομορφής (i) Sai φασματική επιτάχυνση για την ιδιομορφή (i) Οι συντελεστές μετατροπής (c1 και c2) του ισοδύναμου μονοβάθμιου ταλαντωτή (SDOF) δίνονται από την εξίσωση (2) για τις μετακινήσεις και από την εξίσωση (3) για τις φορτίσεις :

2 2 2 2Τxo yo zoο ο

1xo yo

m φ m φ m r φφ Μφ cm m (φ φ )

⋅ + ⋅ + ⋅ ⋅= =

∗ ⋅ + (2)

Τ

xo pox yo poy zo Moο2

po pox poy

φ ψ φ ψ φ ψφ ψ cψ ψ ψ

⋅ + ⋅ + ⋅= =

+ (3)

όπου : m μάζα r ακτίνα αδράνειας

xo

yο

zo

φφ

φ

⎡ ⎤⎢ ⎥⎢ ⎥⎢ ⎥⎣ ⎦

ανηγμένες ιδιομορφικές φασματικές μετακινήσεις

pox

pοy

Mo

ψ

ψ

ψ

⎡ ⎤⎢ ⎥⎢ ⎥⎢ ⎥⎣ ⎦

ανηγμένα ιδιομορφικά φασματικά φορτία

Διαδικασία ανάλυσης Διενεργήθηκαν τρεις στατικές ανελαστικές αναλύσεις (Pushover), δύο συμβατές με το επιταχυνσιογράφημα του EC8 και μία συμβατή με το επιταχυνσιογράφημα Vrancea. Οι παραπάνω αναλύσεις πραγματοποιήθηκαν τόσο κατά τη θετική όσο και κατά την αρνητική φορά κάθε διεύθυνσης (± x , ± y) ώστε να αποτυπωθεί η πιθανή ασυμμετρία στη στατική ανελαστική απόκριση του φορέα. (α) 100 % του EC8 κατά τη διεύθυνση x και 30 % του ΕC8 κατά τη διεύθυνση y ± (EC8 x + 0.3·y)


(β) 100 % του ΕC8 κατά τη διεύθυνση y και 30 % του EC8 κατά τη διεύθυνση x ± (EC8 0.3·x + y) (γ) Συνιστώσα Vrancea NS κατά τη διεύθυνση x και EW κατά τη διεύθυνση y

± (Vrancea x:NS + y:EW) Από την προηγούμενη παράγραφο είναι φανερό ότι προσδιορίστηκαν διαφορετικά στατικά φορτία και συντελεστές μετατροπής για κάθε σεισμική διέγερση. Οι καμπύλες απόκρισης που προέκυψαν από την ανάλυση αρχικά διγραμμικοποιήθηκαν και μετατράπηκαν σε καμπύλες φασματικής επιτάχυνσης-μετακίνησης (ADRS). Στη συνέχεια προσδιορίστηκε η μετακίνηση-στόχος του ισοδύναμου μονοβάθμιου ταλαντωτή με τη χρήση ανελαστικών ικανοτικών φασμάτων σταθερής πλαστιμότητας και τέλος υπολογίστηκε η μετακίνηση-στόχος του πολυβάθμιου ταλαντωτή (MDOF). Η παραπάνω διαδικασία απεικονίζεται σχηματικά στο Σχήμα 9. Τα μεγέθη απόκρισης της στατικής ανελαστικής ανάλυσης, που παρουσιάζονται στην επόμενη ενότητα, αντιστοιχούν στο φορτιστικό βήμα στο οποίο επιτυγχάνεται η μετακίνηση-στόχος του πολυβάθμιου ταλαντωτή.

Σχήμα 9. Προσδιορισμός της μετακίνησης-στόχου του πολυβάθμιου ταλαντωτή με τη χρήση της μεθόδου ADRS.


ΑΝΑΛΥΤΙΚΑ ΑΠΟΤΕΛΕΣΜΑΤΑ Απόκριση σε καθολικό επίπεδο Τα αναλυτικά αποτελέσματα σε καθολικό επίπεδο που παρουσιάζονται είναι τα εξής : (α) Στατικές και δυναμικές καμπύλες αντίστασης (P-δ) (Σχήμα 10) στις οποίες φαίνεται ότι υπάρχει ικανοποιητική σύγκλιση μεταξύ των δύο μεθόδων ανάλυσης (η στατική καμπύλη είναι περιβάλλουσα της δυναμικής) και ότι ο φορέας παρουσιάζει, όπως αναμενόταν, ασύμμετρη απόκριση κατά τη θετική και αρνητική φορά, λόγω της στρεπτικής του ευαισθησίας. Η μέγιστη μετακίνηση οροφής από τη δυναμική ανάλυση για διέγερση συμβατή με τον EC8 είναι 15 cm στη διεύθυνση x και 30 cm στη διεύθυνση y (Πίνακας 3), γεγονός που καθιστά ως ασφαλή εκτίμηση τα 36 cm και 68 cm αντίστοιχα, που προέκυψαν από τον αρχικό ελαστικό σχεδιασμό. (β) Μετακινήσεις ορόφων στις διευθύνσεις x και y για τα όρια της κάτοψης (x : άνω και κάτω όριο, y : αριστερό και δεξιό όριο, βλ. Σχήμα 2) στις οποίες είναι φανερή η στρεπτική απόκριση του φορέα (Σχήμα 11). Στο ίδιο σχήμα παρουσιάζεται και ο υπολογισμός της μετακίνησης-στόχου για τη στατική ανελαστική ανάλυση, σύμφωνα με τη μέθοδο ADRS. Επιπλέον, από τα παραπάνω σχήματα προκύπτουν οι ποσοστιαίες διαφορές των μετακινήσεων μεταξύ δυναμικής και στατικής ανάλυσης (Πίνακας 3) και είναι ενθαρρυντικό ότι οι μεγαλύτερες από αυτές κυμαίνονται μεταξύ 15 % και 20 %, ανάλογα με τη διέγερση και τη θέση μέτρησης. Πίνακας 3. Σύγκριση μετακινήσεων κορυφής στο κέντρο μάζας και στα όρια του φορέα

μεταξύ στατικής και δυναμικής ανελαστικής ανάλυσης

Μετακινήσεις στο κέντρο μάζας (m)

Διέγερση Δυναμική (μέγιστη)

Στατική (στη μετακ. στόχο) Διαφορά

EC8 x + 0.3·y (δx) 0.151 0.141 6.8 % EC8 0.3·x + y (δy) 0.302 0.267 11.8 %

Vrancea x:NS + y:EW (δx) 0.261 0.283 8.4 % Μετακινήσεις στο μέσον του ορίου (m)

Διέγερση Δυναμική (μέγιστη)

Στατική (στη μετακ. στόχο) Διαφορά

EC8 x + 0.3·y (δx) - Άνω 0.123 0.100 18.3 % EC8 x + 0.3·y (δx) – Κάτω 0.179 0.181 1.0 %

EC8 0.3·x + y (δy) - Αριστερά 0.270 0.210 22.2 % EC8 0.3·x + y (δy) - Δεξιά 0.334 0.323 3.4 %

Vrancea x:NS + y:EW (δx) - Άνω 0.222 0.219 1.4 % Vrancea x:NS + y:EW (δx) - Κάτω 0.300 0.347 15.7 %


-60

-40

-20

0

20

40

60

-0.6 -0.4 -0.2 0 0.2 0.4 0.6

Δυναμική

Στατική

Px (MN)

δx (m)

Vrancea x:NS + y:EW Διεύθυνση x

Σχήμα 10. Σύγκριση μεταξύ στατικής και δυναμικής ανελαστικής ανάλυσης για καθολική απόκριση φορέα (P-δ) και για τη διέγερση Vrancea x:NS + y:EW.

0.00

0.05

0.10

0.15

0.20

0.00 0.20 0.40 0.60 0.80 1.00

Sa (g)

Sd (m)

μ = 2 μ = 1

Sd,t = 0.31 → dt = 0.29 m

Vrancea x:NS + y:EW Διεύθυνση x

0

10

20

30

40

50

0.00 0.10 0.20 0.30 0.40

Δυναμική

Στατική σε dt = 0.29 m

h (m)

δx (m)

Άνω Κάτω58.25

Vrancea x:NS + y:EW

Διεύθυνση x

Σχήμα 11. Προσδιορισμός της μετακίνησης-στόχου και προφίλ μετακινήσεων για τη διέγερση Vrancea x:NS + y:EW. Απόκριση σε τοπικό επίπεδο Η αξιολόγηση των αναλυτικών αποτελεσμάτων σε τοπικό επίπεδο έγινε για το γωνιακό υποστύλωμα C13 (άνω-δεξιά γωνία, βλ. Σχήμα 2), ώστε να διερευνηθεί η ανελαστική διαξονική ένταση της διατομής υπό μεταβαλλόμενο αξονικό φορτίο. Στο Σχήμα 12 παρουσιάζονται τα διαγράμματα ροπής – στροφής χορδής για στατική και δυναμική ανάλυση, όπου φαίνεται η μέγιστη απαίτηση της δυναμικής απόκρισης συγκρινόμενη με την αντίστοιχη απαίτηση της στατικής απόκρισης στη μετακίνηση-στόχο. Στο Σχήμα 13 φαίνεται το ιστορικό απόκρισης


των μεγεθών Ν, Μx, My του υποστυλώματος C13 (για στατική ανάλυση συμβατή με τη διέγερση EC8 x + 0.3·y) καθώς και οι αντίστοιχες καμπύλες αλληλεπίδρασης στις θέσεις των μέγιστων τιμών τους. Τα διαγράμματα αυτά δικαιολογούν τη μεγάλη διαφοροποίηση των ροπών αντοχής ανά διεύθυνση που παρατηρείται στο Σχήμα 12, αν και πρόκειται για μια σχεδόν συμμετρική διατομή διαστάσεων 1.20 × 1.30 m. Ωστόσο, θεωρείται ότι η σύγκλιση μεταξύ δυναμικής και στατικής απόκρισης σε τοπικό επίπεδο είναι ενγένει ικανοποιητική, λαμβάνοντας υπόψη τις διαφορετικές διακυμάνσεις του αξονικού φορτίου μεταξύ των δύο μεθόδων.

-12000

-8000

-4000

0

4000

8000

12000

-0.012 -0.008 -0.004 0 0.004 0.008 0.012

Δυναμική Στατική Μετακίνηση-στόχος

ry (rad)

My (kNm)

Υποστύλωμα C13 - Διέυθυνση xEC8 x + 0.3·y

-12000

-8000

-4000

0

4000

8000

12000

-0.012 -0.008 -0.004 0 0.004 0.008 0.012

Δυναμική Στατική Μετακίνηση-στόχος

rx (rad)

Mx (kNm)

Υποστύλωμα C13 - Διέυθυνση yEC8 x + 0.3·y

Σχήμα 12. Σύγκριση μεταξύ στατικής και δυναμικής ανελαστικής ανάλυσης για τοπική απόκριση του υποστυλώματος C13 (M-r) και για τη διέγερση EC8 x + 0.3·y.

ΣΥΜΠΕΡΑΣΜΑΤΑ Από τα αναλυτικά αποτελέσματα που παρουσιάστηκαν στην προηγούμενη ενότητα, μπορούν να εξαχθούν τα ακόλουθα συμπεράσματα σχετικά με την επάρκεια του σχεδιασμού του κτιρίου σύμφωνα με τις απαιτήσεις του EC8, όπως αυτά προέκυψαν από την εφαρμογή των εξελιγμένων μεθόδων ανελαστικής δυναμικής και στατικής ανάλυσης : (α) Το κτίριο διαθέτει την απαιτούμενη υπεραντοχή και πλαστιμότητα, η οποία οδηγεί σε συντελεστή συμπεριφοράς q = 4, όπως θεωρήθηκε και κατά τον ελαστικό σχεδιασμό. (β) Η απόκριση σύμφωνα με το σχεδιασμό κατά ΕC8 (μετακινήσεις και δυνάμεις ορόφων) είναι συντηρητική (κατά περίπου 100 %) σε σχέση με την ανελαστική απόκριση, γεγονός που είναι αναμενόμενο λόγο των προβλεπόμενων κατά ΕC8 συντελεστών ασφάλειας. (γ) Η ανελαστική τέμνουσα αντοχής στη βάση του κτιρίου είναι περίπου 50 ΜΝ για τη διεύθυνση x και 40 ΜΝ για τη διεύθυνση y, ενώ η τέμνουσα βάσης σχεδιασμού είναι περίπου 24 ΜΝ.


(δ) Η τοπική ανελαστική απαίτηση των γωνιακών υποστυλωμάτων, όπως αυτή ελέγχθηκε για το υποστύλωμα C13, αποδεικνύει ότι έχουν σχεδιαστεί επαρκώς.

0

2000

4000

6000

8000

10000Mx (kNm)

0

2000

4000

6000

8000

10000Φορτιστικό βήμα

My (kNm)

0

20000

40000

60000

0 50 100 150

N (kN)

EC8 X+0.3·Y - Υποστύλωμα C13Ιστορικό απόκρισης

Φορτία βαρύτητας

-12000

-8000

-4000

0

4000

8000

12000

-12000 -8000 -4000 0 4000 8000 12000

My

Mx

N = -44385 kN-12000

-8000

-4000

0

4000

8000

12000

-12000 -8000 -4000 0 4000 8000 12000

My

Mx

N = -46397 kN-12000

-8000

-4000

0

4000

8000

12000

-12000 -8000 -4000 0 4000 8000 12000

My

Mx

N = -50388 kN

Μέγιστη τιμή Mx = 9418 kNmΑντίστοιχη τιμή My = 3034 kNmΑντίστοιχη τιμή N = -44385 kN

Μέγιστη τιμή My = 3080 kNmΑντίστοιχη τιμή Mx = 8833 kNmΑντίστοιχη τιμή N = -46397 kN

Μέγιστη τιμή N = -50388 kNmΑντίστοιχη τιμή Mx = 2185 kNmΑντίστοιχη τιμή My = 2103 kN

Σχήμα 13. Ιστορικό απόκρισης του γωνιακού υποστυλώματος C13 για στατική ανελαστική ανάλυση συμβατή με τη διέγερση EC8 x + 0.3·y και αντίστοιχες διαξονικές καμπύλες αλληλεπίδρασης στις θέσεις μέγιστης απόκρισης των μεγεθών Ν, Μx, My. Από τα παραπάνω θεωρείται ασφαλές το συμπέρασμα ότι το παρόν υψηλό και στρεπτικά ευαίσθητο κτίριο Ο/Σ, το οποίο σχεδιάστηκε σύμφωνα με τις διατάξεις


των EC2 και ΕC8 είναι ικανό να ανταποκριθεί, με ένα συντελεστή ασφάλειας περίπου ίσο με 2, στις απαιτήσεις τόσο της σεισμικής διέγερσης συμβατής με το φάσμα του ΕC8 όσο και της διέγερσης Vrancea, η οποία αποτελεί το κρισιμότερο σεισμολογικά γεγονός που έχει καταγραφεί στην περιοχή. Από τη σύγκριση μεταξύ δυναμικής και στατικής ανελαστικής ανάλυσης, είναι προφανές ότι τόσο σε καθολικό όσο σε τοπικό επίπεδο, η λιγότερο χρονοβόρα και απαιτητική υπολογιστικά στατική ανελαστική ανάλυση είναι σε θέση να παράγει αποτελέσματα που συγκλίνουν ικανοποιητικά με τα αντίστοιχα της ακριβέστερης δυναμικής ανάλυσης, καθιστώντας την χρηστική, ακόμα και στις περιπτώσεις στρεπτικά ευαίσθητων κτιρίων. Τέλος, από την προσπάθεια προσομοίωσης του παρόντος σύνθετου φορέα για τη διενέργεια ανελαστικών αναλύσεων, προέκυψαν τα παρακάτω θέματα προς περαιτέρω διερεύνηση : (α) Η λεπτομέρεια στην προσομοίωση των διατομών Ο/Σ παίζει σημαντικό ρόλο όπως π.χ. το μέγεθος και η θέση των οπλισμών, οι απερίσφικτες και περισφιγμένες περιοχές, η προσομοίωση των διατομών τοιχείων και πλακοδοκών κλπ. (β) Η προσομοίωση των πυρήνων Ο/Σ είναι ιδιαίτερα σημαντικό θέμα, το οποίο εξαρτάται άμεσα από τη μέθοδο ανελαστικής ανάλυσης που έχει επιλεγεί. Σε προσομοιώματα ινών κατανεμημένης πλαστικότητας, όπου λαμβάνεται υπόψη η αξονική ανελαστική απόκριση (Ν-δ), οι πυρήνες μπορούν να προσομοιωθούν ως σύνολο μεμονωμένων σκελών, συνδεδεμένων μεταξύ τους με οριζόντια άκαμπτα στοιχεία. Αντίθετα, στην κλασική θεώρηση σημειακών πλαστικών αρθρώσεων, στην οποία οι ροπές δεν μπορούν να λειτουργήσουν ως ζεύγος θλιπτικών-εφελκυστικών αξονικών δυνάμεων (εφόσον δε δίνεται νόμος Ν-δ), ο πυρήνας θα πρέπει να προσομοιώνεται με ένα μεμονωμένο στοιχείο, με διάγραμμα ροπών-στροφών/καμπυλοτήτων που αντιστοιχεί στη σύνθετη διατομή. (γ) Η προσομοίωση των άκαμπτων συνδέσμων (πραγματική προς “αριθμητική” ακαμψία). (δ) Η προσομοίωση των άκαμπτων διαφραγμάτων (π.χ. με τη χρήση κινηματικών δεσμεύσεων ή χιαστί συνδέσμων). (ε) Οι καταστατικοί νόμοι για σκυρόδεμα και χάλυβα, ιδιαίτερα η κράτυνση του χάλυβα συγκρινόμενη με την αναμενόμενη κράτυνση του Ο/Σ ως σύνθετου υλικού. (στ) Η κατανομή των φορτίων βαρύτητας και η κατανομή των μαζών στη δυναμική ανάλυση. (ζ) Ο υπολογισμός και η κατανομή των στατικών οριζόντιων φορτίων και στρεπτικών ροπών, ώστε να αποφευχθούν αριθμητικές αστάθειες.


(η) Η αποδοτική μετεπεξεργασία των αναλυτικών αποτελεσμάτων τα οποία, ιδιαίτερα στη δυναμική ανάλυση, παράγονται σε μεγάλο όγκο. Όλα τα παραπάνω θέματα αποδεικνύουν ότι, επί του παρόντος, η μοναδική πρακτική εναλλακτική λύση στον ελαστικό σχεδιασμό των κτιριακών κατασκευών είναι η απλοποιημένη στατική ανελαστική ανάλυση με θεώρηση σημειακών πλαστικών αρθρώσεων και τη χρήση σχετικών συστάσεων όπως της FEMA 356 (BSSC, 2000), τη στιγμή που οι πλέον εξελιγμένες μέθοδοι εμπεριέχουν πιθανές παγίδες στην εφαρμογή τους. Ωστόσο, θεωρείται πιθανόν ότι η συνεχώς αυξανόμενη διαθέσιμη υπολογιστική ισχύς θα συμβάλλει μεσοπρόθεσμα στην περαιτέρω ανάπτυξη και τελική καθιέρωση της δυναμικής ανελαστικής ανάλυσης ως το ακριβέστερο και πλέον αξιόπιστο αναλυτικό εργαλείο.

ΑΝΑΦΟΡΕΣ CEN (2004a) “Eurocode 2 - Design of Concrete Structures”, European Committee for

Standardization. CEN (2004b) “Eurocode 8 - Design of Structures for Earthquake Resistance”, European

Committee for Standardization. Computers and Structures Inc. (2008) “ETABS reference manual”, Computers and

Structures, Berkeley, California. Elnashai A.S, Papanikolaou V.K and Lee D.H. (2009) “Zeus-NL - A program for

inelastic dynamic analysis of structures”, Mid-America Earthquake Center, University of Illinois at Urbana-Champaign.

Papanikolaou V.K. and Elnashai A.S. (2005) “Evaluation of conventional and adap-tive pushover analysis I : Methodology”, Journal of Earthquake Engineering, Vol. 9, No. 6., pp. 923-941.

Penelis Λογισμικά ΕΠΕ (2008) “eTools – Εγχειρίδιο Χρήσης”, Θεσσαλονίκη. Penelis Gr.G. and Kappos A.J. (2002) “3D pushover analysis: the issue of torsion”,

Proceedings of The Twelfth European Conference on Earthquake Engineer-ing, electronic resource , 10 pp.

Penelis Gr.G. and Kappos A.J. (2005) “Inelastic torsion effects in 3D pushover analy-sis of buildings”, Proceedings of the 4th European Workshop on the Seismic Behaviour of Irregular and Complex Structures, CD ROM, Thessaloniki.

Penelis, G.G., Papanikolaou, V.K. and Pashalidis, K. (2008) “Evaluation of the nonlinear behaviour of a 15-storey torsionally sensitive building designed according to EC2 and EC8”, Fifth European Workshop on the Seismic Behaviour of Irregular and Complex Structures, Catania, Italy, September 2008.

Building Seismic Safety Council (BSSC) (2000) “Prestandard and Commentary for the Seismic Rehabilitation of Buildings, FEMA-356”, Federal Emergency Management Agency, Washington, D.C.


Αποτίμηση ανελαστικής συμπεριφοράς 16- EC8 EC2 Evaluation...

Documents

Transcript of Αποτίμηση ανελαστικής συμπεριφοράς 16- EC8 EC2 Evaluation...