c:/Documents and Settings/dock/Desktop/Matematicka fizika ...rgf.bg.ac.rs/predmet/GO/IV...

18
Glava 1 Teorija polja 1.1 Primeri nekih polja od interesa za fiziku i tehniku Navedimo neke primere potencijalnih polja koja su od posebnog interesa u raznim oblastima fizike i tehnike. Privlaˇ cenje dve taˇ cke u polju gravitacione sile Njutnova sila gravitacije definisana je izrazom F = γ mm 0 R 2 r, (1.1) gde su m i m 0 mase koje se privlaˇ ce, a γ je univerzalna gravitaciona konstanta. U gornjoj relaciji uveli smo slede´ ce oznake (sl. 1.1): r = R R , gde je R = −−−→ M 0 M. (1.2) Slika 1.1: Pitanje je da li postoji potencijal za ovako definisanu silu? Kako su mase m i m 0 , kao i γ konstantne veliˇ cine, to moˇ zemo da ih zamenimo jednom konstantom, recimo c, pa silu moˇ zemo da napiˇ semo u obliku F = c R 2 R R . (1.3) Dalje, ako postoji potencijal U , tada moraju da vaˇ ze slede´ ce relacije, koje povezuju silu F (F = X i+Y j+Z k) 1

Transcript of c:/Documents and Settings/dock/Desktop/Matematicka fizika ...rgf.bg.ac.rs/predmet/GO/IV...

Page 1: c:/Documents and Settings/dock/Desktop/Matematicka fizika ...rgf.bg.ac.rs/predmet/GO/IV semestar/Metode matematicke fizike... · Glava 1 Teorija polja 1.1 Primeri nekih polja od interesa

Glava 1

Teorija polja

1.1 Primeri nekih polja od interesa za fiziku

i tehniku

Navedimo neke primere potencijalnih polja koja su od posebnog interesa u raznim oblastima fizike itehnike.

Privlacenje dve tacke u polju gravitacione sile

Njutnova sila gravitacije definisana je izrazom

F = −γm m0

R2r, (1.1)

gde su m i m0 mase koje se privlace, a γ je univerzalna gravitaciona konstanta.U gornjoj relaciji uveli smo sledece oznake (sl. 1.1):

r =R

R, gde je R =

−−−→M0M. (1.2)

Slika 1.1:

Pitanje je da li postoji potencijal za ovako definisanu silu? Kako su mase m i m0, kao i γ konstantnevelicine, to mozemo da ih zamenimo jednom konstantom, recimo c, pa silu mozemo da napisemo u obliku

F = −c

R2

R

R. (1.3)

Dalje, ako postoji potencijal U , tada moraju da vaze sledece relacije, koje povezuju silu F (F = Xi+Y j+Zk)

1

Page 2: c:/Documents and Settings/dock/Desktop/Matematicka fizika ...rgf.bg.ac.rs/predmet/GO/IV semestar/Metode matematicke fizike... · Glava 1 Teorija polja 1.1 Primeri nekih polja od interesa

i potencijal U (F=gradU):

X =∂U

∂x= −

c

R3(x − x0),

Y =∂U

∂y= −

c

R3(y − y0), (1.4)

Z =∂U

∂z= −

c

R3(z − z0).

Bez gubljenja u opstosti uzmimo da su koordinate tacke M0 (0,0,0), tj.

x0 = y0 = z0 = 0.

Dakle, iz (1.4):∂U

∂x= −

c

R3x ⇒ U = −c

x

R3dx + f(y, z). (1.5)

Dalje, kako je:

R2 = x2 + y2 + z2 i odavde 2xdx = 2R dR, (y i z su fiksirani)

to dobijamo za potencijal U

U =c

R+ f(y, z). (1.6)

Da bismo odredili funkciju f , prvo na -dimo parcijalni izvod U po y:

∂U

∂y=

∂U

∂R

∂R

∂y+

∂f

∂y= −

c

R2

∂R

∂y+

∂f

∂y= −

c

R2

y

R+

∂f

∂y. (1.7)

Da bi bio zadovoljen uslov za Y (1.4) potrebno je da je ∂f∂y = 0, odakle zakljucujemo da f ne zavisi od y,

tj. f = f(z). Na -dimo sada parcijalni izvod ove funkcije po z i iskoristimo (1.4), pa dobijamo:

∂U

∂z= −

c

R3+

df

dz⇒ f = const. (1.8)

Dakle, na osnovu svega, zakljucujemo da gravitaciona sila, definisana sa (1.1), moze da se predstavi sa

F = gradU, (1.9)

tj. sila je potencijalna, a njen potencijal je odre -den izrazom

U =c

R. (1.10)

Ovaj potencijal u literaturi poznat je i kao Njutnov potencijal. Ovde smo uzeli da je f = const. = 0,sto ne umanjuje opstost prethodno izvedenog izraza.

Kako je:∂U

∂x= −

cx

R3⇒ (1.11)

∂2U

∂x2=

= −c∂x∂x (R3) − x

∂(R3)∂x

R6=

= −c

R6

(

R3 − x∂(R3)

∂R·∂R

∂x

)

=

= −c

R6(R3 − 3x2R) ⇒

∂2U

∂x2=

c

R5

(

3x2 − R2)

. (1.12)

2

Page 3: c:/Documents and Settings/dock/Desktop/Matematicka fizika ...rgf.bg.ac.rs/predmet/GO/IV semestar/Metode matematicke fizike... · Glava 1 Teorija polja 1.1 Primeri nekih polja od interesa

Na slican nacin dobijamo i:

∂2U

∂y2=

c

R5(3y2 − R2) odnosno

∂2U

∂z2=

c

R5(3z2 − R2). (1.13)

Iz ovih poslednjih izraza dobijamo:

∂2U

∂x2+

∂2U

∂y2+

∂2U

∂z2=

c

R5(3x2 + 3y2 + 3z2 − 3R2) = 0. (1.14)

Dakle potencijal U zadovoljava Laplasovu jednacinu:

∆U =∂2U

∂x2+

∂2U

∂y2+

∂2U

∂z2= 0. (1.15)

Funkcije koje zadovoljavaju Laplasovu jednacinu zovemo harmonijske funkcije. Dakle Njutnov potencijalje harmonijska funkcija.

Ravanski zadatak. Logaritamski potencijal

Posmatrajmo sada silu u ravni kojom neka tacka O privlaci materijalnu tacku M . Pretpostavimo daje ova sila privlacenja data izrazom

F = −2c

r

r

r(1.16)

gde je r = |r| intenzitet vektora polozaja r.

Slika 1.2:

Ponovo se postavlja pitanje postojanja potencijala U za ovako definisanu silu. Slicno kao i u prethodnomslucaju polazimo od:

X =∂U

∂x= −

2c

r2x,

Y =∂U

∂y= −

2c

r2y.

(1.17)

Odavde dobijamo

U = −2c

x

r2dx + f(y).

Dalje, kako jer2 = x2 + y2, odakle sledi 2xdx = 2r dr,

dobijamo za potencijal U

U = −2c ln r + f(y).

3

Page 4: c:/Documents and Settings/dock/Desktop/Matematicka fizika ...rgf.bg.ac.rs/predmet/GO/IV semestar/Metode matematicke fizike... · Glava 1 Teorija polja 1.1 Primeri nekih polja od interesa

Dalje, na -dimo sada parcijalni izvod po y:

∂U

∂y=

∂y(−2c ln r) +

df

dy=

= −2c∂

∂r(lnr)

∂r

∂y+

df

dy= −2c

y

r2+

df

dy. (1.18)

Dakle, prema (1.17) (uslov za Y ), zakljucujemo da je dfdy = const. I ovde cemo da uzmemo da je ova

konstanta jednaka 0, pa dobijamo za potencijal

U = −2c ln r. (1.19)

Napomenimo da u slucaju kada vise tacaka privlaci jednu tacku, za potencijal imamo:

U = U(x, y, z) =

n∑

i=1

mi

ri=

=n

i=1

mi√

(x − xi)2 + (y − yi)2 + (z − zi)2.

Njutnova sila kojom homogena sferna ljuska privlaci materijalnu tacku

Posmatrajmo sada homogenu sfernu ljusku, poluprecnika R, koja privlaci materijalnu tacku (sl. 1.3).Prvo odredimo silu kojom beskonacno mali deo sfere (dS) privlaci posmatranu tacku

dF =k2m1dm2

r2

r

r. (1.20)

Slika 1.3:

Kako je: dm2 = dV = · d · dS, za jedinicnu debljinu d = 1 dobijamo dm2 = dS. Kroz posmatranutacku M i centar sfere O uvek mozemo da postavimo jednu pravu. Neka u nasem slucaju to bude z – osa,tj. OM = z. Dalje, oznacimo sa r – rastojanje izme -du tacke M i sredista povrsine dS (tacke u kojoj delujesila privlacenja). Ovo rastojanje je, prema slici (primena kosinusne teoreme)

r2 = R2 + z2 − 2 z R cos θ. (1.21)

4

Page 5: c:/Documents and Settings/dock/Desktop/Matematicka fizika ...rgf.bg.ac.rs/predmet/GO/IV semestar/Metode matematicke fizike... · Glava 1 Teorija polja 1.1 Primeri nekih polja od interesa

Sada za intenzitet sile dobijamo

dF =k2m1dS

R2 + z2 − 2 z R cos θ. (1.22)

Projektujuci ovu silu na z – osu, a vodeci racuna o simetriji, zapazamo da se projekcije na normalan pravac,upravno na z, ponistavaju, pa treba voditi racuna samo o Z – projekciji (videti sl. 1.4)

Slika 1.4:

dZ = −dF · cosβ. (1.23)

Dalje, iz relacije: z = r cosβ + R cos θ (vidi sliku 1.4), dobijamo

cosβ =z − R cosθ

r, (1.24)

tako da konacno dobijamo za projekciju elementarne sile dF na z – osu

dZ = −k2m1dS

R2 + z2 − 2 z R cos θ·z − R cos θ

r=

= −k2m1dS(z − R cos θ)

(R2 + z2 − 2z R cos θ)3/2. (1.25)

Sada ukupnu silu dobijamo integraljenjem gornjeg izraza po celoj sferi S

Z = −

S

k2m1(z − R · cos θ)

(R2 + z2 − 2R z cos θ)3/2dS. (1.26)

Izrazimo sada elementarnu povrsinu dS u sfernim koordinatama

dS = R2 sin θ dϕdθ, (1.27)

pa dobijamo izraz

Z = k2m1R2

∫ π

0

∫ 2π

0

(R cos θ − z) · sin θ dϕdθ

(R2 + z2 − 2R z cos θ)3/2. (1.28)

Dalje, kako je: r2 = R2 + z2 − 2R z cos θ, pri cemu su R– kao poluprecnik, i z– kao rastojanje izme -dunepokretne tacke M i centra sfere, konstantna rastojanja, to dobijamo da je: r dr = R z sin θ dθ, odnosno

R cos θ − z = −r2 + z2 − R2

2z, (1.29)

pa je:∫ π

0

(R cos θ − z) sin θ dθ

(R2 + z2 − 2R z cos θ)3/2= −

1

2Rz2

∫ z+R

|z−R|

(

1 +z2 − R2

r2

)

dr = (1.30)

= −1

2Rz2

[

r −z2 − R2

r

]z+R

|z−R|

.

5

Page 6: c:/Documents and Settings/dock/Desktop/Matematicka fizika ...rgf.bg.ac.rs/predmet/GO/IV semestar/Metode matematicke fizike... · Glava 1 Teorija polja 1.1 Primeri nekih polja od interesa

Analizirajmo sada ovaj rezultat. Kako za θ = 0 ⇒ r = |z−R|, to treba razlikovati dva slucaja. Prvi,kada je tacka van sfere tada je z > R, pa je |z − R| = z − R, odakle sledi da je:

Z = −k2 m1 4π R2ρ

z2= −

k2m1m2

z2. (1.31)

Ovde smo iskoristili da je zapremina sfere, jedinicne debljine: V = 4πR2 · d = 4πR2, gde je d – debljinakoja je jednaka 1, pa je masa m2 = V = 4π R2. I drugi slucaj, kada je tacka u unutrasnosti sfere. Uovom slucaju je z < R, pa je |z − R| = R − z, odakle sledi da je Z = 0.

Privlacenje tacke od strane linijskog tela

Odredimo potencijal homogene prave L, koja se poklapa sa z – osom. Kako posmatrana tacka M idata prava L odre -duju jednu ravan, uzmimo da je to xOz ravan, videti sl. 1.5.

Slika 1.5:

Prvo odredimo silu kojom prava privlaci tacku. Elementarna sila je:

|dF| = kdz

r2= k

dz

x2 + z2. (1.32)

Napomenimo da smo posli od uslova da je sila proporcionalna masama, a obrnuto proporcionalanakvadratu rastojanja. Me -dutim, kako je dm = dV = P dz i P = 1, to smo dobili da sila zavisi od dz,dok smo sve ostale konstante oznacili jednim slovom – k. Dalje, slicno kao i u prethodnom primeru, zbogsimetricnosti projekcija, videti sl. 1.5, ponovo je potrebno odrediti samo X – projekciju

dX = dF cosα = −kdz

r2·x

r. (1.33)

Ukupna projekcija je

X = −

∫ +∞

−∞

kx

r·dz

r2. (1.34)

Dalje, kako je

tgα =z

x⇒ z = x tgα ⇒ dz =

x

cos2 αdα, x = const, (1.35)

i

r3 = (x2 + z2)3/2 = x3(1 + tg2α)3/2 =x3

cos3 α, (1.36)

6

Page 7: c:/Documents and Settings/dock/Desktop/Matematicka fizika ...rgf.bg.ac.rs/predmet/GO/IV semestar/Metode matematicke fizike... · Glava 1 Teorija polja 1.1 Primeri nekih polja od interesa

to za X dobijamo

X = −k

∫ −π

2

+ π

2

x·cos3 α

x3

x

cos2 αdα. (1.37)

Kako je x = const, to dobijamo

X = −k1

x

∫ −π

2

+ π

2

cosα dα = −2k

x. (1.38)

Sada mozemo da potrazimo potencijal. Ako postoji potencijal U , on mora da zadovoljava relaciju

X =∂U

∂x= −

2k

x, (1.39)

odakle konacno dobijamo:

U = −2k lnx ili U = 2k ln1

x. (1.40)

Teorema 1 Proizvoljno vektorsko polje F, jednoznacno, neprekidno i ograniceno, moze da se razlozi na

zbir potencijalnog i bezvrtloznog vektorskog polja u obliku:

F = −∇U + ∇× A,

pri cemu je ∇ · A = 0.

Skalarnu funkciju U – nazivamo skalarni potencijal, a vektorsku A – vektorski potencijal vektorskogpolja F.

Ova teorema poznata je u literaturi kao teorema Helmholca. 1 Teoremu smo naveli zbog njene vaznosti.Me -dutim, kako je dokaz relativno slozen, necemo da ga izlozimo, vec citaoca upucujemo na reference: [?](str. 79), [?] (str. 50).

1.2 Generalisane koordinate

Polozaj tacke u trodimenzionom Euklidskom prostoru odre -duje se, u odnosu na jednu unapred odre -denutacku O, koju zovemo pol ili koordinatni pocetak, pomocu vektora polozaja r. U odnosu na Dekartovpravougli sistem koordinata Oxyz, sa pocetkom u polu O, polozaj tacke se odre -duje koordinatama tacke

(x, y, z). Ortogonalne projekcije kraja vektora polozaja na ose ovog koordinatnog sistema poklapaju se sakoordinatama tacke, pa se koordinate vektora polozaja r poklapaju sa ovim koordinatama x, y, z:

r = {x, y, z} ili r = xi + yj + zk. (1.41)

Me -dutim, polozaj tacke u prostoru moze da se odredi i pomocu neka tri me -dusobno nezavisna parametraq1, q2, q3 (ovde 1,2 i 3 nisu eksponenti vec oznake parametara !!!), ili krace qi (i=1,2,3). Kada se parametru

qi daju sve moguce vrednosti, a svakoj tacki prostora odgovara jedan i samo jedan ure -den skup od tri broja(q1, q2, q3) i obrnuto, svakom skupu od tri broja (q1, q2, q3) odgovara jedna i samo jedna tacka u prostoru(trodimenzionom), tada se parametri qi nazivaju opste ili generalisane koordinate tacke. Sada vektorpolozaja mozemo da izrazimo i preko generalisanih koordinata:

r = r(q1, q2, q3), ili krace r = r(qi). (1.42)

Dekartove koordinate ovog vektora mozemo da izrazimo u obliku:

x = x(qi),

y = y(qi), (1.43)

z = z(qi).

1Hermann von Helmholtc (1821–1894), nemacki fizicar. Poznat je po veoma vaznim radovima iz oblasti termodinamike,hidrodinamike i akustike.

7

Page 8: c:/Documents and Settings/dock/Desktop/Matematicka fizika ...rgf.bg.ac.rs/predmet/GO/IV semestar/Metode matematicke fizike... · Glava 1 Teorija polja 1.1 Primeri nekih polja od interesa

Iz zahteva za obostranom jednoznacnoscu izme -du tacaka u prostoru i koordinata qi sledi da svakoj tackisa koordinatama (x, y, z) moraju da odgovaraju tri broja qi, tako da je:

qi = qi(x, y, z),∂(q1, q2, q3)

∂(x, y, z)6= 0. (1.44)

Dakle, jednacine (1.43) uvek zadovoljavaju uslove potrebne da mogu da se rese po qi.Skup jednacina (1.43) i (1.44) predstavljaju jednacine koordinatne transformacije. Ove transformacije

su uzajamno reciprocne – inverzne. Definisimo sada kordinatne linije i koordinatne povrsi.

Definicija.

Kordinatne linije –predstavljaju geometrijsko mesto tacaka koje se dobija akosu dve koordinate konstantne, a treca se menja.

Koordinatne linije mogu da budu prave ili krive i u zavisnosti od toga govorimo o pravolinijskim ilikrivolinijskim koordinatnim sistemima (slika 1.6).

Slika 1.6: Koordinatne linije

Definicija.

Kordinatne povrsi – predstavljaju geometrijsko mesto tacaka koje se dobijaako se menjaju dve koordinate, a treca je konstantna.

Slika 1.7: Koordinatne povrsi

Svakoj tacki prostora odgovaraju tri koordinatne povrsi, tj. tacka predstavlja presek tri koordinatnepovrsi. Koordinatna linija se nalazi u preseku dve koordinatne povrsi (slika 1.7). Kako kroz svaku tackuprolaze po tri ovakve, u opstem slucaju krive linije, to one obrazuju krivolinijski sistem koordinata.

U svakoj tacki prostora mogu da se na koordinatne linije, u toj tacki, povuku tangente, orijentisaneprema smeru u kome rastu vrednosti qi i one cine, u toj tacki, trijedar generalisanih ili krivolinijskihkoordinata.

8

Page 9: c:/Documents and Settings/dock/Desktop/Matematicka fizika ...rgf.bg.ac.rs/predmet/GO/IV semestar/Metode matematicke fizike... · Glava 1 Teorija polja 1.1 Primeri nekih polja od interesa

Bazni vektori ovog sistema koordinata su njihovi tangentni vektori.Koordinatni sistemi, u odnosu na koje se odre -duje polozaj tacaka u prostoru, zovu se i sistemi referencije.

Prema obliku koordinatne linije imamo:

- pravolinijske (slika 1.8) i

- krivolinijske koordinate (slika 1.6).

Ako su koordinatne povrsi, odnosno koordinatne linije me -dusobno upravne u svim tackama prostora, imamoortogonalan (pravougli) sistem koordinata. Specijalno u slucaju Dekartovog koordinatnog sistemapolozaj tacke prostora odre -den je vektorom polozaja

r = xi + yj + zk. (1.45)

Ocigledno je da je:

i =∂r

∂x, j =

∂r

∂y, k =

∂r

∂z(1.46)

jedinicni vektori (ortovi) mogu da dobiju kao parcijalni izvodi vektora polozaja po koordinati tacke.U slucaju generalisanih koordinata polozaj iste tacke je odre -den istim vektorom polozaja r, ali izrazenog

u odnosu na krivolinijske koordinate qi (vidi 1.42). Onda su bazni vektori tog sistema koordinata odre -deni(kao i u slucaju Dekartovih koordinata (1.46)) sa

gi =∂r

∂qi, (i = 1, 2, 3). (1.47)

Dakle, vektori gi su osnovni (bazni) vektori sistema generalisanih koordinata qi i u pravcu su tangentena koordinatne linije. Da su oni bazni vektori zakljucuje se na osnovu njihove linearne nezavisnosti (vidizadatak ???, str. ???)

Uocimo da su indeksi kojima se obelezavaju bazni vektori gi i ei na donjem mestu.Ako koordinate osnovnih vektora izrazimo u odnosu na Dekartove pravougle koordinate, intenziteti ovih

vektora su odre -deni izrazima:

|gi| =

(

∂x

∂qi

)2

+

(

∂y

∂qi

)2

+

(

∂z

∂qi

)2

≡ hi, (i = 1, 2, 3), (1.48)

velicine hi zovemo Lameovi 2 ili metricki koeficijenti.Ako sa ei oznacimo jedinicni vektor u pravcu ovih vektora, imamo:

ei =gi

|gi|⇒ gi = hiei. (1.49)

Slika 1.8: Kosougle koordinatne linije

2Gabriel Lame (1795-1870), francuski matematicar i inzenjer. Dao je veliki doprinos u analitickoj geometriji i analitickojmehanici, kao i u teoriji elasticnosti i teoriji provo -denja toplote.

9

Page 10: c:/Documents and Settings/dock/Desktop/Matematicka fizika ...rgf.bg.ac.rs/predmet/GO/IV semestar/Metode matematicke fizike... · Glava 1 Teorija polja 1.1 Primeri nekih polja od interesa

Pored ovog skupa jedninicnih vektora mozemo da posmatramo i skup vektora gradqi, koji su u pravcunormale na koordinatne povrsi, a definisani su sa

Ei =∇qi

|∇qi|. (1.50)

Vektor Ei je dakle jedinicni vektor normale koordinatne povrsi. Dakle, u svakoj tacki prostora imamodva skupa jedinicnih vektora: ei– jedinicni vektori u pravcu tangente i Ei– jedinicni vektor normale nakoordinatnu povrs (sl. 1.9). Ova dva skupa se poklapaju akko je sistem koordinata ortogonalan. Moze dase pokaze da ovi vektori cine reciprocan sistem vektora.

Slika 1.9:

Posmatrajmo proizvoljan vektor v i predstavimo ga preko ova dva sistema:

v =

= v1e1 + v2e2 + v3e3 =

3∑

i=1

viei (1.51)

= V1E1 + V2E

2 + V3E3 =

3∑

i=1

ViEi.

Tako -de, vektor mozemo da predstavimo i preko osnovnih vektora gi =∂r

∂qi:

v =

=3

i=1

cigi =3

i=1

ci ∂r

∂qi= c1 ∂r

∂q1+ c2 ∂r

∂q2+ c3 ∂r

∂q3(1.52)

=3

i=1

Ci∇qi = C1∇q1 + C2∇q2 + C3∇q3.

Koordinate dobijene razlaganjem vektora na ove nacine zovemo:ci – kontravarijantne iCi – kovarijantne, respektivno.

1.2.1 Element luka i zapremine

Duzinu elementa luka (sl. 1.10) odre -dujemo relacijom:

ds2 = dr · dr. (1.53)

10

Page 11: c:/Documents and Settings/dock/Desktop/Matematicka fizika ...rgf.bg.ac.rs/predmet/GO/IV semestar/Metode matematicke fizike... · Glava 1 Teorija polja 1.1 Primeri nekih polja od interesa

Slika 1.10: Element luka

Kako diferencijal ovog vektora mozemo da izrazimo na jedan od sledecih nacina (u zavisnosti da li suDekartove ili neke druge - generalisane koordinate):

dr =

= dxi + dyj + dzk =

=3

i=1

∂r

∂qidqi =

3∑

n=1

gidqi = (1.54)

=3

n=1

hieidqi.

Ako Dekartove koordinate umesto sa x, y i z oznacimo sa x1, x2 i x3, respektivno, to za element luka(preko Dekartovih koordinata) dobijamo:

ds2 =

3∑

i,j=1

δijdxidxj =

3∑

i=1

dxidxi (1.55)

gde je δij – Kronekerov delta simbol. Dalje, prema (1.43), izracunavamo dxi:

xi = xi(qj) ⇒ dxi =

3∑

j=1

∂xi

∂qjdqj , (1.56)

pa dobijamo za element luka:

ds2 =3

i,j

δijdxidxj =3

i

dxidxi =

=

3∑

i,j=1

3∑

k,l=1

δij∂xi

∂qk·∂xj

∂qldqkdql =

=3

k,l=1

3∑

i,j=1

δij∂xi

∂qk·∂xj

∂ql

dqkdql = (1.57)

=

3∑

k,l=1

gkldqkdql.

11

Page 12: c:/Documents and Settings/dock/Desktop/Matematicka fizika ...rgf.bg.ac.rs/predmet/GO/IV semestar/Metode matematicke fizike... · Glava 1 Teorija polja 1.1 Primeri nekih polja od interesa

Definicija.

Velicine gkl, definisane relacijama:

gkl = gk · gl =

3∑

i,j=1

δij∂xi

∂qk·∂xj

∂ql(1.58)

zovemo osnovni metricki tenzor ili osnovni tenzor prostora.

Napomenimo da je gij = h2i δij za ortogonalni sistem koordinata, naime u tom slucaju imamo

ei · ej = δij , (1.59)

ili u razvijenom oblikue1 · e2 = e2 · e3 = e3 · e1 = 0,

e1 · e1 = e2 · e2 = e3 · e3 = 1,(1.60)

tako da za ortogonalni sistem koordinata, za element luka, dobijamo:

ds2 = dr · dr =

=

3∑

i,j=1

hidqihjdqjei · ej =

3∑

i,j=1

δijhidqihjdqj = (1.61)

=

3∑

i=1

(

hidqi)2

.

Na kraju napisimo izraz za elementarnu zapreminu (sl. 1.11), preko generalisanih ortogonalnih koordi-nata:

dV =∣

(

h1dq1e1

)

·(

h2dq2e2

)

×(

h3dq3e3

)∣

∣ = (1.62)

= h1h2h3dq1dq2dq3 (e1 · e2 × e3) = h1h2h3dq1dq2dq3.

Slika 1.11: Element zapremine

12

Page 13: c:/Documents and Settings/dock/Desktop/Matematicka fizika ...rgf.bg.ac.rs/predmet/GO/IV semestar/Metode matematicke fizike... · Glava 1 Teorija polja 1.1 Primeri nekih polja od interesa

1.2.2 Gradijent, divergencija, rotor i Laplasijan - izrazeni preko generalisanih

koordinata

Gradijent proizvoljne funkcije U , izrazen u generalisanim koordinatama je:

gradU =

∇U =1

h1·

∂U

∂q1e1 +

1

h2·

∂U

∂q2e2 +

1

h3·

∂U

∂q3e3 = (1.63)

=3

i=1

1

hi·∂U

∂qiei.

Divergencija je data izrazom

divA = ∇ · A = (1.64)

=1

h1h2h3

[

∂q1(h2h3A1) +

∂q2(h3h1A2) +

∂q3(h1h2A3)

]

.

Rotor je dat izrazom

rotA = ∇× A =1

h1h2h3·

h1e1 h2e2 h3e3

∂q1

∂q2

∂q3

h1A1 h2A2 h3A3

. (1.65)

Na kraju izrazimo Laplasijan (∆ = ∇2) u odnosu na generalisane koordinate

∆U = (1.66)

=1

h1h2h3

[

∂q1

(

h2h3

h1·

∂U

∂q1

)

+∂

∂q2

(

h3h1

h2·

∂U

∂q2

)

+∂

∂q3

(

h1h2

h3·

∂U

∂q3

)]

.

1.3 Posebni koordinatni sistemi

1. CILINDRICNI (, ϕ, z)

Veza izme -du Dekartovih i cilindricnih koordinata data je sa (sl. 1.12):

13

Page 14: c:/Documents and Settings/dock/Desktop/Matematicka fizika ...rgf.bg.ac.rs/predmet/GO/IV semestar/Metode matematicke fizike... · Glava 1 Teorija polja 1.1 Primeri nekih polja od interesa

Slika 1.12:

x = cosϕ, y = sin ϕ, z = z,

gde je ≥ 0, 0 ≤ ϕ < 2π, −∞ < z < +∞.

Pokazimo sada na primeru ovih koordinata kako se odre -duju Lameovi koeficijenti. Izracunajmo prvodiferencijale dx, dy, dz preko novih koordinata:

dx = cosϕd + (− sinϕ) dϕ,

dy = sin ϕd + cosϕdϕ,

dz = dz.

Zatim izracunamo element luka:

ds2 = dx2 + dy2 + dz2 =

= d2 + 2dϕ2 + dz2, (1.67)

pa poredeci sa (1.57) zakljucujemo da je:

h = 1 hϕ = hz = 1.

2. SFERNI (r, θ, ϕ)

U ovom slucaju veza izme -du ovih koordinata je, prema sl. 1.13:

14

Page 15: c:/Documents and Settings/dock/Desktop/Matematicka fizika ...rgf.bg.ac.rs/predmet/GO/IV semestar/Metode matematicke fizike... · Glava 1 Teorija polja 1.1 Primeri nekih polja od interesa

Slika 1.13:

x = r sin θ cosϕ, y = r sin θ sin ϕ, z = r cos θ,

gde je r ≥ 0, 0 ≤ ϕ < 2π, 0 ≤ θ ≤ π.

Na slican nacin kao i u prethodnom slucaju dobijamo:

hr = 1 hϕ = r sin θ hθ = r.

3. PARABOLICNO-CILINDRICNE (u, v, z)

x =1

2

(

u2 − v2)

, y = uv, z = z

gde je −∞ < u < ∞, v ≥ 0, −∞ < z < ∞

hu = hv =√

u2 + v2, hz = 1,

dok je veza sa cilindricnim kordinatama data relacijama:

u =√

2 cosϕ

2, v =

2 sinϕ

2, z = z.

4. PARABOLOIDNE (u, v, ϕ)

x = u v cosϕ, y = u v sinϕ, z =1

2

(

u2 − v2)

gde je u ≥ 0, v ≥ 0, 0 ≤ ϕ < 2π

i hu = hv =√

u2 + v2, hϕ = u v.

5. ELIPTICKO-CILINDRICNE (u, v, z)

15

Page 16: c:/Documents and Settings/dock/Desktop/Matematicka fizika ...rgf.bg.ac.rs/predmet/GO/IV semestar/Metode matematicke fizike... · Glava 1 Teorija polja 1.1 Primeri nekih polja od interesa

x = a ch u cos v, y = a sh u sin v, z = z

gde je u ≥ 0, 0 ≤ v < 2π, −∞ < z < ∞

i hu = hv = a√

sh2u + sin2 v, hz = 1.

6. SFEROIDNE (ξ, η, ϕ)

a)

x = a sh ξ sin η cosϕ, y = a sh ξ sin η sin ϕ,

z = a ch ξ cos η

gde je ξ ≥ 0, 0 ≤ η ≤ π, 0 ≤ ϕ < 2π

i hξ = hη = a

sh2ξ + sin2 η, hϕ = a shξ sin η.

b)

x = a ch ξ cos η cosϕ, y = a ch ξ cos η sin ϕ,

z = a sh ξ sin η

gde je ξ ≥ 0, −π

2≤ η ≤

π

2, 0 ≤ ϕ < 2π

i hξ = hη = a

sh2ξ + sin2 η, hϕ = a chξ cos η.

7. ELIPSOIDALNE (α, β, ϕ)

x = c sin β cosϕ, y = c sin α sinβ sin ϕ, z = c ch α cosβ

gde su:

0 ≤ α < ∞, 0 ≤ β ≤ π, −π < ϕ ≤ π

a Lameovi koeficijenti su:

hα = hβ = c

sh2 + sin2 β, hϕ = c shα sin β.

8. ELIPSOIDNE (λ, µ, ν)

x2

a2 − λ+

y2

b2 − λ+

z2

c2 − λ= 1 λ < c2 < b2 < a2

x2

a2 − µ+

y2

b2 − µ+

z2

c2 − µ= 1 c2 < µ < b2 < a2

x2

a2 − ν+

y2

b2 − ν+

z2

c2 − ν= 1 c2 < b2 < ν < a2,

16

Page 17: c:/Documents and Settings/dock/Desktop/Matematicka fizika ...rgf.bg.ac.rs/predmet/GO/IV semestar/Metode matematicke fizike... · Glava 1 Teorija polja 1.1 Primeri nekih polja od interesa

pri cemu su:

hλ =1

2

(µ − λ) (ν − λ)

(a2 − λ) (b2 − λ) (c2 − λ)

hµ =1

2

(ν − µ) (λ − µ)

(a2 − µ) (b2 − µ) (c2 − µ)

hν =1

2

(λ − ν) (µ − ν)

(a2 − ν) (b2 − ν) (c2 − ν).

9. BIPOLARNE – ravanske (u, v, z)

x2 + (y − a cotu)2

= a2 csc2u,

(x − a coth v)2

+ y2 = a2 csch2v, z = z

ili x =a sh v

ch v − cos u, y =

a sin u

ch v − cos u, z = z

gde je 0 ≤ u < 2π, −∞ < v < ∞, −∞ < z < ∞.

Za Lameove koeficijente dobijamo:

hu = hv =a

ch v − cos u, hz = 1.

10. BIPOLARNE – bisfericne (u, v, ϕ)

x =c sin u cosϕ

chv − cosu, y =

c sin u sin ϕ

chv − cosu, z =

c shv

chv − cosu

hu = hv =c

chv − cosu, hϕ =

c sin u

chv − cosu.

17

Page 18: c:/Documents and Settings/dock/Desktop/Matematicka fizika ...rgf.bg.ac.rs/predmet/GO/IV semestar/Metode matematicke fizike... · Glava 1 Teorija polja 1.1 Primeri nekih polja od interesa

11. ELIPTICKE (λ, µ, z)

x = c λµ, y = c√

(λ2 − 1) (1 − µ2), z = z

Lameovi (metricki) koeficijenti su:

hλ = c

λ2 − µ2

λ2 − 1, hµ = c

λ2 − µ2

1 − µ2, hz = 1.

12. TOROIDALNE (α, β, ϕ)

x =c shα cosϕ

chα − cosβ, y =

c shα sinϕ

chα − cosβ, z =

c sin β

chα − cosβ

gde je: 0 ≤ α < ∞, −π < β ≤ π, −π < ϕ ≤ π

dok su Lameovi koeficijenti:

hα = hβ =c

chα − cosβ, hϕ =

c shα

chα − cosβ.

14. SFEROIDALNE – a (λ, µ, ϕ)

a)

x = cλµ, y = c√

(λ2 − 1) (1 − µ2) cosϕ, z = c√

(λ2 − 1) (1 − µ2) sinϕ,

gde je: λ ≤ 1, −1 ≤ µ ≤ 1, 0 ≤ ϕ ≤ 2π

gde su Lameovi koeficijenti:

hλ = c

λ2 − µ2

λ2 − 1, hµ = c

λ2 − µ2

1 − µ2, hϕ = c

(λ2 − 1) (1 − µ2).

b)

x = cλµ sin ϕ, y = c√

(λ2 − 1) (1 − µ2), z = cλµ cosϕ,

gde su Lameovi koeficijenti:

hλ = c

λ2 − µ2

λ2 − 1, hµ = c

λ2 − µ2

1 − µ2, hϕ = cλµ.

18