1 ΣΤΟΙΧΕΙΑ ΑΠΟ ΤΗ ΓΡΑΜΜΙΚΗ ΑΛΓΕΒΡΑ ΚΑΙ ΑΝΑΛΥΣΗ

Μηχανική του Συνεχούς Μέσου, Κεφ. 1., Ι. Βαρδουλάκης 2008
5
1 ΣΤΟΙΧΕΙΑ ΑΠΟ ΤΗ ΓΡΑΜΜΙΚΗ ΑΛΓΕΒΡΑ ΚΑΙ
ΑΝΑΛΥΣΗ
1
ΣΤΟΙΧΕΙΑ ΑΠΟ ΤΗ ΓΡΑΜΜΙΚΗ ΑΛΓΕΒΡΑ ΚΑΙ ΑΝΑΛΥΣΗ
1.1
Συµβολισµοί µε τη χρήση δεικτών
1.2
Συµµετρικά και αντισυµµετρικά συστήµατα
1.3
Συµµετρικό και αντισυµµετρικό µέρος συστήµατος 2ας τάξεως
1.4
Το σύµβολο Levi-Civita και τα σύµβολα Kronecker
1.5
Συµµετρίες συστηµάτων 3ης τάξεως
1.6
Ορίζουσες
1.7
Γραµµικά συστήµατα
1.8
Θετικώς ορισµένες τετραγωνικές µορφές
1.9
Η χαρακτηριστική εξίσωση τετραγωνικού πίνακα 3 × 3
1.10 Ο τελεστής " ∇ "
1.10.1 Ο τελεστής βαθµίδας
1.10.2 Ο τελεστής στροβιλισµού
1.10.3 Ο τελεστής αποκλίσεως
1.11 Παράρτηµα I: Χρήσιµες ταυτότητες µεταξύ διαφορικών τελεστών
1.12 Παράρτηµα ΙΙ: Το θεώρηµα αποκλίσεως
5
7
9
10
11
15
16
18
18
20
22
22
23
25
27
28
Στο Κεφάλαιο αυτό εισάγουµε τη σηµειολογία µε χρήση δεικτών και τη χρήση της
αθροιστικής σύµβασης κατά Gauß-Einstein. Ειδικότερα παρουσιάζεται µία σύνοψη
ορισµένων βασικών ορισµών και θεωρηµάτων από την περιοχή της Γραµµικής
Άλγεβρας1,2,3 και της ∆ιανυσµατικής Ανάλυσης4
1
Akivis M.A. and Goldberg V.V. An Introduction to Linear Algebra & Tensors, Dover, 1972.
Pettofrezzo A.J., Matrices and Transformations, Dover, 1966.
3
McConnell A.J., Applications of Tensor Analysis, Dover, 1957.
4
Rutherford D.E., Vector Methods, Dover, 2004.
2
6
Μηχανική του Συνεχούς Μέσου, Κεφ. 1., Ι. Βαρδουλάκης 2008
© 1. ΣΤΟΙΧΕΙΑ ΑΠΟ ΤΗ ΓΡΑΜΜΙΚΗ ΑΛΓΕΒΡΑ, 2008
Ιωάννης Γ. Βαρδουλάκης, Dr-Ing., Καθηγητής της Μηχανικής στο Ε. Μ. Πολυτεχνείο
Τ.Θ. 144, Παιανία 190-02, http://geolab.mechan.ntua.gr/, [email protected]
Μηχανική του Συνεχούς Μέσου, Κεφ. 1., Ι. Βαρδουλάκης 2008
7
1.1 Συµβολισµοί µε τη χρήση δεικτών
Συχνά θα συναντήσουµε στη βιβλιογραφία τη συντοµογραφία µιας αλγεβρικής
ποσότητας, όπως είναι π.χ. οι συντεταγµένες ενός σηµείου, όπου θα γίνεται χρήση του
λεγόµενου βωβού δείκτη. Έτσι όταν θέλουµε να αναφερθούµε στις συντεταγµένες ενός
σηµείου στο χώρο, αντί της πλήρους αναλυτικής αναγραφής αυτών,
x1 , x2 , x3
(1.1)
θα χρησιµοποιήσουµε απλά το συµβολισµό,
xi
(1.2)
αφού προηγουµένως έχουµε καθορίσει ότι ο βωβός (και εµ’ προκειµένω κάτω) δείκτης i
θα παίρνει τις τιµές 1, 2 και 3. Επειδή όµως δεν έχουµε καθορίσει ακόµα τη µαθηµατική
υπόσταση της ποσότητας xi , δεν µπορούµε να χρησιµοποιήσουµε στο σηµείο αυτό
ονοµασίες όπως «διάνυσµα» ή «τανυστής», που έχουµε ακούσει πιθανώς να αναφέρεται
σε άλλα Μαθήµατα Μηχανικής και Μαθηµατικών. Εδώ θα αρκεσθούµε στο να
αποκαλούµε το µέγεθος xi , ως ένα σύστηµα 1ης τάξεως, αφού εµφανίζει ένα µόνο
δείκτη. Επίσης θα παρατηρήσουµε ότι ο δείκτης µπορεί να αναγράφεται είτε κάτω είτε
πάνω, οπότε συχνά θα συναντήσουµε το συµβολισµό
xi
(1.3)
αποσαφηνίζοντας ότι στη περίπτωση αυτή ο άνω δείκτης δεν θα πρέπει να εκληφθεί ως
εκθέτης.
Στο σηµείο αυτό συνηθίζεται να γίνεται αναφορά στις λεγόµενες, α) γραµµικές µορφές,
που είναι αθροίσµατα της µορφής
3
L = ∑ ai xi = a1 x1 + a2 x 2 + a3 x 3
(1.4)
i =1
και β) δι-γραµµικές ή τετραγωνικές µορφές
3
3
Q = ∑∑ aij x i x j = a11 x1 x1 + a12 x1 x 2 + a13 x1 x 3 + "
(1.5)
i =1 j =1
Παράλληλα, για τη συµπύκνωση της γραφής µας, θα κάνουµε χρήση της λεγόµενης
συµβάσεως αθροίσεως πάνω σε επαναλαµβανόµενο δείκτη κατά Gauß-Einstein,
συµφώνως µε την οποία οι παραπάνω γραµµικές ή πολυ-γραµµικές µορφές θα γράφονται
χωρίς το σύµβολο του αθροίσµατος και µε προκαθορισµένο το πεδίο τιµών των δεικτών,
L = ai x i (i, j = 1, 2,3)
Q = aij x i x j (i, j = 1, 2,3)
(1.6)
Παρατηρούµε ότι η µετονοµασία ενός δείκτη π.χ. από i σε k δεν επιφέρει καµία αλλαγή
στο αποτέλεσµα της αθροίσεως,
L = ai x i = ak x k
(1.7)
Μηχανική του Συνεχούς Μέσου, Κεφ. 1., Ι. Βαρδουλάκης 2008
8
Επίσης παρατηρούµε ότι σε κάθε ένα σύστηµα 1ης τάξεως µπορούµε να αντιστοιχίσουµε
αµφιµονοσήµαντα ένα µητρώο-γραµµή ( 1× 3 ) που περιέχει τα στοιχεία του εν λόγω
συστήµατος
ai ↔ {ai } = {a1 , a2 , a3 }
(1.8)
ή ένα µητρώο-στήλη ( 3 × 1 ),
⎧ x1 ⎫
⎪ ⎪
⎪ ⎪
xi ↔ {xi }T = ⎨ x 2 ⎬
⎪ 3⎪
⎪⎩ x ⎪⎭
(1.9)
Η γραµµική µορφή υπολογίζεται από το µητρωικό γινόµενο γραµµής επί στήλη
⎧ x1 ⎫
⎪ ⎪
{a1 , a2 , a3 } ⎨ x 2 ⎬ = {ai xi } ↔ L
⎪ x3 ⎪
⎩ ⎭
(1.10)
Σε κάθε σύστηµα 2ας τάξεως µπορούµε να αντιστοιχήσουµε το µητρώο των στοιχείων
του,
⎡ a11
aij ↔ [aij ] = ⎢⎢ a21
⎢⎣ a31
a12
a22
a32
a13 ⎤
a23 ⎥⎥
a33 ⎥⎦
(1.11)
η δε τετραγωνική µορφή υπολογίζεται από το αντίστοιχο µητρωικό γινόµενο
{x , x , x }
1
2
3
⎡ a11
⎢a
⎢ 21
⎢⎣ a31
a12
a22
a32
a13 ⎤ ⎧ x1 ⎫
⎪ ⎪
a23 ⎥⎥ ⎨ x 2 ⎬ =
a33 ⎥⎦ ⎪⎩ x 3 ⎪⎭
(1.12)
⎧ x1 ⎫
⎪ ⎪
= { x i ai1 , x i ai 2 , x i ai 3 } ⎨ x 2 ⎬ = { x i aij x j } = {aij x i x j } ↔ Q
⎪ x3 ⎪
⎩ ⎭
Άσκηση
Να σχεδιασθούν οι δευτεροβάθµιες επιφάνειες στο χώρο που δίδονται από τις εξισώσεις,
Q = 1 , Εξ. (1.12), για τις κάτωθι περιπτώσεις (Εικ. 1-1):
•
Ελλειψοειδές:
⎡ a −2
⎢
[aij ] = ⎢ 0
⎢
⎢ 0
⎣
0
b −2
0
0⎤
⎥
0⎥ ⇒
⎥
c2 ⎥
⎦
x12
x2
+ 2 +
a 2 b2
x32
c2
=1
Μηχανική του Συνεχούς Μέσου, Κεφ. 1., Ι. Βαρδουλάκης 2008
•
•
Μονόφυλλο υπερβολοειδές:
⎡ a −2
⎢
[aij ] = ⎢ 0
⎢
⎢ 0
⎣
⎡ a −2
⎢
[aij ] = ⎢ 0
⎢
⎢ 0
⎣
∆ίφυλλο υπερβολοειδές:
0
b −2
0
0
−b
−2
0
0 ⎤
⎥
0 ⎥ ⇒
⎥
−c −2 ⎥
⎦
0 ⎤
⎥
0 ⎥ ⇒
⎥
−c −2 ⎥
⎦
9
x2 x2
+ 2 − 3 =1
a2 b2 c2
x12
x12
a2
−
x22
b2
−
x32
c2
=1
Παρατηρούµε ότι για Q = 0 το δεύτερο µητρώο οδηγεί στην εξίσωση ενός πραγµατικού
κώνου.
Εικ. 1-1: (α) Ελλειψοειδές. (β) Μονόφυλλο υπερβολοειδές. (γ) ∆ίφυλλο υπερβολοειδές.
1.2 Συµµετρικά και αντισυµµετρικά συστήµατα
Ένα συµµετρικό σύστηµα 2ας τάξεως aij χαρακτηρίζεται από την ιδιότητα ότι η τιµή του
δεν αλλάζει αν γίνει εναλλαγή των δεικτών,
aij = a ji
(1.13)
Παρατηρούµε ότι στη θεωρούµενη περίπτωση ο πίνακας των στοιχείων του αντίστοιχου
συστήµατος 2ας τάξεως είναι συµµετρικός ως προς την κύρια διαγώνιό του. Π.χ. ένας
συµµετρικός πίνακας 3 × 3 περιέχει µόνο 6 ανεξάρτητα στοιχεία,
⎡ a11
aij = a ji ↔ ⎢⎢ a21
⎢⎣ a31
a12
a22
a32
a13 ⎤
a23 ⎥⎥ =
a33 ⎥⎦
⎡ a11
⎢c
⎢
⎢⎣ b
c
a22
a
b⎤
a ⎥⎥
a33 ⎥⎦
(1.14)
Ως σύµβολο Kronecker ορίζουµε το εξής συµµετρικό σύστηµα 2ας τάξεως δ ij ,
10
Μηχανική του Συνεχούς Μέσου, Κεφ. 1., Ι. Βαρδουλάκης 2008
⎧1 if : i = j
δ ij = ⎨
⇒ δ ij = δ ji
⎩0 if : i ≠ j
(1.15)
Στο σύµβολο Kronecker αντιστοιχεί ο µοναδιαίος πίνακας
⎡δ11 δ12 δ13 ⎤
δ ij = δ ji ↔ [ I ] = ⎢⎢δ 21 δ 22 δ 23 ⎥⎥ =
⎢⎣δ 31 δ 32 δ 33 ⎥⎦
⎡1 0 0 ⎤
⎢0 1 0 ⎥
⎢
⎥
⎢⎣0 0 1 ⎥⎦
(1.16)
Ένα πλήρως συµµετρικό σύστηµα 3ης τάξεως χαρακτηρίζεται από τις εξής σχέσεις,
aijk = akij = a jki = aikj = a jik = akji
(1.17)
Ένα αντισυµµετρικό σύστηµα 2ας τάξεως χαρακτηρίζεται από την ιδιότητα ότι εναλλαγή
των δύο δεικτών επιφέρει αλλαγή του πρόσηµου του,
aij = − a ji
(1.18)
Από τον ορισµό αυτό προκύπτει ότι σε ένα αντισυµµετρικό σύστηµα 2ας τάξεως τα
διαγώνια στοιχεία είναι µηδενικά. Π.χ., σε έναν αντισυµµετρικό πίνακα 3 × 3 έχουµε,
a11 = a22 = a33 = 0
(1.19)
οπότε ο πίνακας αυτός περιέχει µόνο 3 ανεξάρτητα στοιχεία
⎡ a11
aij = −a ji ↔ ⎢⎢ a21
⎢⎣ a31
a12
a22
a32
a13 ⎤
a23 ⎥⎥ =
a33 ⎥⎦
⎡0 −c −b ⎤
⎢c 0 −a ⎥
⎢
⎥
⎢⎣b a
0 ⎥⎦
(1.20)
Αντιστοίχως ένα πλήρως αντισυµµετρικό σύστηµα 3ης τάξεως ικανοποιεί τις παρακάτω
σχέσεις
aijk = akij = a jki = − aikj = − a jik = − akji
(1.21)
1.3 Συµµετρικό και αντισυµµετρικό µέρος συστήµατος 2ας τάξεως
Ως το συµµετρικό και το αντισυµµετρικό µέρος του aij ορίζουµε αντιστοίχως τα
συστήµατα
a(ij ) =
1
(aij + a ji )
2
(1.22)
a[ij ] =
1
(aij − a ji )
2
(1.23)
Με τη χρήση των συµβόλων Kronecker µπορούµε να κατασκευάσουµε δύο συστηµάτων
4ης τάξεως, που έχουν την ιδιότητα να µετατρέπουν ένα τυχόν σύστηµα 2ας τάξεως σε
συµµετρικό και αντισυµµετρικό σύστηµα. Πράγµατι τα συστήµατα
Sijkl =
1
(δ ikδ jl + δ ilδ jk )
2
(1.24)
Μηχανική του Συνεχούς Μέσου, Κεφ. 1., Ι. Βαρδουλάκης 2008
Aijkl =
1 δ ik
2 δ jk
δ il 1
= (δ δ − δ δ )
δ jl 2 ik jl il jk
11
(1.25)
χαρακτηρίζονται αντιστοίχως από τις εξής ιδιότητες,
Sijkl akl = a(ij )
(1.26)
Aijkl akl = a[ij ]
(1.27)
Ασκήσεις
1) Να αποδειχθεί ότι ένα πλήρως αντισυµµετρικό σύστηµα 3ης τάξεως έχει µία µόνο µηµηδενική συνιστώσα,
⎧ a if : (ijk ) = cycl (1, 2,3)
⎪
aijk = ⎨− a if : (ijk ) = cycl (2,1,3)
⎪ 0 else
⎩
(1.28)
2) (Ευθύ) Να αποδειχθεί ότι αν το σύστηµα 2ας τάξεως aij είναι αντισυµµετρικό τότε η
δι-γραµµική µορφή,
aij x i x j = 0
(1.29)
3) (Αντίστροφο) Να αποδειχθεί ότι. αν ισχύει η παραπάνω Εξ. (1.29) για όλες τις τιµές
της µεταβλητής xi , τότε το σύστηµα aij είναι αντισυµµετρικό.
4) Έστω wik και ψ ik η µετατόπιση και στροφή της διατοµής αµφιέρειστης ελαστικής
δοκού ακαµψίας ( EI ) στη θέση xi λόγω µοναδιαίου φορτίου στη θέση xk (Εικ. 1-2).
Να διερευνηθεί και να αιτιολογηθεί αν τα συστήµατα 2ας τάξεως wik και ψ ik είναι
συµµετρικά ή όχι.
Εικ. 1-2: Η επιρροή µοναδιαίου φορτίου στην παραµόρφωση αµφιέρειστης δοκού
1.4 Το σύµβολο Levi-Civita και τα σύµβολα Kronecker
Το λεγόµενο µοναδιαίο αντισυµµετρικό ή αντιµεταθετικό σύστηµα 3ης τάξεως ή σύµβολο
Levi-Civita, ε ijk ορίζεται ως εξής:
Μηχανική του Συνεχούς Μέσου, Κεφ. 1., Ι. Βαρδουλάκης 2008
12
ε ijk
⎧ +1 if : (i, j , k ) = cycl (1, 2,3)
⎪
= ⎨ −1 if : (i, j , k ) = cycl (2,1,3)
⎪ 0 else
⎩
(1.30)
δηλ. ε123 = ε 231 = ε 312 = 1 , ε 213 = ε132 = ε 321 = −1 , ε112 = " = ε 333 = 0 .
Με τη βοήθεια του συµβόλου αυτού το τυχόν αντισυµµετρικό σύστηµα 3ης τάξεως
γράφεται ως εξής (γιατί;)
aijk = a123ε ijk
(1.31)
Όπως ήδη αναφέραµε θα χρειασθεί να διακρίνουµε ανάµεσα σε συστήµατα που έχουν
δείκτες αναγραµµένους είτε «κάτω», όπως το σύστηµα ε ijk , είτε «πάνω», όπως το
σύστηµα ε ijk , το οποίο θα ορισθεί κατ’ αναλογία της (1.30) ως,
ε
ijk
⎧ +1 if : (i, j , k ) = cycl (1, 2,3)
⎪
= ⎨ −1 if : (i, j , k ) = cycl (2,1,3)
⎪ 0 else
⎩
(1.32)
Ένα αντισυµµετρικό σύστηµα 2ας τάξεως προσδιορίζεται πλήρως από ένα σύστηµα 1ης
τάξεως,
⎡0
[ωij ] = ⎢⎢ω3
⎢⎣ω2
−ω3
0
ω1
−ω2 ⎤
−ω1 ⎥⎥ ↔ ωij = −ε ijk ωk
0 ⎥⎦
(1.33)
Στην περίπτωση αυτή έχουµε ότι,
1
2
ωl = − ε mnlωmn
(1.34)
Προφανώς αν ορίσουµε ως ωi′ = −ωi από τις παραπάνω σχέσεις παίρνουµε τις εξής
ισοδύναµες εκφράσεις,
ωij = ε ijk ωk′
1
⇔ ωl′ = ε mnlωmn
2
(1.35)
Το γινόµενο των δύο αντισυµµετρικών συστηµάτων 3ης τάξεως, ορίζει το εξής σύστηµα
6ης τάξεως
rst
ε mnpε rst = δ mnp
(1.36)
Το σύστηµα αυτό καλείται σύµβολο Kronecker 6ης τάξεως. Για τον προσδιορισµό του
rst
συστήµατος δ mnp
ορίζουµε κατ’ αρχήν τα συµµετρικά σύµβολα Kronecker 2ας τάξεως
⎧1 if : i = j
⎩0 if : i ≠ j
δ ij = δ ij = δ i⋅ j = δ ⋅ij = ⎨
Στα σύµβολα Kronecker 2ας τάξεως αντιστοιχεί ο µοναδιαίος πίνακας, π.χ.
(1.37)
Μηχανική του Συνεχούς Μέσου, Κεφ. 1., Ι. Βαρδουλάκης 2008
⎡1 0 0 ⎤
δ ij = δ = δ = δ ↔ [ I ] = ⎢⎢ 0 1 0 ⎥⎥
⎢⎣ 0 0 1 ⎥⎦
ij
⋅j
i
j
⋅i
13
(1.38)
Με τη βοήθεια των συµβόλων Kronecker 2ας τάξεως µπορούµε να ορίσουµε τις πράξεις:
α) της αναβιβάσεως, καταβιβάσεως και εναλλαγής ενός δείκτη
δ ij x j = xi , δ ij x j = xi , δ i j x j = xi , δ ij x j = xi
δ ij x j = xi , δ ij x j = xi , δ i j x j = xi , δ ij x j = xi
(1.39)
δ ik ak = ai , δ ik a k = a i
β) της συστολής5 ενός δείκτη
δ ij a ij = a ii = a11 + a 22 + a 33
δ ij aij = aii
δ i⋅ j a⋅i j = a⋅ii
(1.40)
δ ⋅ij a⋅ji = ai⋅i
Κάνοντας χρήση του ορισµού του συµβόλου Levi-Civita µπορούµε να αποδείξουµε ότι
rst
το σύµβολο Kronecker 6ης τάξεως δ mnp
παίρνει τις εξής τιµές:
rst
δ mnp
= 0 , όταν δύο ή περισσότεροι δείκτες ταυτίζονται,
rst
δ mnp
= +1 , όταν οι δείκτες (r,s,t) και (mnp) διαφέρουν κατά άρτιο αριθµό µεταθέσεων,
rst
δ mnp
= −1 , όταν οι δείκτες (r,s,t) και (mnp) διαφέρουν κατά περιττό αριθµό µεταθέσεων.
Αυτός ο ορισµός µπορεί να συνοψισθεί στην παρακάτω χρήσιµη σχέση, που συνδέει το
rst
µε µία ορίζουσα της οποίας τα στοιχεία είναι σύµβολα Kronecker 2ας
σύµβολο δ mnp
τάξεως,
rst
δ mnp
δ mr δ nr δ pr
= δ ms δ ns δ ps
δ mt δ nt δ pt
(1.41)
rst
Συστολή του συµβόλου δ mnp
δίνει το σύµβολο Kronecker 4ης τάξεως,
rs
rst
rsp
rsp
rs1
rs 2
rs 3
δ mn
= δ mnp
δ tp = δ mnt
δ tp = δ mnp
= δ mn
1 + δ mn 2 + δ mn 3
όπου αποδεικνύεται ότι,
5
Αγγλ. contraction, Γερµ. Verjüngung
(1.42)
Μηχανική του Συνεχούς Μέσου, Κεφ. 1., Ι. Βαρδουλάκης 2008
14
rs
δ mn
=
δ mr δ nr
= δ mr δ ns − δ nrδ ms
δ ms δ ns
(1.43)
οπότε
rst
rs
ε mnpε rsp = ε mnpε rstδ tp = δ mnp
δ tp = δ mn
= δ mr δ ns − δ nrδ ms
(1.44)
Άρα από τις εξ. (1.25), (1.27) και (1.44) παίρνουµε ότι,
a[ij ] =
1
1
(aij − a ji ) = ε ijpε klp akl
2
2
(1.45)
Περαιτέρω συστολή του συµβόλου Kronecker 4ης τάξεως οδηγεί στη σχέση του µε το
σύµβόλο Kronecker 2ας τάξεως,
rn
δ mn
= 2δ mr
(1.46)
Επίσης ισχύουν και οι κάτωθι χρήσιµες ταυτότητες,
ijk
hl
δ mnp
δ rsthlp = δ rstijk δ mn
ijk
δ mnp
δ rsthnp = 2!δ rstijk δ mh
(1.47)
ijk
δ mnp
δ rstmnp = 3!δ rstijk
Ασκήσεις
Nα αποδειχθεί ότι
ε rst =
1
( r − s )( s − t )( t − r ) , r , s, t ∈ {1, 2,3}
2
(1.48)
Να αποδειχθούν οι παρακάτω σχέσεις:
δ rr = 3
rst
δ mst
= 2δ mr
δ rstrst = 3!
(1.49)
rst
rst
rs
δ tmn
= δ ntm
= δ mn
rs mn
δ mn
a = a rs − a sr
rst
δ mnp
a mnp = a rst − a rts + a str − a srt + a trs − a tsr
Μηχανική του Συνεχούς Μέσου, Κεφ. 1., Ι. Βαρδουλάκης 2008
15
1.5 Συµµετρίες συστηµάτων 3ης τάξεως6
Έστω Aijk ένα τυχόν σύστηµα 3ης τάξεως. Με βάση το σύστηµα αυτό µπορούµε να
κατασκευάσουµε ένα πλήρως συµµετρικό σύστηµα 3ης τάξεως,
Sijk =
1
( Aijk + Aikj + Ajki + Ajik + Akij + Akji )
3!
(1.50)
και ένα πλήρως αντισυµµετρικό σύστηµα 3ης τάξεως,
Qijk =
1
( Aijk − Aikj + Ajki − Ajik + Akij − Akji )
3!
(1.51)
Παρατηρούµε όµως ότι το σύστηµα
Rijk = Aijk − ( Sijk + Qijk )
(1.52)
δεν είναι το µηδενικό σύστηµα 3ης τάξεως. Πράγµατι το «υπόλοιπο» Rijk µπορεί να
αναλυθεί ποικιλοτρόπως. Π.χ. το Rijk αναλύεται σε ένα σύστηµα συµµετρικό ως προς
τους δύο πρώτους δείκτες (i, j ) και σε ένα σύστηµα συµµετρικό ως προς τον πρώτο και
τρίτο δείκτη (i, k ) ,
Rijk = Sijk(1) + Sijk(2)
(1.53)
όπου
S
(2)
ijk
(
)
(
)
2
( Aijk − Akji ) + ( Ajik − Akij )
3!
2
=
( Aijk − Ajik ) + ( Akji − Ajki )
3!
Sijk(1) =
(1.54)
Προφανώς τα συστήµατα αυτά µηδενίζονται όταν το σύστηµα Aijk είναι πλήρως
συµµετρικό. Οµοίως µπορούµε να αναλύσουµε το υπόλοιπο Rijk σε ένα σύστηµα
αντισυµµετρικό ως προς τους δείκτες (i, k ) και σε ένα σύστηµα αντισυµµετρικό ως προς
τους δείκτες (i, j )
Rijk = Qijk(1) + Qijk(2)
(1.55)
όπου
(2)
ijk
Q
6
(
)
(
)
2
( Aijk + Ajik ) − ( Akji + Ajki )
3!
2
=
( Aijk + Akji ) − ( Ajik − Akij )
3!
Qijk(1) =
(1.56)
Wade; T. L. and Bruck, R. H. (1944). Types of symmetries. The American Mathematical Monthly, 51,
123-129.
Μηχανική του Συνεχούς Μέσου, Κεφ. 1., Ι. Βαρδουλάκης 2008
16
Τα συστήµατα αυτά µηδενίζονται όταν το σύστηµα Aijk είναι πλήρως αντισυµµετρικό.
Συµπέρασµα των ανωτέρω είναι ότι για να είναι ένα σύστηµα 3ης τάξεως συµµετρικό
πρέπει το αντισυµµετρικό του µέρος Qijk να µηδενίζεται αλλά αυτό δεν αρκεί, διότι
πρέπει επίσης να µηδενίζεται και το υπόλοιπο που συµβολίσαµε µε Rijk .
1.6 Ορίζουσες
Ας θεωρήσουµε την ορίζουσα του πίνακα [asr ]
a11
A = det[asr ] = asr = a12
a13
a12
a22
a23
a31
a32
a33
(1.57)
Παρατηρούµε ότι ισχύουν οι παρακάτω σχέσεις
asr = ε ijk a1i a2j a3k
anm ε rst = ε ijk air a sj akt
(1.58)
anm ε rst = ε ijk ari asj atk
Ας θεωρήσουµε το ανάπτυγµα της ορίζουσας ως προς τα στοιχεία της πρώτης στήλης,
a11
asr = a12
a13
a12
a22
a23
a31
a32 = a11 A11 + a12 A21 + a13 A31
a33
(1.59)
Στο ανάπτυγµα αυτό µε Ar1 έχουµε συµβολίσει τις ελάσσονες ορίζουσες των στοιχείων
a1r , δηλαδή
A11 = (−1)1+1
A = (−1)
1
2
A = (−1)
1
3
2 +1
3+1
a22
a23
a32
= a22 a33 − a23 a32 = ε123 a22 a33 + ε132 a23 a32 = ε1st a2s a3t
a33
a12
a23
a31
= − ( a12 a33 − a23 a31 ) = ε 213 a12 a33 + ε 231a23 a31 = ε 2 st a2s a3t (1.60)
a33
a12
a22
a31
= a12 a32 − a22 a31 = ε 312 a12 a32 + ε 321a22 a31 = ε 3 st a2s a3t
2
a3
Γενικώς ισχύει η σχέση,
ε 123ε rst a2s a3t =
1 1 jk
1 1 jk s t
ε ε rst a sj akt = δ rst
a j ak
2!
2!
οπότε το στοιχείο Ari , παίρνει την εξής µορφή,
(1.61)
Μηχανική του Συνεχούς Μέσου, Κεφ. 1., Ι. Βαρδουλάκης 2008
Ari =
1 ijk s t
δ rst a j ak
2!
17
(1.62)
Με τη βοήθεια αυτής της σχέσεως µπορούµε να υπολογίσουµε το παρακάτω ανάπτυγµα,
1 ijk m s t 1 ijk
δ mst ar a j ak = ε ε mst arm a sj akt
2!
2!
1
1 ijk p
= ε ijk ε rjk aqp = δ rjk
aq
2!
2!
arm Ami =
(1.63)
Οπότε παίρνουµε την παρακάτω έκφραση για το ανάπτυγµα µιας ορίζουσας,
arm Ami = Aδ ri
(1.64)
Οµοίως έχουµε και τη σχέση,
ami Arm = Aδ ri
(1.65)
Όταν η ορίζουσα είναι διάφορη του µηδενός ( A ≠ 0 ) τότε το στοιχείο
α ri =
1 i
Ar
A
(1.66)
είναι το συµπληρωµατικό του στοιχείου air και συµβολίζεται ως,
α ri = co ( air )
(1.67)
Οπότε, έχουµε τις σχέσεις
ami α rm = armα mi = δ ri
(1.68)
ή
⎡ a11
⎢ 2
⎢ a1
⎢ a13
⎣
a12
a22
a23
a31 ⎤ ⎡α11 α 21 α 31 ⎤ ⎡1 0 0 ⎤
⎥⎢
⎥
a32 ⎥ ⎢α12 α 22 α 32 ⎥ = ⎢⎢0 1 0 ⎥⎥
a33 ⎥⎦ ⎢⎣α13 α 23 α 33 ⎥⎦ ⎢⎣0 0 1 ⎥⎦
(1.69)
⎡ a11
⎢ 2
⎢ a1
⎢ a13
⎣
a12
a22
a31 ⎤ ⎡ A11
⎥1⎢
a32 ⎥ ⎢ A12
A
a33 ⎥⎦ ⎢⎣ A13
(1.70)
ή
3
2
a
A21
A22
3
2
A
A31 ⎤ ⎡1 0 0 ⎤
⎥
A32 ⎥ = ⎢⎢0 1 0 ⎥⎥
A33 ⎥⎦ ⎢⎣0 0 1 ⎥⎦
Η παραπάνω σχέσεις ορίζουν τον αντίστροφο ενός τετραγωνικού πίνακα: Έστω,
⎡ a11
[ A] = [asr ] = ⎢⎢ a12
⎢ a13
⎣
a12
a22
3
2
a
a31 ⎤
⎥
a32 ⎥
a33 ⎥⎦
(1.71)
Μηχανική του Συνεχούς Μέσου, Κεφ. 1., Ι. Βαρδουλάκης 2008
18
Αν η ορίζουσα του πίνακα αυτού είναι διάφορη του µηδενός, δηλαδή αν A = det [ A] ≠ 0 ,
τότε ορίζουµε τον αντίστροφο του πίνακα [ A] ως
⎡ A11
1
−1
[ A] = ⎢⎢ A12
A 3
⎢ A1
⎣
A21
A22
3
2
A
A31 ⎤
⎥
A32 ⎥
A33 ⎥⎦
(1.72)
Επίσης µπορούµε να αποδείξουµε τις παρακάτω χρήσιµες σχέσεις,
ε rst Ari = 3ε ijk a sj akt
ε ijk Ari = ε rst a sj akt
(1.73)
∂A
= Ari
r
∂ai
Αντίστοιχο τυπολόγιο µπορεί να αναπτυχθεί και για το σύστηµα amn :
a11
A = det[ars ] = ars = a21
a31
a12
a22
a32
a13
a23
a33
1 ijk rst
ε ε a js akt
2!
amr Ami = arm Aim = Aδ ri
Air =
(1.74)
1 ir
A ( A ≠ 0)
A
= armα im = δ ri
co ( ari ) = α ir =
amrα mi
1.7 Γραµµικά συστήµατα
Η λύση του γραµµικού συστήµατος
amr x m = b r ( A ≠ 0)
(1.75)
1 i m
Amb
A
(1.76)
είναι,
xi =
1.8 Θετικώς ορισµένες τετραγωνικές µορφές
Ορισµός: Η τετραγωνική µορφή, Q = amn x m x n , λέγεται θετικώς ορισµένη, όταν
⎧ 0 if : x r = 0
Q=⎨
r
⎩> 0 if : x ≠ 0
(1.77)
Μηχανική του Συνεχούς Μέσου, Κεφ. 1., Ι. Βαρδουλάκης 2008
19
Θεώρηµα: Θεωρούµε δύο τετραγωνικές µορφές Q = amn x m x n και S = bmn x m x n , όπου τα
συστήµατα amn και bmn είναι συµµετρικά. Επίσης θεωρούµε την ορίζουσα,
Λ ( λ ) = λ amn − bmn
λ a11 − b11 λ a12 − b12 λ a13 − b13
= λ a21 − b21 λ a22 − b22 λ a23 − b23
λ a31 − b31 λ a32 − b32 λ a33 − b33
(1.78)
Αν η τετραγωνική µορφή Q = amn x m x n είναι θετικώς ορισµένη, τότε η εξίσωση,
Λ (λ ) = 0
(1.79)
έχει πραγµατικές ρίζες.
Απόδειξη:
Έστω, ρ = α + i β , ρίζα της Εξ. (1.79). Τότε υπάρχει ένα µη-µηδενικό σύστηµα7,
z n = x n + iy n που να ικανοποιεί το γραµµικό σύστηµα8
[ ρ amn − bmn ]{ z n } = 0
(1.80)
[(α + i β )amn − bmn ]{ x n + iy n } = 0
(1.81)
ή
Εξισώνοντας τα πραγµατικά και φανταστικά µέρη της παραπάνω εξισώσεως παίρνουµε,
α amn x n − β amn y n − bmn x n = 0
(1.82)
α amn y n + β amn x n − bmn y n = 0
(1.83)
Οπότε πολλαπλασιάζοντας την Εξ. (1.82) επί y m και την Εξ. (1.83) επί x m και
αθροίζοντας πάνω στο δείκτη m παίρνουµε,
β (amn x n x m + amn y m y n ) = (bmn − bnm ) x m y n = 0
(1.84)
Τα στοιχεία x r και y r δεν είναι όλα κατ’ ανάγκην µηδέν και η τετραγωνική µορφή
amn x m x n είναι θετικώς ορισµένη, οπότε από την Εξ. (1.84) έπεται ότι β = 0 .
ο.ε.δ.
Αν πολλαπλασιάσουµε την Εξ. (1.82) επί x m και την Εξ. (1.83) επί y m και αθροίσουµε
πάνω στο δείκτη m , τότε παίρνουµε τη σχέση
α (amn x n x m + amn y m y n ) = bmn x m x n + bmn y m y n
7
(1.85)
Ο n είναι εµ’ προκειµένω ένας άνω δείκτης.
Υπενθυµίζουµε ότι ο δείκτης n δεν είναι εκθέτης και ότι πραγµατοποιείται άθροιση πάνω σε κάθε
επαναλαµβανόµενο δείκτη.
8
Μηχανική του Συνεχούς Μέσου, Κεφ. 1., Ι. Βαρδουλάκης 2008
20
Αν τώρα δεχθούµε ότι και η τετραγωνική µορφή bmn x m x n είναι επίσης θετικώς ορισµένη,
τότε από την Εξ. (1.85) παίρνουµε ότι, α > 0 . Άρα αν οι τετραγωνικές µορφές amn x m x n
και bmn x m x n είναι θετικώς ορισµένες και τα συστήµατα amn , bmn είναι συµµετρικά, τότε
οι ρίζες της εξισώσεως, λ amn − bmn = 0 , είναι όλες πραγµατικές και θετικές.
1.9 Η χαρακτηριστική εξίσωση τετραγωνικού πίνακα 3 × 3
Οι ιδιοτιµές ενός τετραγωνικού πίνακα [ Aij ] 3 × 3 ικανοποιούν τη χαρακτηριστική του
εξίσωση (γιατί;)
Aij − αδ ij = 0
(1.86)
ή
α 3 − I Aα 2 + II Aα − III A = 0
(1.87)
όπου I A , II A και III A είναι οι λεγόµενες βασικές αναλλοίωτες του πίνακα [ A]3×3 , που
δίδονται από τις παρακάτω σχέσεις, συναρτήσει των στοιχείων του πίνακα [ A] = [aij ] και
των ιδιοτιµών του α i′ (i′ = 1, 2,3) :
I A = A11 + A22 + A33 = α1 + α 2 + α 3
II A =
A11
A21
A12 A22
+
A22 A32
A23 A11
+
A33 A31
A11
A12
A13
III A = A21
A 31
A22
A32
A23 = α1α 2α 3
A33
(1.88)
A13
= α1α 2 + α 2α 3 + α 3α1
A33
(1.89)
(1.90)
Θεώρηµα
Οι ιδιοτιµές ενός συµµετρικού πραγµατικού πίνακα είναι πραγµατικοί αριθµοί.
Απόδειξη
Έστω
Aij n j = λ n j
,
Aij = Aji ∈ \
(1.91)
και έστω
ni = ai + ibi
, ni = ai − ibi
Από την Εξ. (1.91) παίρνουµε,
(1.92)
Μηχανική του Συνεχούς Μέσου, Κεφ. 1., Ι. Βαρδουλάκης 2008
Aij ni n j =
21
1
( Aij ni n j + Aji n j ni )
2
(
)
1
Aij ( ai + ibi ) ( a j − ib j ) + Aji ( a j + ib j ) ( ai − ibi )
2
1
= ( Aij ai a j − iAij ai b j + iAij bi a j + Aij bi b j + Aji a j ai − iAji a j bi + iAji b j ai + Aji b j bi )
2
1
= ( Aij ai a j − iAij ai b j + iAij bi a j + Aij bi b j + Aij a j ai − iAij a j bi + iAij b j ai + Aij b j bi )
2
= Aij ai a j + Aij bi b j ∈ \
=
αλλά και
1
( λ n j n j + λ ni ni )
2
= λ ( ai + ibi )( ai − ibi )
Aij ni n j =
= λ ( ai ai − iai bi + ibi ai + bi bi )
= λ ( ai ai + bi bi ) ∈ \ ⇒ λ ∈ \
ο.ε.δ.
Ασκήσεις
1. Να γραφεί πρόγραµµα H/Y, που για δεδοµένο συµµετρικό πίνακα να υπολογίζει τις
ιδιοτιµές και τα ιδιο-διανύσµατά του συµπεριλαµβανοµένων και των περιπτώσεων
⎡1 2 0 ⎤
πολλαπλών ιδιοτιµών. Π.χ. ο πίνακας: [ Aij ] = ⎢⎢ 2 2 2 ⎥⎥ έχει τις ιδιοτιµές,
⎢⎣ 0 2 3 ⎥⎦
{α i′ } = {−1, 2,5} .
2. Ποια είναι η γεωµετρική σηµασία των ιδιοτιµών και ιδιοδιανυσµάτων ενός
συµµετρικού πίνακα Aij σε σχέση µε την αντίστοιχη επιφάνεια 2ας τάξεως που
ορίζεται από τη σχέση, {xi }[ Aij ]{x j }T = 1 . Για τη διερεύνηση αυτού του ερωτήµατος
σκόπιµο θα ήταν να σχεδιασθεί η αντίστοιχη επιφάνεια στο σύστηµα κυρίων αξόνων
του πίνακα [ Aij ] .
22
Μηχανική του Συνεχούς Μέσου, Κεφ. 1., Ι. Βαρδουλάκης 2008
1.10 Ο τελεστής " ∇ " 9
1.10.1 Ο τελεστής βαθµίδας
Εικ. 1-3: Η κάθετος σε σηµείο επί επιφανείας στο χώρο
Ας θεωρήσουµε καρτεσιανές συντεταγµένες x1 = x , x2 = y , x3 = z και την εξίσωση,
φ ( xi ) = φ ( x, y, z ) = c = const.
(1.93)
Για διάφορες τιµές της σταθεράς c παίρνουµε µία οικογένεια επιφανειών (Εικ. 1-3). Για
κάποια τιµή της σταθεράς c µία από αυτές τις επιφάνειες διέρχεται από δεδοµένο σηµείο
P ( x, y, z ) . Η κάθετη στην επιφάνεια αυτή στο σηµείο αυτής P ( x, y, z ) δίδεται από τα
συνηµίτονα κατευθύνσεώς της, που όπως γνωρίζουµε από την Αναλυτική Γεωµετρία των
επιφανειών είναι ανάλογα των ποσοτήτων
⎧ ∂φ ⎫
⎧ ∂φ ∂φ ∂φ ⎫
⎨ ⎬ =⎨ , , ⎬
⎩ ∂xi ⎭P ⎩ ∂x ∂y ∂z ⎭ P
(1.94)
Οπότε το διάνυσµα
∂φ G
ei
∂xi P
(1.95)
είναι παράλληλο προς την κάθετο στην επιφάνεια φ ( x, y, z ) = C στο σηµείο P . Έστω
dn το µήκος ενός απειροστικού ευθύγραµµου τµήµατος κατά µήκος της καθέτου στην εν
λόγω επιφάνεια στο σηµείο P . Τα συνηµίτονα κατευθύνσεως της καθέτου αυτής είναι
9
Το κεφάλαιο αυτό είναι βασισµένο στο αντίστοιχο από το βιβλίο: D.E. Rutherford, Vector Methods,
Dover, 2004.
Μηχανική του Συνεχούς Μέσου, Κεφ. 1., Ι. Βαρδουλάκης 2008
⎧ dx dy dz ⎫
⎨ , , ⎬
⎩ dn dn dn ⎭
23
(1.96)
Οπότε το µέγεθος του διανύσµατος που δίδεται από την Εξ. (1.95) είναι
∂φ dx ∂φ dy ∂φ dz ∂φ
+
+
=
∂x dn ∂y dn ∂z dn ∂n
(1.97)
Εισάγοντας τον τελεστή βαθµίδας σε καρτεσιανές συντεταγµένες,
G G ∂
∇ = ei i
∂x
(1.98)
Ορίζουµε την βαθµίδα του βαθµωτού µεγέθους φ ( xi ) ως
G
∂φ G
G
grad φ = ∇φ =
ei = φ,i ei
∂xi
, φ,i =
∂φ
∂xi
(1.99)
1.10.2 Ο τελεστής στροβιλισµού
Εικ. 1-4: Κλειστή όδευση στο επίπεδο
Έστω ένα διανυσµατικό πεδίο,
G
G
G
G
G
v = vi ei = vx ex + v y ey + vz ez
(1.100)
Θεωρούµε µία κλειστή καµπύλη (C ) = ( ABCD) η οποία περικλείει το σηµείο P ( x, y, z )
G G
σε ένα επίπεδο το οποίο είναι κάθετο στο διάνυσµα n = ex (Εικ. 1-4). Επίσης
εξασφαλίζουµε όπως η όδευση της καµπύλης αυτής είναι δεξιόστροφη. Κατά µήκος
αυτής της καµπύλης (C ) µε παράµετρο το µήκος τόξου s
(C ) : x = x( s),
y = y ( s), z = z ( s)
G
ορίζουµε το στοιχειώδες διάνυσµα ds ,
(1.101)
Μηχανική του Συνεχούς Μέσου, Κεφ. 1., Ι. Βαρδουλάκης 2008
24
G ⎛ dx G dy G dz G ⎞
ds = ⎜ ex + ey + ez ⎟ ds
ds
ds ⎠
⎝ ds
(1.102)
Π.χ. για την καµπύλη στην Εικ. 1-4,
G
G ⎧dyey ( AB), (CD)
ds = ⎨ G
⎩ dzez ( BC ), ( DA)
(1.103)
Με βάση τα παραπάνω ορίζουµε την ποσότητα
ΓC =
v∫ v ds
i
(1.104)
i
(Γ)
G
η οποία και καλείται η (δεξιόστροφη) κυκλοφορία κατά Kelvin της ποσότητας v (π.χ.
ταχύτητας) κατά µήκος της καµπύλης C .
G
Ένα διανυσµατικό µέγεθος, το οποίο θα συµβολίσουµε ως rot v , και το οποίο
G
αποκαλούµε στροβιλισµό του διανύσµατος v , έχει συνιστώσα στην κατεύθυνση του
G
διανύσµατος n που ορίζεται ως εξής:
ΓC
S →0 S
(1.105)
lim
όπου S είναι το εµβαδόν της επιφάνειας που περικλείεται από την καµπύλη (C ) .
Αποδεικνύεται ότι όταν το παραπάνω όριο υπάρχει όντως η Εξ. (1.105) ορίζει ένα
διάνυσµα, το οποίο είναι ανεξάρτητο της επιλογής της (C ) .
G
Βάσει της Εικ. 1-4 η αναλυτική έκφραση του διανύσµατος rotv υπολογίζεται ως εξής:
y+β
∫
vi d A i =
∫
vi d A i =
∫
v y dy = 2 β v y ( x, y′, z − γ ) ,
y − β ≤ y′ ≤ y + β
y−β
( AB )
z +γ
∫ v dz = 2γ v ( x, y + β , z′)
z
z
, z − γ ≤ z′ ≤ z + γ
z −γ
( BC )
y−β
∫
vi d A i =
∫
vi d A i =
∫
(1.106)
v y dy = −2 β v y ( x, y′, z + γ ) ,
y − β ≤ y′ ≤ y + β
y+β
( CD )
z −γ
( DA )
∫ v dz = −2γ v ( x, y − β , z′)
z
z
, z − γ ≤ z′ ≤ z + γ
z +γ
Με
S = 4βγ
Από την Εξ. (1.105) και (1.106) παίρνουµε
(1.107)
Μηχανική του Συνεχούς Μέσου, Κεφ. 1., Ι. Βαρδουλάκης 2008
2γ vz ( x, y + β , z′ ) − 2γ vz ( x, y − β , z′ )
β →0 γ → 0
4 βγ
2β v y ( x, y′, z + γ ) − 2β v y ( x, y′, z − γ )
− lim lim
β →0 γ →0
4 βγ
25
lim lim
v ( x, y′, z + γ ) − v y ( x, y′, z − γ )
v ( x, y + β , z ′ ) − v z ( x, y − β , z ′ )
= lim z
− lim y
β →0
γ →0
2β
2γ
∂
∂
= v z ( x, y , z ) − v y ( x, y , z )
∂y
∂z
(1.108)
Οπότε,
G ⎛ ∂v ∂v ⎞ G ⎛ ∂v ∂v
rot v = ⎜ z − y ⎟ ex + ⎜ x − z
⎝ ∂z ∂x
⎝ ∂y ∂z ⎠
⎞ G ⎛ ∂v y ∂vx ⎞ G
−
⎟ ez
⎟ ey + ⎜
∂
∂
x
y
⎠
⎝
⎠
(1.109)
Συµφώνως προς την Εξ. (1.109) ο διαφορικός τελεστής του στροβιλισµού του
G
G
διανύσµατος v ( xi ) ορίζεται ως το εξωτερικό γινόµενο του τελεστή βαθµίδας και του v ,
G G G ⎛G ∂ ⎞
G
rot v = ∇ × v = ⎜ ek
⎟ × ( vl el )
⎝ ∂xk ⎠
Παρατηρούµε ότι ο στροβιλισµός του διανύσµατος
αντισυµµετρικό µέρος της βαθµίδας του.
(1.110)
G
v ( xi ) . ταυτίζεται µε το
1.10.3 Ο τελεστής αποκλίσεως
Εικ. 1-5: Κλειστή επιφάνεια στο χώρο
Κατ’ αναλογία µε τη διαδικασία που εκθέσαµε στη προηγούµενη παράγραφο ορίζουµε
µία κλειστή επιφάνεια ( S ) , που περιβάλει ένα δεδοµένο σηµείο P ( xi ) στο χώρο. Π.χ. η
Μηχανική του Συνεχούς Μέσου, Κεφ. 1., Ι. Βαρδουλάκης 2008
26
κλειστή επιφάνεια µπορεί να αποτελείται από τις πλευρές ενός ορθογωνίου
παραλληλεπιπέδου του οποίου οι πλευρές είναι παράλληλες προς τους άξονες του
καρτεσιανού συστήµατος συντεταγµένων (Εικ. 1-5).
Γενικώς σε κάθε σηµείο της εν λόγω κλειστής επιφάνειας θεωρούµε το εξωτερικό
G
µοναδιαίο διάνυσµα n και αν dS συµβολίζει το εµβαδόν του στοιχείο της επιφάνειας
G
G
τότε ορίζουµε το διάνυσµά ndS . Για ένα δεδοµένο διανυσµατικό µέγεθος v ( xi )
ορίζουµε την ποσότητα,
QS = ∫∫ vi ni dS
(1.111)
(S )
G
καλείται η ροή του v ( xi ) .
G
Ένα βαθµωτό µέγεθος, το οποίο θα συµβολίσουµε ως divv , και το οποίο αποκαλούµε
G
απόκλιση του διανύσµατος v , έχει τιµή που ορίζεται ως εξής:
Q
G
divv = lim S
V →0 V
(1.112)
όπου V είναι ο όγκος που περικλείεται από την επιφάνεια ( S ) .
Άσκηση
Να αποδειχθεί ότι
και ότι,
∂v
G ∂v ∂v ∂v
divv = x + y + z = i = vi ,i
∂x ∂y ∂z ∂xi
(1.113)
G G G
divv = ∇ ⋅ v
(1.114)
Μηχανική του Συνεχούς Μέσου, Κεφ. 1., Ι. Βαρδουλάκης 2008
27
1.11 Παράρτηµα I: Χρήσιµες ταυτότητες µεταξύ διαφορικών τελεστών
Μηχανική του Συνεχούς Μέσου, Κεφ. 1., Ι. Βαρδουλάκης 2008
28
1.12 Παράρτηµα ΙΙ: Το θεώρηµα αποκλίσεως10
Θεωρούµε ένα χωρίο V του ℜ 3 που περιβάλλεται από το σύνορο ∂ V . Στο τυχόν
σηµείο του συνόρου ορίζουµε την στοιχειώδη επιφάνεια d S µε µοναδιαίο εξωτερικό
διάνυσµα n i .
Εικ. 1-6: Κανονικό χωρίο εφαρµογής του θεωρήµατος αποκλίσεως
Έστω στο χωρίο αυτό µία διανυσµατική συνάρτηση qi = qi ( x k ) ( i , k = 1, 2,3) , που
είναι συνεχής και έχει συνεχείς πρώτες παραγώγους. Τότε ισχύει
∂qk
dV = ∫ qk nk dS
∂xk
(V )
( ∂V )
(1.115)
∫ div q dV = ∫
(1.116)
∫
ή
(V )
q ⋅ n dS
( ∂V )
δηλαδή το ολοκλήρωµα της αποκλίσεως ενός διανυσµατικού πεδίου πάνω στο χωρίο V
ισούται µε την συνολική «ροή» του πεδίου δια µέσου του συνόρου αυτού ∂V .
Άσκηση
Να αποδειχθεί ότι11
∂ψ
dV = ∫ niψ dS
∂xk
(V )
( ∂V )
∫
10
(1.117)
Το θεώρηµα αυτό παρουσιάστηκε υπό διαφορετικές µορφές από τους Lagrange (1762), Gauss (1813),
Ostrogradsky (1831) και Green (1828). Καµιά φορά αποκαλείται και Θεώρηµα Αποκλίσεως (Divergence
Theorem). Στη σχετική βιβλιογραφία θα αναζητήσουµε επίσης και τα παρεµφερή Θεωρήµατα Green και
Stokes.
11
Hay, G.E., Vector and Tensor Analysis, Dover, p. 149 ff, 1953.