Κεφάλαιο 1 Βασικές έννοιες

Στο παρόν κεφάλαιο παρουσιάζεται και εξετάζεται το πλαίσιο στο οποίο ορίζεται το πρόβλημα της συνοριακής τιμής στη γραμμική ελαστικότητα.

Αρχικά παρουσιάζεται ο τανυστής των τάσεων (Ενότητα 1.1): Η έννοια της τάσης (Ενότητα 1.1.1), ο τανυστικός χαρακτήρας της τάσης (Ενότητα 1.1.2), οι μετασχηματισμοί των τάσεων (Ενότητα 1.1.3), οι κύριες τάσεις (Ενότητα 1.1.4), o εκτροπικός τανυστής των τάσεων (Ενότητα 1.1.5) και ο κύκλος του Mohr (Ενότητα 1.1.6).

Ακολουθεί η παρουσίαση του τανυστή των παραμορφώσεων, όπου περιγράφονται η έννοια της παραμόρφωσης, η γεωμετρική ερμηνεία της παραμόρφωσης και οι συνθήκες συμβιβαστού (Ενότητα 1.2).

Στη συνέχεια δίνονται οι σχέσεις τάσεων-παραμορφώσεων (Ενότητα 1.3): Εξετάζεται η γενική περίπτωση του Νόμου του Hooke (Ενότητα 1.3.1), ο μετασχηματισμός του τανυστή της ελαστικότητας (Ενότητα 1.3.2) και οι ιδιαίτερες περιπτώσεις ελαστικής ανισοτροπίας που αφορούν εγκαρσίως ισότροπα ελαστικά υλικά (Ενότητα 1.3.3) και ορθότροπα ελαστικά υλικά (Ενότητα 1.3.4). Τέλος παρουσιάζονται τα ισότροπα ελαστικά υλικά (Ενότητα 1.3.5).

Οι περιπτώσεις της επίπεδης συμμετρίας (Ενότητα 1.4) και της αξονοσυμμετρίας (Ενότητα 1.5) εξετάζονται στη συνέχεια. Ειδικά για την περίπτωση της αξονοσυμμετρίας δίνεται το σύστημα συντεταγμένων (Ενότητα 1.5.1), οι μετακινήσεις, οι παραμορφώσεις και τάσεις (Ενότητα 1.5.2), οι κινηματικές εξισώσεις (Ενότητα 1.5.3), οι καταστατικές εξισώσεις (Ενότητα 1.5.4) και οι εξισώσεις ισορροπίας (Ενότητα 1.5.5).

Στο τέλος του κεφαλαίου παρατίθεται η σχετική βιβλιογραφία (Ενότητα 1.6).

1.1 Τάσεις

Ο τανυστής τάσεων του Cauchy, που χρησιμοποιείται για την περιγραφή της εντατικής κατάστασης που αναπτύσσεται σε σώματα που υπόκεινται σε μικρές παραμορφώσεις, περιγράφεται στη συνέχεια.

1.1.1 Ορισμός

Έστω ένα σώμα Ω\Omega που ορίζεται από το κλειστό και φραγμένο σύνορό του Γ\Gamma. Σε ένα τμήμα του συνόρου, που συμβολίζεται με Γu\Gamma_{u}, είναι γνωστές οι μετακινήσεις uu.

Θεωρείται ότι στο σώμα δρουν επιφανειακά φορτία 𝑭𝒕\boldsymbol{F_{t}}, σημειακά φορτία 𝑭𝒑\boldsymbol{F_{p}} και μαζικά φορτία 𝑭𝒃\boldsymbol{F_{b}}, όπως φαίνεται στο Σχήμα 1.1.

Σχήμα 1.1: Ένα στερεό σώμα με φορτία και συνοριακές συνθήκες.
Σχήμα 1.2: Τομή του σώματος κατά μήκος ενός επιπέδου.

Έστω τώρα ένα τυχαίο επίπεδο PP που χωρίζει το σώμα Ω\Omega σε δύο τμήματα, Ω+\Omega^{+} και Ω-\Omega^{-}. Θεωρώντας ένα καρτεσιανό σύστημα συντεταγμένων, η δύναμη FF που ασκείται στην περιοχή ΔA\Delta A του επιπέδου PP μπορεί να οριστεί ως [2]:

𝑭=Fx𝒊+Fy𝒋+Fz𝒌\boldsymbol{F}=F_{x}\boldsymbol{i}+F_{y}\boldsymbol{j}+F_{z}\boldsymbol{k} (1.1)

όπου 𝒊\boldsymbol{i}, 𝒋\boldsymbol{j} και 𝒌\boldsymbol{k} τα μοναδιαία διανύσματα του συστήματος αναφοράς, όπως αυτά περιγράφονται στο Σχήμα 1.2).

Η δύναμη FF μπορεί επίσης να αναλυθεί ως προς ένα διάνυσμα 𝒏\boldsymbol{n} κάθετο στο επίπεδο PP και σε ένα διάνυσμα 𝒔\boldsymbol{s} παράλληλο ως προς αυτό, όπως φαίνεται στο Σχήμα 1.3. Σε αυτή την περίπτωση οι συνιστώσες 𝑭𝒏\boldsymbol{F_{n}} και 𝑭𝒔\boldsymbol{F_{s}} ονομάζονται ορθή και διατμητική δύναμη στην περιοχή ΔA\Delta A αντίστοιχα.

Σχήμα 1.3: Ορθή και διατμητική συνιστώσα της φόρτισης FF.

Το διάνυσμα της τάσης σε ένα σημείο ορίζεται θεωρώντας ότι το ΔA\Delta A γίνεται απειροστό,

𝝈(𝒏):=limΔA0𝑭ΔA,\boldsymbol{\sigma}^{(\boldsymbol{n})}:=\lim_{\Delta A\rightarrow 0}\dfrac{% \boldsymbol{F}}{\Delta A}, (1.2)

ή σε μορφή δεικτών

σi(𝒏):=limΔA0FiΔA.\sigma_{i}^{(\boldsymbol{n})}:=\lim_{\Delta A\rightarrow 0}\dfrac{F_{i}}{% \Delta A}. (1.3)

Αντίστοιχα, η ορθή τάση 𝝈n\boldsymbol{\sigma}_{n} και η διατμητική τάση 𝝉\boldsymbol{\tau} ορίζονται αντίστοιχα ως:

𝝈n\displaystyle\boldsymbol{\sigma}_{n} =limΔA0𝑭nΔA\displaystyle=\lim_{\Delta A\rightarrow 0}\dfrac{\boldsymbol{F}_{n}}{\Delta A} (1.4)
και
𝝉\displaystyle\boldsymbol{\tau} =limΔA0𝑭sΔA.\displaystyle=\lim_{\Delta A\rightarrow 0}\dfrac{\boldsymbol{F}_{s}}{\Delta A}. (1.5)

Σύμφωνα με τη σύμβαση προσήμου που ακολουθείται στην κλασική μηχανική, ως θετική ορθή τάση θεωρείται η εφελκυστική, ενώ η θετική διατμητική τάση έχει θετική φορά στο επίπεδο του οποίου το κάθετο διάνυσμα με κατεύθυνση προς τα έξω είναι επίσης θετικό, όπως αυτές δίνονται στο Σχήμα 1.4.

Σχήμα 1.4: Σύμβαση προσήμου που χρησιμοποιείται στην κλασική μηχανική.

Στην κλασική γεωτεχνική μηχανική ακολουθείται συνήθως η αντίθετη σύμβαση προσήμου, όπως σημειώνεται στο Σχήμα 1.5. Έτσι, οι αρχικές τάσεις, οι οποίες είναι συνήθως θλιπτικές, ορίζονται ως θετικές.

Σχήμα 1.5: Σύμβαση προσήμου που χρησιμοποιείται στη γεωτεχνική μηχανική.

Στο παρόν βιβλίο ακολουθείται η σύμβαση προσήμου της κλασικής μηχανικής για λόγους συμβατότητας με τη βιβλιογραφία και τις τυπικές εφαρμογές της μεθόδου των πεπερασμένων στοιχείων.

1.1.2 Ο τανυστής των τάσεων του Cauchy

Σύμφωνα με τη θεώρηση του Cauchy, οι τάσεις 𝝈(𝐧)\boldsymbol{\sigma}^{(\mathbf{n})} είναι ίδιες για όλες τις επιφάνειες που περνούν από το σημείο PP και έχουν το ίδιο κάθετο μοναδιαίο διάνυσμα 𝒏\boldsymbol{n} στο PP, έχουν δηλαδή κοινή εφαπτομένη στο σημείο PP. Πρακτικά αυτό σημαίνει ότι ο τανυστής των τάσεων είναι συνάρτηση μόνο του κάθετου μοναδιαίου διανύσματος 𝒏\boldsymbol{n} σε ένα σημείο.

Συνέπεια της παραπάνω θεώρησης αποτελεί το γνωστό και ως θεμελιώδες λήμμα του Cauchy, που αναφέρει ότι τα διανύσματα των τάσεων που δρουν στις αντίθετες πλευρές της ίδιας επιφάνειας είναι ίσα όσον αφορά το μέγεθος και αντίθετα όσον αφορά την κατεύθυνση.

Το θεμελιώδες λήμμα του Cauchy που εκφράζεται ως

𝝈(𝒏)=-𝝈(-𝒏)\boldsymbol{\sigma}^{(\boldsymbol{n})}=-\boldsymbol{\sigma}^{(-\boldsymbol{n})} (1.6)

μπορεί να θεωρηθεί και ως μια ισοδύναμη διατύπωση του Τρίτου Νόμου του Νεύτωνα που αναφέρεται στη δράση και την αντίδραση.

Σχήμα 1.6: Τάσεις στον τρισδιάστατο χώρο.

Το θεώρημα των τάσεων του Cauchy ορίζει τον τανυστή των τάσεων του Cauchy σε μια επιφάνεια κάθετη στο 𝒏\boldsymbol{n} ως:

𝝈(𝒏)=𝒏𝝈\boldsymbol{\sigma}^{(\boldsymbol{n})}=\boldsymbol{n}\cdot\boldsymbol{\sigma} (1.7)

ή σε μορφή δεικτών

σin=σijnj.\sigma_{i}^{n}=\sigma_{ij}n_{j}. (1.8)

Θεωρούμε τώρα ένα στοιχειώδες σώμα με επίπεδα κάθετα στο καρτεσιανό σύστημα συντεταγμένων, όπως αυτό περιγράφεται στο Σχήμα 1.6. Τα διανύσματα τάσεων για κάθε επίπεδο (𝝈(𝐞1)\boldsymbol{\sigma}^{(\mathbf{e}_{1})}, 𝝈(𝐞2)\boldsymbol{\sigma}^{(\mathbf{e}_{2})}, 𝝈(𝐞3)\boldsymbol{\sigma}^{(\mathbf{e}_{3})}) μπορούν να αναλυθούν σε μία ορθή και δύο διατμητικές τάσεις κατά τις διευθύνσεις των αξόνων του συστήματος συντεταγμένων. Στην περίπτωση ενός επιπέδου με κάθετο διάνυσμα στην κατεύθυνση του x1x_{1} άξονα, η ορθή τάση ορίζεται ως σ11\sigma_{11} και οι δύο διατμητικές ως σ12\sigma_{12} και σ13\sigma_{13} αντίστοιχα.

Με βάση τα παραπάνω αποδεικνύεται εύκολα ότι:

𝝈(𝐞1)\displaystyle\boldsymbol{\sigma}^{(\mathbf{e}_{1})} =σ1(𝒆1)𝒆1+σ2(𝒆1)𝒆2+σ3(𝒆1)𝒆3=σ11𝒆1+σ12𝒆2+σ13𝒆3\displaystyle=\sigma_{1}^{(\boldsymbol{e}_{1})}\boldsymbol{e}_{1}+\sigma_{2}^{% (\boldsymbol{e}_{1})}\boldsymbol{e}_{2}+\sigma_{3}^{(\boldsymbol{e}_{1})}% \boldsymbol{e}_{3}=\sigma_{11}\boldsymbol{e}_{1}+\sigma_{12}\boldsymbol{e}_{2}% +\sigma_{13}\boldsymbol{e}_{3} (1.9)
𝝈(𝐞2)\displaystyle\boldsymbol{\sigma}^{(\mathbf{e}_{2})} =σ1(𝒆2)𝒆1+σ2(𝒆2)𝒆2+σ3(𝒆2)𝒆3=σ21𝒆1+σ22𝒆2+σ23𝒆3\displaystyle=\sigma_{1}^{(\boldsymbol{e}_{2})}\boldsymbol{e}_{1}+\sigma_{2}^{% (\boldsymbol{e}_{2})}\boldsymbol{e}_{2}+\sigma_{3}^{(\boldsymbol{e}_{2})}% \boldsymbol{e}_{3}=\sigma_{21}\boldsymbol{e}_{1}+\sigma_{22}\boldsymbol{e}_{2}% +\sigma_{23}\boldsymbol{e}_{3} (1.10)
𝝈(𝐞3)\displaystyle\boldsymbol{\sigma}^{(\mathbf{e}_{3})} =σ1(𝒆3)𝒆1+σ2(𝒆3)𝒆2+σ3(𝒆3)𝒆3=σ31𝒆1+σ32𝒆2+σ33𝒆3\displaystyle=\sigma_{1}^{(\boldsymbol{e}_{3})}\boldsymbol{e}_{1}+\sigma_{2}^{% (\boldsymbol{e}_{3})}\boldsymbol{e}_{2}+\sigma_{3}^{(\boldsymbol{e}_{3})}% \boldsymbol{e}_{3}=\sigma_{31}\boldsymbol{e}_{1}+\sigma_{32}\boldsymbol{e}_{2}% +\sigma_{33}\boldsymbol{e}_{3} (1.11)

ή αλλιώς

𝝈(𝒆i)=σj(𝒆i)𝒆j=σij𝒆j\boldsymbol{\sigma}^{(\boldsymbol{e}_{i})}=\sigma_{j}^{(\boldsymbol{e}_{i})}% \boldsymbol{e}_{j}=\sigma_{ij}\boldsymbol{e}_{j} (1.12)

Οι εννέα συνιστώσες σij\sigma_{ij} των διανυσμάτων τάσεων αποτελούν τις συνιστώσες του τανυστή τάσεων του Cauchy, με τον οποίο περιγράφεται πλήρως η εντατική κατάσταση σε ένα σημείο ως:

𝝈=σij=[𝝈(𝐞1)𝝈(𝐞1)𝝈(𝐞1)]=[σ11σ12σ13σ21σ22σ23σ31σ32σ33]\boldsymbol{\sigma}=\sigma_{ij}=\left[\begin{matrix}\boldsymbol{\sigma}^{(% \mathbf{e}_{1})}\\ \boldsymbol{\sigma}^{(\mathbf{e}_{1})}\\ \boldsymbol{\sigma}^{(\mathbf{e}_{1})}\\ \end{matrix}\right]=\left[\begin{matrix}\sigma_{11}&\sigma_{12}&\sigma_{13}\\ \sigma_{21}&\sigma_{22}&\sigma_{23}\\ \sigma_{31}&\sigma_{32}&\sigma_{33}\\ \end{matrix}\right] (1.13)

ή ισοδύναμα στο σύστημα xx-yy-zz

𝝈=σij=[𝝈(𝐞1)𝝈(𝐞1)𝝈(𝐞1)]=[σxxσxyσxzσyxσyyσyzσzxσzyσzz]\boldsymbol{\sigma}=\sigma_{ij}=\left[\begin{matrix}\boldsymbol{\sigma}^{(% \mathbf{e}_{1})}\\ \boldsymbol{\sigma}^{(\mathbf{e}_{1})}\\ \boldsymbol{\sigma}^{(\mathbf{e}_{1})}\\ \end{matrix}\right]=\left[\begin{matrix}\sigma_{xx}&\sigma_{xy}&\sigma_{xz}\\ \sigma_{yx}&\sigma_{yy}&\sigma_{yz}\\ \sigma_{zx}&\sigma_{zy}&\sigma_{zz}\\ \end{matrix}\right] (1.14)

όπου σ11\sigma_{11} (σxx\sigma_{xx}), σ22\sigma_{22} (σyy\sigma_{yy}) και σ33\sigma_{33} (σzz\sigma_{zz}) οι ορθές τάσεις και σ12\sigma_{12} (σxy\sigma_{xy}), σ13\sigma_{13} (σxz\sigma_{xz}), σ21\sigma_{21} (σyx\sigma_{yx}), σ23\sigma_{23} (σyz\sigma_{yz}), σ31\sigma_{31} (σzx\sigma_{zx}) και σ32\sigma_{32} (σzy\sigma_{zy}) οι διατμητικές τάσεις αντίστοιχα.

Λαμβάνοντας υπόψη τη συμμετρία του τανυστή των τάσεων, μπορούμε να τον απεικονίσουμε ως διάνυσμα της μορφής

𝝈=[σ1σ2σ3σ4σ5σ6]T[σ11σ22σ33σ23σ13σ12]T\boldsymbol{\sigma}=\begin{bmatrix}\sigma_{1}&\sigma_{2}&\sigma_{3}&\sigma_{4}% &\sigma_{5}&\sigma_{6}\end{bmatrix}^{T}\equiv\begin{bmatrix}\sigma_{11}&\sigma% _{22}&\sigma_{33}&\sigma_{23}&\sigma_{13}&\sigma_{12}\end{bmatrix}^{T} (1.15)

η οποία ονομάζεται σύμβαση Voigt. Πιο συνηθισμένη είναι η ακόλουθη διατύπωση:

𝝈=[σ1σ2σ3σ4σ5σ6]T[σ11σ22σ33σ12σ21σ31]T\boldsymbol{\sigma}=\begin{bmatrix}\sigma_{1}&\sigma_{2}&\sigma_{3}&\sigma_{4}% &\sigma_{5}&\sigma_{6}\end{bmatrix}^{T}\equiv\begin{bmatrix}\sigma_{11}&\sigma% _{22}&\sigma_{33}&\sigma_{12}&\sigma_{21}&\sigma_{31}\end{bmatrix}^{T} (1.16)

η οποία και θα χρησιμοποιηθεί στη συνέχεια.

Για την περίπτωση δύο διαστάσεων ισχύει

𝝈=[σxxτxyτyxσyy]\boldsymbol{\sigma}=\left[\begin{array}[]{cc}\sigma_{xx}&\tau_{xy}\\ \tau_{yx}&\sigma_{yy}\\ \end{array}\right] (1.17)

ενώ τέλος για την περίπτωση ενός προβλήματος μιας διάστασης γράφουμε απλά

σ=σxx.\sigma=\sigma_{xx}. (1.18)

1.1.3 Μετασχηματισμός των τάσεων

Ο μετασχηματισμός του τανυστή των τάσεων ακολουθεί τους κανόνες μετασχηματισμού των τανυστών δευτέρου βαθμού. Επομένως για τον μετασχηματισμό από ένα σύστημα συντεταγμένων xix_{i} σε ένα σύστημα xix_{i}^{\prime}, οι συνιστώσες σij\sigma_{ij} του αρχικού συστήματος μετασχηματίζονται στις συνιστώσες σij\sigma^{\prime}_{ij} του τελικού συστήματος σύμφωνα με την ακόλουθη σχέση:

σij=aimajnσmn\sigma^{\prime}_{ij}=a_{im}a_{jn}\sigma_{mn} (1.19)

ή σε μητρωική μορφή:

𝝈=𝐀𝝈𝐀T,\boldsymbol{\sigma}^{\prime}=\mathbf{A}\boldsymbol{\sigma}\mathbf{A}^{T}, (1.20)

όπου 𝐀\mathbf{A} είναι ένα μητρώο μετασχηματισμού με συνιστώσες aija_{ij}.

Αναλυτικά, η παραπάνω σχέση δίνεται ως:

[σ11σ12σ13σ21σ22σ23σ31σ32σ33]=[a11a12a13a21a22a23a31a32a33][σ11σ12σ13σ21σ22σ23σ31σ32σ33][a11a21a31a12a22a32a13a23a33].\left[{\begin{matrix}\sigma^{\prime}_{11}&\sigma^{\prime}_{12}&\sigma^{\prime}% _{13}\\ \sigma^{\prime}_{21}&\sigma^{\prime}_{22}&\sigma^{\prime}_{23}\\ \sigma^{\prime}_{31}&\sigma^{\prime}_{32}&\sigma^{\prime}_{33}\\ \end{matrix}}\right]=\left[{\begin{matrix}a_{11}&a_{12}&a_{13}\\ a_{21}&a_{22}&a_{23}\\ a_{31}&a_{32}&a_{33}\\ \end{matrix}}\right]\left[{\begin{matrix}\sigma_{11}&\sigma_{12}&\sigma_{13}\\ \sigma_{21}&\sigma_{22}&\sigma_{23}\\ \sigma_{31}&\sigma_{32}&\sigma_{33}\\ \end{matrix}}\right]\left[{\begin{matrix}a_{11}&a_{21}&a_{31}\\ a_{12}&a_{22}&a_{32}\\ a_{13}&a_{23}&a_{33}\\ \end{matrix}}\right]. (1.21)

1.1.4 Κύριες τάσεις

Κύριες τάσεις ονομάζονται οι ορθές τάσεις οι οποίες δρουν σε επίπεδα στα οποία οι αντίστοιχες διατμητικές τάσεις είναι μηδενικές. Οι τιμές των κυρίων τάσεων είναι ανεξάρτητες του συστήματος συντεταγμένων αναφοράς, οι διευθύνσεις τους όμως εξαρτώνται από αυτό. Σε δισδιάστατα προβλήματα, οι κύριες τάσεις αποτελούν τη συνολικά ελάχιστη και τη μέγιστη κατά μέγεθος τάση που δρουν σε κάθε επίπεδο.

Στη γενική περίπτωση οι κύριες τάσεις προκύπτουν από την επίλυση του ιδιοπροβλήματος [9]:

|σij-λδij|=|σ11-λσ12σ13σ21σ22-λσ23σ31σ32σ33-λ|=0\left|\sigma_{ij}-\lambda\delta_{ij}\right|=\begin{vmatrix}\sigma_{11}-\lambda% &\sigma_{12}&\sigma_{13}\\ \sigma_{21}&\sigma_{22}-\lambda&\sigma_{23}\\ \sigma_{31}&\sigma_{32}&\sigma_{33}-\lambda\\ \end{vmatrix}=0 (1.22)

Από την παραπάνω ορίζουσα προκύπτει η χαρακτηριστική εξίσωση:

|σij-λδij|=-λ3+I1λ2-I2λ+I3=0\left|\sigma_{ij}-\lambda\delta_{ij}\right|=-\lambda^{3}+I_{1}\lambda^{2}-I_{2% }\lambda+I_{3}=0 (1.23)

όπου:

I1\displaystyle I_{1} =σ11+σ22+σ33\displaystyle=\sigma_{11}+\sigma_{22}+\sigma_{33} (1.24)
=σkk\displaystyle=\sigma_{kk} (1.25)
=tr(𝝈)\displaystyle=\text{tr}(\boldsymbol{\sigma}) (1.26)
I2\displaystyle I_{2} =|σ22σ23σ32σ33|+|σ11σ13σ31σ33|+|σ11σ12σ21σ22|\displaystyle=\begin{vmatrix}\sigma_{22}&\sigma_{23}\\ \sigma_{32}&\sigma_{33}\\ \end{vmatrix}+\begin{vmatrix}\sigma_{11}&\sigma_{13}\\ \sigma_{31}&\sigma_{33}\\ \end{vmatrix}+\begin{vmatrix}\sigma_{11}&\sigma_{12}\\ \sigma_{21}&\sigma_{22}\\ \end{vmatrix} (1.27)
=σ11σ22+σ22σ33+σ11σ33-σ122-σ232-σ312\displaystyle=\sigma_{11}\sigma_{22}+\sigma_{22}\sigma_{33}+\sigma_{11}\sigma_% {33}-\sigma_{12}^{2}-\sigma_{23}^{2}-\sigma_{31}^{2} (1.28)
=12(σiiσjj-σijσji)\displaystyle=\dfrac{1}{2}\left(\sigma_{ii}\sigma_{jj}-\sigma_{ij}\sigma_{ji}\right) (1.29)
=12[tr(𝝈)2-tr(𝝈2)]\displaystyle=\dfrac{1}{2}\left[\text{tr}(\boldsymbol{\sigma})^{2}-\text{tr}(% \boldsymbol{\sigma}^{2})\right] (1.30)
I3\displaystyle I_{3} =det(σij)\displaystyle=\det(\sigma_{ij}) (1.31)
=det(𝝈)\displaystyle=\det(\boldsymbol{\sigma}) (1.32)
=σ11σ22σ33+2σ12σ23σ31-σ122σ33-σ232σ11-σ312σ22\displaystyle=\sigma_{11}\sigma_{22}\sigma_{33}+2\sigma_{12}\sigma_{23}\sigma_% {31}-\sigma_{12}^{2}\sigma_{33}-\sigma_{23}^{2}\sigma_{11}-\sigma_{31}^{2}% \sigma_{22} (1.33)

Επειδή ο τανυστής των τάσεων είναι συμμετρικός, η χαρακτηριστική εξίσωση έχει τρεις πραγματικές λύσεις λi\lambda_{i}. Οι κύριες τάσεις ορίζονται ως:

σ1\displaystyle\sigma_{1} =max(λ1,λ2,λ3),\displaystyle=\max\left(\lambda_{1},\lambda_{2},\lambda_{3}\right), (1.35)
σ3\displaystyle\sigma_{3} =min(λ1,λ2,λ3)\displaystyle=\min\left(\lambda_{1},\lambda_{2},\lambda_{3}\right) (1.36)
και
σ2\displaystyle\sigma_{2} =I1-σ1-σ3\displaystyle=I_{1}-\sigma_{1}-\sigma_{3} (1.37)

Από την αναλυτική επίλυση του προβλήματος προκύπτει για τη γενική, τρισδιάστατη περίπτωση:

σ1\displaystyle\sigma_{1} =I13+23(I12-3I2)cosϕ\displaystyle=\dfrac{I_{1}}{3}+\dfrac{2}{3}\left(\sqrt{I^{2}_{1}-3I_{2}}\right% )\cos\phi (1.38)
σ2\displaystyle\sigma_{2} =I13+23(I12-3I2)cos(ϕ+2π3)\displaystyle=\dfrac{I_{1}}{3}+\dfrac{2}{3}\left(\sqrt{I^{2}_{1}-3I_{2}}\right% )\cos\left(\phi+\dfrac{2\pi}{3}\right) (1.39)
σ3\displaystyle\sigma_{3} =I13+23(I12-3I2)cos(ϕ+4π3)\displaystyle=\dfrac{I_{1}}{3}+\dfrac{2}{3}\left(\sqrt{I^{2}_{1}-3I_{2}}\right% )\cos\left(\phi+\dfrac{4\pi}{3}\right) (1.40)

όπου

ϕ=13cos-1(2I13-9I1I2+27I32(I12-3I2)3/2)\phi=\cfrac{1}{3}\cos^{-1}\left(\frac{2I_{1}^{3}-9I_{1}I_{2}+27I_{3}}{2(I_{1}^% {2}-3I_{2})^{3/2}}\right) (1.41)

Αντίστοιχα για την περίπτωση δισδιάστατου προβλήματος, προκύπτει:

σ1\displaystyle\sigma_{1} =σ11+σ222+(σ11-σ222)2+σ122\displaystyle=\frac{\sigma_{11}+\sigma_{22}}{2}+\sqrt{\left(\frac{\sigma_{11}-% \sigma_{22}}{2}\right)^{2}+\sigma_{12}^{2}} (1.42)
σ2\displaystyle\sigma_{2} =σ11+σ222-(σ11-σ222)2+σ122\displaystyle=\frac{\sigma_{11}+\sigma_{22}}{2}-\sqrt{\left(\frac{\sigma_{11}-% \sigma_{22}}{2}\right)^{2}+\sigma_{12}^{2}} (1.43)

1.1.5 Ο εκτροπικός τανυστής των τάσεων

Η χαρακτηριστική εξίσωση για τον εκτροπικό τανυστή των τάσεων αντίστοιχα είναι:

|sij-λδij|=λ3-J1λ2-J2λ-J3=0\left|s_{ij}-\lambda\delta_{ij}\right|=\lambda^{3}-J_{1}\lambda^{2}-J_{2}% \lambda-J_{3}=0 (1.44)

όπου J1J_{1}, J2J_{2} και J3J_{3} είναι η πρώτη, η δεύτερη και η τρίτη αναλλοίωτη αντίστοιχα.

J1\displaystyle J_{1} =skk=0,\displaystyle=s_{kk}=0, (1.45)
J2\displaystyle J_{2} =12sijsji\displaystyle=\textstyle{\dfrac{1}{2}}s_{ij}s_{ji} (1.46)
=12tr(𝒔2)\displaystyle=\dfrac{1}{2}\text{tr}(\boldsymbol{s}^{2}) (1.47)
=12(s12+s22+s32)\displaystyle=\dfrac{1}{2}(s_{1}^{2}+s_{2}^{2}+s_{3}^{2}) (1.48)
=16[(σ11-σ22)2+(σ22-σ33)2+(σ33-σ11)2]+σ122+σ232+σ312\displaystyle=\dfrac{1}{6}\left[(\sigma_{11}-\sigma_{22})^{2}+(\sigma_{22}-% \sigma_{33})^{2}+(\sigma_{33}-\sigma_{11})^{2}\right]+\sigma_{12}^{2}+\sigma_{% 23}^{2}+\sigma_{31}^{2} (1.49)
=16[(σ1-σ2)2+(σ2-σ3)2+(σ3-σ1)2]\displaystyle=\dfrac{1}{6}\left[(\sigma_{1}-\sigma_{2})^{2}+(\sigma_{2}-\sigma% _{3})^{2}+(\sigma_{3}-\sigma_{1})^{2}\right] (1.50)
=13I12-I2\displaystyle=\dfrac{1}{3}I_{1}^{2}-I_{2} (1.51)
=12[tr(𝝈2)-13tr(𝝈)2],\displaystyle=\dfrac{1}{2}\left[\text{tr}(\boldsymbol{\sigma}^{2})-\frac{1}{3}% \text{tr}(\boldsymbol{\sigma})^{2}\right], (1.52)
J3\displaystyle J_{3} =det(sij)\displaystyle=\det(s_{ij}) (1.53)
=13sijsjkski\displaystyle=\dfrac{1}{3}s_{ij}s_{jk}s_{ki} (1.54)
=13tr(𝒔3)\displaystyle=\dfrac{1}{3}\text{tr}(\boldsymbol{s}^{3}) (1.55)
=s1s2s3\displaystyle=s_{1}s_{2}s_{3} (1.56)
=227I13-13I1I2+I3\displaystyle=\dfrac{2}{27}I_{1}^{3}-\dfrac{1}{3}I_{1}I_{2}+I_{3} (1.57)
=13[tr(𝝈3)-tr(𝝈2)tr(𝝈)+29tr(𝝈)3].\displaystyle=\dfrac{1}{3}\left[\text{tr}(\boldsymbol{\sigma}^{3})-\text{tr}(% \boldsymbol{\sigma}^{2})\text{tr}(\boldsymbol{\sigma})+\dfrac{2}{9}\text{tr}(% \boldsymbol{\sigma})^{3}\right]. (1.58)

1.1.6 Ο κύκλος του Mohr

Ο κύκλος του Mohr, όπως ονομάστηκε προς τιμήν του Christian Otto Mohr, αποτελεί μια δισδιάστατη γραφική απεικόνιση του μετασχηματισμού του τανυστή των τάσεων του Cauchy [5] (βλ. π.χ. τον κύκλο του Mohr για μια τρισδιάστατη εντατική κατάσταση στο Σχήμα 1.7).

Σχήμα 1.7: Ο κύκλος του Mohr για μια τρισδιάστατη εντατική κατάσταση.

Από την ισορροπία των δυνάμεων, η ορθή τάση σn\sigma_{\mathrm{n}} και η αντίστοιχη διατμητική τn\tau_{\mathrm{n}} προκύπτουν ως:

σn\displaystyle\sigma_{n} =12(σx+σy)+12(σx-σy)cos2θ+τxysin2θ\displaystyle=\dfrac{1}{2}(\sigma_{x}+\sigma_{y})+\dfrac{1}{2}(\sigma_{x}-% \sigma_{y})\cos 2\theta+\tau_{xy}\sin 2\theta (1.59)
τn\displaystyle\tau_{n} =-12(σx-σy)sin2θ+τxycos2θ\displaystyle=-\frac{1}{2}(\sigma_{x}-\sigma_{y})\sin 2\theta+\tau_{xy}\cos 2\theta (1.60)

1.2 Παραμορφώσεις

Όταν σε ένα σώμα ασκείται μία δύναμη, τότε το σώμα παραμορφώνεται ή μετακινείται. Αν το σώμα δεν είναι επαρκώς περιορισμένο, τότε οι εκτός ισορροπίας δυνάμεις το αναγκάζουν σε περιστροφή ή κίνηση στερεού σώματος. Αν στο σώμα υπάρχουν περιορισμοί στην κίνησή του, τότε το σώμα παραμορφώνεται υπό την επήρεια των δυνάμεων που ασκούνται σε αυτό [11].

Έστω για παράδειγμα ένα εδαφικό δείγμα υπό καθεστώς μονοαξονικής θλίψης. Στην περίπτωση αυτή ένα φορτίο ασκείται στη μία άκρη του δείγματος ενώ στην άλλη δεν επιτρέπονται οι μετακινήσεις. Αν ΔL\Delta L είναι ένα τμήμα του δοκιμίου και Δl\Delta l η αλλαγή του τμήματος, τότε η απειροστική ανηγμένη παραμόρφωση ορίζεται ως:

𝜺=limΔA0δlΔL.\boldsymbol{\varepsilon}=\lim_{\Delta A\rightarrow 0}\dfrac{\delta l}{\Delta L}. (1.61)

Ο παραπάνω ορισμός αντιστοιχεί στην ορθή παραμόρφωση και αναφέρεται στην περίπτωση που αυτή συμβαίνει κατά μιας ευθείας γραμμής. Αντίστοιχα, διατμητική παραμόρφωση παρατηρείται όταν η γωνία μεταξύ δύο αρχικά ευθειών γραμμών αλλάζει.

Σχήμα 1.8: Στοιχειώδες ορθογώνιο ενός επίπεδου σώματος πριν από και μετά την παραμόρφωση.

Στο Σχήμα 1.8 δίνεται ένα στοιχειώδες ορθογώνιο ενός επίπεδου σώματος πριν από και μετά την παραμόρφωση στο οποίο οι μετατοπίσεις uu και vv θεωρούνται συναρτήσεις των συντεταγμένων xx και yy, ισχύει δηλαδή:

u\displaystyle u :=u(x,y)\displaystyle:=u(x,y) (1.62)
v\displaystyle v :=v(x,y)\displaystyle:=v(x,y) (1.63)

Με την παραδοχή των μικρών παραμορφώσεων, η ανοιγμένη παραμόρφωση κατά τη διεύθυνση xx δίνεται από τη σχέση

εx=[dx+(u+uxx)-u]-dxdx=dudx\varepsilon_{x}=\dfrac{\left[dx+\left(u+\dfrac{\partial u}{\partial x}x\right)% -u\right]-dx}{dx}=\dfrac{du}{dx} (1.64)

και αντίστοιχα κατά τη διεύθυνση yy:

εy=dudy\varepsilon_{y}=\dfrac{du}{dy} (1.65)

Η διατμητική ανηγμένη παραμόρφωση ορίζεται ως η αλλαγή της ορθής γωνίας του στοιχειώδους παραλληλογράμμου dxdx-dydy. Με την παραδοχή των μικρών παραμορφώσεων (α1tan(α1)\alpha_{1}\approx\tan(\alpha_{1}), α2tan(α2)\alpha_{2}\approx\tan(\alpha_{2})), προκύπτει:

γxy=α1+α2=(u+uydy-u)dy+(v+vxdy-v)dx=uy+vx\gamma_{xy}=\alpha_{1}+\alpha_{2}=\dfrac{\left(u+\dfrac{\partial u}{\partial y% }dy-u\right)}{dy}+\dfrac{\left(v+\dfrac{\partial v}{\partial x}dy-v\right)}{dx% }=\dfrac{\partial u}{\partial y}+\dfrac{\partial v}{\partial x} (1.66)

Στη γενική περίπτωση ισχύει:

εx\displaystyle\varepsilon_{x} =dudx\displaystyle=\dfrac{du}{dx} (1.67)
εy\displaystyle\varepsilon_{y} =dudy\displaystyle=\dfrac{du}{dy} (1.68)
εz\displaystyle\varepsilon_{z} =dudz\displaystyle=\dfrac{du}{dz} (1.69)
γxy\displaystyle\gamma_{xy} =uy+vx\displaystyle=\dfrac{\partial u}{\partial y}+\dfrac{\partial v}{\partial x} (1.70)
γyz\displaystyle\gamma_{yz} =vz+wy\displaystyle=\dfrac{\partial v}{\partial z}+\dfrac{\partial w}{\partial y} (1.71)
γzx\displaystyle\gamma_{zx} =wx+uz\displaystyle=\dfrac{\partial w}{\partial x}+\dfrac{\partial u}{\partial z} (1.72)

Οι παραπάνω σχέσεις μπορούν να γραφτούν σε μητρωική μορφή ως:

[εxεyεzγxyγyzγzx]=[x000y00zyx00zxz0y][uvw]\left[\begin{array}[]{c}\varepsilon_{x}\\ \varepsilon_{y}\\ \varepsilon_{z}\\ \gamma_{xy}\\ \gamma_{yz}\\ \gamma_{zx}\end{array}\right]=\left[\begin{array}[]{ccc}\dfrac{\partial}{% \partial x}&0&0\\ 0&\dfrac{\partial}{\partial y}&0\\ 0&&\dfrac{\partial}{\partial z}\\ \dfrac{\partial}{\partial y}&\dfrac{\partial}{\partial x}&0\\ 0&\dfrac{\partial}{\partial z}&\dfrac{\partial}{\partial x}\\ \dfrac{\partial}{\partial z}&0&\dfrac{\partial}{\partial y}\\ \end{array}\right]\left[\begin{array}[]{c}u\\ v\\ w\end{array}\right] (1.73)

ή συνεπτυγμένα

𝜺=ε𝒖\boldsymbol{\varepsilon}=\boldsymbol{\partial}_{\varepsilon}\boldsymbol{u} (1.74)

1.3 Σχέσεις τάσεων-ανηγμένων παραμορφώσεων

Οι τάσεις και οι ανηγμένες παραμορφώσεις συνδέονται μεταξύ τους με τις σχέσεις τάσεων-ανηγμένων παραμορφώσεων.

1.3.1 Περιπτώσεις ανισοτροπίας

Στη γενική του μορφή ο Νόμος του Hooke [10] δίνεται ως:

σij=Cijklεkl\sigma_{ij}=C_{ijkl}\varepsilon_{kl} (1.75)

όπου σij\sigma_{ij} και εkl\varepsilon_{kl} είναι τανυστές δεύτερης τάξης και CijklC_{ijkl} ένας τανυστής τέταρτης τάξης.

Λόγω της συμμετρίας του τανυστή των τάσεων του Cauchy,

σij=σji\sigma_{ij}=\sigma_{ji} (1.76)

μπορεί να αποδειχτεί ότι:

Cijkl=Cjikl.C_{ijkl}=C_{jikl}. (1.77)

Παρομοίως, λόγω της συμμετρίας του τανυστή των ανηγμένων παραμορφώσεων,

εkl=εlk\varepsilon_{kl}=\varepsilon_{lk} (1.78)

ισχύει:

Cijkl=Cijlk.C_{ijkl}=C_{ijlk}. (1.79)

Επομένως ο αριθμός των σταθερών που περιγράφουν των τανυστή CijklC_{ijkl} μειώνεται από 81 σε 36. Σε μητρωική μορφή (σύμβαση Voigt) ο Νόμος του Hooke πλέον μπορεί να γραφτεί ως:

𝐂=[C11C12C13C14C15C16C21C22C23C24C25C26C31C32C33C34C35C36C41C42C43C44C45C46C51C52C53C54C55C56C61C62C63C64C65C66]\mathbf{C}=\left[\begin{array}[]{cccccc}C_{11}&C_{12}&C_{13}&C_{14}&C_{15}&C_{% 16}\\ C_{21}&C_{22}&C_{23}&C_{24}&C_{25}&C_{26}\\ C_{31}&C_{32}&C_{33}&C_{34}&C_{35}&C_{36}\\ C_{41}&C_{42}&C_{43}&C_{44}&C_{45}&C_{46}\\ C_{51}&C_{52}&C_{53}&C_{54}&C_{55}&C_{56}\\ C_{61}&C_{62}&C_{63}&C_{64}&C_{65}&C_{66}\end{array}\right] (1.80)

Αν τώρα θεωρήσουμε ένα υπερελαστικό υλικό, τότε οι τάσεις θεωρούνται ότι προκύπτουν από μια συνάρτηση δυναμικού 𝒲\mathcal{W} ως εξής:

σij=𝒲εij\sigma_{ij}=\dfrac{\partial\mathcal{W}}{\partial\varepsilon_{ij}} (1.81)

Η παραπάνω εξίσωση οδηγεί στον ακόλουθο ορισμό για τον τανυστή ελαστικότητας:

Cijkl=2𝒲ϵijεkl.C_{ijkl}=\dfrac{\partial^{2}\mathcal{W}}{\partial\epsilon_{ij}\partial% \varepsilon_{kl}}. (1.82)

Επειδή δεν έχει σημασία η σειρά των παραγωγίσεων, προκύπτει επιπλέον ότι ισχύει:

Cijkl=Clkij.C_{ijkl}=C_{lkij}. (1.83)

Επομένως ο αριθμός των σταθερών που περιγράφουν τον τανυστή CijklC_{ijkl} μειώνεται από 36 σε 21. Σε μητρωική μορφή μπορούμε τώρα να γράψουμε:

𝐂=[C11C12C13C14C15C16C12C22C23C24C25C26C13C23C33C34C35C36C14C24C34C44C45C46C15C25C35C45C55C56C16C26C36C46C56C66]\mathbf{C}=\left[\begin{array}[]{cccccc}C_{11}&C_{12}&C_{13}&C_{14}&C_{15}&C_{% 16}\\ C_{12}&C_{22}&C_{23}&C_{24}&C_{25}&C_{26}\\ C_{13}&C_{23}&C_{33}&C_{34}&C_{35}&C_{36}\\ C_{14}&C_{24}&C_{34}&C_{44}&C_{45}&C_{46}\\ C_{15}&C_{25}&C_{35}&C_{45}&C_{55}&C_{56}\\ C_{16}&C_{26}&C_{36}&C_{46}&C_{56}&C_{66}\end{array}\right] (1.84)

1.3.2 Μετασχηματισμός του ελαστικού τανυστή

Ο μετασχηματισμός του τανυστή CijklC_{ijkl} ακολουθεί τους κανόνες μετασχηματισμού των τανυστών τετάρτης τάξης [1]. Πιο συγκεκριμένα, ορίζεται το μητρώο μετασχηματισμού 𝐓\mathbf{T} ως:

𝐓:=[r112r122r1322r12r132r13r112r11r12r212r222r2322r22r232r23r212r21r22r312r322r3322r32r332r33r312r31r32r21r31r22r32r23r33r22r33+r23r32r23r31+r21r33r21r32+r22r31r31r11r32r12r33r13r32r13+r33r12r33r11+r31r13r31r12+r32r11r11r21r12r22r13r23r12r23+r13r22r13r21+r11r23r11r22+r12r21]\mathbf{T}:=\left[\begin{array}[]{ccc|ccc}r_{11}^{2}&r_{12}^{2}&r_{13}^{2}&2r_% {12}r_{13}&2r_{13}r_{11}&2r_{11}r_{12}\\ r_{21}^{2}&r_{22}^{2}&r_{23}^{2}&2r_{22}r_{23}&2r_{23}r_{21}&2r_{21}r_{22}\\ r_{31}^{2}&r_{32}^{2}&r_{33}^{2}&2r_{32}r_{33}&2r_{33}r_{31}&2r_{31}r_{32}\\ \hline r_{21}r_{31}&r_{22}r_{32}&r_{23}r_{33}&r_{22}r_{33}+r_{23}r_{32}&r_{23}% r_{31}+r_{21}r_{33}&r_{21}r_{32}+r_{22}r_{31}\\ r_{31}r_{11}&r_{32}r_{12}&r_{33}r_{13}&r_{32}r_{13}+r_{33}r_{12}&r_{33}r_{11}+% r_{31}r_{13}&r_{31}r_{12}+r_{32}r_{11}\\ r_{11}r_{21}&r_{12}r_{22}&r_{13}r_{23}&r_{12}r_{23}+r_{13}r_{22}&r_{13}r_{21}+% r_{11}r_{23}&r_{11}r_{22}+r_{12}r_{21}\end{array}\right] (1.85)

Στο παραπάνω πίνακα 𝐓\mathbf{T}, το rijr_{ij} συμβολίζει το συνημίτονο της γωνίας μεταξύ του αρχικού i-i- συστήματος συντεταγμένων και του τελικού j-j- συστήματος συντεταγμένων.

rij=cos(x¯i,xj)r_{ij}=\cos(\bar{x}_{i},x_{j}) (1.86)

Σημειώνεται ότι, λόγω της ορθοκανονικότητας των συνημίτονων κατεύθυνσης, ισχύει:

rijrik=δjkr_{ij}r_{ik}=\delta_{jk} (1.87)

όπου το δij\delta_{ij} ονομάζεται και δέλτα του Kronecker.

Βάσει των παραπάνω θα οριστούν στη συνέχεια τα εγκαρσίως ισότροπα και τα ορθότροπα ελαστικά υλικά, που εμφανίζουν ιδιαίτερο ενδιαφέρον στη γεωτεχνική μηχανική.

1.3.3 Eγκαρσίως ισότροπα ελαστικά υλικά

Μια πρώτη ειδική περίπτωση ενός ανισότροπου υλικού είναι αυτή κατά την οποία το υλικό περιέχει ένα επίπεδο ισοτροπίας, όπως συχνά παρατηρείται στην περίπτωση ιζηματογενών πετρωμάτων. Κάθε στρώση έχει περίπου τις ίδιες ελαστικές ιδιότητες εντός του επιπέδου της αλλά οι ιδιότητες του υλικού διαφοροποιούνται με το βάθος. Το επίπεδο κάθε στρώσης αποτελεί το επίπεδο ισοτροπίας και ο κάθετος σε αυτό άξονας τον αντίστοιχο άξονα συμμετρίας [3]. Άλλες περιπτώσεις εγκαρσίως ισότροπων ελαστικών υλικών δίνονται στο Σχήμα 1.9.

Αυτή η παραδοχή έχει ως αποτέλεσμα να μπορούμε να περιστρέψουμε το στερεό σώμα γύρω από έναν άξονα συμμετρίας, ο οποίος ορίζεται ως ο κάθετος άξονας στο ισότροπο επίπεδο, χωρίς να μεταβάλλεται η ελαστική συμπεριφορά του υλικού.

Σχήμα 1.9: Περιπτώσεις εγκαρσίως ισότροπων ελαστικών υλικών.
Σχήμα 1.10: Εγκαρσίως ισότροπo ελαστικό υλικό, το επίπεδο συμμετρίας (x1x_{1}-x2x_{2}) και ο άξονας συμμετρίας x3x_{3}.

Στη συνέχεια περιγράφεται η διαδικασία εύρεσης του ελαστικού μητρώου για τα εγκαρσίως ισότροπα ελαστικά υλικά, θεωρώντας ότι ο άξονας 𝒆3\boldsymbol{e}_{3} (βλ. Σχήμα 1.10) είναι ο άξονας συμμετρίας, δηλαδή κάθετος στο επίπεδο ισοτροπίας του υλικού.

Περιστροφή κατά 900{}^{0}

Ο πίνακας μετασχηματισμού σε αυτή την περίπτωση δίνεται ως:

𝐓=[010000100000001000000010000-10000000-1]\mathbf{T}=\left[\begin{array}[]{cccccc}0&1&0&0&0&0\\ 1&0&0&0&0&0\\ 0&0&1&0&0&0\\ 0&0&0&0&1&0\\ 0&0&0&-1&0&0\\ 0&0&0&0&0&-1\end{array}\right] (1.88)

Επομένως το μετασχηματισμένο ελαστικό μητρώο προκύπτει σύμφωνα με τον μετασχηματισμό:

𝐂¯=𝐓𝐂𝐓T=[C22C12C23C25-C24-C26C12C11C13C15-C14-C16C23C13C33C35-C34-C36C25C15C35C55-C45-C56-C24-C14-C34-C45C44C46-C26-C16-C36-C56C46C66]\bar{\mathbf{C}}=\mathbf{T}\mathbf{C}\mathbf{T}^{T}=\left[\begin{array}[]{% cccccc}C_{22}&C_{12}&C_{23}&C_{25}&-C_{24}&-C_{26}\\ C_{12}&C_{11}&C_{13}&C_{15}&-C_{14}&-C_{16}\\ C_{23}&C_{13}&C_{33}&C_{35}&-C_{34}&-C_{36}\\ C_{25}&C_{15}&C_{35}&C_{55}&-C_{45}&-C_{56}\\ -C_{24}&-C_{14}&-C_{34}&-C_{45}&C_{44}&C_{46}\\ -C_{26}&-C_{16}&-C_{36}&-C_{56}&C_{46}&C_{66}\end{array}\right] (1.89)

Λόγω των παραδοχών που διέπουν τα εγκαρσίως ισότροπα υλικά, θα πρέπει να ισχύει:

[C11C12C13C14C15C16C12C22C23C24C25C26C13C23C33C34C35C36C14C24C34C44C45C46C15C25C35C45C55C56C16C26C36C46C56C66]=[C22C12C23C25-C24-C26C12C11C13C15-C14-C16C23C13C33C35-C34-C36C25C15C35C55-C45-C56-C24-C14-C34-C45C44C46-C26-C16-C36-C56C46C66]\left[\begin{array}[]{cccccc}C_{11}&C_{12}&C_{13}&C_{14}&C_{15}&C_{16}\\ C_{12}&C_{22}&C_{23}&C_{24}&C_{25}&C_{26}\\ C_{13}&C_{23}&C_{33}&C_{34}&C_{35}&C_{36}\\ C_{14}&C_{24}&C_{34}&C_{44}&C_{45}&C_{46}\\ C_{15}&C_{25}&C_{35}&C_{45}&C_{55}&C_{56}\\ C_{16}&C_{26}&C_{36}&C_{46}&C_{56}&C_{66}\end{array}\right]=\left[\begin{array% }[]{cccccc}C_{22}&C_{12}&C_{23}&C_{25}&-C_{24}&-C_{26}\\ C_{12}&C_{11}&C_{13}&C_{15}&-C_{14}&-C_{16}\\ C_{23}&C_{13}&C_{33}&C_{35}&-C_{34}&-C_{36}\\ C_{25}&C_{15}&C_{35}&C_{55}&-C_{45}&-C_{56}\\ -C_{24}&-C_{14}&-C_{34}&-C_{45}&C_{44}&C_{46}\\ -C_{26}&-C_{16}&-C_{36}&-C_{56}&C_{46}&C_{66}\end{array}\right] (1.90)

Επομένως:

C11\displaystyle C_{11} =C22\displaystyle=C_{22} (1.91)
C25\displaystyle C_{25} =-C14\displaystyle=-C_{14} (1.92)
C56\displaystyle C_{56} =C46\displaystyle=C_{46} (1.93)
C34\displaystyle C_{34} =C35\displaystyle=C_{35} (1.94)
C44\displaystyle C_{44} =C55\displaystyle=C_{55} (1.95)
C16\displaystyle C_{16} =-C26\displaystyle=-C_{26} (1.96)
C46\displaystyle C_{46} =-C56\displaystyle=-C_{56} (1.97)
C23\displaystyle C_{23} =C13\displaystyle=C_{13} (1.98)
C45\displaystyle C_{45} =0\displaystyle=0 (1.99)
C13\displaystyle C_{13} =C23\displaystyle=C_{23} (1.100)
C55\displaystyle C_{55} =C44\displaystyle=C_{44} (1.101)
C15\displaystyle C_{15} =-C24\displaystyle=-C_{24} (1.102)
C35\displaystyle C_{35} =-C34\displaystyle=-C_{34} (1.103)
C22\displaystyle C_{22} =C11\displaystyle=C_{11} (1.104)
C14\displaystyle C_{14} =C25\displaystyle=C_{25} (1.105)
C26\displaystyle C_{26} =-C16\displaystyle=-C_{16} (1.106)
C24\displaystyle C_{24} =C15\displaystyle=C_{15} (1.107)
C36\displaystyle C_{36} =0\displaystyle=0 (1.108)

Σύμφωνα με τα παραπάνω, το μητρώο ελαστικότητας ενός εγκαρσίως ισότροπου υλικού μπορεί να περιγραφεί ως:

𝐂¯=[C11C12C13-C14-C15C16C12C11C13C15-C14-C16C13C13C33-C34-C340-C14C15-C34C440-C46-C15-C14-C340C44C46C16-C160-C46C46C66]\bar{\mathbf{C}}=\left[\begin{array}[]{cccccc}C_{11}&C_{12}&C_{13}&-C_{14}&-C_% {15}&C_{16}\\ C_{12}&C_{11}&C_{13}&C_{15}&-C_{14}&-C_{16}\\ C_{13}&C_{13}&C_{33}&-C_{34}&-C_{34}&0\\ -C_{14}&C_{15}&-C_{34}&C_{44}&0&-C_{46}\\ -C_{15}&-C_{14}&-C_{34}&0&C_{44}&C_{46}\\ C_{16}&-C_{16}&0&-C_{46}&C_{46}&C_{66}\end{array}\right] (1.109)

Περιστροφή κατά 1800{}^{0}

Θεωρείται τώρα μια δεύτερη περιστροφή κατά 1800{}^{0}. Το αντίστοιχο μητρώο μετασχηματισμού σε αυτή την περίπτωση δίνεται ως:

𝐓=[100000010000001000000-1000000-10000001]\mathbf{T}=\left[\begin{array}[]{cccccc}1&0&0&0&0&0\\ 0&1&0&0&0&0\\ 0&0&1&0&0&0\\ 0&0&0&-1&0&0\\ 0&0&0&0&-1&0\\ 0&0&0&0&0&1\end{array}\right] (1.110)

Ο μετασχηματισμένος τανυστής της ελαστικότητας δίνεται τώρα σύμφωνα με τον μετασχηματισμό:

𝐂¯¯=𝐓𝐂¯𝐓T=[C11C12C13C14C15C16C12C11C13-C15C14-C16C13C13C33C34C340C14-C15C34C440C46C15C14C340C44-C46C16-C160C46-C46C66]\bar{\bar{\mathbf{C}}}=\mathbf{T}\bar{\mathbf{C}}\mathbf{T}^{T}=\left[\begin{% array}[]{cccccc}C_{11}&C_{12}&C_{13}&C_{14}&C_{15}&C_{16}\\ C_{12}&C_{11}&C_{13}&-C_{15}&C_{14}&-C_{16}\\ C_{13}&C_{13}&C_{33}&C_{34}&C_{34}&0\\ C_{14}&-C_{15}&C_{34}&C_{44}&0&C_{46}\\ C_{15}&C_{14}&C_{34}&0&C_{44}&-C_{46}\\ C_{16}&-C_{16}&0&C_{46}&-C_{46}&C_{66}\end{array}\right] (1.111)

Και τώρα θα πρέπει να ισχύει:

[C11C12C13-C14-C15C16C12C11C13C15-C14-C16C13C13C33-C34-C340-C14C15-C34C440-C46-C15-C14-C340C44C46C16-C160-C46C46C66]=[C11C12C13C14C15C16C12C11C13-C15C14-C16C13C13C33C34C340C14-C15C34C440C46C15C14C340C44-C46C16-C160C46-C46C66]\left[\begin{array}[]{cccccc}C_{11}&C_{12}&C_{13}&-C_{14}&-C_{15}&C_{16}\\ C_{12}&C_{11}&C_{13}&C_{15}&-C_{14}&-C_{16}\\ C_{13}&C_{13}&C_{33}&-C_{34}&-C_{34}&0\\ -C_{14}&C_{15}&-C_{34}&C_{44}&0&-C_{46}\\ -C_{15}&-C_{14}&-C_{34}&0&C_{44}&C_{46}\\ C_{16}&-C_{16}&0&-C_{46}&C_{46}&C_{66}\end{array}\right]=\left[\begin{array}[]% {cccccc}C_{11}&C_{12}&C_{13}&C_{14}&C_{15}&C_{16}\\ C_{12}&C_{11}&C_{13}&-C_{15}&C_{14}&-C_{16}\\ C_{13}&C_{13}&C_{33}&C_{34}&C_{34}&0\\ C_{14}&-C_{15}&C_{34}&C_{44}&0&C_{46}\\ C_{15}&C_{14}&C_{34}&0&C_{44}&-C_{46}\\ C_{16}&-C_{16}&0&C_{46}&-C_{46}&C_{66}\end{array}\right] (1.112)

Επομένως:

C46\displaystyle C_{46} =0\displaystyle=0 (1.113)
C14\displaystyle C_{14} =0\displaystyle=0 (1.114)
C15\displaystyle C_{15} =0\displaystyle=0 (1.115)
C34\displaystyle C_{34} =0\displaystyle=0 (1.116)

Σύμφωνα με τα παραπάνω, το μητρώο ελαστικότητας ενός εγκαρσίως ισότροπου υλικού μπορεί να περιγραφεί ως:

𝐂¯¯=[C11C12C13-C14-C15C16C12C11C13C15-C14-C16C13C13C33-C34-C340-C14C15-C34C440-C46-C15-C14-C340C44C46C16-C160-C46C46C66]\bar{\bar{\mathbf{C}}}=\left[\begin{array}[]{cccccc}C_{11}&C_{12}&C_{13}&-C_{1% 4}&-C_{15}&C_{16}\\ C_{12}&C_{11}&C_{13}&C_{15}&-C_{14}&-C_{16}\\ C_{13}&C_{13}&C_{33}&-C_{34}&-C_{34}&0\\ -C_{14}&C_{15}&-C_{34}&C_{44}&0&-C_{46}\\ -C_{15}&-C_{14}&-C_{34}&0&C_{44}&C_{46}\\ C_{16}&-C_{16}&0&-C_{46}&C_{46}&C_{66}\end{array}\right] (1.117)

Περιστροφή κατά 450{}^{0}

Τέλος θεωρείται μια περιστροφή κατά 450{}^{0}. Το αντίστοιχο μητρώο μετασχηματισμού σε αυτή την περίπτωση δίνεται ως:

𝐓=[121200011212000-1001000000122-12200001221220-12120000]\mathbf{T}=\left[\begin{array}[]{cccccc}\frac{1}{2}&\frac{1}{2}&0&0&0&1\\ \frac{1}{2}&\frac{1}{2}&0&0&0&-1\\ 0&0&1&0&0&0\\ 0&0&0&\frac{1}{2}\sqrt{2}&-\frac{1}{2}\sqrt{2}&0\\ 0&0&0&\frac{1}{2}\sqrt{2}&\frac{1}{2}\sqrt{2}&0\\ -\frac{1}{2}&\frac{1}{2}&0&0&0&0\end{array}\right] (1.118)

Ο μετασχηματισμένος τανυστής της ελαστικότητας δίνεται τώρα σύμφωνα με τον μετασχηματισμό:

𝐂¯¯¯=𝐓𝐂¯¯𝐓T=[12C11+12C12+C6612C11+12C12-C66C1300-C1612C11+12C12-C6612C11+12C12+C66C1300C16C13C13C33000000C44000000C440-C16C1600012C11-12C12]\bar{\bar{\bar{\mathbf{C}}}}=\mathbf{T}\bar{\bar{\mathbf{C}}}\mathbf{T}^{T}=% \left[\begin{array}[]{cccccc}\frac{1}{2}C_{11}+\frac{1}{2}C_{12}+C_{66}&\frac{% 1}{2}C_{11}+\frac{1}{2}C_{12}-C_{66}&C_{13}&0&0&-C_{16}\\ \frac{1}{2}C_{11}+\frac{1}{2}C_{12}-C_{66}&\frac{1}{2}C_{11}+\frac{1}{2}C_{12}% +C_{66}&C_{13}&0&0&C_{16}\\ C_{13}&C_{13}&C_{33}&0&0&0\\ 0&0&0&C_{44}&0&0\\ 0&0&0&0&C_{44}&0\\ -C_{16}&C_{16}&0&0&0&\frac{1}{2}C_{11}-\frac{1}{2}C_{12}\end{array}\right] (1.119)

Θα πρέπει επίσης να ισχύει:

[C11C12C1300C16C12C11C1300-C16C13C13C33000000C44000000C440C16-C16000C66]=[12C11+12C12+C6612C11+12C12-C66C1300-C1612C11+12C12-C6612C11+12C12+C66C1300C16C13C13C33000000C44000000C440-C16C1600012C11-12C12]\left[\begin{array}[]{cccccc}C_{11}&C_{12}&C_{13}&0&0&C_{16}\\ C_{12}&C_{11}&C_{13}&0&0&-C_{16}\\ C_{13}&C_{13}&C_{33}&0&0&0\\ 0&0&0&C_{44}&0&0\\ 0&0&0&0&C_{44}&0\\ C_{16}&-C_{16}&0&0&0&C_{66}\end{array}\right]=\\ \left[\begin{array}[]{cccccc}\frac{1}{2}C_{11}+\frac{1}{2}C_{12}+C_{66}&\frac{% 1}{2}C_{11}+\frac{1}{2}C_{12}-C_{66}&C_{13}&0&0&-C_{16}\\ \frac{1}{2}C_{11}+\frac{1}{2}C_{12}-C_{66}&\frac{1}{2}C_{11}+\frac{1}{2}C_{12}% +C_{66}&C_{13}&0&0&C_{16}\\ C_{13}&C_{13}&C_{33}&0&0&0\\ 0&0&0&C_{44}&0&0\\ 0&0&0&0&C_{44}&0\\ -C_{16}&C_{16}&0&0&0&\frac{1}{2}C_{11}-\frac{1}{2}C_{12}\end{array}\right] (1.120)

Και επομένως:

C11\displaystyle C_{11} =C12+2C66\displaystyle=C_{12}+2C_{66} (1.121)
C12\displaystyle C_{12} =C11-2C66\displaystyle=C_{11}-2C_{66} (1.122)
C16\displaystyle C_{16} =0\displaystyle=0 (1.123)
C66\displaystyle C_{66} =12C11-12C12\displaystyle=\frac{1}{2}C_{11}-\frac{1}{2}C_{12} (1.124)

Σύμφωνα με τα παραπάνω, το μητρώο ελαστικότητας ενός εγκαρσίως ισότροπου υλικού μπορεί τελικά να περιγραφεί ως:

𝐂¯¯¯=[C11C12C13000C12C11C13000C13C13C33000000C44000000C4400000012C11-12C12]\bar{\bar{\bar{\mathbf{C}}}}=\left[\begin{array}[]{cccccc}C_{11}&C_{12}&C_{13}% &0&0&0\\ C_{12}&C_{11}&C_{13}&0&0&0\\ C_{13}&C_{13}&C_{33}&0&0&0\\ 0&0&0&C_{44}&0&0\\ 0&0&0&0&C_{44}&0\\ 0&0&0&0&0&\frac{1}{2}C_{11}-\frac{1}{2}C_{12}\end{array}\right] (1.126)

Περιστροφή κατά 1350{}^{0}

Μπορεί να αποδειχτεί ότι κάθε επιπλέον περιστροφή αφήνει το μητρώο ελαστικότητας ανεπηρέαστο. Π.χ. το μητρώο μετασχηματισμού για μια γωνία ίση με 1350{}^{0} δίνεται ως:

𝐓=[1212000-1121200010010000001221220000-122122012-120000]\mathbf{T}=\left[\begin{array}[]{cccccc}\frac{1}{2}&\frac{1}{2}&0&0&0&-1\\ \frac{1}{2}&\frac{1}{2}&0&0&0&1\\ 0&0&1&0&0&0\\ 0&0&0&\frac{1}{2}\sqrt{2}&\frac{1}{2}\sqrt{2}&0\\ 0&0&0&-\frac{1}{2}\sqrt{2}&\frac{1}{2}\sqrt{2}&0\\ \frac{1}{2}&-\frac{1}{2}&0&0&0&0\end{array}\right] (1.127)

Εύκολα διακρίνεται ότι ο μετασχηματισμένος τανυστής της ελαστικότητας,

𝐂¯¯¯¯=𝐓𝐂¯¯¯𝐓T=[C11C12C13000C12C11C13000C13C13C33000000C44000000C4400000012C11-12C12]\bar{\bar{\bar{\bar{\mathbf{C}}}}}=\mathbf{T}\bar{\bar{\bar{\mathbf{C}}}}% \mathbf{T}^{T}=\left[\begin{array}[]{cccccc}C_{11}&C_{12}&C_{13}&0&0&0\\ C_{12}&C_{11}&C_{13}&0&0&0\\ C_{13}&C_{13}&C_{33}&0&0&0\\ 0&0&0&C_{44}&0&0\\ 0&0&0&0&C_{44}&0\\ 0&0&0&0&0&\frac{1}{2}C_{11}-\frac{1}{2}C_{12}\end{array}\right] (1.128)

παραμένει ανεπηρέαστος.

Τελική διατύπωση

Από τα παραπάνω καταλήγουμε στο συμπέρασμα ότι ο ελαστικός τανυστής για την περίπτωση εγκαρσίως ισότροπων υλικών δίνεται ως:

𝐂transverse=[C11C12C13000C12C11C13000C13C13C33000000C44000000C4400000012C11-12C12]\mathbf{C}^{transverse}=\left[\begin{array}[]{cccccc}C_{11}&C_{12}&C_{13}&0&0&% 0\\ C_{12}&C_{11}&C_{13}&0&0&0\\ C_{13}&C_{13}&C_{33}&0&0&0\\ 0&0&0&C_{44}&0&0\\ 0&0&0&0&C_{44}&0\\ 0&0&0&0&0&\frac{1}{2}C_{11}-\frac{1}{2}C_{12}\end{array}\right] (1.129)

Για τον ορισμό του απαιτούνται πέντε ανεξάρτητες σταθερές, οι C11C_{11}, C33C_{33}, C12C_{12}, C13C_{13} και C44C_{44} αντίστοιχα.

1.3.4 Ορθότροπα υλικά

Ένα ορθότροπο υλικό έχει τουλάχιστον δύο επίπεδα συμμετρίας ορθογώνια μεταξύ τους, όπου οι ιδιότητες του υλικού είναι ανεξάρτητες από τη διεύθυνση εντός των επιπέδων αυτών [6]. Η παραδοχή αυτή χρησιμοποιείται στη συνέχεια για τον προσδιορισμό του αντίστοιχου μητρώου ελαστικότητας. Οι άξονες του συστήματος συντεταγμένων που χρησιμοποιείται στη συνέχεια θεωρείται ότι ταυτίζονται με τα τρία επίπεδα συμμετρίας του υλικού.

Αντανάκλαση ως προς το επίπεδο 2-3

Ας θεωρήσουμε αρχικά μια αντανάκλαση ως προς το επίπεδο 2-3. Το αντίστοιχο μητρώο μετασχηματισμού σε αυτή την περίπτωση δίνεται από την ακόλουθη εξίσωση:

𝐓=[1000000100000010000001000000-1000000-1]\mathbf{T}=\left[\begin{array}[]{cccccc}1&0&0&0&0&0\\ 0&1&0&0&0&0\\ 0&0&1&0&0&0\\ 0&0&0&1&0&0\\ 0&0&0&0&-1&0\\ 0&0&0&0&0&-1\end{array}\right] (1.130)

Ο μετασχηματισμένος τανυστής της ελαστικότητας δίνεται τώρα σύμφωνα με τη σχέση:

𝐂¯=[C11C12C13C14-C15-C16C12C22C23C24-C25-C26C13C23C33C34-C35-C36C14C24C34C44-C45-C46-C15-C25-C35-C45C55C56-C16-C26-C36-C46C56C66]\bar{\mathbf{C}}=\left[\begin{array}[]{cccccc}C_{11}&C_{12}&C_{13}&C_{14}&-C_{% 15}&-C_{16}\\ C_{12}&C_{22}&C_{23}&C_{24}&-C_{25}&-C_{26}\\ C_{13}&C_{23}&C_{33}&C_{34}&-C_{35}&-C_{36}\\ C_{14}&C_{24}&C_{34}&C_{44}&-C_{45}&-C_{46}\\ -C_{15}&-C_{25}&-C_{35}&-C_{45}&C_{55}&C_{56}\\ -C_{16}&-C_{26}&-C_{36}&-C_{46}&C_{56}&C_{66}\end{array}\right] (1.131)

Βάσει των παραδοχών για τα ορθότροπα ελαστικά υλικά, ισχύει:

𝐂=𝐂¯\mathbf{C}=\bar{\mathbf{C}} (1.132)

Επομένως:

[C11C12C13C14C15C16C12C22C23C24C25C26C13C23C33C34C35C36C14C24C34C44C45C46C15C25C35C45C55C56C16C26C36C46C56C66]=[C11C12C13C14-C15-C16C12C22C23C24-C25-C26C13C23C33C34-C35-C36C14C24C34C44-C45-C46-C15-C25-C35-C45C55C56-C16-C26-C36-C46C56C66]\left[\begin{array}[]{cccccc}C_{11}&C_{12}&C_{13}&C_{14}&C_{15}&C_{16}\\ C_{12}&C_{22}&C_{23}&C_{24}&C_{25}&C_{26}\\ C_{13}&C_{23}&C_{33}&C_{34}&C_{35}&C_{36}\\ C_{14}&C_{24}&C_{34}&C_{44}&C_{45}&C_{46}\\ C_{15}&C_{25}&C_{35}&C_{45}&C_{55}&C_{56}\\ C_{16}&C_{26}&C_{36}&C_{46}&C_{56}&C_{66}\end{array}\right]=\left[\begin{array% }[]{cccccc}C_{11}&C_{12}&C_{13}&C_{14}&-C_{15}&-C_{16}\\ C_{12}&C_{22}&C_{23}&C_{24}&-C_{25}&-C_{26}\\ C_{13}&C_{23}&C_{33}&C_{34}&-C_{35}&-C_{36}\\ C_{14}&C_{24}&C_{34}&C_{44}&-C_{45}&-C_{46}\\ -C_{15}&-C_{25}&-C_{35}&-C_{45}&C_{55}&C_{56}\\ -C_{16}&-C_{26}&-C_{36}&-C_{46}&C_{56}&C_{66}\end{array}\right] (1.133)

Η παραπάνω απαίτηση μπορεί να ικανοποιηθεί μόνο όταν:

C25\displaystyle C_{25} =0\displaystyle=0 (1.134)
C16\displaystyle C_{16} =0\displaystyle=0 (1.135)
C46\displaystyle C_{46} =0\displaystyle=0 (1.136)
C45\displaystyle C_{45} =0\displaystyle=0 (1.137)
C35\displaystyle C_{35} =0\displaystyle=0 (1.138)
C15\displaystyle C_{15} =0\displaystyle=0 (1.139)
C36\displaystyle C_{36} =0\displaystyle=0 (1.140)
C26\displaystyle C_{26} =0\displaystyle=0 (1.141)

Επομένως, ο τανυστής της ελαστικότητας ισούται με:

𝐂¯=[C11C12C13C1400C12C22C23C2400C13C23C33C3400C14C24C34C44000000C55C560000C56C66]\bar{\mathbf{C}}=\left[\begin{array}[]{cccccc}C_{11}&C_{12}&C_{13}&C_{14}&0&0% \\ C_{12}&C_{22}&C_{23}&C_{24}&0&0\\ C_{13}&C_{23}&C_{33}&C_{34}&0&0\\ C_{14}&C_{24}&C_{34}&C_{44}&0&0\\ 0&0&0&0&C_{55}&C_{56}\\ 0&0&0&0&C_{56}&C_{66}\end{array}\right] (1.142)

Αντανάκλαση ως προς το επίπεδο 3-1

Ας θεωρήσουμε τώρα μια αντανάκλαση ως προς το επίπεδο 3-1. Το αντίστοιχο μητρώο μετασχηματισμού σε αυτή την περίπτωση δίνεται από την ακόλουθη εξίσωση:

𝐓=[100000010000001000000-10000001000000-1]\mathbf{T}=\left[\begin{array}[]{cccccc}1&0&0&0&0&0\\ 0&1&0&0&0&0\\ 0&0&1&0&0&0\\ 0&0&0&-1&0&0\\ 0&0&0&0&1&0\\ 0&0&0&0&0&-1\end{array}\right] (1.143)

Ο μετασχηματισμένος τανυστής της ελαστικότητας δίνεται τώρα σύμφωνα με τη σχέση:

𝐂¯¯=[C11C12C13-C1400C12C22C23-C2400C13C23C33-C3400-C14-C24-C34C44000000C55-C560000-C56C66]\bar{\bar{\mathbf{C}}}=\left[\begin{array}[]{cccccc}C_{11}&C_{12}&C_{13}&-C_{1% 4}&0&0\\ C_{12}&C_{22}&C_{23}&-C_{24}&0&0\\ C_{13}&C_{23}&C_{33}&-C_{34}&0&0\\ -C_{14}&-C_{24}&-C_{34}&C_{44}&0&0\\ 0&0&0&0&C_{55}&-C_{56}\\ 0&0&0&0&-C_{56}&C_{66}\end{array}\right] (1.144)

Βάσει των παραδοχών για τα ορθότροπα ελαστικά υλικά, ισχύει:

𝐂¯=𝐂¯¯\bar{\mathbf{C}}=\bar{\bar{\mathbf{C}}} (1.145)

Επομένως:

[C11C12C13C1400C12C22C23C2400C13C23C33C3400C14C24C34C44000000C55C560000C56C66]=[C11C12C13-C1400C12C22C23-C2400C13C23C33-C3400-C14-C24-C34C44000000C55-C560000-C56C66]\left[\begin{array}[]{cccccc}C_{11}&C_{12}&C_{13}&C_{14}&0&0\\ C_{12}&C_{22}&C_{23}&C_{24}&0&0\\ C_{13}&C_{23}&C_{33}&C_{34}&0&0\\ C_{14}&C_{24}&C_{34}&C_{44}&0&0\\ 0&0&0&0&C_{55}&C_{56}\\ 0&0&0&0&C_{56}&C_{66}\end{array}\right]=\left[\begin{array}[]{cccccc}C_{11}&C_% {12}&C_{13}&-C_{14}&0&0\\ C_{12}&C_{22}&C_{23}&-C_{24}&0&0\\ C_{13}&C_{23}&C_{33}&-C_{34}&0&0\\ -C_{14}&-C_{24}&-C_{34}&C_{44}&0&0\\ 0&0&0&0&C_{55}&-C_{56}\\ 0&0&0&0&-C_{56}&C_{66}\end{array}\right] (1.146)

Η παραπάνω απαίτηση μπορεί να ικανοποιηθεί μόνο όταν:

C14\displaystyle C_{14} =0\displaystyle=0 (1.147)
C56\displaystyle C_{56} =0\displaystyle=0 (1.148)
C34\displaystyle C_{34} =0\displaystyle=0 (1.149)
C24\displaystyle C_{24} =0\displaystyle=0 (1.150)

Επομένως, ο τανυστής της ελαστικότητας μπορεί να γραφτεί ως:

𝐂¯¯=[C11C12C13000C12C22C23000C13C23C33000000C44000000C55000000C66]\bar{{\bar{\mathbf{C}}}}=\left[\begin{array}[]{cccccc}C_{11}&C_{12}&C_{13}&0&0% &0\\ C_{12}&C_{22}&C_{23}&0&0&0\\ C_{13}&C_{23}&C_{33}&0&0&0\\ 0&0&0&C_{44}&0&0\\ 0&0&0&0&C_{55}&0\\ 0&0&0&0&0&C_{66}\end{array}\right] (1.151)

Αντανάκλαση ως προς το επίπεδο 1-2

Μπορεί να αποδειχτεί ότι κάθε επιπλέον μετασχηματισμός αφήνει το μητρώο ελαστικότητας ανεπηρέαστο. Έστω π.χ. μια αντανάκλαση ως προς το επίπεδο 1-2. Τότε το μητρώο μετασχηματισμού γράφεται ως:

𝐓=[100000010000001000000-1000000-10000001]\mathbf{T}=\left[\begin{array}[]{cccccc}1&0&0&0&0&0\\ 0&1&0&0&0&0\\ 0&0&1&0&0&0\\ 0&0&0&-1&0&0\\ 0&0&0&0&-1&0\\ 0&0&0&0&0&1\end{array}\right] (1.152)

Ο μετασχηματισμένος τανυστής της ελαστικότητας δίνεται τώρα σύμφωνα με τη σχέση:

𝐂¯¯=[C11C12C13000C12C22C23000C13C23C33000000C44000000C55000000C66]\bar{\bar{\mathbf{C}}}=\left[\begin{array}[]{cccccc}C_{11}&C_{12}&C_{13}&0&0&0% \\ C_{12}&C_{22}&C_{23}&0&0&0\\ C_{13}&C_{23}&C_{33}&0&0&0\\ 0&0&0&C_{44}&0&0\\ 0&0&0&0&C_{55}&0\\ 0&0&0&0&0&C_{66}\end{array}\right] (1.153)

ο οποίος είναι ίδιος με τον προηγούμενο.

Τελική διατύπωση

Από τα παραπάνω καταλήγουμε στο συμπέρασμα ότι ο ελαστικός τανυστής για την περίπτωση ορθότροπων υλικών δίνεται ως:

𝐂orthotropic=[C11C12C13000C12C22C23000C13C23C33000000C44000000C55000000C66]\mathbf{C}^{orthotropic}=\left[\begin{array}[]{cccccc}C_{11}&C_{12}&C_{13}&0&0% &0\\ C_{12}&C_{22}&C_{23}&0&0&0\\ C_{13}&C_{23}&C_{33}&0&0&0\\ 0&0&0&C_{44}&0&0\\ 0&0&0&0&C_{55}&0\\ 0&0&0&0&0&C_{66}\end{array}\right] (1.154)

και για τον ορισμό του απαιτούνται εννέα ανεξάρτητες ελαστικές σταθερές, δηλαδή οι C11C_{11}, C22C_{22}, C33C_{33}, C12C_{12}, C13C_{13} C23C_{23}, C44C_{44}, C55C_{55} και C66C_{66}.

1.3.5 Ισότροπα υλικά

Ισότροπα υλικά θεωρούνται τα υλικά εκείνα των οποίων οι μηχανικές τους ιδιότητες σε κάθε σημείο τους είναι ίδιες ως προς κάθε διεύθυνση [8]. Επομένως το μητρώο ελαστικότητας των ισότροπων υλικών θα πρέπει να παραμένει ανεπηρέαστο, ανεξαρτήτως του μετασχηματισμού που μπορεί να εφαρμόσουμε σε αυτό.

Για τη μελέτη των ισότροπων υλικών αρχικά ορίζουμε τον τανυστή των παραμορφώσεων ως:

εij=13εkkδij+(εij-13εkkδij)\varepsilon_{ij}=\dfrac{1}{3}\varepsilon_{kk}\delta_{ij}+\left(\varepsilon_{ij% }-\dfrac{1}{3}\varepsilon_{kk}\delta_{ij}\right) (1.155)

όπου δij\delta_{ij} το δ\delta του Dirac, το οποίο ορίζεται ως:

δij={1ανij0αλλιώς.\delta_{ij}=\begin{cases}1&\text{αν}\,i\equiv j\\ 0&\text{αλλιώς.}\end{cases} (1.156)

Ο τανυστής των παραμορφώσεων (1.155) επομένως ορίζεται ως το άθροισμα,

  • ενός σταθερού τανυστή, 13εkkδij\dfrac{1}{3}\varepsilon_{kk}\delta_{ij}, ο οποίος και ονομάζεται και ογκομετρικός τανυστής των παραμορφώσεων, ή vol(𝜺)\text{vol}(\boldsymbol{\varepsilon}),

  • και ενός συμμετρικού τανυστή με μηδενικό ίχνος, εij-13εkkδij\varepsilon_{ij}-\dfrac{1}{3}\varepsilon_{kk}\delta_{ij}, ο οποίος και ονομάζεται και αποκλίνων τανυστής των παραμορφώσεων, ή dev(𝜺)\text{dev}(\boldsymbol{\varepsilon}).

Η (1.155) τώρα μπορεί να γραφτεί και ως:

𝜺\displaystyle\boldsymbol{\varepsilon} =vol(𝜺)+dev(𝜺)\displaystyle=\mathrm{vol}(\boldsymbol{\varepsilon})+\mathrm{dev}(\boldsymbol{% \varepsilon}) (1.157)
=13tr(𝜺)𝐈+𝜺-(vol(𝜺))\displaystyle=\dfrac{1}{3}\mathrm{tr}(\boldsymbol{\varepsilon})\mathbf{I}+% \boldsymbol{\varepsilon}-\left(\mathrm{vol}(\boldsymbol{\varepsilon})\right) (1.158)
=[13tr(𝜺)00013tr(𝜺)00013tr(𝜺)]+[ε11-13tr(𝜺)ε12ε13ε21ε22-13tr(𝜺)ε23ε31ε32ε33-13tr(𝜺)]\displaystyle=\left[\begin{array}[]{ccc}\frac{1}{3}\mathrm{tr}(\boldsymbol{% \varepsilon})&0&0\\ 0&\frac{1}{3}\mathrm{tr}(\boldsymbol{\varepsilon})&0\\ 0&0&\frac{1}{3}\mathrm{tr}(\boldsymbol{\varepsilon})\\ \end{array}\right]+\left[\begin{array}[]{ccc}\varepsilon_{11}-\frac{1}{3}% \mathrm{tr}(\boldsymbol{\varepsilon})&\varepsilon_{12}&\varepsilon_{13}\\ \varepsilon_{21}&\varepsilon_{22}-\frac{1}{3}\mathrm{tr}(\boldsymbol{% \varepsilon})&\varepsilon_{23}\\ \varepsilon_{31}&\varepsilon_{32}&\varepsilon_{33}-\frac{1}{3}\mathrm{tr}(% \boldsymbol{\varepsilon})\\ \end{array}\right] (1.159)

Στην πιο γενική του μορφή ο Νόμος του Hooke ορίζει για τα ελαστικά, ισότροπα υλικά τη σχέση των τάσεων παραμορφώσεων ως έναν γραμμικό συνδυασμό του ογκομετρικού και του αποκλίνοντα τανυστή των παραμορφώσεων, σύμφωνα με τη σχέση:

σij=3K(13εkkδij)+2G(εij-13εkkδij)\sigma_{ij}=3K\left(\dfrac{1}{3}\varepsilon_{kk}\delta_{ij}\right)+2G\left(% \varepsilon_{ij}-\dfrac{1}{3}\varepsilon_{kk}\delta_{ij}\right) (1.160)

όπου:

  • KK

    είναι το μέτρο διόγκωσης και

  • GG

    το μέτρο διάτμησης.

Χρησιμοποιώντας τις σχέσεις μεταξύ των ελαστικών σταθερών, όπως αυτές παρουσιάζονται στον Πίνακα 1.1, η (1.160) μπορεί να γραφτεί με διάφορους τρόπους.

Χρησιμοποιώντας τις σταθερές του Lamé, λ\lambda και μ\mu, η (1.160) δίνεται ως:

𝝈=λtr(𝜺)𝐈+2μ𝜺𝝈𝐂𝜺\boldsymbol{\sigma}=\lambda\mathrm{tr}(\boldsymbol{\varepsilon})\mathbf{I}+2% \mu\boldsymbol{\varepsilon}\quad\boldsymbol{\sigma}\mathbf{C}\boldsymbol{\varepsilon} (1.161)

όπου 𝐈\mathbf{I} το μοναδιαίο μητρώο (3×3)(3\times 3) και 𝐂\mathbf{C} ο τανυστής της ελαστικότητας,

𝐂=λ𝐈𝐈+2μ𝖨.\mathbf{C}=\lambda\mathbf{I}\otimes\mathbf{I}+2\mu\mathsf{I}. (1.162)

Στην παραπάνω σχέση το 𝖨\mathsf{I} είναι το συμμετρικό τμήμα του μοναδιαίου τανυστή τετάρτης τάξης. Σε μορφή δεικτών η (1.161) δίνεται ως:

σij\displaystyle\sigma_{ij} =λεkkδij+2μεij\displaystyle=\lambda\varepsilon_{kk}\delta_{ij}+2\mu\varepsilon_{ij}
=Cijklεkl\displaystyle=C_{ijkl}~{}\varepsilon_{kl} (1.163)

όπου

Cijkl=λδijδkl+μ(δikδjl+δilδjk)C_{ijkl}=\lambda\delta_{ij}\delta_{kl}+\mu(\delta_{ik}\delta_{jl}+\delta_{il}~% {}\delta_{jk}) (1.164)

Χρησιμοποιώντας το μέτρο του Young EE και τον λόγο του Poisson ν\nu, ο Νόμος του Hooke μπορεί να γραφτεί σε μορφή δεικτών:

εij=1E(σij-ν[σkkδij-σij])\varepsilon_{ij}=\dfrac{1}{E}\left(\sigma_{ij}-\nu\left[\sigma_{kk}\delta_{ij}% -\sigma_{ij}\right]\right) (1.165)

η οποία μπορεί να αναπτυχθεί ως:

ε11\displaystyle\varepsilon_{11} =1E[σ11-ν(σ22+σ33)]\displaystyle=\dfrac{1}{E}\left[\sigma_{11}-\nu(\sigma_{22}+\sigma_{33})\right] (1.166)
ε22\displaystyle\varepsilon_{22} =1E[σ22-ν(σ11+σ33)]\displaystyle=\dfrac{1}{E}\left[\sigma_{22}-\nu(\sigma_{11}+\sigma_{33})\right] (1.167)
ε33\displaystyle\varepsilon_{33} =1E[σ33-ν(σ11+σ22)]\displaystyle=\dfrac{1}{E}\left[\sigma_{33}-\nu(\sigma_{11}+\sigma_{22})\right] (1.168)
ε12\displaystyle\varepsilon_{12} =12Gσ12\displaystyle=\dfrac{1}{2G}\sigma_{12} (1.169)
ε13\displaystyle\varepsilon_{13} =12Gσ13\displaystyle=\dfrac{1}{2G}\sigma_{13} (1.170)
ε23\displaystyle\varepsilon_{23} =12Gσ23\displaystyle=\dfrac{1}{2G}\sigma_{23} (1.171)

Σε μητρωική μορφή, ο Νόμος του Hooke για ισότροπα υλικά δίνεται ως:

[ε11ε22ε332ε232ε132ε12]=1E[1-ν-ν000-ν1-ν000-ν-ν10000002(1+ν)0000002(1+ν)0000002(1+ν)][σ11σ22σ33σ23σ13σ12]\left[\begin{array}[]{c}\varepsilon_{11}\\ \varepsilon_{22}\\ \varepsilon_{33}\\ 2\varepsilon_{23}\\ 2\varepsilon_{13}\\ 2\varepsilon_{12}\\ \end{array}\right]=\dfrac{1}{E}\left[\begin{array}[]{cccccc}1&-\nu&-\nu&0&0&0% \\ -\nu&1&-\nu&0&0&0\\ -\nu&-\nu&1&0&0&0\\ 0&0&0&2(1+\nu)&0&0\\ 0&0&0&0&2(1+\nu)&0\\ 0&0&0&0&0&2(1+\nu)\\ \end{array}\right]\left[\begin{array}[]{c}\sigma_{11}\\ \sigma_{22}\\ \sigma_{33}\\ \sigma_{23}\\ \sigma_{13}\\ \sigma_{12}\\ \end{array}\right] (1.172)

όπου, σύμφωνα και με όσα προαναφέρθηκαν:

γij:=2εij\gamma_{ij}:=2\varepsilon_{ij} (1.173)

Η αντίστροφη σχέση τώρα δίνεται ως:

[σ11σ22σ33σ23σ13σ12]=E(1+ν)(1-2ν)[1-ννν000ν1-νν000νν1-ν000000(1-2ν)/2000000(1-2ν)/2000000(1-2ν)/2][ε11ε22ε332ε232ε132ε12]\left[\begin{array}[]{c}\sigma_{11}\\ \sigma_{22}\\ \sigma_{33}\\ \sigma_{23}\\ \sigma_{13}\\ \sigma_{12}\\ \end{array}\right]=\dfrac{E}{(1+\nu)(1-2\nu)}\left[\begin{array}[]{cccccc}1-% \nu&\nu&\nu&0&0&0\\ \nu&1-\nu&\nu&0&0&0\\ \nu&\nu&1-\nu&0&0&0\\ 0&0&0&(1-2\nu)/2&0&0\\ 0&0&0&0&(1-2\nu)/2&0\\ 0&0&0&0&0&(1-2\nu)/2\\ \end{array}\right]\left[\begin{array}[]{c}\varepsilon_{11}\\ \varepsilon_{22}\\ \varepsilon_{33}\\ 2\varepsilon_{23}\\ 2\varepsilon_{13}\\ 2\varepsilon_{12}\\ \end{array}\right] (1.174)

ή με τη χρήση των σταθερών του Lamé:

[σ11σ22σ33σ23σ13σ12]=1E[2μ+λλλ000λ2μ+λλ000λλ2μ+λ000000μ000000μ000000μ][ε11ε22ε332ε232ε132ε12]\left[\begin{array}[]{c}\sigma_{11}\\ \sigma_{22}\\ \sigma_{33}\\ \sigma_{23}\\ \sigma_{13}\\ \sigma_{12}\\ \end{array}\right]=\dfrac{1}{E}\left[\begin{array}[]{cccccc}2\mu+\lambda&% \lambda&\lambda&0&0&0\\ \lambda&2\mu+\lambda&\lambda&0&0&0\\ \lambda&\lambda&2\mu+\lambda&0&0&0\\ 0&0&0&\mu&0&0\\ 0&0&0&0&\mu&0\\ 0&0&0&0&0&\mu\\ \end{array}\right]\left[\begin{array}[]{c}\varepsilon_{11}\\ \varepsilon_{22}\\ \varepsilon_{33}\\ 2\varepsilon_{23}\\ 2\varepsilon_{13}\\ 2\varepsilon_{12}\\ \end{array}\right] (1.175)
EE ν\nu λ\lambda GG (μ\mu) KK
EE, ν\nu Eν(1+ν)(1-2ν)\dfrac{E\nu}{(1+\nu)(1-2\nu)} E2(1+ν)\dfrac{E}{2(1+\nu)} E3(1-2ν)\dfrac{E}{3(1-2\nu)}
EE, GG E2G-1\dfrac{E}{2G}-1 G(E-2G)3G-E\dfrac{G(E-2G)}{3G-E} EG3(3G-E)\dfrac{EG}{3(3G-E)}
EE, KK 3K-E6K\dfrac{3K-E}{6K} 3K(3K-E)9K-E\dfrac{3K(3K-E)}{9K-E} 3KE9K-E\dfrac{3KE}{9K-E}
ν\nu, GG 2G(1+ν)2G(1+\nu) 2Gν1-2ν\dfrac{2G\nu}{1-2\nu} 2G(1+ν)3(1-2ν)\dfrac{2G(1+\nu)}{3(1-2\nu)}
ν\nu, KK 3K(1-2ν)3K(1-2\nu) 3Kν1+ν\dfrac{3K\nu}{1+\nu} 3K(1-2ν)2(1+ν)\dfrac{3K(1-2\nu)}{2(1+\nu)}
λ\lambda, μ\mu λ2(λ+μ)\dfrac{\lambda}{2(\lambda+\mu)} λ+2μ\lambda+2\mu λ+2μ3\lambda+\dfrac{2\mu}{3}
Πίνακας 1.1: Κύριες σχέσεις μεταξύ των ελαστικών σταθερών.

1.4 Επίπεδη παραμόρφωση

H εντατική κατάσταση της επίπεδης παραμόρφωσης αποτελεί μια αρκετά συνηθισμένη εντατική κατάσταση σε προβλήματα γεωτεχνικής μηχανικής. Αναφέρεται σε περιπτώσεις κατά τις οποίες οι παραμορφώσεις οι σχετικές με ένα επίπεδο, κατά σύμβαση το επίπεδο κάθετα στον άξονα zz, είναι μηδενικές ή σχεδόν μηδενικές [7], ισχύει δηλαδή:

εzzγyzγzx0\varepsilon_{zz}\equiv\gamma_{yz}\equiv\gamma_{zx}\approx 0 (1.176)

Συνθήκες επίπεδης παραμόρφωσης θεωρούμε ότι επικρατούν σε σώματα στα οποία δεν μεταβάλλεται ούτε η γεωμετρία αλλά ούτε και η φόρτιση κατά μήκος του άξονα zz, όπως φαίνεται στο Σχήμα 1.11, και:

  1. 1.

    Υπάρχουν οι κατάλληλοι περιορισμοί που επιβάλλουν τη συνθήκη που περιγράφεται από την (1.176).

  2. 2.

    Το σώμα είναι αρκετά επίμηκες (LzL_{z}\rightarrow\infty) ούτως ώστε να θεωρήσουμε ότι ισχύει η (1.176).

Σχήμα 1.11: Περίπτωση επίπεδης παραμόρφωσης.

Στη δεύτερη περίπτωση ανήκουν πολλά προβλήματα της γεωτεχνικής μηχανικής, όπως φράγματα, επιχώματα, ορύγματα, τοίχοι αντιστήριξης κ.ά. (βλ. και Σχήμα 1.12).

Σχήμα 1.12: Τυπικά προβλήματα γεωτεχνικής μηχανικής στα οποία επικρατούν συνθήκες επίπεδης παραμόρφωσης.

Οι σχέσεις τάσεων παραμορφώσεων προκύπτουν από το τρισδιάστατο μητρώο ελαστικότητας μηδενίζοντας τους αντίστοιχους όρους που ορίζονται στην (1.176). Για την περίπτωση της ισότροπης ελαστικότητας π.χ. ισχύει:

[σ11σ22σ33σ23σ13σ12]=E(1+ν)(1-2ν)[1-ννν000ν1-νν000νν1-ν000000(1-2ν)/2000000(1-2ν)/2000000(1-2ν)/2][ε11ε22ε332ε232ε132ε12]\left[\begin{array}[]{c}\sigma_{11}\\ \sigma_{22}\\ \sigma_{33}\\ \sigma_{23}\\ \sigma_{13}\\ \sigma_{12}\\ \end{array}\right]=\dfrac{E}{(1+\nu)(1-2\nu)}\left[\begin{array}[]{cccccc}1-% \nu&\nu&\nu&0&0&0\\ \nu&1-\nu&\nu&0&0&0\\ \nu&\nu&1-\nu&0&0&0\\ 0&0&0&(1-2\nu)/2&0&0\\ 0&0&0&0&(1-2\nu)/2&0\\ 0&0&0&0&0&(1-2\nu)/2\\ \end{array}\right]\left[\begin{array}[]{c}\varepsilon_{11}\\ \varepsilon_{22}\\ \varepsilon_{33}\\ 2\varepsilon_{23}\\ 2\varepsilon_{13}\\ 2\varepsilon_{12}\\ \end{array}\right] (1.177)

και αντικαθιστώντας την (1.176) προκύπτει:

[σ11σ22σ33σ23σ13σ12]=E(1+ν)(1-2ν)[1-ννν000ν1-νν000νν1-ν000000(1-2ν)/2000000(1-2ν)/2000000(1-2ν)/2][ε11ε220002ε12]\left[\begin{array}[]{c}\sigma_{11}\\ \sigma_{22}\\ \sigma_{33}\\ \sigma_{23}\\ \sigma_{13}\\ \sigma_{12}\\ \end{array}\right]=\dfrac{E}{(1+\nu)(1-2\nu)}\left[\begin{array}[]{cccccc}1-% \nu&\nu&\nu&0&0&0\\ \nu&1-\nu&\nu&0&0&0\\ \nu&\nu&1-\nu&0&0&0\\ 0&0&0&(1-2\nu)/2&0&0\\ 0&0&0&0&(1-2\nu)/2&0\\ 0&0&0&0&0&(1-2\nu)/2\\ \end{array}\right]\left[\begin{array}[]{c}\varepsilon_{11}\\ \varepsilon_{22}\\ 0\\ 0\\ 0\\ 2\varepsilon_{12}\\ \end{array}\right] (1.178)

Από την παραπάνω σχέση εύκολα προκύπτει ότι

σ23\displaystyle\sigma_{23} =0\displaystyle=0 (1.179)
σ13\displaystyle\sigma_{13} =0\displaystyle=0 (1.180)

και

σ33=ν(ε11+ε22).\sigma_{33}=\nu(\varepsilon_{11}+\varepsilon_{22}). (1.181)

Επομένως οι σχέσεις τάσεων παραμορφώσεων μπορούν πλέον να περιγραφούν από την ακόλουθη σχέση:

[σ11σ22σ12]=E(1+ν)(1-2ν)[1-νν0ν1-ν000(1-2ν)/2][ε11ε222ε12]\left[\begin{array}[]{c}\sigma_{11}\\ \sigma_{22}\\ \sigma_{12}\\ \end{array}\right]=\dfrac{E}{(1+\nu)(1-2\nu)}\left[\begin{array}[]{cccccc}1-% \nu&\nu&0\\ \nu&1-\nu&0\\ 0&0&(1-2\nu)/2\\ \end{array}\right]\left[\begin{array}[]{c}\varepsilon_{11}\\ \varepsilon_{22}\\ 2\varepsilon_{12}\\ \end{array}\right] (1.182)

1.5 Αξονοσυμμετρία

Αντίστοιχη εντατική κατάσταση είναι και η κατάσταση της αξονοσυμμετρίας που αφορά σώματα που εμφανίζουν περιστροφική συμμετρία, τόσο όσον αφορά τη γεωμετρία αλλά και όσον αφορά και τη φόρτιση γύρω από έναν άξονα [4] (βλ. π.χ. Σχήμα 1.13.)

Σχήμα 1.13: Περίπτωση αξονοσυμμετρίας.

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

1.5.1 Το σύστημα συντεταγμένων

Για τον προσδιορισμό των καταστατικών εξισώσεων του αξονοσυμμετρικού προβλήματος, χρησιμοποιείται ένα κυλινδρικό σύστημα συντεταγμένων (r,z,θ)(r,z,\theta), όπου:

  • rr

    η ακτινική συντεταγμένη,

  • zz

    η συντεταγμένη κατά τον άξονα συμμετρίας και

  • θ\theta

    η περιμετρική συντεταγμένη.

Σημειώνεται ότι το (r,z)(r,z) αποτελεί ένα δεξιόστροφο καρτεσιανό σύστημα συντεταγμένων για τα επίπεδα όπου θ=σταθερό\theta=\text{σταθερό}, ενώ για τα επίπεδα όπου ισχύει z=constz=const σχηματίζεται το πολικό σύστημα συντεταγμένων (r,θ)(r,\theta).

Το παραπάνω σύστημα συντεταγμένων δίνεται στο Σχήμα 1.5.

1.5.2 Μετακινήσεις, παραμορφώσεις και τάσεις

Το πεδίο μετακινήσεων είναι μια συνάρτηση μόνο της ακτινικής συντεταγμένης rr και της συντεταγμένης κατά τον άξονα συμμετρίας zz και ορίζεται ως:

𝒖(r,z)=[ur(r,z)uz(r,z)]\boldsymbol{u}(r,z)=\left[\begin{array}[]{c}u_{r}(r,z)\\ u_{z}(r,z)\end{array}\right] (1.183)

Ο τανυστής των απειροστών παραμορφώσεων στο κυλινδρικό σύστημα συντεταγμένων δίνεται από το 3×33\times 3:

𝜺(r,z)=[εrrεrzεrθεzrεzzεzθεrθεzθεθθ]\boldsymbol{\varepsilon}(r,z)=\left[\begin{array}[]{ccc}\varepsilon_{rr}&% \varepsilon_{rz}&\varepsilon_{r\theta}\\ \varepsilon_{zr}&\varepsilon_{zz}&\varepsilon_{z\theta}\\ \varepsilon_{r\theta}&\varepsilon_{z\theta}&\varepsilon_{\theta\theta}\\ \end{array}\right] (1.184)

Το παραπάνω μητρώο είναι συμμετρικό, ισχύει δηλαδή ότι:

εij=εjiγιαi,j=1,2,3.\varepsilon_{ij}=\varepsilon_{ji}\quad\text{για}\quad i,j=1,2,3. (1.185)

Λόγω των παραδοχών της αξονοσυμμετρικής κατάστασης, ισχύει:

εrθ\displaystyle\varepsilon_{r\theta} =0\displaystyle=0 (1.186)
εzθ\displaystyle\varepsilon_{z\theta} =0\displaystyle=0 (1.187)

και επομένως ο τανυστής των παραμορφώσεων ξαναγράφεται ως:

𝜺(r,z)=[εrrεrz0εzrεzzεzθ0εzθεθθ]\boldsymbol{\varepsilon}(r,z)=\left[\begin{array}[]{ccc}\varepsilon_{rr}&% \varepsilon_{rz}&0\\ \varepsilon_{zr}&\varepsilon_{zz}&\varepsilon_{z\theta}\\ 0&\varepsilon_{z\theta}&\varepsilon_{\theta\theta}\\ \end{array}\right] (1.188)

Εύκολα διαπιστώνεται ότι οι όροι του μητρώου 𝜺(r,z)\boldsymbol{\varepsilon}(r,z) είναι συναρτήσεις μόνο των συντεταγμένων rr και zz.

Όπως και στις προηγούμενες ενότητες, παραλείποντας τους μηδενικούς όρους το παραπάνω μητρώο επαναδιατυπώνεται ως ένα διάνυσμα 4×14\times 1:

𝜺(r,z)=[εrrεzzεθθεrz]\boldsymbol{\varepsilon}(r,z)=\left[\begin{array}[]{c}\varepsilon_{rr}\\ \varepsilon_{zz}\\ \varepsilon_{\theta\theta}\\ \varepsilon_{rz}\\ \end{array}\right] (1.189)

Συχνά στη βιβλιογραφία συναντούμε το παραπάνω διάνυσμα και ως:

𝜺(r,z)=[εrrεzzεθθγrz]\boldsymbol{\varepsilon}(r,z)=\left[\begin{array}[]{c}\varepsilon_{rr}\\ \varepsilon_{zz}\\ \varepsilon_{\theta\theta}\\ \gamma_{rz}\\ \end{array}\right] (1.190)

όπου

γrz=εrz+εzr=2εrz.\gamma_{rz}=\varepsilon_{rz}+\varepsilon_{zr}=2\varepsilon_{rz}. (1.191)

Στην περίπτωση αυτή χρήζει ιδιαίτερης προσοχής ο ορισμός του μητρώου ελαστικότητας.

Ο τανυστής των τάσεων σε κυλινδρικές συντεταγμένες δίνεται αντίστοιχα από το συμμετρικό μητρώο:

𝝈(r,z)=[σrrσrzσrθσzrσzzσzθσrθσzθσθθ]\boldsymbol{\sigma}(r,z)=\left[\begin{array}[]{ccc}\sigma_{rr}&\sigma_{rz}&% \sigma_{r\theta}\\ \sigma_{zr}&\sigma_{zz}&\sigma_{z\theta}\\ \sigma_{r\theta}&\sigma_{z\theta}&\sigma_{\theta\theta}\\ \end{array}\right] (1.192)

Και σε αυτή την περίπτωση, λόγω των παραδοχών της αξονοσυμμετρικής κατάστασης, ισχύει:

σrθ\displaystyle\sigma_{r\theta} =0\displaystyle=0 (1.193)
σzθ\displaystyle\sigma_{z\theta} =0\displaystyle=0 (1.194)

και επομένως ο τανυστής των τάσεων ξαναγράφεται ως:

𝝈(r,z)=[σrrσrz0σzrσzzσzθ0σzθσθθ]\boldsymbol{\sigma}(r,z)=\left[\begin{array}[]{ccc}\sigma_{rr}&\sigma_{rz}&0\\ \sigma_{zr}&\sigma_{zz}&\sigma_{z\theta}\\ 0&\sigma_{z\theta}&\sigma_{\theta\theta}\\ \end{array}\right] (1.195)

Αντίστοιχα με τον τανυστή των παραμορφώσεων, ο τανυστής των τάσεων επαναδιατυπώνεται ως ένα διάνυσμα 4×14\times 1, του οποίου οι όροι είναι συναρτήσεις μόνο των συντεταγμένων rr και zz:

𝝈(r,z)=[σrrσzzσθθσrz]\boldsymbol{\sigma}(r,z)=\left[\begin{array}[]{c}\sigma_{rr}\\ \sigma_{zz}\\ \sigma_{\theta\theta}\\ \sigma_{rz}\\ \end{array}\right] (1.196)

1.5.3 Κινηματικές εξισώσεις

Συνεχίζουμε με τον προσδιορισμό των κινηματικών εξισώσεων την αξονοσυμμετρική κατάσταση. Οι σχέσεις μετακινήσεων-παραμορφώσεων για το πρόβλημα δίνονται ως:

εrr\displaystyle\varepsilon_{rr} =urr\displaystyle=\dfrac{\partial u_{r}}{\partial r} (1.197)
εzz\displaystyle\varepsilon_{zz} =uzz\displaystyle=\dfrac{\partial u_{z}}{\partial z} (1.198)
εθθ\displaystyle\varepsilon_{\theta\theta} =urr\displaystyle=\dfrac{u_{r}}{r} (1.199)
εrz\displaystyle\varepsilon_{rz} =12(urz+uzr)\displaystyle=\dfrac{1}{2}\left(\dfrac{\partial u_{r}}{\partial z}+\dfrac{% \partial u_{z}}{\partial r}\right) (1.200)

Οι παραπάνω εξισώσεις δίνονται σε μητρωική μορφή ως:

𝜺𝒓𝒓=[εrrεzzεθθεrz]=[r00z1r0zr][εrεz]\boldsymbol{\varepsilon_{rr}}=\left[\begin{array}[]{c}\varepsilon_{rr}\\ \varepsilon_{zz}\\ \varepsilon_{\theta\theta}\\ \varepsilon_{rz}\\ \end{array}\right]=\left[\begin{array}[]{cc}\dfrac{\partial}{\partial r}&0\\ 0&\dfrac{\partial}{\partial z}\\ \dfrac{1}{r}&0\\ \dfrac{\partial}{\partial z}&\dfrac{\partial}{\partial r}\\ \end{array}\right]\left[\begin{array}[]{c}\varepsilon_{r}\\ \varepsilon_{z}\end{array}\right] (1.201)

ή αλλιώς σε συνεπτυγμένη μορφή,

𝜺=𝐁𝒖\boldsymbol{\varepsilon}=\mathbf{B}\boldsymbol{u} (1.202)

όπου 𝐁\mathbf{B} το μητρώο 4×24\times 2 που ορίζεται ως

𝐁=[r00z1r0zr][εrεz]\mathbf{B}=\left[\begin{array}[]{cc}\dfrac{\partial}{\partial r}&0\\ 0&\dfrac{\partial}{\partial z}\\ \dfrac{1}{r}&0\\ \dfrac{\partial}{\partial z}&\dfrac{\partial}{\partial r}\\ \end{array}\right]\left[\begin{array}[]{c}\varepsilon_{r}\\ \varepsilon_{z}\end{array}\right] (1.203)

είναι ένας διαφορικός τελεστής που συνδέει τις μετακινήσεις με τις παραμορφώσεις.

1.5.4 Καταστατικές εξισώσεις

Οι καταστατικές εξισώσεις, οι σχέσεις τάσεων-παραμορφώσεων δηλαδή, δίνονται από τη σχέση

𝜺𝒓𝒓=[εrrεzzεθθεrz]=[C11C21C31C41C12C22C32C42C13C23C33C43C14C24C34C44][εrrεzzεθθεrz]\boldsymbol{\varepsilon_{rr}}=\left[\begin{array}[]{c}\varepsilon_{rr}\\ \varepsilon_{zz}\\ \varepsilon_{\theta\theta}\\ \varepsilon_{rz}\\ \end{array}\right]=\left[\begin{array}[]{cccc}C_{11}&C_{21}&C_{31}&C_{41}\\ C_{12}&C_{22}&C_{32}&C_{42}\\ C_{13}&C_{23}&C_{33}&C_{43}\\ C_{14}&C_{24}&C_{34}&C_{44}\\ \end{array}\right]\left[\begin{array}[]{c}\varepsilon_{rr}\\ \varepsilon_{zz}\\ \varepsilon_{\theta\theta}\\ \varepsilon_{rz}\\ \end{array}\right] (1.204)

ή αλλιώς σε συνεπτυγμένη μορφή,

𝜺=𝐂𝝈\boldsymbol{\varepsilon}=\mathbf{C}\boldsymbol{\sigma} (1.205)

όπου 𝐂\mathbf{C} το μητρώο του τανυστή της ελαστικότητας

𝐂=[C11C21C31C41C12C22C32C42C13C23C33C43C14C24C34C44].\mathbf{C}=\left[\begin{array}[]{cccc}C_{11}&C_{21}&C_{31}&C_{41}\\ C_{12}&C_{22}&C_{32}&C_{42}\\ C_{13}&C_{23}&C_{33}&C_{43}\\ C_{14}&C_{24}&C_{34}&C_{44}\\ \end{array}\right]. (1.206)

Λόγω της παραδοχής της αξονοσυμμετρίας, οι περιμετρικές τάσεις σθθ\sigma_{\theta\theta} δεν θα πρέπει να συνδέονται με τις διατμητικές παραμορφώσεις εrz\varepsilon_{rz}. Θα πρέπει δηλαδή να ισχύει

σθθ=C13εrr+C23εzz+C33εθθ\sigma_{\theta\theta}=C_{13}\varepsilon_{rr}+C_{23}\varepsilon_{zz}+C_{33}% \varepsilon_{\theta\theta} (1.207)

μόνο.

Επομένως μπορούμε να ξαναγράψουμε το μητρώο του τανυστή της ελαστικότητας ως:

𝐂=[C11C21C31C41C12C22C32C42C13C23C330C14C240C44]\mathbf{C}=\left[\begin{array}[]{cccc}C_{11}&C_{21}&C_{31}&C_{41}\\ C_{12}&C_{22}&C_{32}&C_{42}\\ C_{13}&C_{23}&C_{33}&0\\ C_{14}&C_{24}&0&C_{44}\\ \end{array}\right] (1.208)

γιατί εξαιτίας της συμμετρίας θα πρέπει να ισχύει:

C43=C34=0C_{43}=C_{34}=0 (1.209)

Θεωρώντας τώρα ένα ομογενές και ισότροπο υλικό με μέτρο του Young EE και λόγο του Poisson ν\nu, το μητρώο του τανυστή της ελαστικότητας για την περίπτωση της αξονοσυμμετρίας δίνεται ως:

𝐂=E(1+ν)(1-ν)[1-ννν0ν1-νν0νν1-ν000012(1-2ν)]\mathbf{C}=\dfrac{E}{(1+\nu)(1-\nu)}\left[\begin{array}[]{cccc}1-\nu&\nu&\nu&0% \\ \nu&1-\nu&\nu&0\\ \nu&\nu&1-\nu&0\\ 0&0&0&\dfrac{1}{2}(1-2\nu)\\ \end{array}\right] (1.210)

1.5.5 Εξισώσεις ισορροπίας

Τέλος δίνονται οι γενικές (τρισδιάστατες) μερικές διαφορικές εξισώσεις ισορροπίας σε κυλινδρικές συντεταγμένες, οι οποίες ορίζονται ως:

1rrσrr+1rθσrθ+z(rσz)-σθθr+br\displaystyle\dfrac{1}{r}\dfrac{\partial}{\partial r}\sigma_{rr}+\dfrac{1}{r}% \dfrac{\partial}{\partial\theta}\sigma_{r\theta}+\dfrac{\partial}{\partial z}% \left(r\sigma_{z}\right)-\dfrac{\sigma_{\theta\theta}}{r}+b_{r} =0\displaystyle=0 (1.211)
1rr(rσzr)+1rσzθ+zσzz+bz\displaystyle\dfrac{1}{r}\dfrac{\partial}{\partial r}\left(r\sigma_{zr}\right)% +\dfrac{1}{r}\sigma_{z\theta}+\dfrac{\partial}{\partial z}\sigma_{zz}+b_{z} =0\displaystyle=0 (1.212)
1r2r(r2σθr)+1rzσθθ+zσθz+bθ\displaystyle\dfrac{1}{r^{2}}\dfrac{\partial}{\partial r}\left(r^{2}\sigma_{% \theta r}\right)+\dfrac{1}{r}\dfrac{\partial}{\partial z}\sigma_{\theta\theta}% +\dfrac{\partial}{\partial z}\sigma_{\theta z}+b_{\theta} =0\displaystyle=0 (1.213)

όπου brb_{r}, bzb_{z}, bθb_{\theta}, τα στοιχεία του διανύσματος 𝒃\boldsymbol{b} των μαζικών δυνάμεων

𝒃=[brbzbθ]\boldsymbol{b}=\left[\begin{array}[]{c}b_{r}\\ b_{z}\\ b_{\theta}\\ \end{array}\right] (1.214)

κατά τις διευθύνσεις rr, zz και θ\theta αντίστοιχα.

Για το αξονοσυμμετρικό πρόβλημα οι παραπάνω εξισώσεις (1.211) γράφονται ως:

1rrσrr+z(rσz)-σθθr+br\displaystyle\dfrac{1}{r}\dfrac{\partial}{\partial r}\sigma_{rr}+\dfrac{% \partial}{\partial z}\left(r\sigma_{z}\right)-\dfrac{\sigma_{\theta\theta}}{r}% +b_{r} =0\displaystyle=0 (1.215)
1rr(rσzr)+zσzz+bz\displaystyle\dfrac{1}{r}\dfrac{\partial}{\partial r}\left(r\sigma_{zr}\right)% +\dfrac{\partial}{\partial z}\sigma_{zz}+b_{z} =0\displaystyle=0 (1.216)

Η τρίτη εξίσωση από τη σχέση (1.211) ισχύει ταυτοτικά αν bθ=0b_{\theta}=0. Αν bθ0b_{\theta}\neq 0, το πρόβλημα δεν μπορεί να θεωρηθεί αξονοσυμμετρικό.

1.6 Βιβλιογραφία

  • 1 R. D. Cook, D. S. Malkus, M. E. Plesha and R. J. Witt. Concepts and applications of finite element analysis. John Wiley & Sons, 1989.
  • 2 Debasis Deb. Finite element methods: concepts and applications in geomechanics. PHI Learning Pvt. Ltd., 2006.
  • 3 Haojiang Ding, Weiqiu Chen and Liangchi Zhang. Elasticity of transversely isotropic materials, volume 126. Springer Science & Business Media, 2006.
  • 4 Carlos A. Felippa. Lecture notes in advanced finite element methods, 2013.
  • 5 J. C. Jaeger, N. G. W. Cook and R. W. Zimmerman. Fundamentals of rock mechanics. Wiley-Blackwell, 2007.
  • 6 Robert M. Jones. Mechanics of composite materials, volume 1. McGraw-Hill New York, 1975.
  • 7 Marc A. Meyers and Krishan Kumar Chawla. Mechanical behavior of materials, volume 2. Cambridge University Press, 2009.
  • 8 Graeme W. Milton. The theory of composites, volume 6. Cambridge University Press, 2002.
  • 9 Eugenio Oñate and Roger Owen. Computational plasticity, volume 7. Springer, 2007.
  • 10 Ivan Stephen Sokolnikoff and Robert Dickerson Specht. Mathematical theory of elasticity. McGraw-Hill New York, 1956.
  • 11 Μ. Παπαδρακάκης. Ανάλυση φορέων με τη μέθοδο των πεπερασμένων στοιχείων. Παπασωτηρίου, 2001.