Probabilistic finite element analysis of structures using the multiplicative dimensional reduction method

Abstract

It is widely accepted that uncertainty may be present for many engineering problems, such as in input variables (loading, material properties, etc.), in response variables (displacements, stresses, etc.) and in the relationships between them. Reliability analysis is capable of dealing with all these uncertainties providing the engineers with accurate predictions of the probability of a structure performing adequately during its lifetime.In probabilistic finite element analysis (FEA), approximate methods such as Taylor series methods are used in order to compute the mean and the variance of the response, while the distribution of the response is usually approximated based on the Monte Carlo simulation (MCS) method. This study advances probabilistic FEA by combining it with the multiplicative form of dimensional reduction method (M-DRM). This combination allows fairly accurate estimations of both the statistical moments and the probability distribution of the response of interest. The response probability distribution is obtained using the fractional moments, which are calculated from M-DRM, together with the maximum entropy (MaxEnt) principle. In addition, the global variance-based sensitivity coefficients are also obtained as a by-product of the previous analysis. Therefore, no extra analytical work is required for sensitivity analysis. The proposed approach is integrated with the OpenSees FEA software using Tcl programing and with the ABAQUS FEA software using Python programing. OpenSees is used to analyze structures under seismic loading, where both pushover analysis and dynamic analysis is performed. ABAQUS is used to analyze structures under static loading, where the concrete damage plasticity model is used for the modeling of concrete. Thus, the efficient applicability of the proposed method is illustrated and its numerical accuracy is examined, through several examples of nonlinear FEA of structures. This research shows that the proposed method, which is based on a small number of finite element analyses, is robust, computational effective and easily applicable, providing a feasible alternative for finite element reliability and sensitivity analysis of practical and real life problems. The results of such work have significance in future studies for the estimation of the probability of the response exceeding a safety limit and for establishing safety factors related to acceptable probabilities of structural failures.

Είναι ευρέως αποδεκτό ότι αβεβαιότητα μπορεί να υπάρξει σε πολλά προβλήματα μηχανικής, όπως στις εισαγόμενες παραμέτρους (φορτίο, γεωμετρία, ιδιότητες υλικών, κλπ.), στις εξαγόμενες παραμέτρους (μετατοπίσεις, τάσεις, κλπ.) και στην μεταξύ τους συσχέτιση. Η ανάλυση αξιοπιστίας είναι σε θέση να αντιμετωπίσει όλες αυτές τις αβεβαιότητες παρέχοντας στους μηχανικούς ακριβείς προβλέψεις σχετικά με την πιθανότητα μιας κατασκευής να συμπεριφέρεται επαρκώς (και με ασφάλεια) στις εξωτερικές φορτίσεις κατά τη διάρκεια της ζωής της.Στην πιθανοτική ανάλυση με τη χρήση των πεπερασμένων στοιχείων (FEA), χρησιμοποιούνται προσεγγιστικές μέθοδοι (όπως οι μέθοδοι σειράς Taylor) προκειμένου να υπολογιστεί ο μέσος όρος και η διακύμανση της δομικής απόκρισης, ενώ η πιθανοτική κατανομή της απόκρισης συνήθως προσεγγίζεται με βάση τη μέθοδο προσομοίωσης Monte Carlo (MCS). Αυτή η διατριβή προάγει την πιθανοτική ανάλυση με τη χρήση των πεπερασμένων στοιχείων FEA, δηλαδή μέσω μη γραμμικής στατικής και δυναμικής ανάλυσης της κατασκευής με τη χρήση αριθμητικών μεθόδων ανάλυσης, συνδυάζοντάς την με την πολλαπλασιαστική μορφή της μεθόδου μείωσης των διαστάσεων (M-DRM). Αυτός ο συνδυασμός επιτρέπει αρκετά ακριβείς εκτιμήσεις τόσο των στατιστικών παραμέτρων (μέση τιμή, τυπική απόκλιση, κλπ) όσο και την εκτίμηση της πιθανοτικής κατανομής της δομικής απόκρισης (τάση, μετατόπιση, κλπ.).Η προτεινόμενη προσέγγιση ενσωματώνεται σε δύο ευρέως γνωστά λογισμικά ανάλυσης πεπερασμένων στοιχείων, δηλαδή στο λογισμικό OpenSees χρησιμοποιώντας γλώσσα προγραμματισμού Tcl και στο λογισμικό ABAQUS χρησιμοποιώντας γλώσσα προγραμματισμού Python. Εν συνεχεία, το OpenSees χρησιμοποιείται για την ανάλυση δομών υπό σεισμική φόρτιση, όπου εκτελούνται ανελαστική στατική ανάλυση (pushover) και δυναμική ανάλυση χρησιμοποιώντας επεξεργασμένα επιταχυνσιογραφήματα από προηγούμενους σεισμούς και λαμβάνοντας υπόψη και επαναλαμβανόμενους σεισμούς. Επίσης, το ABAQUS χρησιμοποιείται για την ανάλυση δομών υπό μη γραμμική στατική φόρτωση, για συνδέσεις πλακών-υποστυλώματος από οπλισμένο σκυρόδεμα και για τοιχία πυρηνικών αντιδραστήρων από προεντεταμένο σκυρόδεμα. Αυτή η έρευνα δείχνει ότι η προτεινόμενη μέθοδος, η οποία βασίζεται σε μικρό αριθμό αναλύσεων πεπερασμένων στοιχείων, είναι αποτελεσματική, υπολογιστικά εφικτή και εύκολα εφαρμόσιμη, παρέχοντας μια εφικτή εναλλακτική λύση για την ανάλυση αξιοπιστίας (και ανάλυση ευαισθησίας) με την χρήση των πεπερασμένων στοιχείων για πραγματικά κατασκευαστικά προβλήματα (κατασκευές τρισδιάστατες και μεγάλης κλίμακας). Τα αποτελέσματα μιας τέτοιας εργασίας έχουν σημασία σε μελλοντικές μελέτες για την εκτίμηση της πιθανότητας της δομικής απόκρισης να υπερβεί ένα όριο ασφαλείας και για τον καθορισμό παραγόντων ασφάλειας που σχετίζονται με αποδεκτές πιθανότητες δομικών βλαβών.