vektorok, tenzorok - Műszaki Mechanikai Tanszék - Budapesti ...
vektorok, tenzorok - Műszaki Mechanikai Tanszék - Budapesti ...
vektorok, tenzorok - Műszaki Mechanikai Tanszék - Budapesti ...
You also want an ePaper? Increase the reach of your titles
YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.
BUDAPESTI MŰSZAKI ÉS GAZDASÁGTUDOMÁNYI EGYETEM<br />
Műszaki <strong>Mechanikai</strong> Tanszék<br />
_____________________________________________________________________<br />
VEKTOROK ÉS TENZOROK<br />
OKTATÁSI SEGÉDLET
2<br />
___________________________________________________________________________<br />
Összeállította:<br />
dr. Vörös Gábor, egyetemi docens, 2003. szeptember<br />
módosítva: 2005. január
3<br />
___________________________________________________________________________<br />
TARTALOM<br />
1. Az összegzési konvenció........................................................................................................ 4<br />
1.1. Vektorok szorzása ........................................................................................................... 5<br />
2. Másodrendű tenzor, a lineáris vektor-vektor függvény ......................................................... 7<br />
2.1. Tenzorok szorzása......................................................................................................... 10<br />
3. Transzformációk................................................................................................................... 11<br />
4. Tenzor sajátértékei, a skalár invariánsok ............................................................................. 13<br />
5. Görbevonalú koordináta rendszerek..................................................................................... 14<br />
5.1. A kovariáns derivált ...................................................................................................... 16<br />
5.2. Differenciálási szabályok .............................................................................................. 17<br />
5.3. Hengerkoordináta rendszer ........................................................................................... 19<br />
6. Ajánlott irodalom ................................................................................................................. 22
4<br />
___________________________________________________________________________<br />
1. Az összegzési konvenció<br />
1<br />
e1<br />
3<br />
e3<br />
e2<br />
R<br />
2<br />
1. ábra<br />
Az 1. ábra szerinti derékszögű – Descartes féle – koordináta rendszer három tengelyének<br />
irányát az e1, e2 és e3 egység<strong>vektorok</strong>, más szóval a bázis<strong>vektorok</strong> jelölik ki. Ha a koordináta<br />
tengelyeken mért ívhosszak x1, x2, és x3, akkor a P pont R helyvektora<br />
P<br />
F<br />
3<br />
∑<br />
k=<br />
1<br />
R = x e + x e + x e = x e ,<br />
1<br />
és egy, a P ponthoz kötött F vektor pedig a következő módon irható fel:<br />
1<br />
2<br />
2<br />
3<br />
3<br />
3<br />
∑<br />
k=<br />
1<br />
F = F e + F e + F e = F e ,<br />
1<br />
1<br />
2<br />
2<br />
ahol az Fi skalár mennyiségek az F vektor koordinátái. A továbbiakban az ilyen<br />
kifejezésekben az összegzés Σ jelét elhagyjuk, és megállapodunk abban, hogy a két azonos<br />
latin betű index – a néma indexpár – egy szorzat kifejezésen belül összegzési utasítást jelent.<br />
Az összegző indexek értékének felső határa térbeli feladatoknál 3. Természetesen, mivel nem<br />
az indexként felhasznált betű konkrét formája, hanem az indexpár megjelenése utal az<br />
összegzésre, bármely kis vagy nagy latin betű felhasználható. Ezek szerint, az összegzési<br />
konvenció alkalmazásával az F vektor röviden az<br />
alakban is felírható.<br />
k<br />
3<br />
k<br />
3<br />
p<br />
p<br />
k<br />
k<br />
k<br />
k<br />
F = F e = F e<br />
(1)<br />
Ha egy szorzat kifejezésben több néma indexpár is szerepel, akkor az többszörös összegzést<br />
jelent. Például három néma indexpárral egy huszonhét tagú összeg irható le röviden a<br />
következő formában:<br />
3<br />
∑∑∑<br />
a b c d = a b c d .<br />
ijk<br />
i<br />
j<br />
k<br />
3<br />
3<br />
i=<br />
1 j=<br />
1 k=<br />
1<br />
ijk<br />
i<br />
j<br />
k
5<br />
___________________________________________________________________________<br />
Ilyen esetekben ügyelni kell arra, hogy a néma indexekhez már felhasznált betűjelet ne<br />
használjuk még egyszer, mert akkor az összegzési konvenció elveszíti az egyértelműségét.<br />
Ha azt akarjuk, hogy egy kifejezésben a néma indexpár ne jelentsen összegzést, akkor egyik<br />
indexet egy megkülönböztető jellel, aláhúzással, „kizárjuk”. Ezek szerint az a ke k nem az (1)<br />
szerinti a vektor, hanem annak csak a k-adik vektorkomponense. A bázisvektor k indexe<br />
ilyenkor szabad index, a koordináta k indexe pedig zárt.<br />
A most bevezetett indexes jelölés nagy előnye a tömörebb leírás mellett az, hogy formálisan<br />
szétválasztottuk a vektort két részre, a skalár koordináták és a bázis egység<strong>vektorok</strong><br />
szorzatára. Ez a vektor és tenzor algebra és analízis körében elvégzendő számításokat<br />
jelentősen leegyszerűsíti, lerövidíti.<br />
A kontinuumechanika körében gyakran használják együtt az alsó és felső (kovariáns és<br />
kontravariáns) indexet. Erre akkor van szükség, ha a koordináta rendszerünk nem mindig<br />
derékszögű. A továbbiakban kizárólag a mérnöki gyakorlatban fontos ortogonális<br />
rendszereket használunk.<br />
1.1. Vektorok szorzása<br />
A következőkben röviden összefoglaljuk a <strong>vektorok</strong> körében ismert szorzási szabályokat. Ha<br />
felhasználjuk az (1) összegzési konvenciót, akkor a <strong>vektorok</strong> különböző szorzatait csak az<br />
ortogonális bázis egység<strong>vektorok</strong> között kell értelmezni, mivel a skalár koordináták szorzata<br />
egyértelmű.<br />
A. Skalár szorzat: Két bázis egységvektor szorzata a kétindexes Kronecker szimbólum:<br />
δ<br />
ij<br />
= 1<br />
ha<br />
e<br />
i<br />
⋅ e<br />
j<br />
i = j<br />
= δ<br />
és<br />
ij<br />
,<br />
δ<br />
ij<br />
= 0<br />
ha i ≠<br />
Ezzel két tetszőleges a és b vektor skalár szorzata a következő módon irható fel:<br />
a ⋅ b = (a e ) ⋅ (b e ) = a b ( e ⋅e<br />
) = a b δ = a b = a b + a b + a b .<br />
i<br />
i<br />
j<br />
j<br />
i<br />
j<br />
i<br />
j<br />
A kilenctagú összegnek csak a három nem feltétlenül zérus tagja van. A továbbiakban<br />
gyakran felhasználjuk a Kronecker szimbólumnak a definíciójából következő indexcserélő<br />
szerepét:<br />
j<br />
i<br />
ij<br />
j<br />
ij<br />
i<br />
i<br />
i<br />
1<br />
1<br />
j.<br />
2<br />
2<br />
3<br />
3<br />
(2)<br />
b δ = b . (3)<br />
Ha kikötjük, hogy a Kronecker szimbólum két indexe legyen mindig azonos, akkor az<br />
összegzési konvenció alkalmazásával belátható, hogy<br />
δii 11 22 33<br />
= δ + δ + δ = 3 . (4)
6<br />
___________________________________________________________________________<br />
B. Vektor szorzat: Két bázis egységvektor vektor szorzatát a permutációs szimbólum<br />
segítségével írjuk fel:<br />
e × = . (5a)<br />
p eq<br />
e pqr e r<br />
A vektor szorzat tulajdonságaiból következik, hogy az epqr ferdén szimmetrikus, azaz két<br />
szomszédos index felcserélése előjelváltást jelent és értéke zérus, ha bármely két, vagy mind a<br />
három index azonos. A zérustól különböző hat érték:<br />
Az a és b <strong>vektorok</strong> vektorszorzata:<br />
e123 312 231<br />
213 321 132<br />
= e = e = + 1 , e = e = e = −1<br />
.<br />
(5b)<br />
c = a× b = (a e ) × (b e ) = ab( e × e ) = abe e =<br />
i i j j i j i j i j ijk k<br />
= (a b − a b ) e + (a b − a b ) e + (a b − a b ) e = c ek<br />
.<br />
2 3 3 2 1 3 1 1 3 2 1 2 2 1 3 k<br />
A huszonhét tagú összeg többi huszonegy tagja a permutációs szimbólum (5b) definíciójából<br />
következően zérus.<br />
Az eddigiek felhasználásával számítsuk ki a bázis egység<strong>vektorok</strong> kettős vektor szorzatát:<br />
Másrészről, a vektoralgebrából ismert<br />
kifejtési tétel alapján:<br />
n<br />
n<br />
( e j × e k ) = e n × ( e jkt e t ) e jkte<br />
ntp e p<br />
e × = .<br />
( b × c)<br />
= b(<br />
a ⋅ c)<br />
− c(<br />
a b)<br />
a ×<br />
⋅<br />
( e j × e k ) = e j(<br />
e n ⋅e<br />
k ) − e k ( e n ⋅e<br />
j ) = e jδ<br />
nk e kδ<br />
nj<br />
e × − .<br />
Mivel a kétféle módon kiszámolt eredmény azonos,<br />
e e e = e δ − e δ .<br />
jkt<br />
ntp<br />
p<br />
Fontos, hogy az egyenlőség mindkét oldalán ugyanazok a j, k, n szabad indexek szerepelnek.<br />
Mivel ezek lehetséges értéke 1, 2 és 3, ez az egy sorban felirt egyenlőség elvileg huszonhét<br />
egyenlőséget fejez ki, továbbá a bal oldalon - a t és p néma indexpárok lehetséges értékeit<br />
számba véve – kilenc, illetve az index érték azonosságát kizárva, hat tagú összegek<br />
szerepelnek.<br />
Szorozzuk meg a fenti egyenlőség mindkét oldalát skalárisan az es bázis egységvektorral:<br />
jkt<br />
ntp<br />
ps<br />
j<br />
js<br />
nk<br />
e e δ = δ δ − δ δ .<br />
Az (3) és a permutációs szimbólum indexsorrendjére vonatkozó (5b) szabály alapján az<br />
egyenlőség bal oldala átalakítható:<br />
jkt<br />
snt<br />
nk<br />
js<br />
nk<br />
k<br />
ks<br />
nj<br />
nj<br />
e e = δ δ − δ δ . (6a)<br />
ks<br />
nj
7<br />
___________________________________________________________________________<br />
Ha kikötjük, hogy ebben az egyenlőségben a k és n szabad indexek legyenek azonosak, akkor<br />
a (3) és (4) felhasználásával:<br />
továbbá, a j és s indexek azonosságából<br />
következik.<br />
Három vektor vegyes szorzata:<br />
e e = δ δ − δ δ = 3δ<br />
− δ = 2δ<br />
, (6b)<br />
jkt<br />
skt<br />
js<br />
kk<br />
ks<br />
kj<br />
js<br />
js<br />
e jkt e jkt = 2δ<br />
jj = 6<br />
(6c)<br />
( b × ) = a ib<br />
jc<br />
ke<br />
ijk<br />
a ⋅ c<br />
. (7)<br />
Ha ez a vegyes szorzat nem zérus, akkor a három vektor nincs egy síkban, más szóval<br />
lineárisan függetlenek. A vegyes szorzat eredménye geometriailag úgy is értelmezhető, mint a<br />
három vektor által meghatározott párhuzamos oldalélű test (paralel epipedon) térfogata. Az<br />
ortogonális bázis<strong>vektorok</strong> – egység<strong>vektorok</strong> - vegyes szorzata az (5) permutációs szimbólum<br />
( p × eq<br />
) ⋅e<br />
r = e pqr<br />
e .<br />
C. Diadikus szorzat: Két vektor (a b) diadikus – vagy általános – szorzata definíció szerint<br />
legyen a következő tulajdonságú:<br />
( ab) ∗ c = a(<br />
b ∗ c)<br />
, c ∗ ( ab)<br />
= ( c ∗ a)b<br />
js<br />
, (8)<br />
ahol a * jel helyére bármelyik szorzás – a skalár vagy a vektor szorzás - jelét írhatjuk. Például<br />
( a b)<br />
⋅ c = a ( b ⋅ c)<br />
, c ⋅ ( ab)<br />
= ( c ⋅ a)<br />
b vagy ( ab)<br />
× c = a(<br />
b × c)<br />
általában (a b)≠ (b a).<br />
2. Másodrendű tenzor, a lineáris vektor-vektor függvény<br />
. Ez alapján nyilvánvaló, hogy<br />
A (8) definíció alapján számítsuk ki az a1 és b1 <strong>vektorok</strong> diádjának és az r helyvektornak a<br />
jobb oldali skaláris szorzatát:<br />
( a1b<br />
1 ) ⋅ r = a1(<br />
b1<br />
⋅ r)<br />
a1b1kx<br />
k<br />
R = = .<br />
Látszik, hogy itt az R vektor, az r irányától függetlenül, mindig párhuzamos az a1 vektorral.<br />
Ez tehát egy olyan lineáris homogén függvénykapcsolat az R és r <strong>vektorok</strong> között, ami a b1<br />
vektorra merőleges sík pontjaihoz – ahol a b1 és az r skalár szorzata azonos - az a1 irányú<br />
egyenesen ugyanazt a pontot rendeli. Az ilyen függvény elfajuló (nem invertálható), abban az<br />
értelemben, hogy a háromméretű teret egyméretű altérre – egy egyenesre – képezi le.
8<br />
___________________________________________________________________________<br />
Bővítsük fenti szorzatot még egy diáddal úgy, hogy az a1 és a2, valamint a b1 és b2 <strong>vektorok</strong><br />
ne legyenek párhuzamosak, azaz a1×a2 ≠ 0, b1×b2 ≠ 0:<br />
[ ( a b ) + ( a b ) ] ⋅ r = a ( b ⋅ r)<br />
+ a ( b r)<br />
R ⋅<br />
= 1 1 2 2<br />
1 1 2 2<br />
Most az R vektor az a1 és a2 <strong>vektorok</strong> által meghatározott síkban van. Ez a kibővített lineáris<br />
függvény a b1 és b2 <strong>vektorok</strong>ra merőleges síkok metszésvonalán lévő pontokhoz az a1, a2<br />
<strong>vektorok</strong> síkjában egy pontot rendel, vagyis a háromméretű teret egy kétméretű altérre, síkra<br />
képezi le.<br />
Ezek után már könnyen belátható, hogy a legáltalánosabb lineáris homogén vektor-vektor<br />
függvény három diád összegeként adható meg:<br />
[ ( a b ) + ( a b ) + ( a b ) ] ⋅r<br />
= a ( b ⋅r<br />
) + a ( b ⋅r<br />
) + a ( b r)<br />
R ⋅<br />
= 1 1 2 2 3 3<br />
1 1 2 2<br />
3 3 ,<br />
feltéve, hogy az ak és bk <strong>vektorok</strong> lineárisan függetlenek, más szóval az (7) vegyes szorzatuk<br />
nem zérus: a1.(a2×a3) ≠0 és b1.(b2×b3) ≠0.<br />
A továbbiakban, ha a bk <strong>vektorok</strong> helyett az ek bázis egység<strong>vektorok</strong>at alkalmazzuk, az (1)<br />
összegzési konvenció és az ak = apkep azonosság felhasználásával:<br />
ahol<br />
[ ( a e ) + ( a e ) + ( a e ) ] ⋅r<br />
= ( a e ) ⋅r<br />
= ( e e ) ⋅r<br />
= A r<br />
R = 1 1 2 2 3 3<br />
k k a pk p k ⋅ , (9a)<br />
A = a e e<br />
(9b)<br />
egy másodrendű tenzor és apk az A tenzor 3×3 méretű mátrixa, p a sor és k az oszlop<br />
sorszáma. Az ak <strong>vektorok</strong> az A tenzor mátrixának oszlopvektorai:<br />
[ ]<br />
⎣ pk⎦ 1 2 3 ⎢ 21 22 23⎥<br />
⎢⎣ a31 a32 a ⎥ 33 ⎦<br />
pk<br />
p<br />
k<br />
⎡a11 a12 a13<br />
⎤<br />
⎡a ⎤ = a a a =<br />
⎢<br />
a a a<br />
⎥<br />
. (9c)<br />
Egy tenzor által leirt leképzés, a lineáris vektor – vektor függvény – tehát akkor nem elfajuló,<br />
ha az ap oszlop<strong>vektorok</strong> lineárisan függetlenek, más szóval az (7) vegyes szorzatuk nem<br />
zérus. Az oszlop<strong>vektorok</strong> vegyes szorzata a tenzor determinánsa:<br />
det<br />
( A)<br />
= a ⋅ ( a × a ) = a a a e ⋅ ( e × e ) = a a a e ⋅ ( e e )<br />
1<br />
= a<br />
k1<br />
a<br />
2<br />
p2<br />
a<br />
q3<br />
3<br />
e<br />
pqs<br />
δ<br />
k1<br />
ks<br />
p2<br />
= a<br />
k1<br />
q3<br />
a<br />
k<br />
p2<br />
a<br />
q3<br />
p<br />
e<br />
pqk<br />
q<br />
= a<br />
k1<br />
a<br />
k1<br />
p2<br />
a<br />
p2<br />
q3<br />
e<br />
q3<br />
.<br />
kpq<br />
k<br />
.<br />
pqs<br />
s<br />
=<br />
(10a)<br />
A tenzor oszlopvektor sorrendjének cseréjével a determináns értéke nem változik, csak<br />
esetleg az előjele. Például<br />
( ) 1 ( 2 3) 1 ( 3 2)<br />
det A = a ⋅ a × a = −a ⋅ a × a .
9<br />
___________________________________________________________________________<br />
A hat lehetséges oszlopvektor sorrend eredményét egyben, az (5b) felhasználásával felírva<br />
( )<br />
e det A = a a a e . (10b)<br />
rst kr ps qt kpq<br />
Ha ennek az egyenletnek mindkét oldalát megszorozzuk az erst szimbólummal, akkor a (6c)<br />
eredmény helyettesítése után az A tenzor determinánsának kiszámítására alkalmas<br />
összefüggést kapunk:<br />
Elfajuló leképzésnél det(A) = 0.<br />
1<br />
det ( A ) = akrapsaqtekpqerst . (10c)<br />
6<br />
Ha apk az A tenzor mátrixa, akkor akp az A T transzponált mátrixa. Az A szimmetrikus, ha<br />
és ferdén szimmetrikus (aszimmetrikus), ha<br />
pk<br />
kp<br />
T<br />
a = a , A = A , (11a)<br />
pk<br />
kp<br />
T<br />
a = −a<br />
, A = −A<br />
. (11b)<br />
Minden másodrendű tenzor felbontható egy szimmetrikus és egy ferdén szimmetrikus részre:<br />
T T<br />
( ) ( )<br />
1 1<br />
A = A+ A + A− A = E+ Ω<br />
2 2<br />
( ) ( )<br />
1 T T 1<br />
T T<br />
E = A+ A = E , Ω = A− A = −Ω,<br />
2 2<br />
1 1<br />
E pq = (a pq + a qp) , Ω pq = (a pq −a qp) = −Ωqp,<br />
2 2<br />
ahol az E szimmetrikus az Ω pedig aszimmetrikus <strong>tenzorok</strong>. Mivel az aszimmetrikus tenzor<br />
főátlójában lévő elemek értéke mindig zérus, annak csak három zérustól különböző eleme<br />
lehet. Ha ezt a három elemet egy ω vektor három koordinátájának tekintjük, akkor a<br />
permutációs szimbólum (5b) és az aszimmetrikus tenzor (11b) tulajdonságainak hasonlósága<br />
alapján belátható, hogy minden aszimmetrikus tenzor felírható a következő alakban is:<br />
⎡ 0 −ω3 ω2<br />
⎤<br />
Ωpq = ωe i iqp = −ωie ipq , ⎡Ω ⎤ pq =<br />
⎢<br />
3 0<br />
⎥<br />
⎣ ⎦ ⎢<br />
ω −ω1⎥<br />
⎢−ω ω 0 ⎥<br />
⎣ 2 1 ⎦<br />
Ebből az egyenlőségből kifejezhetjük az ωi vektort, ha mindkét oldalt megszorozzuk az ejpq<br />
permutációs szimbólummal. A (6b) felhasználásával<br />
amiből<br />
Ω e = −ω<br />
e e = −ω<br />
2δ ,<br />
pq<br />
jpq<br />
i<br />
ipq<br />
jpq<br />
i<br />
ij<br />
,<br />
.<br />
(12)
10<br />
___________________________________________________________________________<br />
1<br />
ω Ω e<br />
2<br />
j =− pq jpq<br />
, (13)<br />
következik. Az ω vektor az Ω tenzor – illetve a (12) kapcsolat alapján az A tenzor – vektor<br />
invariánsa:<br />
1<br />
ω j =− Apqe jpq ,<br />
2<br />
1 1 1<br />
ω 1 = ( A32 −A 23) , ω 2 = ( A31 −A 13) , ω 3 = ( A21 −A<br />
12)<br />
.<br />
2 2 2<br />
A szimmetrikus E tenzor vektor invariánsa zérus, mivel egyszerűen belátható, hogy<br />
E pq e jpq<br />
=<br />
0 .<br />
Az I egységtenzor mátrixa a (2) Kronecker delta:<br />
( δ e e ) ⋅ ( e ) = δ a δ e = a e a<br />
I a = a = . (14)<br />
⋅ kt k t s s kt s ts k k k<br />
A (9b) definícióban az A másodrendű tenzor, ahol a másodrendű szó a bázis<strong>vektorok</strong> számára<br />
utal. Ebből kiindulva kézenfekvő az általánosítás, például egy negyedrendű tenzor<br />
C = C e e e e<br />
alakban irható fel. Ennek a negyedrendű tenzornak 81 skalár koordinátája van.<br />
2.1. Tenzorok szorzása<br />
pqrs<br />
A <strong>vektorok</strong>nál használatos szorzási műveleteket a (8) definíció és az összegzési konvenció<br />
alapján most már a <strong>tenzorok</strong>ra is kiterjeszthetjük. Ezek közül nézzünk néhány, a továbbiakban<br />
is gyakran előforduló esetet.<br />
Egy tenzor és egy vektor jobb oldali skalár szorzatának eredménye egy vektor:<br />
( a e e ) ⋅ ( b e ) = a b e ( e ⋅ e )<br />
c = A ⋅ b =<br />
= a b δ e = a b e<br />
pk<br />
p<br />
k<br />
s<br />
s<br />
c<br />
p<br />
pk<br />
s<br />
= a<br />
ahol felhasználtuk a Kronecker szimbólum (2)-(3) tulajdonságait. A bal oldali skalár szorzat<br />
d = b ⋅ A = a b e , d = a b .<br />
Ha az A szimmetrikus, azaz asp = aps , akkor c = d.<br />
Egy tenzor és egy vektor jobb oldali vektor szorzatának eredménye egy tenzor:<br />
sp<br />
( a e e ) × ( b e ) = a b e ( e × e )<br />
C = A × b =<br />
= a b e e e<br />
c<br />
s<br />
= a<br />
p<br />
p<br />
ps<br />
p<br />
b<br />
q<br />
s<br />
k<br />
b e<br />
ahol megjelent az (5) szerinti permutációs szimbólum.<br />
pk<br />
p<br />
k<br />
s<br />
pt<br />
s<br />
pk<br />
pk<br />
s<br />
s<br />
,<br />
r<br />
kst<br />
p<br />
p<br />
,<br />
s<br />
s<br />
k<br />
sp<br />
pk<br />
s<br />
s<br />
s<br />
ks<br />
pk<br />
p<br />
s<br />
kst<br />
ps<br />
p<br />
s<br />
t<br />
,<br />
p<br />
,<br />
,
11<br />
___________________________________________________________________________<br />
Két tenzor skalár szorzata egy ugyanolyan méretű tenzor:<br />
( a e e ) ⋅ ( b e e ) = a b e ( e ⋅e<br />
)<br />
C = A ⋅ B =<br />
e = a b δ e e = a b e e<br />
pk<br />
p<br />
k<br />
tm<br />
t<br />
m<br />
pk<br />
c<br />
pm<br />
tm<br />
p<br />
= a<br />
Két tenzor kettős skalár szorzata egy skalár szám:<br />
( a e e ) ⋅⋅(<br />
b e e ) = a b ( e ⋅e<br />
)( e ⋅ )<br />
c = A ⋅⋅B<br />
=<br />
e<br />
pk<br />
p<br />
k<br />
tm<br />
t<br />
m<br />
pk<br />
tm<br />
c = a<br />
Két tenzor vektor szorzata egy harmadrendű tenzor:<br />
( a e e ) × ( b e e ) = a b e ( e × e )<br />
c<br />
= a<br />
pk<br />
pk<br />
b<br />
k<br />
b<br />
b<br />
k<br />
km<br />
kp<br />
e<br />
t<br />
.<br />
.<br />
t<br />
.<br />
m<br />
p<br />
m<br />
pk<br />
tm<br />
= a<br />
pk<br />
kt<br />
b<br />
p<br />
tm<br />
δ<br />
m<br />
kt<br />
δ<br />
pm<br />
pk<br />
= a<br />
C = A × B =<br />
e = a b e e e e<br />
pk<br />
p<br />
A harmadrendű tenzornak 27 skalár koordinátája van.<br />
k<br />
tm<br />
Egy tenzor és két vektor jobb és baloldali skalár szorzata egy skalár szám:<br />
t<br />
m<br />
psm<br />
( ) ( ) ( ) ( ) ( )<br />
pk<br />
pk<br />
c = ⋅ ⋅ = b ⋅ a ⋅ d = b a d ⋅ ⋅ = b a d δ δ ,<br />
b A d rer pkepek ses r pk s er ep ek es r pk s rp ks<br />
3. Transzformációk<br />
E1<br />
E3<br />
E2<br />
tm<br />
tm<br />
p<br />
kts<br />
c = b a d .<br />
a<br />
r rs s<br />
2. ábra<br />
Gyakran szükség van arra, hogy egy vektor vagy egy tenzor koordinátáit különböző<br />
ortogonális koordináta rendszerekben kell megadni. A koordináták átszámítása, vagy<br />
transzformációja a bázis egység<strong>vektorok</strong> kapcsolatából következik.<br />
A 2. ábra szerint jelölje EQ a kiinduló, vagy eredeti, és ep az új rendszer bázis egységvektorait.<br />
Az egyszerűbb megkülönböztetés miatt az eredeti rendszerre vonatkozó koordinátáknál és<br />
indexeknél használjuk a nagy betűs, míg az új rendszerben a kisbetűs jelöléseket.<br />
Természetesen, ez eddig használt összegzési konvenció mindkét rendszerben változatlan<br />
formában érvényes. Induljunk ki abból, hogy ismerjük az új bázis egység<strong>vektorok</strong>at, mint az<br />
eredeti bázis egység<strong>vektorok</strong> lineáris kombinációit:<br />
e1<br />
k<br />
e2<br />
e3<br />
t<br />
m<br />
e = . (15)<br />
p t pQ EQ<br />
pk<br />
tm<br />
kts<br />
p<br />
s<br />
km<br />
pk<br />
m<br />
b<br />
,<br />
p<br />
kp<br />
,<br />
m<br />
,
12<br />
___________________________________________________________________________<br />
Ha ennek az egyenletnek mindkét oldalát skalárisan megszorozzuk az EK vektorral,<br />
figyelembe véve, hogy E Q ⋅ E K = δQK<br />
a (12) szerinti Kronecker szimbólum, a transzformáció<br />
mátrixának kiszámítására alkalmas összefüggést kapunk:<br />
t<br />
pK<br />
⎡e1<br />
⋅ E1<br />
e1<br />
⋅ E 2 e1<br />
⋅ E3<br />
⎤<br />
= e ⋅ [ ] =<br />
⎢<br />
⎥<br />
p E K , t pK ⎢<br />
e 2 ⋅ E1<br />
e 2 ⋅ E 2 e 2 ⋅ E3<br />
⎥<br />
. (16)<br />
⎢⎣<br />
e ⋅ 1 ⋅ 2 ⋅ 3 ⎥<br />
3 E e3<br />
E e3<br />
E ⎦<br />
Hasonló módon, felírhatjuk a bázis egység<strong>vektorok</strong> inverz kapcsolatát is az<br />
−1<br />
( ) q<br />
E K = t e<br />
formában, amit most skalárisan szorozva az ep bázis vektorral, a<br />
eredményt kapjuk. Mivel<br />
Kq<br />
−1<br />
( t ) Kp = EK ⋅e<br />
p<br />
−1<br />
T<br />
T<br />
( t ) = K ⋅e<br />
p = ( e p ⋅ E K ) = ( t pK ) = t Kp<br />
Kp<br />
E ,<br />
megállapíthatjuk, hogy a transzformáció inverze megegyezik a transzponáltjával. Ezért, ha a<br />
transzformációs mátrixot megszorozzuk a transzponáltjával, az eredmény az egységmátrix, a<br />
Kronecker delta lesz:<br />
-1<br />
-1<br />
( t ) t t = δ , ( t ) t = t t = δ .<br />
t qK Kp qK Kp qp<br />
Kp pQ Kp pQ KQ<br />
Ezzel a bázis egység<strong>vektorok</strong> inverz kapcsolata<br />
= (17)<br />
E = . (18)<br />
K t Kq eq<br />
A (16) alapján belátható, hogy a [tKq] mátrix oszlopvektorai az eq bázis egység<strong>vektorok</strong>, ezért<br />
det(tKq) = 1.<br />
A vektor és tenzor mennyiségeket a koordináta rendszer megválasztása nem befolyásolhatja,<br />
más szóval, azok a koordináta transzformációval szemben invariánsak. A transzformáció<br />
során csak a <strong>vektorok</strong> vagy <strong>tenzorok</strong> koordinátái változhatnak:<br />
u = U E = u e , A = A E E = a e e .<br />
P<br />
P<br />
m<br />
m<br />
Az u vektor koordinátáinak kapcsolata a (16) vagy a (18) behelyettesítésével:<br />
u = U E = U t e = u e , u = U t ,<br />
P P P Pm m m m<br />
m P Pm<br />
u = u e = u t E = U E , U = u t ,<br />
m m m mP P P P<br />
P m mP<br />
és hasonló módon az A tenzor koordinátáinak kapcsolata:<br />
( )<br />
PK<br />
a = A t t = t A t<br />
mn PK Pm Kn mP PK Kn<br />
( )<br />
A = a t t = t a t<br />
PK mn mP nK Pm mn nK<br />
T<br />
T<br />
P<br />
K<br />
mn<br />
m<br />
n<br />
(19a)<br />
. (19b)
13<br />
___________________________________________________________________________<br />
4. Tenzor sajátértékei, a skalár invariánsok<br />
A másodrendű tenzor a (9a) definíciója szerint R = A ⋅r<br />
, ahol r a lineáris függvényben a<br />
független változó. Vizsgáljuk meg, hogy mi a feltétele annak, hogy az R és r <strong>vektorok</strong><br />
párhuzamosak legyenek, azaz r = n és R = λr = λI n<br />
( ) ( )<br />
A⋅r− R = A−λ I ⋅ n = 0 ,<br />
(20)<br />
ahol λ egy skalár szám, 0 a zérus vektor és n egyenlőre ismeretlen egységvektor. A (15)<br />
szerint – tekintettel a másodrendű tenzor (9a) alakjára – az n vektor az (A – λI) tenzor<br />
oszlopvektoraira merőleges, vagyis mind a három oszlopvektor az n irányára merőleges<br />
síkban van, pontosabban, egyiknek sincs n irányú vetülete. Tehát az (A – λI) tenzor biztosan<br />
elfajuló és ezért a (10) szerinti determinánsa, ami a három oszlopvektor vegyes szorzata,<br />
zérus:<br />
( A I)<br />
( ) ( ) ( )<br />
det − λ = 0= a −λδ a −λδ a −λδ<br />
e<br />
( ) (<br />
k1 k1 p2 p2 q3 q3 kpq<br />
=− λ e δ δ δ + λ e δ δ a + δ a δ + a δ δ −<br />
3 2<br />
kpq k1 p2 q3 kpq k1 p2 q3 k1 p2 q3 k1 p2 q3<br />
( )<br />
−λ e a a δ + a δ a + δ a a + e a a a .<br />
kpq k1 p2 q3 k1 p2 q3 k1 p2 q3 kpq k1 p2 q3<br />
A (2), (5b) és (10) definíciókból következő egyszerűsítések után az A tenzor karakterisztikus<br />
egyenlete<br />
2<br />
− A λ + A λ − A = 0<br />
(21)<br />
3<br />
λ I II III<br />
ahol a polinom együtthatói az A tenzor skalár invariánsai. Az első invariáns a főátlóban lévő<br />
elemek összege, a második a főátlóhoz tartozó aldeterminánsok összege, a harmadik pedig az<br />
A determinánsa:<br />
A<br />
A<br />
A<br />
1<br />
II<br />
III<br />
= a<br />
= e<br />
=<br />
kk<br />
(a<br />
= e<br />
1pq<br />
= a<br />
22<br />
kpq<br />
a<br />
a<br />
a<br />
p2<br />
33<br />
k1<br />
11<br />
a<br />
- a<br />
a<br />
+ a<br />
q3<br />
p2<br />
+ e<br />
32<br />
a<br />
22<br />
a<br />
q3<br />
+ a<br />
k2q<br />
23<br />
a<br />
33<br />
k1<br />
a<br />
,<br />
) + (a<br />
q3<br />
11<br />
a<br />
= det(<br />
A).<br />
+ e<br />
33<br />
3kp<br />
- a<br />
a<br />
31<br />
k1<br />
a<br />
a<br />
13<br />
p2<br />
=<br />
) + (a<br />
11<br />
a<br />
22<br />
- a<br />
21<br />
a<br />
12<br />
) ,<br />
)<br />
=<br />
(22a)
14<br />
___________________________________________________________________________<br />
Az invariáns megnevezés arra utal, hogy egy tenzornak ezek a skalár adatai az (19)<br />
transzformáció során nem változnak. Például az első skalár invariáns:<br />
( ) T<br />
A = a = A t t = A t t = A δ = A K ,<br />
I mm PK Pm Km PK mP Km PK PK K<br />
mivel a transzformációs mátrix transzponáltja megegyezik az inverzével. Hasonló módon<br />
igazolható az AII és AIII és a (13) szerinti ω vektor invariáns tulajdonsága.<br />
A (21) polinom λ1, λ2, és λ3 gyökei az A tenzor sajátértékei. (Tételezzük fel, hogy az A olyan<br />
tulajdonságú, hogy ezek valós számok.) Ha a sajátértékeket rendre visszahelyettesítjük az (20)<br />
egyenletbe, akkor az<br />
( A − λ I)<br />
⋅n<br />
= 0 , i 1, 2, 3,<br />
i =<br />
homogén lineáris egyenletek megoldásával meghatározhatjuk a három n1, n2, n3 sajátvektort.<br />
Igazolható, hogy szimmetrikus <strong>tenzorok</strong> sajátvektorai egymásra merőlegesek. Az ortogonális<br />
saját<strong>vektorok</strong> (egység<strong>vektorok</strong>) lehetnek egy koordináta rendszer bázis vektorai. Ha a<br />
szimmetrikus A tenzort a (19) formulák szerint ebbe a rendszerbe transzformáljuk, akkor ott a<br />
mátrixa<br />
A = a n n<br />
pq<br />
p<br />
q<br />
,<br />
[ a ]<br />
pq<br />
⎡λ1<br />
=<br />
⎢<br />
⎢<br />
0<br />
⎢⎣<br />
0<br />
diagonál alakú lesz. A szimmetrikus tenzor (22a) skalár invariánsai a sajátértékekkel kifejezve<br />
a következők:<br />
5. Görbevonalú koordináta rendszerek<br />
A<br />
A<br />
A<br />
1<br />
II<br />
III<br />
= a<br />
kk<br />
1<br />
= λ<br />
2<br />
1<br />
3<br />
+ λ<br />
2<br />
0<br />
λ<br />
0<br />
= λ 2λ<br />
3 + λ 2λ1<br />
+ λ1λ<br />
= λ λ λ = det(<br />
A).<br />
2<br />
+ λ<br />
3<br />
0 ⎤<br />
0<br />
⎥<br />
⎥<br />
λ ⎥ 3 ⎦<br />
,<br />
2<br />
,<br />
(22b)<br />
A tér egy pontjának helyét az xp koordináták mellett még más, célszerűen megválasztott<br />
számhármassal is megadhatjuk. Például, ha hengerkoordináták használunk: q1 a sugár, q2 a<br />
szög és q3 az alkotó mentén mért hosszúság. Tételezzük fel, hogy ismerjük az xp Descartes<br />
féle és a qk általános koordináták kapcsolatát leíró xp(qk) függvényeket és ezek qk(xp)<br />
inverzeit is. Egy tetszőleges, de rögzített q1 = a, q2 = b, q3 = c koordinátájú P ponton áthaladó<br />
koordináta görbék egy paraméteres vektoregyenletei a következők:
15<br />
___________________________________________________________________________<br />
E1<br />
r(q<br />
, q<br />
r(q<br />
r(q<br />
1<br />
1<br />
1<br />
=<br />
=<br />
E2<br />
2<br />
=<br />
a, q<br />
a, q<br />
E3 r<br />
b, q<br />
2<br />
2<br />
, q<br />
=<br />
3<br />
3<br />
= c) = r<br />
= c) = r<br />
b, q<br />
P<br />
3<br />
q1<br />
1<br />
2<br />
( q )<br />
,<br />
( q )<br />
,<br />
) = r ( q ).<br />
e1<br />
3<br />
e3<br />
1<br />
2<br />
e2<br />
3<br />
q3<br />
q2<br />
3. ábra Görbevonalú koordináták<br />
A P pontbeli lokális görbe vonalú koordináta rendszer bázisvektorai legyenek a P ponton<br />
átmenő koordináta vonalak érintő egységvektorai. A továbbiakban csak olyan görbe vonalú<br />
rendszereket használunk, amelyeknek az ily módon megszerkesztett bázisvektorai<br />
ortogonálisak. A p-edik (p = 1, 2, 3) koordináta görbe érintő vektora<br />
g<br />
∂r<br />
∂x<br />
= = e , (23)<br />
p ∂q<br />
p<br />
k<br />
∂q<br />
p<br />
aminek a hossza, az abszolút értéke a Hp Lamé tényező:<br />
p<br />
p<br />
p<br />
k<br />
g = + g ⋅ g = H . (24)<br />
Az index aláhúzása itt arra utal, hogy a gyökjel alatt nem egy összeg szerepel, a két p index<br />
nem egy néma indexpár, mert az egyiket ”kizártuk”. Ezzel az ortogonális görbe vonalú<br />
rendszer bázis egységvektorai a következő módon határozhatók meg.<br />
p<br />
1<br />
e p = g , p = 1, 2, 3.<br />
(25)<br />
p H<br />
p<br />
A bázis egység<strong>vektorok</strong> iránya pontról pontra változhat. Ha valamelyik qp általános<br />
koordináta a koordináta görbe mentén mért ívhossz, akkor a megfelelő érintő vektor<br />
egységvektor, azaz gp = ep és Hp = 1.
16<br />
___________________________________________________________________________<br />
5.1. A kovariáns derivált<br />
Tenzor és vektormennyiségek használata során gyakran meg kell határozni ezeknek a<br />
helykoordináták szerinti megváltozását, differenciálját. Lényeges, hogy görbe vonalú<br />
koordináta rendszerek használata esetén a vektor vagy tenzor differenciálja nem csak a<br />
koordináták, hanem a bázis egységvektor irányok változását is tartalmazza.<br />
A bázis egység<strong>vektorok</strong>nak a qk koordináták szerinti deriváltjait, amik szintén <strong>vektorok</strong><br />
lesznek, a következő alakban írjuk fel:<br />
∂e<br />
∂q<br />
p<br />
k<br />
⎛ s ⎞<br />
= ⎜ ⎟ e<br />
⎝p<br />
k⎠<br />
s<br />
,<br />
⎛ s ⎞ ⎛ ∂e<br />
⎜ ⎟ =<br />
p k ⎜<br />
⎝ ⎠ ⎝ ∂q<br />
p<br />
k<br />
⎞<br />
⎟ ⋅ e<br />
⎠<br />
s<br />
, (26)<br />
ahol az alsó két index jelöli, hogy honnan származik a derivált egységvektor egy felső index<br />
irányú koordinátája.<br />
Az r(qk) helyvektor dr megváltozása vagy differenciálja<br />
∂r<br />
r g e , (27)<br />
d dqs dq s s Hs<br />
s s<br />
qs<br />
q d<br />
= = =<br />
∂<br />
ahol felhasználtuk a (23) és a (25) összefüggéseket is.<br />
Ezek alapján számítsuk ki egy F(r) = F(qk) skalár függvény megváltozását, amit formális<br />
átalakítások után két vektor skalár szorzataként írhatunk fel:<br />
ahol<br />
d F ∂F<br />
∂F<br />
H s ⎛ ⎞<br />
d d r d<br />
d ⎜ 1 ∂F<br />
F = ⋅ = q =<br />
= e ⎟<br />
k δ ks q s<br />
k ⋅ s s s F ⋅<br />
d r ∂q<br />
∂q<br />
H ⎜ ⎟<br />
k<br />
k k ⎝<br />
H k ∂q<br />
k ⎠<br />
d F<br />
=<br />
d r<br />
( ∇F)<br />
= grad(F)<br />
( H e d q ) = ( ∇ ) d r<br />
az F függvény gradiens vektora és ennek megfelelően a Hamilton féle, vektor értékű<br />
differenciál operátor definíciója a következő:<br />
d<br />
∇ = = e<br />
d r<br />
k<br />
1<br />
H<br />
k<br />
∂<br />
∂q<br />
k<br />
. (28)<br />
A ∇ szimbólum neve nabla. A következőkben határozzuk meg egy u(r) = u(qk) vektor mező<br />
differenciálját. Az eddigiek alapján:<br />
d u ∂u<br />
d u = ⋅ d r = d q<br />
d r ∂q<br />
k<br />
k<br />
∂u<br />
= δ<br />
∂q<br />
=<br />
H<br />
H<br />
d q<br />
( u∇)<br />
⋅ d r = D ⋅ d r<br />
⎛<br />
= ⎜<br />
⎜<br />
⎝<br />
1<br />
H<br />
∂u<br />
∂q<br />
e<br />
⎞<br />
⎟ ⋅<br />
⎟<br />
⎠<br />
( H e dq<br />
)<br />
ahol D az u vektor derivált tenzora, vagy gradiense, amit tovább részletezve:<br />
k<br />
ks<br />
s<br />
k<br />
s<br />
k<br />
k<br />
k<br />
s<br />
s<br />
s<br />
=<br />
,<br />
, (29)
17<br />
___________________________________________________________________________<br />
D = u∇<br />
=<br />
=<br />
1<br />
H<br />
k<br />
1<br />
H<br />
k<br />
∂<br />
⎛ ∂u<br />
⎜<br />
⎝ ∂q<br />
( u e )<br />
s<br />
k<br />
∂q<br />
δ<br />
s<br />
st<br />
k<br />
s<br />
e<br />
t<br />
e<br />
k<br />
=<br />
1<br />
H<br />
⎛ ∂u<br />
⎜<br />
⎝ ∂q<br />
⎛ t ⎞ ⎞<br />
+ u s ⎜ ⎟e<br />
⎟ t e<br />
s k ⎟<br />
⎝ ⎠ ⎠<br />
k<br />
s<br />
k<br />
k<br />
e<br />
s<br />
=<br />
+ u<br />
1<br />
H<br />
k<br />
s<br />
∂e<br />
∂q<br />
s<br />
k<br />
⎛ ∂u<br />
⎜<br />
⎝ ∂q<br />
⎞<br />
⎟ e<br />
⎠<br />
t<br />
k<br />
k<br />
=<br />
1<br />
H<br />
⎛ ∂u<br />
⎜<br />
⎝ ∂q<br />
⎛ t ⎞⎞<br />
+ u s ⎜ ⎟⎟e<br />
t e<br />
s k ⎟<br />
⎝ ⎠⎠<br />
k<br />
s<br />
k<br />
k<br />
e<br />
s<br />
⎛ t ⎞ ⎞<br />
+ u s ⎜ ⎟e<br />
⎟ t e<br />
s k ⎟<br />
⎝ ⎠ ⎠<br />
A parciális deriváltra és a derivált tenzor koordinátáira gyakran használjuk a következő<br />
egyszerű jelöléseket:<br />
∂u<br />
∂q<br />
t =<br />
k<br />
u<br />
t , k<br />
,<br />
= D<br />
ahol ut,k a parciális derivált (a koordináta és a deriválási index között vessző van),<br />
tk<br />
e<br />
t<br />
e<br />
k<br />
.<br />
k<br />
=<br />
(30)<br />
1<br />
t<br />
D tk u t ; k u t , k u s ,<br />
H k s k ⎟ ⎛ ⎛ ⎞⎞<br />
= = ⎜<br />
+ ⎜ ⎟<br />
(31)<br />
⎝ ⎝ ⎠⎠<br />
az ut;k pedig a kovariáns derivált (a koordináta és a deriválási index között pontosvessző van).<br />
Látszik, hogy a parciális derivált és a kovariáns derivált közötti különbség a bázis<br />
egység<strong>vektorok</strong> hely szerinti megváltozásából származik. Egyenes vonalú derékszögű<br />
koordináta rendszerben nyilvánvalóan ut,k = ut;k .<br />
Két Hamilton operátor skaláris szorzata a Laplace operátor:<br />
∆ = ∇ ⋅ ∇<br />
5.2. Differenciálási szabályok<br />
⎡ 1 ∂ ⎛ ⎞⎤<br />
⎢ ⎜ 1 ∂<br />
= e ⎟<br />
k ⋅ e t ⎥ . (32)<br />
⎢H<br />
∂ ⎜ ⎟<br />
⎣ k q k ⎝<br />
H t ∂q<br />
t ⎠⎥⎦<br />
A Hamilton és a Laplace féle differenciál operátorok (28), (32) alakjainál ismét<br />
megmutatkozik az indexes jelölés használatának előnye, mivel formálisan szétválik a<br />
koordináta – a műveleti utasítás – és a bázisvektor. A <strong>vektorok</strong> és <strong>tenzorok</strong> körében<br />
használatos szorzási műveleteket az összegzési konvenció alapján most már a deriválási<br />
műveletekre is kiterjeszthetjük. Ezek közül nézzünk meg részletesebben néhány, a<br />
továbbiakban gyakran előforduló esetet.<br />
Az F(qk) skalár függvény gradiens vektora<br />
1 ∂F<br />
1 1 1<br />
G = grad + e . (33)<br />
F = ∇F<br />
= e k = e1<br />
F,<br />
1 + e 2 F,<br />
2 3 F,<br />
3<br />
H k ∂q<br />
k H 1 H 2 H 3<br />
Az u(qk) jobb oldali gradiense – a vektor és a differenciál operátor (8) diadikus szorzata - a<br />
derivált vagy gradiens tenzor, ami a (30) alapján:
18<br />
___________________________________________________________________________<br />
D = u∇<br />
=<br />
1<br />
H<br />
k<br />
∂<br />
( u e )<br />
∂q<br />
D<br />
s<br />
tk<br />
k<br />
s<br />
=<br />
e<br />
k<br />
1<br />
H<br />
k<br />
=<br />
⎛<br />
⎜<br />
u<br />
⎝<br />
1<br />
H<br />
k<br />
t , k<br />
⎛ ∂u<br />
⎜<br />
⎝ ∂q<br />
s<br />
k<br />
e<br />
s<br />
+ u<br />
∂q<br />
⎛ t ⎞⎞<br />
+ u s ⎜ ⎟⎟<br />
= u t<br />
s k ⎟<br />
⎝ ⎠⎠<br />
s<br />
∂e<br />
s<br />
k<br />
; k<br />
⎞<br />
⎟ e<br />
⎠<br />
.<br />
k<br />
= D<br />
A vektor bal oldali gradiense a derivált tenzor transzponáltja: ∇u = D T .<br />
Az u vektor divergenciája a vektor és a differenciál operátor skalár szorzata:<br />
F = divu<br />
= ∇ ⋅ u =<br />
=<br />
1<br />
H<br />
k<br />
1<br />
H<br />
k<br />
∂<br />
⎛ ∂u<br />
⎜<br />
⎝ ∂q<br />
s<br />
k<br />
( u e )<br />
s<br />
∂q<br />
δ<br />
sk<br />
k<br />
s<br />
⋅ e<br />
k<br />
=<br />
1<br />
H<br />
⎛ ∂u<br />
⎜<br />
⎝ ∂q<br />
⎛ t ⎞ ⎞<br />
+ u s ⎜ ⎟δ<br />
⎟ tk ,<br />
s k ⎟<br />
⎝ ⎠ ⎠<br />
k<br />
s<br />
k<br />
e<br />
s<br />
tk<br />
e<br />
t<br />
e<br />
⎛ t ⎞ ⎞<br />
+ u s ⎜ ⎟e<br />
⎟ t ⋅<br />
s k ⎟<br />
e<br />
⎝ ⎠ ⎠<br />
k<br />
,<br />
k<br />
=<br />
(34)<br />
1 ⎛ ⎛ k ⎞⎞<br />
= ⎜u<br />
u ⎟<br />
k , k s = u k ; .<br />
H ⎜<br />
+ ⎜ ⎟<br />
k s k ⎟<br />
(35)<br />
⎝ ⎝ ⎠⎠<br />
F k<br />
Az u vektor rotációja a vektor és a differenciál operátor vektor szorzata:<br />
=<br />
c = rot u = u × ∇ =<br />
1<br />
H<br />
k<br />
⎛ ∂u<br />
⎜<br />
⎝ ∂q<br />
s<br />
k<br />
e<br />
skp<br />
e<br />
p<br />
1<br />
H<br />
k<br />
∂<br />
( u e )<br />
∂q<br />
k<br />
⎛ t ⎞<br />
+ u s ⎜ ⎟e<br />
⎝s<br />
k⎠<br />
s<br />
s<br />
tkp<br />
× e<br />
e<br />
p<br />
k<br />
⎞<br />
⎟<br />
=<br />
⎠<br />
=<br />
1<br />
H<br />
1<br />
H<br />
k<br />
k<br />
⎛ ∂u<br />
⎜<br />
⎝ ∂q<br />
⎛ ∂u<br />
⎜<br />
⎝ ∂q<br />
t<br />
c p<br />
u s e tkp<br />
H ⎜<br />
=<br />
k ∂q<br />
k s k ⎟<br />
t<br />
k<br />
s<br />
k<br />
e<br />
s<br />
⎛ t ⎞ ⎞<br />
+ u s ⎜ ⎟e<br />
⎟ t × e<br />
s k ⎟<br />
⎝ ⎠ ⎠<br />
⎛ t ⎞⎞<br />
+ u ⎟<br />
s ⎜ ⎟ e<br />
s k ⎟<br />
⎝ ⎠⎠<br />
tkp<br />
e<br />
p<br />
k<br />
= c e<br />
1 ⎛ ∂u<br />
⎛ t ⎞⎞<br />
= ⎜ + ⎜ ⎟⎟<br />
u t ; ke<br />
tkp . (36)<br />
⎝ ⎝ ⎠⎠<br />
Az A tenzor jobb oldali divergenciája az operátor és a tenzor skalár szorzata:<br />
c = A ⋅ ∇ =<br />
1<br />
H<br />
=<br />
k<br />
∂<br />
1<br />
H<br />
( a e e )<br />
k<br />
sq<br />
∂q<br />
s<br />
k<br />
⎛ ∂a<br />
⎜<br />
⎝ ∂q<br />
sq<br />
k<br />
q<br />
e δ<br />
s<br />
⋅e<br />
qk<br />
k<br />
=<br />
+ a<br />
1<br />
H<br />
sq<br />
k<br />
⎛ ∂a<br />
⎜<br />
⎝ ∂q<br />
sq<br />
k<br />
e<br />
⎛ t ⎞<br />
⎜ ⎟e<br />
tδ<br />
⎝s<br />
k⎠<br />
s<br />
qk<br />
e<br />
q<br />
+ a<br />
+ a<br />
sq<br />
sq<br />
⎛ t ⎞<br />
⎜ ⎟e<br />
t e<br />
⎝s<br />
k⎠<br />
⎛ t ⎞<br />
⎜ ⎟esδ<br />
⎝q<br />
k⎠<br />
tk<br />
q<br />
+ a<br />
sq<br />
⎞<br />
⎟<br />
= c t e<br />
⎠<br />
t<br />
,<br />
p<br />
=<br />
⎛ t ⎞ ⎞<br />
⎜ ⎟es<br />
e ⎟ t ⋅e<br />
q k ⎟<br />
⎝ ⎠ ⎠<br />
1 ⎛ ⎛ t ⎞ ⎛ k ⎞⎞<br />
c t = ⎜a<br />
tk , k a sk a ⎟<br />
tq = a tk ; k<br />
H ⎜<br />
+ ⎜ ⎟ + ⎜ ⎟<br />
k<br />
s k q k ⎟<br />
. (37)<br />
⎝ ⎝ ⎠ ⎝ ⎠⎠<br />
p<br />
,<br />
k<br />
=
19<br />
___________________________________________________________________________<br />
Az F(qk) skalár függvény (32) szerinti második deriváltja:<br />
∆F<br />
= ∇ ⋅ ∇F<br />
= div<br />
=<br />
1<br />
H<br />
s<br />
∂<br />
∂q<br />
5.3. Hengerkoordináta rendszer<br />
s<br />
⎛<br />
⎜<br />
⎜<br />
⎝<br />
( grad F)<br />
1<br />
H<br />
s<br />
=<br />
∂F<br />
⎞<br />
⎟ +<br />
∂q<br />
⎟<br />
s ⎠<br />
1<br />
H<br />
1<br />
H<br />
k<br />
s<br />
∂<br />
∂q<br />
1<br />
H<br />
s<br />
s<br />
⎛<br />
⎜<br />
⎜<br />
⎝<br />
1<br />
H<br />
k<br />
∂F<br />
e<br />
∂q<br />
k<br />
∂F<br />
⎛ s ⎞<br />
⎜ ⎟<br />
∂q<br />
s ⎝k<br />
s⎠<br />
k<br />
⎞<br />
⎟ ⋅ e<br />
⎟<br />
⎠<br />
s<br />
=<br />
. (38)<br />
Az előző két fejezetben tárgyalt összefüggéseket írjuk fel a gyakorlatban legtöbbször használt<br />
hengerkoordináta rendszerben. Az általános koordináták legyenek a 4. ábra szerinti r sugár, a<br />
φ szög és a z alkotó menti magasság.<br />
A Descartes féle és a görbe vonalú koordináták xp(qk) kapcsolata:<br />
∂x<br />
=<br />
∂q<br />
= q cos q = r cos ϕ , x = q sin q = r sin ϕ , x = q = z .<br />
x1 1 2<br />
2 1 2<br />
3 3<br />
A (23) érintő <strong>vektorok</strong>:<br />
∂x<br />
g =<br />
1 ∂q<br />
g<br />
2<br />
k<br />
2<br />
E<br />
k<br />
∂x<br />
k<br />
= E<br />
∂ϕ<br />
k<br />
k<br />
1<br />
E<br />
k<br />
E3<br />
E1<br />
r<br />
E2<br />
φ<br />
e3<br />
z<br />
P<br />
e2<br />
e1<br />
4. ábra Henger koordináták<br />
∂x<br />
=<br />
∂r<br />
k<br />
= −r<br />
sinϕ<br />
E<br />
E<br />
1<br />
k<br />
= cosϕ<br />
E<br />
+ r cosϕ<br />
E<br />
és a gk <strong>vektorok</strong> abszolút értékei, a (24) Lamé tényezők:<br />
H<br />
1<br />
2<br />
1<br />
+ sinϕ<br />
E<br />
,<br />
g<br />
3<br />
2<br />
,<br />
∂x<br />
=<br />
∂q<br />
= g<br />
= 1,<br />
H = g = r , H = g<br />
1<br />
2<br />
2<br />
3<br />
3<br />
k<br />
3<br />
E<br />
k<br />
= 1 .<br />
∂x<br />
k<br />
= E<br />
∂z<br />
k<br />
= E<br />
3<br />
,
20<br />
___________________________________________________________________________<br />
A (25) bázis egység<strong>vektorok</strong>:<br />
1<br />
1<br />
1<br />
e 1 = g = cos ϕ E1<br />
+ sin ϕ E2<br />
, e 2 = g = − sin ϕ E1<br />
+ cos ϕ E 2 , e<br />
1<br />
2<br />
3 = g = E<br />
3<br />
H<br />
H<br />
H<br />
1<br />
2<br />
A bázis egység<strong>vektorok</strong> (26) deriváltjai, figyelembe véve, hogy azok most csak a q2 = φ<br />
koordináta függvényei:<br />
∂e<br />
∂q<br />
∂e<br />
∂q<br />
1<br />
2<br />
2<br />
2<br />
⎛ k ⎞<br />
= ⎜ ⎟ e<br />
⎝12⎠<br />
k<br />
⎛ k ⎞<br />
= ⎜ ⎟ e<br />
⎝2<br />
2⎠<br />
k<br />
= − sin ϕ E<br />
1<br />
= − cos ϕ E<br />
+ cos ϕ E<br />
1<br />
2<br />
− sin ϕ E<br />
2<br />
= e<br />
2<br />
= −e<br />
Tehát a derivált bázisvektor koordináták zérustól különböző elemei:<br />
⎛ 2 ⎞ ⎛ 1 ⎞<br />
⎜ ⎟ = 1 , ⎜ ⎟ = −1<br />
.<br />
⎝12⎠<br />
⎝2<br />
2⎠<br />
A koordináta rendszer jellemzőinek ismeretében fel lehet írni a különböző mennyiségek<br />
deriváltjait.<br />
Az F(r, φ, z) skalár függvény (33) gradiens vektora<br />
1 ∂F<br />
1 ∂F<br />
∂F<br />
1<br />
∇F = e k F,<br />
k = e1<br />
+ e 2 + e3<br />
= e1F,<br />
r + e 2 F,<br />
ϕ + e3F,<br />
z .<br />
H ∂r<br />
r ∂ϕ<br />
∂z<br />
r<br />
A (34) szerinti D = u∇<br />
derivált tenzornak a<br />
k<br />
1 ⎛ ⎛ t ⎞⎞<br />
D tk = ⎜u<br />
t k u ⎟<br />
, s = u t ; k<br />
H ⎜<br />
+ ⎜ ⎟<br />
k s k ⎟<br />
⎝ ⎝ ⎠⎠<br />
mátrixából részletezzük például az első sor második elem számítását:<br />
1 ⎛ ⎛ 1 ⎞⎞<br />
1 ⎛ ∂u1<br />
⎞<br />
D ⎜<br />
⎟<br />
12 =<br />
⎜<br />
u1<br />
, 2 + u s ⎜ ⎟<br />
⎟<br />
= ⎜ − u 2 ⎟<br />
H 2 ⎝ ⎝s<br />
2⎠⎠<br />
r ⎝ ∂ϕ<br />
⎠<br />
.<br />
Hasonló módon meghatározható a gradiens tenzor mátrixának többi eleme is:<br />
[ D ]<br />
tk<br />
⎡ ∂u<br />
⎢<br />
⎢<br />
∂r<br />
⎢∂u<br />
= ⎢<br />
⎢<br />
∂r<br />
⎢∂u<br />
⎢<br />
⎣<br />
∂r<br />
1<br />
2<br />
3<br />
1 ⎛ ∂u1<br />
⎞<br />
⎜ − u 2 ⎟<br />
r ⎝ ∂ϕ<br />
⎠<br />
1 ⎛ ∂u<br />
2 ⎞<br />
⎜ + u1<br />
⎟<br />
r ⎝ ∂ϕ<br />
⎠<br />
1 ∂u<br />
3<br />
r ∂ϕ<br />
∂u<br />
⎤ 1<br />
⎥<br />
∂z<br />
⎥<br />
∂u<br />
2 ⎥<br />
∂z<br />
⎥<br />
⎥<br />
∂u<br />
3 ⎥<br />
∂z<br />
⎥<br />
⎦<br />
,<br />
1<br />
.<br />
3<br />
3<br />
,
21<br />
___________________________________________________________________________<br />
Az u vektor (35) divergenciája:<br />
1 ⎛ ⎛ k ⎞⎞<br />
∇ ⋅ u = ⎜u<br />
k k u ⎟<br />
s = u k<br />
H ⎜ , + ⎜ ⎟<br />
k s k ⎟ ;<br />
⎝ ⎝ ⎠⎠<br />
k<br />
∂u<br />
=<br />
∂r<br />
1<br />
1 ∂u<br />
2 ∂u<br />
3 1<br />
+ + + u1<br />
r ∂ϕ<br />
∂z<br />
r<br />
a gradiens tenzor mátrixából a főátlóban lévő elemek összege. Ez a tenzornak a (22a) szerint<br />
az első skalár invariánsa.<br />
Az A tenzor c = A ⋅∇<br />
jobb oldali divergenciájának (37) koordinátái:<br />
c<br />
=<br />
1<br />
=<br />
1<br />
H<br />
k<br />
1<br />
H<br />
k<br />
⎛<br />
⎜<br />
a<br />
⎝<br />
⎛<br />
⎜<br />
a<br />
⎝<br />
1k , k<br />
1k , k<br />
+ a<br />
+ a<br />
1k<br />
1 ⎛ ⎛ t ⎞ ⎛ k ⎞⎞<br />
c t = ⎜a<br />
tk , k a sk a ⎟<br />
tq = a tk ; k<br />
H ⎜<br />
+ ⎜ ⎟ + ⎜ ⎟<br />
k<br />
s k q k ⎟<br />
,<br />
⎝ ⎝ ⎠ ⎝ ⎠⎠<br />
sk<br />
⎛ 1 ⎞<br />
⎜ ⎟ + a<br />
⎝s<br />
k ⎠<br />
⎛ 1 ⎞<br />
⎜ ⎟ + a<br />
⎝1<br />
k ⎠<br />
c<br />
c<br />
2<br />
3<br />
2k<br />
1q<br />
∂a<br />
=<br />
∂r<br />
∂a<br />
=<br />
∂r<br />
⎛ k ⎞⎞<br />
⎜ ⎟⎟<br />
=<br />
q k ⎟<br />
⎝ ⎠⎠<br />
⎛ 1 ⎞<br />
⎜ ⎟ + a<br />
⎝2<br />
k ⎠<br />
11<br />
21<br />
∂a<br />
=<br />
∂r<br />
31<br />
3k<br />
⎛ 1 ⎞<br />
⎜ ⎟ + a<br />
⎝3<br />
k ⎠<br />
1 ∂a12<br />
∂a<br />
+ +<br />
r ∂ϕ<br />
∂z<br />
1 ∂a32<br />
∂a<br />
+ +<br />
r ∂ϕ<br />
∂z<br />
13<br />
1 ∂a22<br />
∂a<br />
+ +<br />
r ∂ϕ<br />
∂z<br />
23<br />
33<br />
11<br />
1<br />
− a<br />
r<br />
1<br />
+ a<br />
r<br />
1<br />
+ a<br />
r<br />
Az F(r, φ, z) skalár függvény (38) szerinti második deriváltja:<br />
⎛ k ⎞<br />
⎜ ⎟ + a<br />
⎝1<br />
k ⎠<br />
22<br />
12<br />
31<br />
1<br />
+ a<br />
r<br />
1<br />
+ a<br />
r<br />
.<br />
11<br />
21<br />
12<br />
,<br />
,<br />
⎛ k ⎞<br />
⎜ ⎟ + a<br />
⎝2<br />
k ⎠<br />
13<br />
⎛ k ⎞⎞<br />
⎜ ⎟⎟<br />
=<br />
3 k ⎟<br />
⎝ ⎠⎠<br />
2<br />
2 2<br />
1 ∂ ⎛ 1 F ⎞ 1 1 F s F 1 F F 1 F<br />
F ⎜ ∂<br />
∂ ⎛ ⎞ ∂ ∂ ∂ ∂<br />
∆ =<br />
⎟ +<br />
⎜ ⎟ = + + + .<br />
2 2 2 2<br />
H s ∂q<br />
⎜<br />
s H s q ⎟<br />
⎝<br />
∂ s ⎠<br />
H k H s ∂q<br />
s ⎝k<br />
s⎠<br />
∂r<br />
r ∂ϕ<br />
∂z<br />
r ∂r
22<br />
___________________________________________________________________________<br />
6. Ajánlott irodalom<br />
Béda Gyula: Kontinuum mechanika I. Tankönyvkiadó, 1980.<br />
Uj József: Kontinuum mechanika példatár I. Tankönyvkiadó 1987.<br />
Elter Pálné - Vörös Gábor: Alkalmazott mechanika I. Tankönyvkiadó 1980.<br />
Lánczos Kornél: A geometriai térfogalom fejlődése. Gondolat Könyvkiadó 1976.<br />
Rózsa Pál: Lineáris algebra és alkalmazásai. Műszaki Könyvkiadó 1976.<br />
Simmonds, J.G: Tenzoranalízis dióhéjban. Műszaki Könyvkiadó 1985.