Θεωρία ελαστικότητας
Η Θεωρία ελαστικότητας είναι κλάδος της μηχανικής και των μαθηματικών που αφορά στη μοντελοποίηση της συμπεριφοράς στερεού σώματος, υπό την επίδραση εξωτερικών δράσεων (δυνάμεων και ροπών, μετατοπίσεων, θερμοκρασιακών διαφορών) όταν αυτό θεωρηθεί ως συνεχές, ομογενές, ισότροπο και ελαστικό.
Πιο συγκεκριμένα, η Θεωρία ελαστικότητας είναι κλάδος της Μηχανικής παραμορφωσίμων σωμάτων, η οποία με τη σειρά της κλάδος της Μηχανικής του συνεχούς μέσου.
Η Μηχανική του συνεχούς μέσου άλλωστε περιλαμβάνει τη Μηχανική παραμορφωσίμων σωμάτων και τη Μηχανική των ρευστών.
Στη γενική περίπτωση ένα στερεό σώμα μπορεί υπό την επίδραση δυνάμεων ή άλλων δράσεων να επιταχύνεται και/ή να παραμορφώνεται.
Ωστόσο αν το μοντέλο που επιλέξουμε είναι γραμμικό τότε τα δύο μπορούν να εξεταστούν ανεξάρτητα.
Άρα λοιπόν στη Θεωρία ελαστικότητας σπάνια εξετάζουμε μη-γραμμικά μοντέλα.
Επίσης σπάνια εξετάζουμε μοντέλα που περιλαμβάνουν δυναμική συμπεριφορά (επιταχύνσεις).
Με άλλα λόγια μια εισαγωγή στη θεωρία ελαστικότητας καταρχήν εξετάζει τις παραμορφώσεις γραμμικού ελαστικού σώματος υπό την επίδραση δυνάμεων και ροπών και εξωτερικά επιβαλλόμενων μετατοπίσεων.
Η Μηχανική παραμορφωσίμων σωμάτων έχει λοιπόν αρκετούς άλλους κλάδους όπως την Αντοχή των υλικών, τη Θεωρία πεπερασμένων στοιχείων, τη μοντελοποίηση δυναμικών αποκρίσεων, τη μοντελοποίηση παραμορφώσεων υπό την επίδραση θερμικής διαστολής/συστολής, τη θεωρία γραμμικών τέλεια πλαστικών σωμάτων και άλλους.
Ο τανυστής των τάσεων
thumb|right|260px thumb|left|280px
Οι εξισώσεις των τάσεων προκύπτουν απ'το δεύτερο νόμο του Νεύτωνα που είναι σε μορφή τανυστή:
\\ \\sigma_{ij,j}+ F_i = \\rho \\partial_{tt} u_i
Για μηδενικές επιταχύνσεις έχουμε:
\\frac{\\partial \\sigma_x}{\\partial x} + \\frac{\\partial \\tau_{yx}}{\\partial y} + \\frac{\\partial \\tau_{zx}}{\\partial z} + F_x = 0
\\frac{\\partial \\tau_{xy}}{\\partial x} + \\frac{\\partial \\sigma_y}{\\partial y} + \\frac{\\partial \\tau_{zy}}{\\partial z} + F_y = 0
\\frac{\\partial \\tau_{xz}}{\\partial x} + \\frac{\\partial \\tau_{yz}}{\\partial y} + \\frac{\\partial \\sigma_z}{\\partial z} + F_z = 0
Πρόκειται για 3 γραμμικά ανεξάρτητες εξισώσεις με 6 αγνώστους.
Η παραπάνω μορφή προκύπτει για ορθογώνιο σύστημα συντεταγμένων.
Ο τανυστής των τροπών
thumb|right|350px
Οι εξισώσεις των τροπών, δηλαδή των ανηγμένων παραμορφώσεων, προκύπτουν με χρήση διαφορικής ανάλυσης για τον στοιχειώδη κύβο όταν αυτός παραμορφωθεί.
Αφορούν δηλαδή απλά τη γεωμετρία του παραμορφωμένου κύβου όταν αυτός θεωρηθεί πολύ μικρός.
Ανηγμένη παραμόρφωση στη μία διάσταση σημαίνει παραμόρφωση διά το αρχικό μήκος (αδιάστατο δηλαδή μέγεθος).
Στις τρεις διαστάσεις τις ονομάζουμε τροπές και γενικά λαμβάνουν πιο πολύπλοκη έκφραση:
\\epsilon_x=\\frac{\\partial u_x}{\\partial x}
\\gamma_{xy}=\\frac{\\partial u_x}{\\partial y}+\\frac{\\partial u_y}{\\partial x}
\\epsilon_y=\\frac{\\partial u_y}{\\partial y}
\\gamma_{yz}=\\frac{\\partial u_y}{\\partial z}+\\frac{\\partial u_z}{\\partial y}
\\epsilon_z=\\frac{\\partial u_z}{\\partial z}
\\gamma_{zx}=\\frac{\\partial u_z}{\\partial x}+\\frac{\\partial u_x}{\\partial z}
Οι οποίες είναι 6 εξισώσεις με 9 αγνώστους.
Σε μορφή τανυστή μπορούν να γραφούν:
\\varepsilon_{ij} =\\frac{1}{2} (u_{j,i}+u_{i,j})
Οι εξισώσεις συμβιβαστότητας
Οι 6 εξισώσεις των τροπών δεν είναι ανεξάρτητες μεταξύ τους.
Απαλείφοντας τις μετατοπίσεις u_x \\,, u_y \\,, u_z \\, προκύπτουν οι 6 εξισώσεις συμβιβαστότητας:
\\frac{\\partial^2 \\epsilon_x}{\\partial y^2} + \\frac{\\partial^2 \\epsilon_y}{\\partial x^2} = \\frac{\\partial^2 \\gamma_{xy}}{\\partial x \\partial y}
\\frac{\\partial^2 \\epsilon_y}{\\partial z^2} + \\frac{\\partial^2 \\epsilon_z}{\\partial y^2} = \\frac{\\partial^2 \\gamma_{yz}}{\\partial y \\partial z}
\\frac{\\partial^2 \\epsilon_z}{\\partial x^2} + \\frac{\\partial^2 \\epsilon_x}{\\partial z^2} = \\frac{\\partial^2 \\gamma_{zx}}{\\partial x \\partial z}
2 \\frac{\\partial^2 \\epsilon_x}{\\partial y \\partial z} = \\frac{\\partial}{\\partial x} ( - \\frac{\\partial \\gamma_{yz}}{\\partial x} + \\frac{\\partial \\gamma_{xz}}{\\partial y} + \\frac{\\partial \\gamma_{xy}}{\\partial z})
2 \\frac{\\partial^2 \\epsilon_y}{\\partial x \\partial z} = \\frac{\\partial}{\\partial y} ( - \\frac{\\partial \\gamma_{xz}}{\\partial y} + \\frac{\\partial \\gamma_{xy}}{\\partial z} + \\frac{\\partial \\gamma_{yz}}{\\partial x})
2 \\frac{\\partial^2 \\epsilon_z}{\\partial x \\partial y} = \\frac{\\partial}{\\partial z} ( - \\frac{\\partial \\gamma_{xy}}{\\partial z} + \\frac{\\partial \\gamma_{yz}}{\\partial x} + \\frac{\\partial \\gamma_{xz}}{\\partial y})
Όταν ικανοποιούνται οι συνθήκες συμβιβαστότητας των παραμορφώσεων εξασφαλίζεται η ύπαρξη μονοσήμαντου και συνεχούς πεδίου τιμών για την ανηγμένη παραμόρφωση.
Γραμμική ελαστικότητα
Για τον προσδιορισμό των έξι τάσεων διαθέτουμε τρεις διαφορικές εξισώσεις.
Άρα το πρόβλημα είναι τρεις φορές αόριστο.
Για τον προσδιορισμό των έξι ανηγμένων παραμορφώσεων και των τριών μετατοπίσεων διαθέτουμε έξι διαφορικές εξισώσεις.
Πάλι μας λείπουν τρεις σχέσεις.
Συνολικά δεκαπέντε άγνωστοι και εννέα εξισώσεις.
Η αοριστία μπορεί να αρθεί αν εισάγουμε τις λεγόμενες καταστατικές εξισώσεις για το υλικό.
Στην περίπτωση γραμμικά ελαστικού, ομογενούς και ισότροπου υλικού ο νόμος του Hook στις τρεις διαστάσεις γράφεται:
\\epsilon_x = \\frac{1}{E} [\\sigma_x - \\nu(\\sigma_y + \\sigma_z)] + \\alpha \\Delta T
\\epsilon_y = \\frac{1}{E} [\\sigma_y - \\nu(\\sigma_x + \\sigma_z)] + \\alpha \\Delta T
\\epsilon_z = \\frac{1}{E} [\\sigma_z - \\nu(\\sigma_x + \\sigma_y)] + \\alpha \\Delta T
\\gamma_{xy} = \\frac{1}{G} \\tau_{xy}
\\gamma_{xz} = \\frac{1}{G} \\tau_{xz}
\\gamma_{yz} = \\frac{1}{G} \\tau_{yz}
Σύμβολα που χρησιμοποιήθηκαν
Αριθμητική επίλυση
Η μορφή που παίρνει το μοντέλο ώστε να είναι επιλύσιμο με τη μέθοδο των πεπερασμένων στοιχείων αλλάζει (μετασχηματίζεται) σε σχέση με το αρχικό μοντέλο.
Πολλές αριθμητικές μέθοδοι μπορούν να χρησιμοποιηθούν και η επιλογή μπορεί να εξαρτάται απ'τις συνοριακές συνθήκες, τη γεωμετρία του συγκεκριμένου προβλήματος που εξετάζεται, απ'το αν υπάρχει ήδη διαθέσιμος κώδικας και από άλλους τεχνικούς παράγοντες.
Ακριβής επίλυση
Ακριβής επίλυση είναι δυνατή μόνο σε πολύ συγκεκριμένες περιπτώσεις προβλημάτων, κατά κανόνα με μείωση των διαστάσεων σε μία, και συγκεκριμένες συνοριακές συνθήκες.
Ιστορικά η ακριβής επίλυση προηγήθηκε και αποτελεί τον κλάδο της Αντοχής των Υλικών.
Βιβλιογραφία - Παραπομπές
