20.01.2013 Views

kart_kreiv_ir_pav_integralai.pdf

kart_kreiv_ir_pav_integralai.pdf

kart_kreiv_ir_pav_integralai.pdf

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

V skyrius<br />

KARTOTINIAI, KREIVINIAI IR PAVIRŠINIAI<br />

INTEGRALAI<br />

§ 1. DVILYPIS INTEGRALAS IR CILINDROIDO TŪRIS ....................................................83<br />

§ 2. DVILYPIO INTEGRALO APSKAIČIAVIMAS STAČIAKAMPĖJE KOORDINAČIŲ<br />

SISTEMOJE ...............................................................................................................................84<br />

§ 3. DVILYPIS INTEGRALAS POLINĖJE KOORDINAČIŲ SISTEMOJE IR JO<br />

APSKAIČIAVIMAS..................................................................................................................86<br />

§ 4. TRILYPIS INTEGRALAS, JO SAVYBĖS IR APSKAIČIAVIMAS STAČIAKAMPĖJE<br />

KOORDINAČIŲ SISTEMOJE..................................................................................................88<br />

§ 5. TRILYPIO INTEGRALO APSKAIČIAVIMAS CILINDRINĖJE IR SFERINĖJE<br />

KOORDINAČIŲ SISTEMOSE .................................................................................................90<br />

§ 6. KARTOTINIŲ INTEGRALŲ TAIKYMAS MECHANIKOJE .........................................91<br />

§ 7. SKALIARINIAI IR VEKTORINIAI LAUKAI..................................................................94<br />

§ 8. LYGIO LINIJOS IR PAVIRŠIAI .......................................................................................94<br />

§ 9. GRADIENTAS IR JO SAVYBĖS......................................................................................95<br />

§ 10. PIRMOJO IR ANTROJO TIPŲ KREIVINIAI INTEGRALAI........................................97<br />

§ 11. KREIVINIŲ INTEGRALŲ APSKAIČIAVIMAS ...........................................................98<br />

§ 12. PLOTŲ APSKAIČIAVIMAS KREIVINIU INTEGRALU .............................................99<br />

§ 12. GRINO FORMULĖ ........................................................................................................100<br />

§ 13. SĄLYGA, KADA KREIVINIS INTEGRALAS NEPRIKLAUSO NUO<br />

INTEGRAVIMO KELIO.........................................................................................................101<br />

§ 14. VIENPUSIAI IR DVIPUSIAI PAVIRŠIAI. PAVIRŠIAUS ORIENTACIJA ...............103<br />

§ 15. PAVIRŠINIAI INTEGRALAI........................................................................................104<br />

§ 16. PAVIRŠINIO INTEGRALO APSKAIČIAVIMAS .......................................................105<br />

§ 17. PAVIRŠINIO IR KREIVINIO INTEGRALŲ RYŠYS (STOKSO FORMULĖ)...........106<br />

§ 18. GAUSO – OSTROGRADSKIO FORMULĖ .................................................................109<br />

© A. Laurutis, D.Šiaučiūnas Paskaitų konspektas 82


§ 1. DVILYPIS INTEGRALAS IR CILINDROIDO TŪRIS<br />

Imkime xOy plokštumos uždarą sritį D , apribotą linija L . Šiai sričiai priklauso ne<br />

tik vidiniai srities D taškai, bet <strong>ir</strong> linijosLtaškai.Tegul šioje srityje D duota tolydi funkcija<br />

z = f ( x,<br />

y)<br />

.<br />

y<br />

Sritį D bet kokiomis linijomis išdalinkime į n dalių<br />

D<br />

∆ s1 , ∆s2,<br />

∆s3,...,<br />

∆sn<br />

, kurias vadinsime sritelėmis.<br />

Nenorėdami įvesti naujų simbolių, sritelių<br />

Pi<br />

∆s<br />

i<br />

∆ s1 , ∆s2,<br />

∆s3,...,<br />

∆sn<br />

∆ s1 , ∆s2,<br />

∆s3,...,<br />

∆sn<br />

.<br />

plotą žymėkime taip pat<br />

Kiekvienos ∆ si<br />

sritelės viduje imkime taškus P i .<br />

L Turime n<br />

apskaičiuojame<br />

taškų P 1 , P2<br />

, P3<br />

,..., Pn<br />

,<br />

funkcijos z = f ( x,<br />

y)<br />

kuriuose<br />

reikšmes<br />

f P f P , f P ,..., f P .<br />

Šią sandaugų ( i ) i<br />

f ( x,<br />

y)<br />

integraline suma srityje D .<br />

Jei f ( x,<br />

y)<br />

≥ 0 , tai kiekvieną iš sandaugų ( Pi<br />

) si<br />

( 1 ) , ( 2 ) ( 3 ) ( n )<br />

Sudarykime n sandaugų ( Pi<br />

) si<br />

f P ∆ s sumą vadinsime funkcijos<br />

f ∆ galime<br />

įsivaizduoti kaip mažo cilindro, turinčio pagrindą ∆ si<br />

<strong>ir</strong><br />

aukštinę f ( Pi<br />

) , tūrį. Tada integralinė suma V n reikš tam<br />

tikro laiptuoto cilindrinio kūno tūrį.<br />

Sritį D galėjome kitaip skaidyti į sriteles, galėjome<br />

kitaip pas<strong>ir</strong>inkti taškus P i jose. Todėl galėjome gauti<br />

n<br />

∑<br />

i=<br />

1<br />

daugybę integralinių sumų, V = f ( P )<br />

n<br />

i<br />

V<br />

∆s<br />

i<br />

n<br />

f ∆ <strong>ir</strong> šių sandaugų sumą<br />

( P ) ∆s<br />

+ f ( P ) ∆s<br />

+ f ( P ) ∆s<br />

= f ( P )<br />

© A. Laurutis, D.Šiaučiūnas Paskaitų konspektas 83<br />

=<br />

f<br />

1<br />

1<br />

kurios visos<br />

sudaro integralinių sumų seką V n , V , ,..., ,...<br />

1 n V 2 n V<br />

3 nk<br />

( k ∈ 1,<br />

m ). Apskaičiuokime tokios sekos ribą, kai<br />

maksimalios sritelės plotas ∆ si<br />

artėja į nulį<br />

( max s → 0)<br />

arba n → ∞ .<br />

∆ i<br />

Teorema. Jei funkcija ( x y)<br />

V<br />

n<br />

=<br />

n<br />

∑<br />

i=<br />

1<br />

f<br />

( P )<br />

i<br />

∆s<br />

i<br />

2<br />

2<br />

n<br />

n<br />

n<br />

∑<br />

i=<br />

1<br />

i<br />

∆s<br />

,<br />

f , yra tolydi tam tikroje uždaroje srityje D , tai integralinių sumų<br />

seka turi ribą, kai maksimalios sritelės plotas artėja į nulį ( → ∞)<br />

n . Ši riba<br />

nepriklauso nei nuo srities D skaidymo į aikšteles būdo, nei taško P i pas<strong>ir</strong>inkimo aikštelėje ∆ Si<br />

būdo. Ši riba vadinama funkcijos ( x y)<br />

n<br />

∑<br />

lim<br />

∆Si<br />

→0<br />

i=<br />

1<br />

Sritis D vadinama integravimo sritimi.<br />

f<br />

x<br />

f , dvilypiu integralu srityje D <strong>ir</strong> žymima<br />

( P ) ∆S<br />

= f ( x,<br />

y)<br />

ds = f ( x,<br />

y)<br />

i<br />

i<br />

∫∫<br />

D<br />

∫∫<br />

D<br />

x<br />

z<br />

dxdy .<br />

f(<br />

iP )<br />

iP<br />

D<br />

S<br />

∆s<br />

i<br />

i<br />

( z y)<br />

z = f ,<br />

y


Jei f ( x,<br />

y)<br />

≥ 0 , tai dvilypis integralas srityje D reiškia cilindroido, sudaryto tarp plokštumos<br />

z = 0 , <strong>pav</strong><strong>ir</strong>šiaus z = f ( x,<br />

y)<br />

bei cilindrinio <strong>pav</strong><strong>ir</strong>šiaus, kurio vedančioji linija yra L , tūriui.<br />

Pastebėkime, kad srities D plotas yra lygus srities D plotui, todėl<br />

∫∫<br />

D<br />

( x y)<br />

ds = f ( x y)<br />

© A. Laurutis, D.Šiaučiūnas Paskaitų konspektas 84<br />

∫∫<br />

f , , dxdy .<br />

Dvilypio integralo savybes nusakykime šiomis teoremomis.<br />

I teorema. Dviejų funkcijų algebrinės sumos f ( x,<br />

y)<br />

ϕ(<br />

x,<br />

y)<br />

+ dvilypis integralas srityje D yra<br />

lygus sumai dviejų dvilypių integralų, apskaičiuotų kiekvienai funkcijai atsk<strong>ir</strong>ai, t.y.<br />

∫∫<br />

D<br />

[ f ( x,<br />

y)<br />

+ ( x,<br />

y)<br />

] ds = f ( x,<br />

y)<br />

ds + ϕ ( x,<br />

y)<br />

D<br />

∫∫<br />

D<br />

∫∫<br />

ϕ ds .<br />

II teorema. Pastovų daugiklį galime iškelti prieš integralo ženklą, t.y., jei a = const , tai<br />

af x,<br />

y ds = a f x,<br />

y ds .<br />

∫∫<br />

D<br />

( ) ( )<br />

Šias teoremas galime įrodyti analogiškai, kaip atitinkamas apibrėžtinio integralo teoremas.<br />

III teorema. Jeigu sritis D susideda iš dviejų sričių 1 D <strong>ir</strong> D 2 , neturinčių bendrų vidinių taškų <strong>ir</strong><br />

x y<br />

f x,<br />

y ds = f x,<br />

y ds + f x,<br />

y ds .<br />

f ( , ) yra tolydi visoje D tai ( ) ( ) ( )<br />

> Srities D integralinę sumą galima užrašyti taip<br />

∑<br />

D<br />

∫∫<br />

D<br />

f<br />

∫∫<br />

D1<br />

∫∫<br />

D<br />

∫∫<br />

D2<br />

( P ) ∆s<br />

= f ( P ) ∆s<br />

+ f ( P )<br />

i<br />

Apskaičiuokime šios lygybės ribą, kai max s → 0 , tada gausime<br />

∫∫<br />

D<br />

i<br />

∑<br />

D1<br />

i<br />

∆ i<br />

i<br />

∑<br />

D2<br />

( x,<br />

y)<br />

ds = f ( x,<br />

y)<br />

ds + f ( x,<br />

y)<br />

∫∫<br />

D1<br />

∫∫<br />

D2<br />

D<br />

i<br />

∆s<br />

f ds <<br />

Ši teorema galioja bet kuriam srities D dalių skaičiui.<br />

§ 2. DVILYPIO INTEGRALO APSKAIČIAVIMAS STAČIAKAMPĖJE KOORDINAČIŲ<br />

SISTEMOJE<br />

y<br />

O<br />

A<br />

M 1<br />

y = ϕ2<br />

y = ϕ1<br />

D<br />

( x)<br />

( x)<br />

N1<br />

N 2<br />

M 2<br />

a b<br />

L<br />

B<br />

x<br />

i<br />

Tegul sritis D , esanti xOy plokštumoje<br />

yra tokia, kad kiekviena tiesė, lygiagreti<br />

vienai koordinačių ašiai (pvz.: Oy ) <strong>ir</strong><br />

einanti per bet kurį vidinį srities D tašką<br />

kerta šios srities ribą dviejuose taškuose 1 N<br />

<strong>ir</strong> 2<br />

N . Tokią sritį vadiname teisinga Oy<br />

ašies kryptimi. Analogiškai apibrėžiama <strong>ir</strong><br />

teisinga sritis ašies Ox kryptimi. Sritį,<br />

teisingą Ox ašies <strong>ir</strong> Oy ašies kryptimi<br />

vadiname teisinga sritimi. Brėžinyje<br />

<strong>pav</strong>aizduota teisinga sritis D .


Matome, kad nagrinėjamoji sritis D , ( x a,<br />

b)<br />

atitinkamai, linijomis y = ϕ1(<br />

x)<br />

<strong>ir</strong> y ϕ 2(<br />

x)<br />

ϕ1( x) ≤ ϕ2(<br />

x)<br />

<strong>ir</strong> a < b , funkcijos ϕ 1(<br />

x) <strong>ir</strong> ϕ2<br />

( x)<br />

yra tolydžios atkarpoje [<br />

D yra žinoma tolydi funkcija z = f ( x,<br />

y)<br />

. Tada reiškinys<br />

b]<br />

yra funkcijos ( x y)<br />

I =<br />

∈ iš apačios <strong>ir</strong> iš v<strong>ir</strong>šaus yra ribojamos,<br />

= , kurios abi sudaro liniją L . Čia, suprantama,<br />

a, . Tegul, be to srityje<br />

2<br />

( x,<br />

y)<br />

dxdy = ⎜ f ( x,<br />

y)<br />

© A. Laurutis, D.Šiaučiūnas Paskaitų konspektas 85<br />

⎛<br />

⎜<br />

⎝<br />

ϕ ( x)<br />

∫∫ ∫ ∫ ⎟ ⎟<br />

f<br />

dy<br />

D<br />

f , dvilypis integralas srityje D . Jis apskaičiuojamas dvigubu integralu:<br />

pradžioje integralas, esantis skliaustuose, integruojant jį pagal kintamąjį y (kai x = const ). Šį<br />

suintegravus, gaunama tolydi kintamojo x funkcija<br />

Φ<br />

b<br />

a<br />

ϕ ( x)<br />

ϕ ( x)<br />

2<br />

( x)<br />

= f ( x,<br />

y)<br />

∫<br />

ϕ ( x )<br />

1<br />

1<br />

⎞<br />

dx<br />

⎠<br />

(1)<br />

dy . (2)<br />

Funkciją Φ ( x)<br />

toliau integruojame išoriniu integralu pagal kintamąjį x atkarpoje [ a, b]<br />

b<br />

∫ Φ = I dx . (3)<br />

a<br />

Gautas rezultatas (3) yra pastovus skaičius <strong>ir</strong> nusako dvilypio integralo reikšmę<br />

stačiakampėje koordinačių sistemoje.<br />

ϕ x = ϕ x )<br />

Gali atsitikti, kad sritis D yra tokia, kad vienos iš funkcijų y = 1(<br />

) (arba y 2(<br />

)<br />

atkarpoje [ a, b]<br />

negalime aprašyti viena analitine išraiška. Tegul, <strong>pav</strong>yzdžiui<br />

⎧<br />

D1<br />

: ⎨<br />

⎩<br />

⎧ψ<br />

( x)<br />

, kai<br />

1(<br />

x)<br />

= ⎨<br />

⎩χ<br />

( x)<br />

, kai<br />

x ∈ ( a,<br />

c)<br />

,<br />

y ∈ ( ψ ( x)<br />

, ϕ ( x)<br />

)<br />

( x)<br />

a ≤ x ≤ c<br />

ϕ , a < c < b , tada sritį D tenka išdalinti į dvi sritis<br />

c ≤ x ≤ b<br />

2<br />

<strong>ir</strong><br />

D<br />

2<br />

⎧x<br />

∈<br />

: ⎨<br />

⎩y<br />

∈<br />

( c,<br />

b)<br />

,<br />

( χ(<br />

x)<br />

, ϕ ( x)<br />

)<br />

2<br />

, todėl dvilypis integralas<br />

apskaičiuojamas pasinaudojant trečiąja savybe<br />

( )<br />

( )<br />

∫∫ ( ) ∫ ∫ ( ) ∫ ∫ ( )<br />

⎟<br />

( )<br />

( )<br />

⎟<br />

c ϕ ⎛ 2 x<br />

b ϕ ⎞ ⎛ 2 x ⎞<br />

f x,<br />

y dxdy = ⎜ f x,<br />

y dy⎟dx<br />

+ ⎜ f x,<br />

y dy dx . (4)<br />

⎜<br />

⎟ ⎜<br />

D<br />

a ⎝ ψ x ⎠ c ⎝ χ x ⎠<br />

d<br />

c<br />

y<br />

= ψ<br />

x 1<br />

( y)<br />

D<br />

L<br />

M<br />

N<br />

x = ψ 2<br />

( y)<br />

x<br />

ψ ( y )<br />

Integravimo sritį D galime aprašyti <strong>ir</strong><br />

kita tvarka. Dabar nagrinėjamoji sritis D ,<br />

( y ∈ c,<br />

d ) iš ka<strong>ir</strong>ės <strong>ir</strong> iš dešinės yra<br />

ribojamos, atitinkamai, linijomis x = ψ 1(<br />

y)<br />

<strong>ir</strong><br />

= ψ y , kurios abi sudaro liniją L .<br />

x 2<br />

( )<br />

( )<br />

( ( ) ( y)<br />

) .<br />

⎧y<br />

∈ c,<br />

d ,<br />

Todėl D : ⎨<br />

⎩x<br />

∈ ψ 1 y , ψ 2<br />

Taip aprašytai sričiai ∫∫ ( ) ∫ ∫ ( )<br />

⎟<br />

( )<br />

⎟<br />

d ⎛ 2 ⎞<br />

I = f x,<br />

y dxdy = ⎜ f x,<br />

y dx dy . (5)<br />

⎜<br />

D<br />

c ⎝ψ<br />

1 y ⎠<br />

Apskaičiuota dvilypio integralo reikšmė nepriklauso nuo pas<strong>ir</strong>inktos integravimo tvarkos.


I =<br />

⎛<br />

⎜<br />

⎝<br />

ϕ ( x)<br />

2<br />

( x,<br />

y)<br />

dxdy = ⎜ f ( x,<br />

y)<br />

b<br />

ϕ ( x )<br />

ψ ( y )<br />

ψ ( y )<br />

( x,<br />

y)<br />

∫∫ ∫ ∫ ∫ ∫ ⎟ ⎟<br />

f<br />

dy⎟dx<br />

= ⎜ f dx<br />

D<br />

a<br />

1<br />

⎞<br />

dy . (6)<br />

⎠<br />

Taigi, dvilypį integralą galima apskaičiuoti, pa<strong>kart</strong>otinai integruojant pagal vieną, o paskui<br />

pagal kitą kintamąjį. Kiekvienu konkrečiu atveju integravimo tvarka pas<strong>ir</strong>enkama, priklausomai<br />

z = f x,<br />

y <strong>ir</strong> nuo integravimo srities D .<br />

nuo integruojamos funkcijos ( )<br />

Pavyzdys. Pakeisti integravimo tvarką<br />

integrale ∫ ∫ ( )<br />

⎟ 1 x ⎛ ⎞<br />

I = ⎜ f x,<br />

y dy dx .<br />

⎜<br />

0 ⎝ x ⎠<br />

⎧x<br />

∈(<br />

0,<br />

1 ) ,<br />

> Duotajame integrale sritis yra D : ⎨<br />

⎩y<br />

∈(<br />

x,<br />

x ) .<br />

Pakeiskime integravimo tvarką, tada gauname<br />

⎧y<br />

∈ ( 0,<br />

1)<br />

,<br />

D : ⎨ 2<br />

⎩x<br />

∈ ( y , y)<br />

.<br />

∫ ∫ ( ) ∫ ∫ ( )<br />

⎟ ⎟<br />

1 ⎛ x ⎞ 1 ⎛ y ⎞<br />

I = ⎜ f x,<br />

y dy⎟dx<br />

= ⎜ f x,<br />

y dx dy .<br />

⎜ ⎟ ⎜ 2<br />

0 ⎝ x ⎠ 0 ⎝ y ⎠<br />

§ 3. DVILYPIS INTEGRALAS POLINĖJE KOORDINAČIŲ SISTEMOJE IR JO<br />

APSKAIČIAVIMAS<br />

A<br />

L<br />

D<br />

ρ = Φ1<br />

( θ )<br />

∆Si<br />

ρ = Φ2<br />

© A. Laurutis, D.Šiaučiūnas Paskaitų konspektas 86<br />

⎞<br />

⎟<br />

⎠<br />

d<br />

c<br />

⎛<br />

⎜<br />

⎝<br />

2<br />

1<br />

Tegul polinėje koordinačių sistemoje θ , ρ<br />

duota tokia sritis D , kad kiekviena tiesė,<br />

einanti iš poliaus P , <strong>ir</strong> kiekvienas<br />

apskritimas, kurio centras yra poliuje,<br />

kerta sritį nedaugiau kaip dviejuose<br />

taškuose. Iš brėžinio matome, kad sritis D<br />

yra apribota <strong>kreiv</strong>ėmis ρ 1 ( θ ) Φ = ,<br />

ρ 2(<br />

θ ) Φ = <strong>ir</strong> tiesėmis θ = α <strong>ir</strong> θ = β , čia<br />

Φ ≤ Φ θ <strong>ir</strong> α < β .<br />

1<br />

( θ ) 2(<br />

)<br />

⎧θ<br />

∈ ( α,<br />

β ) ,<br />

D : ⎨<br />

ρ ∈ ( ρ = Φ ( θ ) , ρ = Φ ( θ ) )<br />

P<br />

β<br />

∆ ϕi<br />

α<br />

ρ = ρi<br />

B ρ = ρi−1<br />

p<br />

Taigi<br />

⎩<br />

1<br />

2<br />

Tokia sritis D vadinama teisinga sritimi<br />

polinėje koordinačių sistemoje.<br />

Tegul šioje srityje D yra žinoma tolydi<br />

funkcija z = F(<br />

θ, ρ)<br />

polinėje koordinatėse<br />

Išdalinkime sritį D į sriteles ∆ s1 , ∆s2,<br />

∆s3,...,<br />

∆si,...,<br />

∆sn<br />

. Šiose sritelėse laisvai parenkame taškus<br />

z = F θ , ρ = F P reikšmes.<br />

P i , o juose apskaičiuojame funkcijos ( ) ( i )<br />

Sudarome integralinę sumą V =<br />

n<br />

f ( P ) ∆s<br />

n<br />

∑<br />

i=<br />

1<br />

Pagal dvilypio integralo prasmę, jei = F(<br />

θ , ρ ) ≥ 0<br />

∆sk → 0 ( n → ∞)<br />

V = F(<br />

ρ ) ds<br />

( θ )<br />

θ = θn<br />

∆θn<br />

max yra cilindroido tūris<br />

θ<br />

i<br />

= θn−1<br />

∆ρi<br />

∫∫<br />

D<br />

i<br />

y<br />

1<br />

O<br />

, čia P i yra bet kuris sritelės ∆ si<br />

taškas.<br />

D<br />

z , šios integralinės sumos riba, kai<br />

y =<br />

y = x,<br />

x = y<br />

θ, . (1)<br />

1<br />

x<br />

2<br />

x = y , 0 ≤ y ≤1<br />

x


x<br />

Kaip jau suprantame, integralinės sumos riba nepriklauso nuo srities D dalijimo į sriteles ∆ si<br />

būdo, todėl šią sritį galime išdalinti patogiausiu būdu: tiesėmis θ = θ0<br />

, θ = θ1,<br />

θ = θ2<br />

, ..., θ = θn<br />

,<br />

čiaθ 0 = α,<br />

θn<br />

= β , θ0<br />

< θ1<br />

< θ2<br />

< ... < θn<br />

<strong>ir</strong> apskritimais ρ = ρ0<br />

, ρ = ρ1,<br />

ρ = ρ2<br />

,..., ρ = ρn<br />

, čia<br />

ρ 0 < ρ1<br />

< ρ2<br />

< ... < ρn<br />

, o ρ 0 - minimali funkcijos Φ 1(<br />

θ, ρ)<br />

reikšmė <strong>ir</strong> ρ n - maksimali funkcijos<br />

Φ 2(<br />

θ, ρ ) reikšmė. Dėl šitokio srities D dalijimo į sriteles būdo, sritelės plotas<br />

∆ si ≈ ρi ⋅ ∆ρi<br />

⋅ ∆θi<br />

, kuris riboje, kai max ∆sk → 0 ( n → ∞)<br />

bus lygus ds = ρ ⋅ dρ<br />

⋅ dθ<br />

, todėl<br />

dvilypis integralas polinėje koordinačių sistemoje apskaičiuojamas <strong>kart</strong>otiniu integralu<br />

P<br />

V<br />

=<br />

2<br />

( ρ ) ds = ⎜ F(<br />

θ , ρ )<br />

© A. Laurutis, D.Šiaučiūnas Paskaitų konspektas 87<br />

( θ )<br />

β Φ ⎛<br />

⎞<br />

θ , dθ<br />

. (2)<br />

⎜<br />

α ⎝ Φ1(<br />

θ ) ⎠<br />

∫∫ ∫ ∫ ⎟ ⎟<br />

F<br />

ρdρ<br />

D<br />

Galime pas<strong>ir</strong>inkti <strong>ir</strong> kitą integravimo tvarką dvilypiame integrale polinėje koordinačių<br />

sistemoje. Kai pradžioje integruosime pagal θ , o paskui<br />

pagal ρ , galime naudoti kitą formulę:<br />

θ 2<br />

= ω<br />

( ρ )<br />

θ 1<br />

= ω<br />

( ρ )<br />

ρ1 ρ2<br />

O<br />

z<br />

p<br />

2 2 2<br />

x + y + z =<br />

2a<br />

V<br />

=<br />

ϖ ( θ )<br />

ρ ⎛<br />

⎞<br />

θ ρdρ<br />

. (3)<br />

⎜<br />

ρ1⎝<br />

ϖ1<br />

( θ ) ⎠<br />

2 2<br />

( , ρ ) ds = ⎜ F(<br />

θ , ρ )<br />

∫∫ ∫ ∫ ⎟ ⎟<br />

F<br />

dθ<br />

D<br />

Brėžinys ka<strong>ir</strong>ėje atitinka šią formulę. Galime nustatyti ryšį<br />

tarp dvilypio integralo polinėje <strong>ir</strong> stačiakampėje<br />

koordinačių sistemose. Kaip žinome<br />

⎧<br />

2 2<br />

⎧x<br />

= ρ cosθ<br />

,<br />

ρ = x + y ,<br />

⎪<br />

⎨<br />

<strong>ir</strong> ⎨<br />

⎩y<br />

= ρ sinθ<br />

,<br />

y<br />

⎪θ<br />

= arctg .<br />

⎩ x<br />

Be to ( x,<br />

y)<br />

f ( ρ cosθ<br />

, ρ sinθ<br />

)<br />

2<br />

x = ± 2ay − y<br />

Brėžinyje nupiešta<br />

pusė apskaičiuojamo<br />

kūno tūrio (kita<br />

pusė yra simetriška šiam kūnui xOy plokštumos atžvilgiu).<br />

Kūno tūris yra<br />

2<br />

2 2<br />

2 2 2<br />

2∫<br />

∫ 4<br />

.<br />

⎟<br />

0<br />

2<br />

2<br />

⎟<br />

2 2<br />

x + y −2ay=<br />

0<br />

y<br />

a⎛<br />

ay−<br />

y<br />

⎞<br />

V =<br />

⎜<br />

a − x − y dx dy<br />

⎜<br />

⎝ − ay−<br />

y<br />

⎠<br />

Gautąjį integralą lengviau apskaičiuosime polinėje<br />

4a<br />

2<br />

f = , todėl<br />

( θ )<br />

b<br />

∫∫ f ( x y)<br />

dxdy ∫ ∫ f ( ρ ϕ ρ ϕ)<br />

ρdρ<br />

dθ<br />

⎟<br />

D<br />

a ( θ )<br />

⎟<br />

Φ ⎛ 2<br />

⎞<br />

, = ⎜ cos , sin . (4)<br />

⎜<br />

⎝ Φ1<br />

⎠<br />

Pavyzdys. Apskaičiuoti tūrį kūno, gauto, susik<strong>ir</strong>tus<br />

2 2 2 2<br />

sferai x + y + z = 4a <strong>ir</strong> cilindrui<br />

2 2<br />

x + y − 2ay<br />

= 0 .<br />

> Sritis D yra tarp taškų y = 0 <strong>ir</strong> y = 2a<br />

, bei linijų<br />

P<br />

θ<br />

x, p<br />

a<br />

ρ<br />

2 a<br />

D<br />

1<br />

ρ = 2a sinθ<br />

y


koordinačių sistemoje. Žinome ryšio lygtis<br />

⎧<br />

2 2<br />

⎧x<br />

= ρ cosθ<br />

,<br />

ρ = x + y ,<br />

⎪<br />

⎨<br />

<strong>ir</strong> ⎨<br />

⎩y<br />

= ρ sinθ<br />

,<br />

y<br />

⎪θ<br />

= arctg .<br />

⎩ x<br />

2 2<br />

Cilindro x + y − 2ay<br />

= 0 lygtis polinėje koordinačių sistemoje yra ρ = 2a sinθ<br />

.<br />

Sritis D yra<br />

π<br />

2<br />

⎧ ⎛ π ⎞<br />

⎪θ<br />

∈⎜<br />

0,<br />

⎟,<br />

D : ⎨ ⎝ 2 ⎠ pointegralinė funkcija<br />

⎪<br />

⎩ρ<br />

∈ ( 0,<br />

2a<br />

sinθ<br />

) ,<br />

2 1<br />

⎛<br />

⎜<br />

⎝<br />

= ∫ ∫<br />

∫ ∫<br />

⎞<br />

⎟<br />

⎠<br />

2 2<br />

2 2<br />

todėl V 4 4a<br />

− ρ ρdρ<br />

dθ<br />

= −2<br />

4a<br />

− ρ ( − 2ρ<br />

) dρ<br />

dθ<br />

=<br />

4<br />

= −<br />

3<br />

π<br />

2<br />

∫<br />

0<br />

32a<br />

=<br />

3<br />

⎡<br />

⎢<br />

⎣<br />

π<br />

3<br />

0<br />

2a<br />

sinθ<br />

0<br />

π<br />

2<br />

0<br />

⎛<br />

⎜<br />

⎝<br />

2a<br />

sinθ<br />

3<br />

2<br />

3<br />

2 2 ⎤ 4 ⎡ 2 2 2<br />

2<br />

( 4a<br />

− ρ ) 2<br />

( 4 4 sin ) 2<br />

⎥ dθ<br />

= − ∫ ⎢ a − a θ − ( 4a<br />

)<br />

0<br />

2 2 2<br />

2 2<br />

z = 4a − x − y = 4a<br />

− ρ ,<br />

2<br />

3 2<br />

3 2<br />

3<br />

3<br />

3 32a<br />

32a<br />

3 16a<br />

π 32a<br />

2!<br />

! 16 3<br />

∫ ( 1−<br />

cos θ ) dθ<br />

= ∫ dθ<br />

− ∫cos<br />

θdθ<br />

= − = a ( 3π<br />

− 4).<br />

0<br />

⎦<br />

2a<br />

sinθ<br />

0<br />

3<br />

π<br />

0<br />

3<br />

π<br />

0<br />

⎣<br />

3<br />

π<br />

0<br />

© A. Laurutis, D.Šiaučiūnas Paskaitų konspektas 88<br />

3<br />

3<br />

2<br />

3<br />

⎞<br />

⎟<br />

⎠<br />

⎤<br />

⎥dθ<br />

=<br />

⎦<br />

§ 4. TRILYPIS INTEGRALAS, JO SAVYBĖS IR APSKAIČIAVIMAS STAČIAKAMPĖJE<br />

KOORDINAČIŲ SISTEMOJE<br />

Tegul trimatėje erdvėje duota tam tikra sritis V (erdvinio kūnas) apribotas uždaru<br />

<strong>pav</strong><strong>ir</strong>šiumi S (erdvinis kūno <strong>pav</strong><strong>ir</strong>šius). Tegul šioje srityje V <strong>ir</strong> jo <strong>pav</strong><strong>ir</strong>šiuje S žinoma tam tikra<br />

tolydi funkcija w = f ( x,<br />

y,<br />

z)<br />

, čia x , y,<br />

z stačiakampės srities V koordinatės. Tuo atveju, kai<br />

f ( x,<br />

y,<br />

z)<br />

≥ 0 , šią funkciją galime įsivaizduoti erdvinio kūno tankio funkcija.<br />

Išdalinkime V į tam tikras daleles ∆ Vi<br />

, pažymėdami ∆ Vi<br />

ne tik atitinkamą srities V dalelę, bet<br />

<strong>ir</strong> jos tūrį. Kiekvienoje ∆Vi dalelėje paimkime tašką P i <strong>ir</strong> šiame taške apskaičiuokime funkcijos<br />

= f ( x,<br />

y z)<br />

reikšmę, t. y. f ( x y,<br />

z)<br />

= f ( P )<br />

w ,<br />

P<br />

i<br />

i<br />

3!<br />

!<br />

, . Sudarykime integralinę sumą<br />

∑ f ( Pi<br />

) ∆Vi<br />

i=1<br />

<strong>ir</strong> apskaičiuokime jos ribą, kai ∆V → 0,<br />

( n → ∞)<br />

diam i<br />

n<br />

max . Ši riba vadinama trilypiu integralu<br />

n<br />

∑<br />

max diam∆Vi<br />

→0<br />

i=<br />

1<br />

( P ) ∆V<br />

= f ( P)<br />

i<br />

i<br />

∫∫∫<br />

lim f<br />

dv arba, (1)<br />

kai skaidymas į sriteles vykdomas plokštumomis, lygiagrečiomis koordinačių plokštumoms<br />

lim<br />

n<br />

∑<br />

max diam∆Vi<br />

→0<br />

i=<br />

1<br />

f<br />

( P ) ∆V<br />

= f ( x,<br />

y,<br />

z)<br />

Jeigu ( x,<br />

y,<br />

z)<br />

≥ 0<br />

kūno masę M . Atsk<strong>ir</strong>u atveju, kai f ( x y,<br />

z)<br />

= const = ρ<br />

i<br />

i<br />

∫∫∫<br />

V<br />

V<br />

dxdydz<br />

f yra erdvinio kūno tankio funkcija, tai integralas (1) apskaičiuoja erdvinio<br />

, , tai<br />

9<br />

(2)


∫∫∫<br />

M = ρ dxdydz . (3)<br />

V<br />

Tegul trimatis erdvinis kūnas V yra apribotas uždaru <strong>pav</strong><strong>ir</strong>šiumi S , pasižyminčiu savybėmis:<br />

1) kiekviena tiesė, lygiagreti Oz ašiai, <strong>ir</strong> einanti per srities V vidinį tašką, kerta <strong>pav</strong><strong>ir</strong>šių S<br />

dviejuose taškuose,<br />

2) visa sritis V , suprojektuota į xOy plokštumą, reiškia teisingą dvimatę sritį D ,<br />

3) kiekviena srities V dalis, atk<strong>ir</strong>sta plokštumomis, lygiagrečiomis koordinačių plokštumoms<br />

xOy , xOz <strong>ir</strong> yOz , taip pat pasižymi 1 <strong>ir</strong> 2 savybėmis.<br />

Tokia sritis V vadinama teisinga trimate sritimi.<br />

b<br />

x<br />

a<br />

z<br />

y = ϕ1<br />

Funkcijos w f ( x,<br />

y,<br />

z)<br />

Aprašykime šią sritį V . Iš brėžinio<br />

matome, kad sritį V iš apačios <strong>ir</strong> iš<br />

v<strong>ir</strong>šaus riboja du sk<strong>ir</strong>tingi <strong>pav</strong><strong>ir</strong>šiai<br />

z = ψ 1 ( x , y ) <strong>ir</strong> z = ψ 2(<br />

x,<br />

y)<br />

,<br />

atitinkamai. Plokščią sritį D iš ka<strong>ir</strong>ės <strong>ir</strong> iš<br />

dešinės riboja linijos y = ϕ1(<br />

x)<br />

<strong>ir</strong><br />

= ϕ x , atitinkamai, kai x ∈ a,<br />

b .Tada<br />

y 2(<br />

)<br />

⎧x<br />

∈ ( a,<br />

b)<br />

,<br />

⎪<br />

V : ⎨y<br />

∈ ( ϕ 1(<br />

x)<br />

, ϕ 2 ( x)<br />

) ,<br />

⎪<br />

⎩z<br />

∈ ( z = ψ 1(<br />

x,<br />

y)<br />

, z = ψ 2(<br />

x,<br />

y)<br />

) ,<br />

<strong>ir</strong> funkcijos w f ( x,<br />

y,<br />

z)<br />

=<br />

integralas srityje V yra<br />

trigubas<br />

ϕ ( x)<br />

ψ ( x,<br />

y )<br />

ϕ ( x)<br />

ψ ( x,<br />

y )<br />

( x,<br />

y,<br />

z)<br />

∫ ∫ ∫ ⎟ ⎟<br />

⎜ ⎜ f dz⎟dy<br />

© A. Laurutis, D.Šiaučiūnas Paskaitų konspektas 89<br />

I<br />

V<br />

=<br />

b<br />

a<br />

⎛<br />

⎜<br />

⎝<br />

2<br />

1<br />

⎛<br />

⎜<br />

⎝<br />

2<br />

1<br />

⎞<br />

⎟<br />

⎠<br />

⎞<br />

dx .<br />

⎠<br />

= trilypis integralas teisingoje srityje V yra lygus trigubam integralui<br />

toje pačioje srityje<br />

( ) ( )<br />

∫∫∫ ( ) ∫ ∫ ∫ ( )<br />

⎟<br />

( ) ( )<br />

⎟<br />

b ⎛ϕ<br />

2 x ψ ⎛ 2 x,<br />

y ⎞ ⎞<br />

I =<br />

= ⎜ ⎜<br />

⎟<br />

V f x,<br />

y,<br />

z dv<br />

f x,<br />

y,<br />

z dz dy dx . (4)<br />

⎜ ⎜<br />

⎟<br />

V<br />

a ⎝ ϕ1<br />

x ⎝ψ<br />

1 x,<br />

y ⎠ ⎠<br />

Galime nustatyti tokias trilypio integralo savybes:<br />

1. Jeigu sritį V perk<strong>ir</strong>sime plokštuma, lygiagrečia kuriai nors iš koordinačių plokštumų, į dvi<br />

sritis 1 V <strong>ir</strong> 2 V , tai trilypis integralas srityje V yra lygus trilypių integralų srityse 1 V <strong>ir</strong> V 2 sumai.<br />

Šią savybę galima išplėsti baigtiniam sričių V 1 ,..., Vn<br />

skaičiui.<br />

I = I + I + ... + I<br />

V<br />

V<br />

1<br />

2. (Trilypio integralo įvertinimas). Jei m <strong>ir</strong> M yra, atitinkamai, mažiausia <strong>ir</strong> didžiausia funkcijos<br />

w = f x,<br />

y,<br />

z reikšmės srityje V , tai mV ≤ IV<br />

≤ MV , čia V nagrinėjamos srities V tūris.<br />

( )<br />

( x)<br />

V<br />

D<br />

y = ϕ2<br />

( x)<br />

( x y)<br />

z = ψ ,<br />

2<br />

y<br />

( x y)<br />

z =<br />

ψ ,<br />

1<br />

3. (Vidurkio savybė). Trigubas integralas srityje V yra lygus šios srities tūriui V , padaugintam iš<br />

funkcijos w = f ( x,<br />

y,<br />

z)<br />

reikšmės, apskaičiuotos tam tikrame srities taške P .<br />

IV = f ( P)V<br />

.<br />

V<br />

2<br />

V<br />

n


§ 5. TRILYPIO INTEGRALO APSKAIČIAVIMAS CILINDRINĖJE IR SFERINĖJE<br />

KOORDINAČIŲ SISTEMOSE<br />

Jau žinome, kad funkcijos w f ( x,<br />

y,<br />

z)<br />

= trilypis integralas srityje V , aprašytoje<br />

( a,<br />

b)<br />

,<br />

( ϕ 1(<br />

x)<br />

, ϕ 2 ( x)<br />

) ,<br />

( z = ψ ( x,<br />

y)<br />

, z = ψ ( x,<br />

y)<br />

) ,<br />

⎧x<br />

∈<br />

⎪<br />

stačiakampėje koordinačių sistemoje reiškiniu V : ⎨y<br />

∈<br />

⎪<br />

⎩z<br />

∈<br />

apskaičiuojamas, pakeičiant jį trigubu integralu<br />

1<br />

2<br />

( ) ( )<br />

∫∫∫ ( ) ∫ ∫ ∫ ( )<br />

⎟<br />

( ) ( )<br />

⎟<br />

b ⎛ ϕ2<br />

x ψ ⎛ 2 x,<br />

y ⎞ ⎞<br />

I =<br />

= ⎜ ⎜<br />

⎟<br />

V f x,<br />

y,<br />

z dv<br />

f x,<br />

y,<br />

z dz dy dx .<br />

⎜ ⎜<br />

⎟<br />

V<br />

a ⎝ ϕ1<br />

x ⎝ψ<br />

1 x,<br />

y ⎠ ⎠<br />

Būtina pabrėžti, kad integravimo tvarka gali būti <strong>ir</strong> kitokia, priklausomai nuo konkretaus<br />

sprendžiamo uždavinio.<br />

1. Cilindrinėje koordinačių sistemoje<br />

taško M padėtis erdvėje nusakoma<br />

koordinatėmis θ , ρ,<br />

z , t. y. M ( θ , ρ,<br />

z)<br />

, čia θ<br />

<strong>ir</strong> ρ polinės taško A koordinatės<br />

plokštumoje, o z taško M aplikatė.<br />

x, p<br />

z<br />

O, P<br />

θ<br />

ρ<br />

Ryšio lygtys tarp cilindrinių <strong>ir</strong> stačiakampių<br />

⎧x<br />

= ρ cosθ<br />

,<br />

⎪<br />

koordinačių yra ⎨y<br />

= ρ sinθ<br />

,<br />

⎪<br />

⎩z<br />

= z.<br />

Cilindrinė koordinačių sistemoje dažniausiai naudojama tais atvejais, kai funkcija f ( x,<br />

y,<br />

z)<br />

arba sritis V yra cilindriniai bei sferiniai <strong>pav</strong><strong>ir</strong>šiai.<br />

Kadangi ∆v = ρ ∆θ∆ρ∆z<br />

, tai dv = ρdθdρdz<br />

<strong>ir</strong><br />

∫∫∫<br />

čia f ( x,<br />

y,<br />

z)<br />

f ( ρ cosθ<br />

, ρ sinθ<br />

, z)<br />

= F(<br />

θ , ρ,<br />

z)<br />

V<br />

( x y,<br />

z)<br />

dv = F(<br />

θ , ρ,<br />

z)<br />

= .<br />

© A. Laurutis, D.Šiaučiūnas Paskaitų konspektas 90<br />

∫∫∫<br />

f , ρdθdρdz<br />

, (1)<br />

Pavyzdys. Apskaičiuoti pusrutulio masę M , jei rutulio centras yra koordinačių pradžioje,<br />

spindulys – R , o tankio funkcija – F = kz .<br />

> V<strong>ir</strong>šutinės sferos dalies <strong>pav</strong><strong>ir</strong>šiaus lygtis yra<br />

Cilindrinėje koordinačių sistemoje<br />

z<br />

( θ,<br />

ρ z)<br />

M ,<br />

N ( θ, ρ )<br />

y<br />

V<br />

x, p<br />

z −<br />

2 2 2<br />

= R − x y .<br />

2 2<br />

2 2<br />

z = R − ρ , nes = x + y<br />

z<br />

O, P<br />

∆θ<br />

∆ρ<br />

ρ∆θ<br />

∆z<br />

ρ , todėl<br />

∆V = ρ<br />

⋅ ∆θ<br />

⋅ ∆ρ<br />

⋅ ∆z<br />

y


M<br />

=<br />

=<br />

⎛<br />

⎜<br />

⎝<br />

2π<br />

R<br />

∫∫<br />

0 0<br />

∫∫∫<br />

V<br />

k<br />

2<br />

kzρdθdρdz<br />

=<br />

2π<br />

R<br />

0 0 0<br />

⎞<br />

2π<br />

R<br />

2π<br />

2 4<br />

2π<br />

4 4<br />

4<br />

4<br />

2 2 k 2 ρ ρ k ⎛ R R ⎞ k R kπR<br />

( R − ρ ) ρdρ<br />

⎟dθ<br />

⎜ R<br />

⎟<br />

⎟<br />

= ∫ − dθ<br />

= ∫ ⎜ − ⎟<br />

⎟dθ<br />

= 2π<br />

= . <<br />

⎠<br />

⎛<br />

⎜<br />

⎜<br />

⎝<br />

⎛<br />

⎜<br />

⎜<br />

⎝<br />

2 2<br />

R −ρ<br />

⎞ ⎞<br />

kzdz<br />

⎟ ⎟<br />

ρdρ<br />

dθ<br />

=<br />

⎟ ⎟<br />

⎠ ⎠<br />

⎛<br />

⎜<br />

⎜<br />

⎝<br />

⎛ R R ⎞<br />

2 ⎜<br />

0 ⎝<br />

2 0 4 ⎟<br />

0 ⎠<br />

∫∫ ∫<br />

© A. Laurutis, D.Šiaučiūnas Paskaitų konspektas 91<br />

∫∫<br />

kz<br />

2<br />

2<br />

0 0 0<br />

2<br />

2<br />

R −<br />

0<br />

⎝<br />

ρ<br />

2<br />

2<br />

⎞<br />

⎟<br />

ρdρ<br />

⎟dθ<br />

=<br />

⎟<br />

⎠<br />

2. Sferinėje koordinačių sistemoje taško M padėtį erdvėje nusakome koordinatėmis θ ,r , ϕ ,<br />

0 ≤ r < ∞ , ϕ – kampas, kurį vektorius<br />

0 ≤ ϕ ≤ π , θ – kampas, kurį vektoriaus OM r<br />

projekcija į xOy plokštumą sudaro su<br />

z<br />

teigiama Ox ašies<br />

čia r – taško atstumas nuo koordinačių pradžios taško ( )<br />

→<br />

OM sudaro su teigiama Oz ašies kryptimi ( )<br />

x<br />

z = r cosϕ<br />

x = rsinϕ<br />

sinθ<br />

O<br />

θ<br />

r<br />

ϕ<br />

( θ , ρ z)<br />

M ,<br />

Bet kuriam erdvės taškui M ( x,<br />

y,<br />

z)<br />

⎧x<br />

= r sinϕ<br />

cosθ<br />

,<br />

⎪<br />

⎨y<br />

= r sinϕ<br />

sinθ<br />

,<br />

⎪<br />

⎩z<br />

= r cosϕ.<br />

2<br />

Galime rasti <strong>ir</strong> dv = r sinϕdrdθdϕ<br />

, tada<br />

∫∫∫<br />

V<br />

y = rsinϕ<br />

sinθ<br />

ρ = rsinϕ<br />

A<br />

( θ ,r sinϕ<br />

)<br />

( x y,<br />

z)<br />

dv = F(<br />

r,<br />

θ , ϕ )<br />

∫∫∫<br />

V<br />

4<br />

⎠<br />

2<br />

4<br />

kryptimi ( 0 θ ≤ 2π<br />

)<br />

≤ .<br />

2<br />

f , r drdθdϕ<br />

. (2)<br />

§ 6. KARTOTINIŲ INTEGRALŲ TAIKYMAS MECHANIKOJE<br />

1. Dvilypio integralo taikymas<br />

a) plokštelės inercijos momentas<br />

Materialaus taško M , kurio masė yra m , inercijos momentas koordinačių pradžios taško O<br />

2<br />

atžvilgiu yra I = mr , čia r = OM . Tada materialiųjų taškų M i , turinčių mases m i , sistemos<br />

inercijos momentą koordinačių pradžios taško O atžvilgiu apskaičiuosime sumuodami<br />

I<br />

n<br />

=<br />

n<br />

∑<br />

i=<br />

1<br />

2<br />

m r . (1)<br />

Dabar apskaičiuokime plokštelės inercijos momentą taško O atžvilgiu<br />

y<br />

x<br />

i i<br />

∆ϕ<br />

O<br />

θ<br />

z<br />

ϕ<br />

∆θ<br />

r<br />

∆ r<br />

4<br />

2<br />

∆V = r sinϕ∆r∆θ∆ϕ<br />

y


Plokštelę, kurios forma atitinka sritį D , o tankis γ = const = 1 , išdaliname į elementarias sriteles<br />

y<br />

D<br />

∆ Si<br />

( i = 1,<br />

2,...,<br />

n)<br />

. Kiekvienoje ∆ Si<br />

sritelėje<br />

paimame tašką Pi ( ξ i,<br />

ηi<br />

) . Tokios elementarios<br />

sritelės inercijos momentas yra<br />

2<br />

2 2<br />

∆Ii<br />

= ri<br />

∆Si<br />

= ( ξ i + ηi<br />

) ∆Si<br />

.<br />

Pi<br />

Susumuokime visus šiuos elementarių sritelių<br />

ηi<br />

∆Si<br />

i ∈ 1,<br />

n momentus<br />

I<br />

0<br />

2 2 ( x y )<br />

= ∫∫ +<br />

D<br />

( )<br />

2 2 ( ξ + η )<br />

© A. Laurutis, D.Šiaučiūnas Paskaitų konspektas 92<br />

n<br />

∑<br />

i=<br />

1<br />

i<br />

i<br />

∆S<br />

<strong>ir</strong> apskaičiuokime šios sumos ribą, kai<br />

max diam ∆Si → 0 ( n → ∞)<br />

. Ši sumos riba yra<br />

plokštelės inercijos momentas koordinačių<br />

pradžios taško O atžvilgiu.<br />

I<br />

2 2 ( + )<br />

n<br />

0 = lim ∑ i ηi<br />

max diam<br />

∆S<br />

→ i=<br />

1<br />

i 0<br />

i<br />

ξ ∆S<br />

arba<br />

dxdy . (2)<br />

Šios plokštelės inercijos momentus ašių Ox <strong>ir</strong> Oy atžvilgiu galime apskaičiuoti analogiškai.<br />

I xx<br />

2<br />

= ∫∫ y dxdy , I yy = ∫∫<br />

D<br />

D<br />

2<br />

x dxdy . (3),(4)<br />

Pavyzdys. Apskaičiuokite skritulio, kurio centras yra koordinačių pradžioje, o spindulys – R ,<br />

inercijos momentą jo centro atžvilgiu.<br />

2 2<br />

> Kadangi I0<br />

= ∫∫ ( x + y ) dxdy , tikslinga integruoti polinėje koordinačių sistemoje, tada<br />

( 0,<br />

2π<br />

) ,<br />

( 0,<br />

R)<br />

,<br />

D<br />

2π<br />

R<br />

4<br />

⎧θ<br />

∈<br />

⎛ ⎞ 2 πR<br />

D : ⎨ <strong>ir</strong> I0<br />

= ρ ρdρ<br />

dϕ<br />

=<br />

⎩ρ<br />

∈ ∫∫ ⎜ ⎟<br />

⎜ ⎟<br />

.<br />

2<br />

0⎝0 ⎠<br />

Išvada: „vilkelio“ sukimosi apie ašį inercijos momentas yra proporcingas<br />

Pastaba. Jei plokštelės tankis nėra lygus 1, bet kinta, t.y. γ ( x, y)<br />

I<br />

0<br />

=<br />

ri<br />

O ξi<br />

∫∫<br />

D<br />

2 2<br />

( x,<br />

y)(<br />

x + y )<br />

2<br />

γ dxdy , I = ( x,<br />

y)<br />

y dxdy<br />

b) plokštelės masės centras<br />

Materialių taškų P , P2<br />

,..., Pi<br />

,..., Pn<br />

koordinates apskaičiuojame pagal formules<br />

xx<br />

∫∫<br />

D<br />

γ = , tai<br />

∫∫<br />

γ , I = ( x,<br />

y)<br />

1 , turinčių mases m , m2,...,<br />

mi<br />

,..., mn<br />

∑<br />

∑<br />

x<br />

i<br />

∑<br />

∑<br />

yy<br />

D<br />

4<br />

R .<<br />

2<br />

γ x dxdy . (5),(6),(7)<br />

1 sistemos masės centro<br />

ximi<br />

yimi<br />

x c = , y c = . (8)<br />

m<br />

m<br />

Dabar apskaičiuokime plokštelės masės centro koordinates. Analogiškai aukščiau užrašytoms<br />

formulėms (8) galime gauti<br />

i<br />

i


x<br />

n<br />

∑<br />

i=<br />

1<br />

c ≈ n<br />

ξi∆Si<br />

, y<br />

∆S<br />

∑<br />

i=<br />

1<br />

i<br />

© A. Laurutis, D.Šiaučiūnas Paskaitų konspektas 93<br />

n<br />

∑<br />

i=<br />

1<br />

c ≈ n<br />

η ∆S<br />

Šių integralinių sumų ribos, kai max ∆Si<br />

artėja į nulį, atitinka dvilypius integralus skaitikliuose <strong>ir</strong><br />

vardikliuose. Tada<br />

∫∫<br />

D<br />

∫∫<br />

D<br />

xdxdy<br />

∑<br />

i=<br />

1<br />

∫∫<br />

D<br />

∫∫<br />

i<br />

∆S<br />

i<br />

i<br />

.<br />

x c = , y c = . (9)<br />

dxdy dxdy<br />

D<br />

ydxdy<br />

Tuo atveju, kai plokštelės medžiagos tankis yra nusakomas funkcija γ ( x, y)<br />

∫∫<br />

x<br />

c<br />

=<br />

∫∫<br />

D<br />

∫∫<br />

D<br />

xγ<br />

γ<br />

( x,<br />

y)<br />

( x,<br />

y)<br />

dxdy<br />

dxdy<br />

čia išraiškos x ( x,<br />

y)<br />

dxdy <strong>ir</strong> yγ<br />

( x,<br />

y)<br />

D<br />

∫∫<br />

D<br />

,<br />

y<br />

c<br />

=<br />

∫∫<br />

D<br />

∫∫<br />

D<br />

yγ<br />

γ<br />

( x,<br />

y)<br />

( x,<br />

y)<br />

dxdy<br />

dxdy<br />

γ = , tai<br />

, (10)<br />

γ dxdy vadinamos plokščios figūros D statiniais<br />

momentais (sukimo momentas apie ašis y <strong>ir</strong> x<br />

γ x, y dxdy - plokštelės D masė.<br />

∫∫<br />

D<br />

atitinkamai), o ( )<br />

Pavyzdys. Apskaičiuoti ketv<strong>ir</strong>tadalio elipsės <strong>pav</strong>idalo<br />

plokštelės masės centro koordinates, kai γ = 1.<br />

> x<br />

c<br />

=<br />

a<br />

⎛<br />

⎜<br />

⎜<br />

⎜<br />

⎝<br />

⎛<br />

⎜<br />

⎜<br />

⎜<br />

⎝<br />

∫ ∫<br />

a<br />

b 2<br />

a<br />

a<br />

0 0<br />

b<br />

a<br />

a<br />

∫ ∫<br />

0 0<br />

2 2 ( a − x )<br />

3<br />

2<br />

2<br />

−x<br />

⎞<br />

⎟<br />

xdy⎟dx<br />

⎟<br />

⎠<br />

=<br />

2 2<br />

−x<br />

⎞<br />

⎟<br />

dy⎟dx<br />

⎟<br />

⎠<br />

a<br />

b<br />

a<br />

− 4<br />

0 4a<br />

= , y<br />

2<br />

3πa<br />

3π<br />

c<br />

a<br />

∫<br />

0<br />

=<br />

a<br />

2<br />

− x<br />

1<br />

πab<br />

4<br />

a<br />

⎛<br />

⎜<br />

⎜<br />

⎜<br />

⎝<br />

b<br />

a<br />

∫ ∫<br />

0 0<br />

2) Trilypio integralo taikymas<br />

a) kūno inercijos momentas<br />

Taško M ( x y,<br />

z)<br />

2<br />

xdx<br />

=<br />

2 2<br />

a −x<br />

⎞<br />

⎟<br />

ydy⎟dx<br />

⎟<br />

⎠ 4b<br />

= .<<br />

1<br />

πab<br />

3π<br />

4<br />

, , kurio masė yra m , inercijos momentas koordinačių ašių Ox, Oy <strong>ir</strong> Oz<br />

2 2<br />

2 2<br />

2 2<br />

atžvilgiu apskaičiuojamas pagal formules ( ) , ( ) , ( ) .<br />

I xx = y + z m I yy = x + z m I zz = x<br />

Kūno inercijos momentai šių ašių atžvilgiu apskaičiuojami trilypiais <strong>integralai</strong>s<br />

b<br />

O<br />

y<br />

b<br />

y =<br />

a −<br />

a<br />

2 2<br />

x y<br />

+ = 1<br />

2 2<br />

a b<br />

2 2<br />

x<br />

a<br />

+ y<br />

x<br />

m


I<br />

I<br />

I<br />

xx<br />

yy<br />

zz<br />

=<br />

=<br />

=<br />

2 2<br />

∫∫∫(<br />

y + z ) γ ( x,<br />

y,<br />

z)<br />

V<br />

2 2<br />

∫∫∫(<br />

x + z ) γ ( x,<br />

y,<br />

z)<br />

V<br />

2 2<br />

∫∫∫(<br />

x + y ) γ ( x,<br />

y,<br />

z)<br />

V<br />

čia ( x , y,<br />

z)<br />

dxdydz<br />

dxdydz , (11)<br />

dxdydz<br />

γ medžiagos tankio funkcija.<br />

b) erdvinio kūno masės centro<br />

koordinatės<br />

Analogiškai plokštelės masės centro<br />

koordinačių formulėms, galime užrašyti<br />

kūno masės centro koordinačių formules<br />

x<br />

c<br />

∫∫∫<br />

=<br />

∫∫∫<br />

V<br />

∫∫∫<br />

V<br />

xγ<br />

γ<br />

( x,<br />

y,<br />

z)<br />

( x,<br />

y,<br />

z)<br />

dxdydz<br />

dxdydz<br />

∫∫∫<br />

,<br />

y<br />

c<br />

=<br />

∫∫∫<br />

V<br />

∫∫∫<br />

V<br />

yγ<br />

γ<br />

( x,<br />

y,<br />

z)<br />

( x,<br />

y,<br />

z)<br />

dxdydz<br />

dxdydz<br />

čia x γ ( x,<br />

y,<br />

z)<br />

dxdydz,<br />

yγ<br />

( x,<br />

y,<br />

z)<br />

dxdydz,<br />

z ( x,<br />

y,<br />

z)<br />

V<br />

momentai apie ašis Ox Oy,<br />

Oz<br />

V<br />

∫∫∫<br />

V<br />

, , o ( x , y,<br />

z)<br />

2<br />

y +<br />

© A. Laurutis, D.Šiaučiūnas Paskaitų konspektas 94<br />

∫∫∫<br />

V<br />

,<br />

z<br />

c<br />

=<br />

∫∫∫<br />

V<br />

∫∫∫<br />

V<br />

zγ<br />

γ<br />

( x,<br />

y,<br />

z)<br />

( x,<br />

y,<br />

z)<br />

dxdydz<br />

dxdydz<br />

, (12)<br />

γ dxdydz yra statiniai (sukimo)<br />

γ dxdydz – kūno masė, atitinkamai.<br />

§ 7. SKALIARINIAI IR VEKTORINIAI LAUKAI<br />

Sakoma, kad erdvėje yra žinomas laukas, jei kiekviename erdvės taške žinoma šio lauko<br />

reikšmė. Tokio lauko <strong>pav</strong>yzdžiu gali būti temperatūros laukas (kiekviename taške temperatūra<br />

apibrėžta) arba krūvio sudaromas elektrinis laukas. Laukai gali būti skaliariniai <strong>ir</strong> vektoriniai,<br />

priklausomai nuo nagrinėjamo dydžio charakterio (temperatūros laukas yra skaliarinis laukas, o<br />

elektrinis – vektorinis). Be to, laukas gali būti stacionarus (nusistovėjęs), jei kiekviename erdvės<br />

taške lauko reikšmė nekinta laike <strong>ir</strong> nestacionarus (nenusistovėjęs), jei lauko reikšmė priklauso<br />

nuo laiko momento.<br />

Imkime skaliarinį lauką, kurio dydis kiekviename erdvės taške M yra u . Galime užrašyti<br />

funkcinį ryšį u = f ( M ) . Tai, suprantama, kitokio pobūdžio funkcija, lyginant ją su anksčiau<br />

nagrinėtomis, kadangi taškas nėra dydis. Bet juk funkcinį ryšį galime nagrinėti platesne prasme.<br />

Tai vienos rūšies „objektų“ priklausomybė nuo kitos rūšies objektų. Tuo atveju, kai laukas<br />

nestacionarus u = f ( M , t)<br />

.<br />

Jei taškas M = M ( x,<br />

y,<br />

z)<br />

yra trimatės erdvės taškas, o u = f ( x,<br />

y,<br />

z)<br />

priklauso nuo visų šios<br />

erdvės koordinačių tai šis laukas yra erdvinis <strong>ir</strong> jį tenka nagrinėti stačiakampėje arba kurioje nors<br />

kitoje koordinačių sistemoje. Jeigu nagrinėjamas laukas nepriklauso nuo vienos iš koordinačių, jį<br />

galime įsivaizduoti plokščiu lauku, t.y. lauku plokštumoje.<br />

§ 8. LYGIO LINIJOS IR PAVIRŠIAI<br />

= , tai sakoma, kad<br />

šioje srityje yra žinomas skaliarinis laukas (temperatūrų skaliarinis laukas, slėgių skaliarinis<br />

laukas <strong>ir</strong> pan.).<br />

Jeigu trimatėje erdvėje yra sritis V , kurioje žinoma funkcija u f ( x,<br />

y,<br />

z)<br />

x<br />

z<br />

2<br />

z<br />

y<br />

z<br />

O<br />

x<br />

2<br />

x +<br />

( x,<br />

y z)<br />

M ,<br />

y<br />

2<br />

x +<br />

2<br />

y<br />

z<br />

2


Imkime srities V taškus, kuriuose funkcija u f ( x,<br />

y,<br />

z)<br />

= turi pastovią reikšmę c .<br />

( x,<br />

y z)<br />

c<br />

u = f , =<br />

Tokių taškų aibė sudaro <strong>pav</strong><strong>ir</strong>šių f ( x,<br />

y,<br />

z)<br />

= c . Kai paimsime kitą c reikšmę, gausime kitą<br />

<strong>pav</strong><strong>ir</strong>šių. Šie <strong>pav</strong><strong>ir</strong>šiai vadinami lygio <strong>pav</strong><strong>ir</strong>šiais.<br />

2 2 2<br />

x y z<br />

Pavyzdys. Apskaičiuokime skaliarinį lauko f ( x,<br />

y,<br />

z)<br />

= + + lygio <strong>pav</strong><strong>ir</strong>šius.<br />

4 9 16<br />

2 2 2<br />

x y z<br />

> Lygio <strong>pav</strong><strong>ir</strong>šiai yra + + = c , t.y. elipsoidai.<<br />

4 9 16<br />

x<br />

u = z<br />

O<br />

x<br />

( x y)<br />

u = z = f ,<br />

( x,<br />

y)<br />

c<br />

f =<br />

§ 9. GRADIENTAS IR JO SAVYBĖS<br />

Tegul erdvėje yra duotas stacionarus laukas bus u u(<br />

x,<br />

y,<br />

z)<br />

= yra<br />

plokščias laukas, o lauko u reikšmes<br />

(lauko intensyvumas) atidėsime z<br />

ašyje, tai laukas u = z = f ( x,<br />

y)<br />

yra<br />

dviejų kintamųjų funkcija <strong>ir</strong> galime<br />

kalbėti apie tokio lauko lygio linijas<br />

f ( x,<br />

y)<br />

= c .<br />

Šiuo atveju lygio linijos bus z = c ,<br />

t.y. <strong>pav</strong><strong>ir</strong>šiaus z = f ( x,<br />

y)<br />

susik<strong>ir</strong>timo<br />

linijos su plokštumomis z = c . Žinant<br />

lygio linijas, patogu nagrinėti lauko<br />

charakterį.<br />

Tuo atveju, kai u f ( x,<br />

y,<br />

z)<br />

= . Šioje erdvėje yra <strong>ir</strong> tam tikra<br />

linija L , einanti per tašką M , kuriame linijos L liestinės kryptis yra →<br />

l . Lauko u išvestinė pagal<br />

kryptį →<br />

l reiškia šio lauko kitimo greitį duotąja kryptimi<br />

∂u<br />

u(<br />

N ) − u(<br />

M )<br />

= lim . (1)<br />

∂l<br />

∆s→0 ∆s<br />

Kai linija L duota parametriniu <strong>pav</strong>idalu r r () s<br />

→ →<br />

l<br />

= , tai<br />

laukas u išilgai linijos L yra sudėtinė kelio s funkcija<br />

τ<br />

N L<br />

u ( s)<br />

= u(<br />

x(<br />

s)<br />

, y(<br />

s)<br />

, z(<br />

s)<br />

) .<br />

∆<br />

s<br />

→<br />

du<br />

Lauko u išvestinę pagal kryptį l galime<br />

ds<br />

M<br />

apskaičiuoti, kaip sudėtinės funkcijos išvestinę<br />

du ∂u<br />

dx ∂u<br />

dy ∂u<br />

dz<br />

= + + ,<br />

ds ∂x<br />

ds ∂y<br />

ds ∂z<br />

ds<br />

čia dešinėje pusėje galime rasti dviejų vektorių skaliarinę sandaugą, t.y.<br />

© A. Laurutis, D.Šiaučiūnas Paskaitų konspektas 95<br />

y<br />

y


du ⎛ ∂u<br />

→ ∂u<br />

→ ∂u<br />

→⎞<br />

⎛ dx →<br />

= ⎜ i + j+<br />

k ⎟ ⋅ ⎜ i +<br />

dl ⎝ ∂x<br />

∂y<br />

∂z<br />

⎠ ⎝ dS<br />

P<strong>ir</strong>masis iš šių vektorių yra vadinamas lauko u gradientu <strong>ir</strong> žymimas<br />

© A. Laurutis, D.Šiaučiūnas Paskaitų konspektas 96<br />

dy<br />

dS<br />

Jo fizikinę prasmę išsiaiškinsime vėliau. Antrasis vektorius<br />

yra krypties →<br />

l vienetinis vektorius. Todėl<br />

Taigi skaliarinio lauko u u(<br />

x,<br />

y,<br />

z)<br />

→<br />

j+<br />

dz<br />

dS<br />

→⎞<br />

k ⎟<br />

⎠<br />

∂u<br />

→ ∂u<br />

→ ∂u<br />

→<br />

grad u = i + j+<br />

k . (2)<br />

∂x<br />

∂y<br />

∂z<br />

→ → →<br />

d<br />

⎛<br />

⎞ →<br />

→ → → ⎜ x i + y j+<br />

z k ⎟<br />

dx dy dz<br />

→<br />

+ + =<br />

⎝<br />

⎠ d r<br />

i j k<br />

= = τ<br />

ds ds ds ds ds<br />

du<br />

ds<br />

→<br />

= grad u ⋅τ<br />

(4)<br />

= išvestinė →<br />

l kryptimi yra lygi dviejų vektorių grad u <strong>ir</strong><br />

→<br />

τ skaliarinei sandaugai. Pasinaudokime žinoma vektorinės algebros formule, tada<br />

du<br />

→<br />

= grad u ⋅τ<br />

= prτ<br />

( grad u)<br />

= gradlu<br />

(5)<br />

ds<br />

yra gradiento projekcija(skaliaras) į vektorių →<br />

l .<br />

∂u<br />

∂u<br />

∂u<br />

Pažymėkime, kad išvestinės , , yra gradiento projekcijos.<br />

∂x<br />

∂y<br />

∂z<br />

u = u x,<br />

y,<br />

z pilnąjį diferencialą<br />

Apskaičiuokime funkcijos ( )<br />

∂u<br />

∂u<br />

∂u<br />

⎛ ∂u<br />

→ ∂u<br />

→ ∂u<br />

→⎞<br />

→ → →<br />

du = dx + dy + dz =<br />

⎛<br />

⎞<br />

⎜ i + j+<br />

k ⎟⎜dx<br />

i + dy j+<br />

dz k ⎟ =<br />

∂x<br />

∂y<br />

∂z<br />

⎝ ∂x<br />

∂y<br />

∂z<br />

⎠⎝<br />

⎠<br />

→ → →<br />

→<br />

= grad u ⋅ d<br />

⎛<br />

⎞<br />

⎜ x i + y j+<br />

z k ⎟ = grad u ⋅ d r<br />

⎝<br />

⎠<br />

Dabar išsiaiškinkime grad u fizikinę prasmę. Iš<br />

du<br />

(5) formulės matome, kad yra maksimalus,<br />

ds<br />

kai grad u <strong>ir</strong> →<br />

τ kryptys sutampa, nes tada<br />

projekcija yra didžiausia. Taigi grad u reikšmė<br />

taške M nusako absoliutinį lauko u = u(<br />

x,<br />

y,<br />

z)<br />

kitimo greitį, o jo kryptis šio lauko kitimo<br />

(augimo) kryptį.<br />

Jei tam tikroje aplinkoje yra „krosnis“, tai visur<br />

temperatūros lauko gradientas nukreiptas į<br />

„krosnį“. Kuo arčiau „krosnies“, tuo<br />

temperatūros lauko modulis didesnis, nes ten<br />

temperatūros lauko kitimas yra didesnis.<br />

(3)


taško i<br />

P iki i+<br />

1<br />

→ →<br />

i = Fi<br />

⋅∆<br />

si<br />

§ 10. PIRMOJO IR ANTROJO TIPŲ KREIVINIAI INTEGRALAI<br />

P , t.y.<br />

→ →<br />

∆ i = Pi<br />

Pi<br />

+ 1<br />

A (skaliarinė sandauga).<br />

o<br />

y<br />

y i<br />

O<br />

i+<br />

1<br />

y<br />

∆<br />

yi<br />

P 2<br />

P 1<br />

M<br />

→<br />

Fi<br />

xi<br />

P i<br />

∆ xi<br />

Tegul turime tam tikrą liniją L , kurioje<br />

tarp taškų M <strong>ir</strong> N , jėgos →<br />

F veikiamas, juda<br />

P x,<br />

y . Jėga →<br />

F gali keisti tiek savo<br />

taškas ( )<br />

→<br />

kryptį, tiek <strong>ir</strong> dydį, t.y. F<br />

→<br />

F(<br />

P)<br />

= .<br />

Apskaičiuokime jėgos →<br />

F atliekamą<br />

darbą A , kai taškas P nueina kelią nuo taško<br />

M iki taško N .<br />

Kreivės L dalį nuo taško M iki taško<br />

N taškaisP i išdaliname į n dalių,<br />

pažymėdami →<br />

∆s i vektorių<br />

→<br />

P iPi<br />

+1 . Lankus,<br />

esančius tarp taškų P i <strong>ir</strong>P i+<br />

1 , pakeičiame<br />

styga ∆ si<br />

, tiksliau vektoriumi, einančiu nuo<br />

s . Tada elementarioje atkarpoje atliekamas darbas yra<br />

Jei ( ) ( ) →<br />

→<br />

→<br />

F = X x,<br />

y i + Y x,<br />

y j , čia X ( x,<br />

y)<br />

<strong>ir</strong> ( x y)<br />

→ → →<br />

si = ∆xi<br />

i + ∆yi<br />

P<br />

∆ j , tai<br />

x<br />

i+<br />

1<br />

i+<br />

1<br />

A ≈<br />

∆ si<br />

i ∆ si<br />

→ →<br />

Y , jėgos →<br />

F projekcijos,<br />

( xi<br />

yi<br />

) ∆xi<br />

+ Y ( xi<br />

, yi<br />

) yi<br />

F = X , ∆ <strong>ir</strong><br />

n<br />

∑<br />

i=<br />

1<br />

→<br />

→<br />

F ∆ s<br />

i<br />

i<br />

=<br />

n<br />

∑<br />

i=<br />

1<br />

[ X ( x y ) ∆x<br />

+ Y ( x , y ) ∆y<br />

]<br />

© A. Laurutis, D.Šiaučiūnas Paskaitų konspektas 97<br />

i<br />

, (1)<br />

Apskaičiuojame reiškinio (1) ribą, kai max s → 0 ( x → 0 , y → 0)<br />

<strong>ir</strong> gauname<br />

∫<br />

n<br />

∑ ∆xi→0<br />

∆y<br />

→ i=<br />

1<br />

i 0<br />

∆ i<br />

i<br />

i<br />

∆ i<br />

[ X ( x , y ) ∆x<br />

+ Y ( x , y ) ∆y<br />

]<br />

A = lim .<br />

i<br />

i<br />

i<br />

( N )<br />

i<br />

i<br />

i<br />

i<br />

i<br />

i<br />

∆ i<br />

( x y)<br />

dx + Y ( x,<br />

y)<br />

dy = X ( x,<br />

y)<br />

dx + Y ( x y)<br />

∫<br />

A = X , , dy . (2)<br />

L<br />

Ši integralinės sumos riba vadinama funkcijų X ( x,<br />

y)<br />

<strong>ir</strong> ( x y)<br />

Y , antrojo tipo <strong>kreiv</strong>iniu<br />

integralu <strong>kreiv</strong>e L .<br />

Pastaba. Integrale (2) M <strong>ir</strong> N yra skliaustuose, nes tai ne skaičiai, o linijos, kuria vyksta<br />

integravimas, pradžios <strong>ir</strong> galo taškai.<br />

Tuo atveju, kai <strong>kreiv</strong>ė L yra erdvinė <strong>kreiv</strong>ė, tada trijų funkcijų X ( x,<br />

y,<br />

z)<br />

, Y ( x,<br />

y,<br />

z)<br />

, Z(<br />

x,<br />

y,<br />

z)<br />

<strong>kreiv</strong>inis integralas apskaičiuojamas analogiškai<br />

∫<br />

( M )<br />

( x y,<br />

z)<br />

dx + Y ( x,<br />

y,<br />

z)<br />

dy + Z(<br />

x,<br />

y z)<br />

A = X , , dz , (3)<br />

L<br />

N<br />

L<br />

x


čia L - erdvinė integravimo <strong>kreiv</strong>ė.<br />

Atsk<strong>ir</strong>u atveju, kai F ( P)<br />

yra skaliarinė funkcija (pvz.: linijos tankis), tada galime gauti<br />

p<strong>ir</strong>mojo tipo <strong>kreiv</strong>inį integralą, kuris reiškia linijos masę<br />

⎧x<br />

= ϕ<br />

⎨<br />

⎩y<br />

= ψ<br />

M<br />

=<br />

n<br />

∑<br />

lim<br />

max ∆si<br />

→0<br />

i=<br />

1<br />

F ∆s<br />

( x,<br />

y)<br />

§ 11. KREIVINIŲ INTEGRALŲ APSKAIČIAVIMAS<br />

Apskaičiuokime p<strong>ir</strong>mojo tipo <strong>kreiv</strong>inį integralą įva<strong>ir</strong>ias atvejais.<br />

1. <strong>kreiv</strong>ė L duodama parametrinėmis lygtimis:<br />

() t ,<br />

() t ,<br />

tada<br />

⎧dx<br />

= ϕ<br />

⎨<br />

⎩dy<br />

= ψ<br />

′ () t<br />

′ () t<br />

© A. Laurutis, D.Šiaučiūnas Paskaitų konspektas 98<br />

i<br />

i<br />

=<br />

∫<br />

L<br />

F<br />

ds<br />

. (4)<br />

dt,<br />

2 2<br />

2<br />

2<br />

, čia t1 ≤ t ≤ t2<br />

. Kadangi ds = dx + dy = [ ϕ ′ () t ] + [ ψ ′ () t ] dt ,<br />

dt,<br />

2<br />

2<br />

2<br />

tai F(<br />

x,<br />

y)<br />

ds = F(<br />

() t , ψ () t ) ϕ ′ ( t)<br />

+ ψ ′ ( t)dt<br />

⎧x<br />

= ϕ<br />

⎪<br />

⎨y<br />

= ψ<br />

⎪<br />

⎩z<br />

= η<br />

∫<br />

L<br />

2. <strong>kreiv</strong>ė L duota kryptimi f ( x)<br />

∫<br />

L<br />

F<br />

t<br />

∫<br />

t1<br />

ϕ . (1)<br />

2<br />

y = , a ≤ x ≤ b , tada ds = 1 + f ′ ( x)dx<br />

<strong>ir</strong><br />

b<br />

2<br />

( x y)<br />

ds = F(<br />

x,<br />

f ( x)<br />

) 1+<br />

f ′ ( x)<br />

∫<br />

, dx . (2)<br />

3. <strong>kreiv</strong>ė L yra erdvinė <strong>ir</strong> duota parametrinėmis lygtimis:<br />

() t ,<br />

() t ,<br />

() t ,<br />

čia t1 t ≤ t2<br />

2 2 2<br />

≤ , tada ds () t () t ()dt t<br />

∫<br />

L<br />

F<br />

a<br />

= ϕ ′ + ψ ′ + η′<br />

<strong>ir</strong><br />

t<br />

2<br />

2<br />

2 2<br />

( x,<br />

y,<br />

z)<br />

ds = F(<br />

() t , ψ () t , η()<br />

t ) ϕ ′ () t + ψ ′ () t + η′<br />

()dt t<br />

∫<br />

t1<br />

ϕ . (3)<br />

Antrojo tipo <strong>kreiv</strong>inius integralus apskaičiuojame analogiškai. Tegul plokščia <strong>kreiv</strong>ė L yra<br />

duota parametrinėmis lygtimis:<br />

⎧x<br />

= ϕ<br />

⎨<br />

⎩y<br />

= ψ<br />

() t ,<br />

() t ,<br />

čia t1 t ≤ t2<br />

∫<br />

L<br />

≤ , tada<br />

⎧dx<br />

= ϕ<br />

⎨<br />

⎩dy<br />

= ψ<br />

′ ( t)<br />

′ () t<br />

t<br />

dt,<br />

, todėl<br />

dt,<br />

2<br />

( x,<br />

y)<br />

dx + Y ( x,<br />

y)<br />

dy = ( X ( ( t)<br />

, ψ ( t)<br />

) ϕ ′ () t + Y ( ϕ () t , ψ () t ) ′ ( t)<br />

)<br />

∫<br />

X ϕ ψ dt . (4)<br />

Pavyzdys. Apskaičiuoti integralą ∫ ydx − xdy v<strong>ir</strong>šutine elipsės dalimi<br />

> Elipsės parametrinės lygtys yra<br />

L<br />

t1


⎧x<br />

= a cost,<br />

⎨ jos v<strong>ir</strong>šutinė dalis, kai 0 ≤ t ≤ π . Apskaičiuojame ⎨<br />

⎩y<br />

= bcost,<br />

⎩ ⎧dx<br />

= −a<br />

sin tdt,<br />

tada<br />

dy = bcostdt,<br />

π<br />

∫<br />

0<br />

2 2<br />

( asin<br />

t)<br />

dt − a costbcostdt<br />

= −ab<br />

( sin t + cos t)<br />

dt = − ab<br />

bsin t − ∫<br />

π . <<br />

π<br />

0<br />

§ 12. PLOTŲ APSKAIČIAVIMAS KREIVINIU INTEGRALU<br />

Tegul plokštumoje turime sritį D , apribotą<br />

uždara linija L . Sritį D aprašome taip:<br />

⎧x<br />

∈(<br />

a,<br />

b)<br />

,<br />

⎨<br />

Tokios srities plotas yra<br />

⎩y<br />

∈(<br />

y1(<br />

x)<br />

, y2<br />

( x)<br />

) .<br />

∫∫<br />

b<br />

∫<br />

y2<br />

( x)<br />

∫<br />

( x)<br />

∫<br />

( x)<br />

dx − y ( x)<br />

Q = dxdy = dx dy = y2<br />

1<br />

D<br />

a<br />

y1<br />

Paskutiniuosius du apibrėžtinius integralus<br />

galime pakeisti <strong>kreiv</strong>iniais <strong>integralai</strong>s. P<strong>ir</strong>majame<br />

iš šių <strong>kreiv</strong>inių integralų integruojame <strong>kreiv</strong>e,<br />

priešinga <strong>kreiv</strong>ei l 2 , nes y = y2(<br />

x)<br />

, a ≤ x ≤ b yra<br />

šios <strong>kreiv</strong>ės − l2<br />

lygtis, todėl<br />

b<br />

∫<br />

a<br />

2<br />

( x)<br />

dx = ∫ ydx = −∫<br />

y ydx . Tuo tarpu, antrajame<br />

−l2 l2<br />

b<br />

a<br />

integrale integruojame <strong>kreiv</strong>e l 1 <strong>ir</strong> jis yra ∫ y 1(<br />

x)<br />

dx = ∫ ydx , todėl<br />

Q = −<br />

b<br />

∫<br />

a<br />

b<br />

a<br />

∫<br />

l2<br />

ydx −<br />

© A. Laurutis, D.Šiaučiūnas Paskaitų konspektas 99<br />

dx<br />

∫<br />

l1<br />

l1<br />

ydx = −<br />

∫<br />

L<br />

ydx . (1)<br />

Analogiškai galime įrodyti, kad Q = ∫ xdy . (2)<br />

Sudedame (1) <strong>ir</strong> (2), o po to padaliname iš 2 <strong>ir</strong> gauname<br />

y<br />

Q<br />

φ<br />

M<br />

x<br />

L<br />

O<br />

y<br />

1<br />

Q = ∫ xdy − ydx . (3)<br />

2<br />

L<br />

Pavyzdys. Apskaičiuoti De<strong>kart</strong>o lapo kilpos<br />

plotą, jei šios <strong>kreiv</strong>ės parametrinės lygtys yra<br />

⎧ 3at<br />

⎪<br />

x =<br />

1+<br />

t<br />

⎨<br />

⎪ 3at<br />

⎪<br />

y =<br />

⎩ 1+<br />

t<br />

> Apskaičiuojame<br />

3<br />

2<br />

3<br />

a b<br />

,<br />

.<br />

l 1<br />

D<br />

L<br />

y1 ( x)<br />

L<br />

y2 ( x)<br />

l 2<br />

x


3<br />

⎧ 1−<br />

2t<br />

⎪dx<br />

= 3a<br />

dt,<br />

3 2<br />

⎪ ( 1+<br />

t )<br />

⎨<br />

4<br />

⎪ 2t<br />

− t<br />

dy = 3a<br />

dt.<br />

2<br />

⎪<br />

3<br />

⎩ ( 1+<br />

t )<br />

Kadangi<br />

y<br />

π<br />

= t = tgϕ<br />

, o 0 ≤ ϕ ≤ , tai 0 ≤ t ≤ ∞ .<br />

x<br />

2<br />

Todėl pagal (3) formulę<br />

4 ( 2t<br />

− t )<br />

3 3<br />

( 1+<br />

t )<br />

2 3 ( 1−<br />

2t<br />

)<br />

3 3<br />

( 1+<br />

t )<br />

2 3<br />

2 ∞ 2<br />

( 1+<br />

t ) 9a<br />

t<br />

dt =<br />

3 3<br />

3<br />

( 1+<br />

t ) 2 0 ( 1+<br />

t )<br />

2 ∞ ⎡<br />

⎤<br />

2 ∞<br />

9a<br />

t t<br />

9a<br />

t<br />

Q = ∫ ⎢ − ⎥dt<br />

=<br />

2<br />

∫ ∫<br />

0 ⎢⎣<br />

⎥⎦<br />

2 0<br />

3a<br />

=<br />

2<br />

2<br />

N<br />

2<br />

2<br />

2<br />

3t<br />

3a<br />

⎛ 1 ⎞ 3a<br />

lim dt = lim ⎜−<br />

⎟ =<br />

N →∞ ∫ .<<br />

0<br />

3 2 ( 1 ) 2 N →∞<br />

3<br />

+ t<br />

⎝ 1 + t ⎠ 2 0<br />

N<br />

§ 12. GRINO FORMULĖ<br />

Nustatykime ryšį tarp dvilypio integralo tam tikroje srityje D <strong>ir</strong> <strong>kreiv</strong>inio integralo linija L ,<br />

ribojančia šią sritį D .<br />

Tegul turime teisingą sritį D <strong>ir</strong> tolydžias šioje<br />

y<br />

M<br />

kiekvieną iš jų.<br />

D<br />

y 2 = f ( x)<br />

b<br />

∫<br />

a<br />

dt =<br />

srityje funkcijas X ( x,<br />

y)<br />

<strong>ir</strong> ( x y)<br />

turi <strong>ir</strong> tolydžias dalines išvestines.<br />

Panagrinėkime integralą<br />

( x,<br />

y)<br />

y2<br />

( x)<br />

( x,<br />

y)<br />

y ( x)<br />

© A. Laurutis, D.Šiaučiūnas Paskaitų konspektas 100<br />

∫∫<br />

D<br />

b<br />

∫<br />

a<br />

∂X<br />

∂y<br />

X<br />

(1)<br />

( x y ( x)<br />

) dx = X ( x y)<br />

X ,<br />

, 2<br />

čia linijos MPN parametrinės lygtys yra<br />

Q<br />

P<br />

( )<br />

x<br />

f<br />

y =<br />

1<br />

a b<br />

b<br />

∫<br />

a<br />

⎧ = ,<br />

⎨<br />

⎩ = y y<br />

x x<br />

1<br />

dxdy =<br />

dx =<br />

b<br />

∫<br />

⎛<br />

⎜<br />

⎜<br />

⎝<br />

y2<br />

∫ ∫<br />

b<br />

a<br />

a<br />

y1<br />

2<br />

( x)<br />

( x)<br />

Y , , kurios, be to<br />

∂X<br />

⎞<br />

dy⎟dx<br />

=<br />

∂y<br />

⎟<br />

⎠<br />

[ X ( x,<br />

y ( x)<br />

) − X ( x,<br />

y ( x)<br />

) ] dx.<br />

Pastebėkime, kad integralo (1) dešinioji pusė<br />

susideda iš dviejų integralų. Išnagrinėkime<br />

∫<br />

MPN<br />

( ) ,<br />

2 x<br />

( x y ( x)<br />

) dx = X ( x y)<br />

∫<br />

MQN<br />

Integralus (2) <strong>ir</strong> (3) įrašome į išraišką (1) <strong>ir</strong> gauname<br />

L<br />

N<br />

x<br />

2<br />

dx , (2)<br />

x – parametras. Analogiškai<br />

X , 1 , dx<br />

(3)<br />

1


∫∫<br />

D<br />

∂X<br />

dxdy =<br />

∂y<br />

∫<br />

X<br />

MPNQM<br />

∫<br />

MPN<br />

X<br />

( x,<br />

y)<br />

dx − X ( x,<br />

y)<br />

dx = X ( x,<br />

y)<br />

dx + X ( x,<br />

y)<br />

( x,<br />

y)<br />

dx = X ( x,<br />

y)<br />

∫<br />

L pagal laikr.rod<br />

∫<br />

MQN<br />

dx.<br />

D<br />

© A. Laurutis, D.Šiaučiūnas Paskaitų konspektas 101<br />

∫<br />

MPN<br />

∂Y<br />

Galime taip pat apskaičiuoti integralą ∫∫ dxdy <strong>ir</strong> gauti<br />

∂x<br />

∫∫<br />

D<br />

∂Y<br />

dxdy = −<br />

∂x<br />

∫<br />

Y<br />

( x,<br />

y)<br />

Lpagal<br />

laikr.rod.<br />

Kai iš (4) integralo atimame (5), randame taip vadinamą Grino formulę<br />

∫∫<br />

D<br />

⎛ ∂X<br />

∂Y<br />

⎞<br />

⎜ − ⎟dxdy<br />

=<br />

⎝ ∂y<br />

∂x<br />

⎠<br />

∫<br />

dy<br />

L pagal laikr.rod.<br />

∫<br />

NQM<br />

Xdx + Ydy<br />

kuri atvejui, kai kontūras L apeinamas prieš laikrodžio rodyklę yra kitokia<br />

∫∫<br />

D<br />

⎛ ∂Y<br />

∂X<br />

⎞<br />

⎜ − ⎟dxdy<br />

=<br />

⎝ ∂x<br />

∂y<br />

⎠<br />

∫<br />

L<br />

dx =<br />

(4)<br />

(5)<br />

, 6)<br />

Xdx + Ydy . (7)<br />

Pastaba. Kai kontūras L apeinamas prieš laikrodžio rodyklę, t.y. teigiama kontūro apėjimo<br />

kryptimi, prie L nieko specialiai nežymime, bet apėjimas pagal laikrodžio rodyklę visada<br />

pažymimas.<br />

§ 13. SĄLYGA, KADA KREIVINIS INTEGRALAS NEPRIKLAUSO NUO INTEGRAVIMO<br />

KELIO<br />

Tegul turime antrojo tipo <strong>kreiv</strong>inį integralą <strong>kreiv</strong>e L , esančia plokštumoje<br />

( N )<br />

∫<br />

( M )<br />

( x y)<br />

dx + Y ( x y)<br />

X , , dy ,<br />

čia M <strong>ir</strong> N yra <strong>kreiv</strong>ės L pradžios <strong>ir</strong> galo taškai, o funkcijos X ( x,<br />

y)<br />

<strong>ir</strong> ( x y)<br />

Y , yra tolydžios <strong>ir</strong><br />

turi tolydžias dalines išvestines tam tikroje srityje D , talpinančioje <strong>ir</strong> <strong>kreiv</strong>ę L . Išsiaiškinkime,<br />

kokią sąlygą turi tenkinti <strong>kreiv</strong>inis integralas, kad jo reikšmė nepriklausytų nuo integravimo kelio,<br />

o tik nuo taškų M <strong>ir</strong> N padėties.<br />

M N<br />

P<br />

Q<br />

L<br />

L<br />

Tegul srityje D yra dvi <strong>kreiv</strong>ės,<br />

jungiančios taškus M <strong>ir</strong> N , t.y. MPN <strong>ir</strong><br />

MQN . Jei norime, kad <strong>kreiv</strong>inis<br />

nepriklausytų nuo integravimo kelio, tai<br />

turime pareikalauti, kad būtų patenkinta<br />

sąlyga<br />

Xdx + Ydy = Xdx + Ydy arba<br />

∫<br />

MPN<br />

∫<br />

MPN<br />

∫<br />

MQN<br />

Xdx + Ydy − ∫ Xdx + Ydy = 0.<br />

MQN


Pakeiskime integravimo kryptį, tada<br />

∫<br />

NQM<br />

∫<br />

MPN<br />

Xdx + Ydy = − Xdx + Ydy <strong>ir</strong><br />

∫<br />

MQN<br />

Xdx + Ydy + ∫ Xdx + Ydy = 0 arba<br />

NQM<br />

∫ L<br />

Xdx + Ydy = 0 . (1)<br />

Paskutiniojo integralo integravimo <strong>kreiv</strong>ė yra L = MPNQM . Matome, kad <strong>kreiv</strong>inio<br />

integralo reikšmė nepriklauso nuo integravimo kelio, o tik nuo taškų M <strong>ir</strong> N padėties, kai<br />

antrojo tipo <strong>kreiv</strong>inis integralas <strong>kreiv</strong>e L = MPNQM yra lygus nuliui. Ir atv<strong>ir</strong>kščias tv<strong>ir</strong>tinimas<br />

teisingas. Jei <strong>kreiv</strong>inis integralas uždaru kontūru yra lygus nuliui, tai antrojo tipo <strong>kreiv</strong>inis<br />

integralas bet kuria šio uždaro kontūro dalies <strong>kreiv</strong>e, jungiančia bet kuriuos šio kontūro taškus,<br />

priklauso tik nuo šių taškų padėties.<br />

Teorema. Jei tam tikroje srityje D funkcijos X ( x,<br />

y)<br />

<strong>ir</strong> ( x y)<br />

∂ X ( x,<br />

y)<br />

∂ Y ( x,<br />

y)<br />

Y , yra tolydžios <strong>ir</strong> turi tolydžias<br />

dalines išvestines<br />

∂y<br />

<strong>ir</strong><br />

∂x<br />

, tai <strong>kreiv</strong>inis integralas uždaru kontūru L , esančiu šioje<br />

∂ X ( x,<br />

y)<br />

∂Y<br />

( x,<br />

y)<br />

srityje D , yra lygus nuliui, t. y. ∫ Xdx + Ydy = 0 , kai = (būtina <strong>ir</strong> pakankama<br />

∂y<br />

∂x<br />

L<br />

sąlyga).<br />

> Jau žinome, kad ∫ Xdx + Ydy = ∫ Xdx + Ydy , kai Xdx + Ydy = 0 . (1)<br />

MPN<br />

MQN<br />

© A. Laurutis, D.Šiaučiūnas Paskaitų konspektas 102<br />

∫ L<br />

( x,<br />

y)<br />

∂X<br />

( x y)<br />

⎛ ∂Y<br />

, ⎞<br />

Taip pat žinome Grino formulę ∫∫ ⎜ − ⎟dxdy<br />

= ∫ Xdx + Ydy .<br />

⎝ ∂x<br />

∂y<br />

D<br />

⎠<br />

L<br />

∂Y<br />

( x,<br />

y)<br />

∂X<br />

( x,<br />

y)<br />

Iš šios formulės seka, kad sąlyga − = 0 garantuoja, kad <strong>ir</strong> + = 0<br />

∂x<br />

∂y<br />

∫ Xdx Ydy ,<br />

L<br />

∂ Y(<br />

x,<br />

y)<br />

∂X<br />

( x,<br />

y)<br />

todėl = .<<br />

∂x<br />

∂y<br />

∂ Y ( x,<br />

y)<br />

∂X<br />

( x,<br />

y)<br />

Įdomu, ką gali reikšti sąlyga = . Galime įrodyti, kad ji yra tapatinga kitai<br />

∂x<br />

∂y<br />

sąlygai du ( x,<br />

y)<br />

= X ( x,<br />

y)<br />

dx + Y ( x,<br />

y)dy<br />

, t. y. <strong>kreiv</strong>inio integralo pointegralinis reiškinys yra tam<br />

tikros funkcijos u ( x,<br />

y)<br />

pilnasis diferencialas. Iš tiesų<br />

∂u<br />

( )<br />

( x,<br />

y)<br />

∂u(<br />

x,<br />

y)<br />

du x,<br />

y = dx + dy = X ( x,<br />

y)<br />

dx + Y ( x,<br />

y)dy<br />

, jei (2)<br />

∂x<br />

∂<br />

⎧∂u(<br />

x,<br />

y)<br />

⎪<br />

= X ( x,<br />

y)<br />

,<br />

∂x<br />

⎨<br />

(3)<br />

⎪<br />

∂u(<br />

x,<br />

y)<br />

= Y ( x,<br />

y)<br />

.<br />

⎪⎩<br />

∂y<br />

2 ⎧∂<br />

u(<br />

x,<br />

y)<br />

∂X<br />

( x,<br />

y)<br />

⎪ = ,<br />

⎪ ∂x∂y<br />

∂y<br />

Iš (3) sąlygos galime gauti ⎨<br />

(4)<br />

2<br />

⎪∂<br />

u(<br />

x,<br />

y)<br />

∂Y<br />

( x,<br />

y)<br />

= .<br />

⎪⎩<br />

∂y∂x<br />

∂x


Kaip žinome, visada<br />

tik patenkinta sąlyga (2).<br />

2 ( x,<br />

y)<br />

∂ u(<br />

x y)<br />

2<br />

∂ u<br />

,<br />

=<br />

∂x∂y<br />

∂y∂x<br />

( N )<br />

, o todėl<br />

Išvada. Kreivino integralo X ( x y)<br />

dx + Y ( x,<br />

y)<br />

∫<br />

( M )<br />

( x,<br />

y)<br />

∂Y<br />

( x y)<br />

∂ X<br />

,<br />

=<br />

∂y<br />

∂x<br />

visada patenkinta, jei<br />

, dy reikšmė nepriklauso nuo integravimo<br />

kelio (<strong>kreiv</strong>iniai <strong>integralai</strong> <strong>kreiv</strong>ėmis MPN <strong>ir</strong> MQN sutampa), o tik nuo taškų M <strong>ir</strong> N padėties,<br />

kai<br />

∂ X ( x,<br />

y)<br />

∂Y<br />

( x,<br />

y)<br />

= , o tai yra tas pats, kad <strong>kreiv</strong>inio integralo pointegralinis reiškinys<br />

∂y<br />

∂x<br />

X ( x,<br />

y)<br />

dx + Y ( x,<br />

y)dy<br />

yra tam tikros funkcijos u ( x,<br />

y)<br />

pilnasis diferencialas<br />

du x,<br />

y = X x,<br />

y dx + Y x,<br />

y .<br />

( ) ( ) ( )dy<br />

§ 14. VIENPUSIAI IR DVIPUSIAI PAVIRŠIAI. PAVIRŠIAUS ORIENTACIJA<br />

Bet kuriame <strong>pav</strong><strong>ir</strong>šiuje parinkime tašką <strong>ir</strong> fiksuokime normalės kryptį šiame taške. Slinkime<br />

pas<strong>ir</strong>inktą <strong>pav</strong><strong>ir</strong>šiumi bet kuria<br />

A<br />

C kryptimi, neperžengdami<br />

B<br />

D<br />

<strong>pav</strong><strong>ir</strong>šiaus krašto. Jei pas<strong>ir</strong>odys,<br />

kad taškui slenkant šiuo <strong>pav</strong><strong>ir</strong>šiumi<br />

tolydžiai, galime vėl grįžti į pradinį<br />

tašką su ta pačia normalės kryptimi, tai tokį<br />

<strong>pav</strong><strong>ir</strong>šių vadiname dvipusiu. Vienpusio <strong>pav</strong><strong>ir</strong>šiaus<br />

atveju į pradinį tašką grįžtame su priešinga<br />

normalės kryptimi. Paprasčiausiu vienpusio<br />

<strong>pav</strong><strong>ir</strong>šiaus <strong>pav</strong>yzdžiu gali būti Mebijuso lapas. Jis<br />

gaunamas iš popieriaus juostelės, kai ją<br />

suklijuosime persukę, kaip parodyta <strong>pav</strong>yzdyje.<br />

Pradėję šį <strong>pav</strong><strong>ir</strong>šių dažyti, nudažysime jį iš<br />

abiejų šonų, neperžengdami <strong>pav</strong><strong>ir</strong>šiaus kraštų.<br />

Dvipusio <strong>pav</strong><strong>ir</strong>šiaus taške parinkta normalės<br />

A B<br />

D C<br />

kryptis nustato šio <strong>pav</strong><strong>ir</strong>šiaus pusę arba, kaip sakoma, šį <strong>pav</strong><strong>ir</strong>šių orientuoja. Pav<strong>ir</strong>šiaus orientacija<br />

parodo, kuria jo puse yra atliekamas vienas<br />

teigiamos<br />

orientacijos<br />

x<br />

z<br />

N<br />

S<br />

n r<br />

neigiamos orientacijos<br />

y<br />

ar kitas procesas (pvz.: <strong>pav</strong><strong>ir</strong>šinio integralo<br />

apskaičiavimas).<br />

Pav<strong>ir</strong>šiuje S paimkime tašką N <strong>ir</strong> apie jį<br />

apibrėžkime pakankamai mažo spindulio<br />

apskritimu. Šio taško teigiama apėjimo<br />

kryptimi laikysime tą, kai, einant apskritimo<br />

linija, jo aprėžiamas plotas lieka ka<strong>ir</strong>ėje<br />

pusėje. Išveskime šiame taške normalę.<br />

Pav<strong>ir</strong>šiaus orientaciją laikysime<br />

teigiama, jei, žiūrint iš normalės galo į<br />

<strong>pav</strong><strong>ir</strong>šiaus tašką, pastarojo apėjimo krytis bus<br />

prieš laikrodžio rodyklę.<br />

Pav<strong>ir</strong>šiaus S v<strong>ir</strong>šutinis šonas yra teigiamos<br />

orientacijos, o apatinis šonas – neigiamos.<br />

Teigiamos orientacijos atveju <strong>pav</strong><strong>ir</strong>šiaus<br />

© A. Laurutis, D.Šiaučiūnas Paskaitų konspektas 103


normalė →<br />

n su teigiama z ašies kryptimi (vektoriumi →<br />

k ) sudaro smailą kampą, tuo neigiamos<br />

orientacijos atveju – buką.<br />

∧<br />

⎛ → →⎞<br />

Jei <strong>pav</strong><strong>ir</strong>šiaus lygtis yra z = f ( x,<br />

y)<br />

, tai kampo γ =<br />

⎜<br />

n, k<br />

⎟<br />

kryptį nusakysime krypties kosinusu<br />

⎜ ⎟<br />

⎝ ⎠<br />

1<br />

cosγ<br />

= ±<br />

, čia ženklas „+“ taikomas v<strong>ir</strong>šutinei <strong>pav</strong><strong>ir</strong>šiaus pusei, o „–“<br />

2<br />

2<br />

⎛ ∂z<br />

⎞ ⎛ ∂z<br />

⎞<br />

1+<br />

⎜ ⎟ + ⎜ ⎟<br />

⎝ ∂x<br />

⎠ ⎝ ∂y<br />

⎠<br />

apatinei pusei.<br />

§ 15. PAVIRŠINIAI INTEGRALAI<br />

Stačiakampėje koordinačių sistemoje duota tam tikra sritis V . Tegul šioje srityje yra <strong>ir</strong> tam<br />

tikras <strong>pav</strong><strong>ir</strong>šiaus σ , apribotas erdvine linija λ . Pav<strong>ir</strong>šiaus σ orientacija žinoma, nes<br />

→<br />

kiekviename jo taške P ( x,<br />

y,<br />

z)<br />

yra žinoma jo normalės – vienetinio vektoriaus n(<br />

P)<br />

teigiama<br />

kryptis, o šio vektoriaus krypties kosinusai cos α, cos β,<br />

cosγ<br />

yra tolydžios funkcijos visame<br />

<strong>pav</strong><strong>ir</strong>šiuje σ .<br />

Tegul bet kuriame <strong>pav</strong><strong>ir</strong>šiaus σ taške yra žinoma funkcija<br />

→<br />

( ) ( ) ( ) →<br />

→<br />

→<br />

x y,<br />

z i + Y x,<br />

y,<br />

z j+<br />

Z x,<br />

y z<br />

F = X , , k , (1)<br />

čia X ( x,<br />

y,<br />

z)<br />

, Y ( x,<br />

y,<br />

z)<br />

, Z(<br />

x,<br />

y,<br />

z)<br />

yra šios funkcijos →<br />

F koordinatės.<br />

Pav<strong>ir</strong>šių σ išskaidykime į n bet kokių <strong>pav</strong><strong>ir</strong>šiaus dalių, kurių kiekvienos plotas yra ∆ σ i<br />

( i = 1,<br />

n)<br />

<strong>ir</strong> kiekvienoje iš jų laisvai paimkime taškus P i , kuriuose apskaičiuokime funkcijos (1)<br />

reikšmę. Sudarykime integralinę sumą<br />

čia ( P )<br />

→<br />

→<br />

F →<br />

i<br />

– vektoriaus →<br />

F reikšmė taške i<br />

F ⋅ n – skaliarinė sandauga.<br />

→ →<br />

∧<br />

⎛ r r ⎞<br />

Fi<br />

⋅ni<br />

= Fi<br />

cos⎜<br />

ni,<br />

Fi<br />

⎟<br />

⎝ ⎠<br />

ni r<br />

Pi<br />

∆σi<br />

∧<br />

⎛ ⎞<br />

⎜ni<br />

Fi<br />

⎟<br />

⎝ ⎠<br />

r r<br />

,<br />

→<br />

Fi<br />

n →<br />

∑ F<br />

i=<br />

1<br />

P , ( P )<br />

( P ) n(<br />

P )<br />

© A. Laurutis, D.Šiaučiūnas Paskaitų konspektas 104<br />

→<br />

n<br />

i<br />

i<br />

→<br />

i<br />

∆σ<br />

, (2)<br />

i<br />

– <strong>pav</strong><strong>ir</strong>šiaus normalės →<br />

n reikšmė taške P i ,<br />

∆ i<br />

Integralinės sumos (2) riba, kai max σ → 0<br />

vadinamas antrojo tipo <strong>pav</strong><strong>ir</strong>šiniu integralu<br />

<strong>pav</strong><strong>ir</strong>šiumi σ <strong>ir</strong> žymimas ∫∫ → →<br />

F n dσ<br />

, jei tik ši riba<br />

egzistuoja. Taigi<br />

σ<br />

n<br />

∑<br />

max ∆σi<br />

→0<br />

i=<br />

1<br />

→ →<br />

∫∫<br />

→ →<br />

lim F n ∆σ<br />

= F n dσ<br />

.<br />

i<br />

i<br />

i<br />

σ<br />

(3)<br />

Kiekvienas sumos (2) narys yra<br />

∧<br />

→ →<br />

⎛ → → ⎞<br />

F ∆ = ∆<br />

⎜ ⎟<br />

i ni<br />

σ i Fi<br />

σ i cos F<br />

⎜ i,<br />

ni<br />

, (4)<br />

⎟<br />

⎝ ⎠


eiškia tūrį cilindro, kurio pagrindas ∆ σ i , o aukštinė –<br />

∧<br />

⎛ → → ⎞ →<br />

→<br />

F<br />

⎜ ⎟<br />

i cos F<br />

⎜ i,<br />

ni<br />

, nes F<br />

⎟<br />

i = Fi<br />

, n i = 1.<br />

⎝ ⎠<br />

Jeigu vektorius →<br />

F yra skysčio, pratekančio pro <strong>pav</strong><strong>ir</strong>šių σ greitis, tai (4) reiškia kiekį<br />

skysčio, pratekančio →<br />

n i kryptimi pro <strong>pav</strong><strong>ir</strong>šiaus elementą ∆ σ i per laiko vienetą. Tada ∫∫ → →<br />

F n dσ<br />

apskaičiuoja bendrą skysčio kiekį, kuris prateka pro σ teigiama kryptimi per laiko vienetą.<br />

→<br />

→<br />

Kadangi = i cosα<br />

+ j cos β + k cosγ<br />

⋅ → →<br />

→<br />

→<br />

n , t.y. = { cosα<br />

, cosβ<br />

, cosγ<br />

}<br />

→<br />

→<br />

n , o<br />

{ X ( x , y , z ) , Y ( x , y , z ) , Z ( x , y z ) }<br />

F =<br />

,<br />

© A. Laurutis, D.Šiaučiūnas Paskaitų konspektas 105<br />

, tai<br />

( x,<br />

y,<br />

z)<br />

cosα<br />

+ Y ( x,<br />

y,<br />

z)<br />

cos β Z(<br />

x,<br />

y,<br />

) γ<br />

F n = X cosα + Y cos β + Z cosγ<br />

= X<br />

+ z cos<br />

∆σ<br />

xz<br />

x<br />

∫∫<br />

σ<br />

→ →<br />

∫∫<br />

σ<br />

→ →<br />

∫∫<br />

( X α + Y cos β + Z γ )<br />

F n dσ<br />

= cos cos dσ<br />

(5)<br />

∫∫<br />

σ<br />

σ<br />

, todėl<br />

Kaip matome iš brėžinio, sandaugos<br />

⎧∆σ<br />

cosα<br />

= ∆σ<br />

yz ,<br />

⎪<br />

⎨∆σ<br />

cos β = ∆σ<br />

xz ,<br />

⎪<br />

⎩∆σ<br />

cosγ<br />

= ∆σ<br />

xy<br />

( X α + Y cos β + Z γ ) dσ<br />

= ∫∫<br />

yra <strong>pav</strong><strong>ir</strong>šiaus elemento ∆ σ projekcijos į<br />

atitinkamas koordinačių plokštumas.<br />

Būtent todėl <strong>pav</strong><strong>ir</strong>šinį integralą (5) galime<br />

perrašyti į <strong>pav</strong>idalą<br />

F n dσ<br />

= cos cos Xdydz + Ydxdz + Zdxdy . (6)<br />

σ<br />

§ 16. PAVIRŠINIO INTEGRALO APSKAIČIAVIMAS<br />

Pav<strong>ir</strong>šinį integralą bet kuriuo <strong>pav</strong><strong>ir</strong>šiumi galime pakeisti dvilypiais <strong>integralai</strong>s plokštumose<br />

∫∫<br />

σ<br />

=<br />

→ →<br />

F n dσ<br />

=<br />

∫∫<br />

σ<br />

∫∫<br />

σ<br />

( X cosα<br />

+ Y cos β + Z cosγ<br />

)<br />

X cosαdσ<br />

+ Y cos βdσ<br />

+ Z cosγdσ<br />

=<br />

∫∫<br />

σ<br />

dσ<br />

=<br />

σ<br />

Xdydz + Ydxdz + Zdxdy.<br />

Savo ruožtu, norint apskaičiuoti integralą ∫∫ Xdydz + Ydxdz + Zdxdy , pakanka mokėti<br />

apskaičiuoti atitinkamus integralus ∫∫ X cos αdσ , ∫∫Y<br />

cos βdσ<br />

, ∫∫ Z cosγdσ<br />

arba, kas yra tas<br />

pats, integralus ∫∫ Xdydz,<br />

∫∫Ydxdz,<br />

∫∫ Zdxdy .<br />

σ<br />

yz<br />

z<br />

∆σ<br />

xy<br />

n r<br />

σ<br />

xz<br />

∆σ<br />

σ<br />

∆σ<br />

yz<br />

σ<br />

xy<br />

y<br />

σ<br />

σ<br />

σ


Apskaičiavimą pademonstruokime integralu ∫∫<br />

σ<br />

© A. Laurutis, D.Šiaučiūnas Paskaitų konspektas 106<br />

γdσ<br />

Z cos , čia<br />

σ yra <strong>pav</strong><strong>ir</strong>šius z f ( x,<br />

y)<br />

= .<br />

Jei D <strong>pav</strong>adinsime <strong>pav</strong><strong>ir</strong>šiaus σ projekciją į xOy plokštumą, tai nurodytą <strong>pav</strong><strong>ir</strong>šinį integralą<br />

galėsime pakeisti dvilypiu integralu<br />

Z cosγ<br />

dσ<br />

= ± Z x,<br />

y,<br />

f x,<br />

y dxdy ,<br />

∫∫<br />

σ<br />

∫∫<br />

D<br />

( ( ) )<br />

čia ženklą „+“ arba „–“ imsime priklausomai nuo cos γ ženklo. Jei cos γ > 0 , imamas „+“, jei<br />

cosγ ≤ 0 , – „–“.<br />

Pavyzdys. Apskaičiuokime ∫∫ x yzdxdy<br />

S<br />

3<br />

spindulio r sferos išoriniu <strong>pav</strong><strong>ir</strong>šiumi, patenkančiu į<br />

p<strong>ir</strong>mąjį oktantą.<br />

2 2 2<br />

> Integravimo <strong>pav</strong><strong>ir</strong>šius S yra sferos išorė z = r − x − y , x > 0,<br />

y > 0 .<br />

Pav<strong>ir</strong>šius yra teigiamai orientuotas, nes jo normalė su Oz ašimi sudaro smailą kampą. Tada<br />

I<br />

=<br />

∫∫<br />

S<br />

3<br />

x yzdxdy =<br />

∫∫<br />

p<strong>ir</strong>mąjį ketv<strong>ir</strong>tį ( > 0 , y > 0)<br />

r<br />

2 2<br />

r −x<br />

D<br />

x<br />

3<br />

y<br />

x , t.y.<br />

2 2 2<br />

r − x − y dxdy , čia D yra skritulio x ≤<br />

⎪⎧<br />

x ∈<br />

D : ⎨<br />

⎪⎩ y ∈<br />

( 0,<br />

r)<br />

,<br />

2 2 ( 0,<br />

r − x ) .<br />

2 2<br />

r −x<br />

3<br />

2 2 2<br />

3 ⎛ 1 ⎞<br />

I = ∫ x dx ∫ y r − x − y dy = ∫ x dx⎜−<br />

⎟ ∫ ⎝ 2 ⎠<br />

0<br />

1<br />

= −<br />

2<br />

r<br />

∫<br />

0<br />

x<br />

3<br />

0<br />

r<br />

0<br />

0<br />

2 2 2<br />

+ y r dalis, patenkanti į<br />

1<br />

2 2 2 2 2 2<br />

( r − x − y ) 2 d(<br />

r − x − y )=<br />

3<br />

2 2 2<br />

r<br />

3<br />

r<br />

( r − x − y ) 2 2<br />

1 3⎛<br />

2 2 ⎞ 1 3 2 2<br />

dx = − x 0 ( r x ) 2 ∫ ⎜ − − ⎟dx<br />

= − ∫ x ( r − x )<br />

2 2 2 2 2<br />

⎡r<br />

− x = u , x = r − u<br />

⎢<br />

= ⎢−<br />

xdx = udu,<br />

x = r,<br />

u = 0,<br />

⎢<br />

⎣x<br />

= 0,<br />

u = r<br />

3<br />

0<br />

2 2<br />

r −x<br />

2<br />

⎤<br />

⎥ 1<br />

⎥ = −<br />

3<br />

⎥<br />

⎦<br />

0<br />

∫<br />

r<br />

3<br />

u<br />

4<br />

0<br />

⎝<br />

7 2 5<br />

7 7 7<br />

2 2 1 ⎛ u r u ⎞ 1 ⎛ r r ⎞ 2r<br />

( r − u ) du = ⎜ − ⎟ = ⎜−<br />

+ ⎟ = . <<br />

3 ⎜<br />

⎝<br />

7<br />

⎠<br />

5<br />

3<br />

0<br />

⎟<br />

⎠<br />

0<br />

r<br />

3⎜<br />

⎝<br />

7<br />

3<br />

2<br />

dx =<br />

§ 17. PAVIRŠINIO IR KREIVINIO INTEGRALŲ RYŠYS (STOKSO FORMULĖ)<br />

Tegul erdvėje yra <strong>pav</strong><strong>ir</strong>šiusσ :<br />

z = f ( x,<br />

y)<br />

toks, kad kiekviena tiesė,<br />

lygiagreti Oz ašiai, kerta šį <strong>pav</strong><strong>ir</strong>šių<br />

viename taške. Pav<strong>ir</strong>šių σ riboja<br />

linija λ . Teigiamą normalės →<br />

n kryptį<br />

imkime tokią, kad ji su Oz ašimi<br />

sudarytų smailą kampą. Tada krypties<br />

kosinusai yra<br />

x<br />

z<br />

O<br />

λ<br />

D<br />

L<br />

n r<br />

σ<br />

5<br />

⎟<br />

⎠<br />

akis<br />

y<br />

105<br />

prieš<br />

laikrodžio<br />

rodyklę


⎧<br />

⎪<br />

⎪cosα<br />

=<br />

⎪<br />

⎪<br />

⎪<br />

⎪<br />

⎪<br />

⎨cos<br />

β =<br />

⎪<br />

⎪<br />

⎪<br />

⎪<br />

⎪cosγ<br />

=<br />

⎪<br />

⎪<br />

⎪<br />

⎩<br />

∂f<br />

−<br />

∂x<br />

,<br />

2<br />

2<br />

⎛ ∂f<br />

⎞ ⎛ ∂f<br />

⎞<br />

1+<br />

⎜ ⎟ + ⎜ ⎟<br />

⎝ ∂x<br />

⎠ ⎝ ∂y<br />

⎠<br />

∂f<br />

−<br />

∂y<br />

,<br />

2<br />

2<br />

⎛ ∂f<br />

⎞ ⎛ ∂f<br />

⎞<br />

1+<br />

⎜ ⎟ + ⎜ ⎟<br />

⎝ ∂x<br />

⎠ ⎝ ∂y<br />

⎠<br />

1<br />

.<br />

2<br />

2<br />

⎛ ∂f<br />

⎞ ⎛ ∂f<br />

⎞<br />

1+<br />

⎜ ⎟ + ⎜ ⎟<br />

⎝ ∂x<br />

⎠ ⎝ ∂y<br />

⎠<br />

(1)<br />

Tegul visas <strong>pav</strong><strong>ir</strong>šius σ yra srityje V , kurioje duota funkcija X ( x y,<br />

z)<br />

, , tolydi <strong>kart</strong>u su jos<br />

p<strong>ir</strong>mos eilės dalinėmis išvestinėmis. Išnagrinėkime <strong>kreiv</strong>inį integralą <strong>kreiv</strong>e λ<br />

Ir linijoje λ galioja lygybė z f ( x,<br />

y)<br />

∫ λ<br />

( x y z)<br />

X , , dx .<br />

= , čia x <strong>ir</strong> y linijos L koordinatės. Galime užrašyti<br />

∫<br />

λ<br />

( x,<br />

y,<br />

z)<br />

dx = X ( x,<br />

y,<br />

f ( x y)<br />

)<br />

X ,<br />

Paskutinysis integralas (2) yra <strong>kreiv</strong>inis integralas <strong>kreiv</strong>e L . Pertvarkykime jį taip, kad<br />

galėtume pritaikyti Grino formulę<br />

∫<br />

L<br />

Xdx<br />

© A. Laurutis, D.Šiaučiūnas Paskaitų konspektas 107<br />

∫<br />

⎛ ∂Y<br />

∂X<br />

⎞<br />

+ Ydy<br />

= ∫∫⎜<br />

− ⎟dxdy .<br />

⎝ ∂X<br />

∂Y<br />

⎠<br />

Mūsų atveju, X ( x,<br />

y)<br />

= X ( x,<br />

y,<br />

f ( x,<br />

y)<br />

) , ( x,<br />

y)<br />

= 0<br />

rodyklę. Gauname tokį rezultatą<br />

−<br />

∫∫<br />

D<br />

∂X<br />

( x,<br />

y,<br />

f ( x,<br />

y)<br />

)<br />

∂y<br />

D<br />

L<br />

dx<br />

(2)<br />

Y , L apėjimo kryptis yra prieš laikrodžio<br />

dxdy =<br />

∫<br />

L<br />

X<br />

( x,<br />

y,<br />

f ( x,<br />

y)<br />

)<br />

Apskaičiuojame sudėtinės funkcijos X ( x,<br />

y,<br />

z)<br />

, čia z f ( x,<br />

y)<br />

kintamąjį y<br />

( x,<br />

y,<br />

f ( x,<br />

y)<br />

) ∂X<br />

( x,<br />

y,<br />

z)<br />

∂f<br />

( x,<br />

y)<br />

∂X<br />

( x,<br />

y z)<br />

∂ X<br />

,<br />

=<br />

+<br />

∂y<br />

∂z<br />

∂y<br />

∂y<br />

Reiškinį (4) įrašome į (3) <strong>ir</strong> gauname<br />

dx . (3)<br />

= dalinę išvestinę pagal<br />

⎛ ∂X<br />

∂X<br />

∂z<br />

⎞ ∂X<br />

∂X<br />

∂z<br />

∫ X ( x,<br />

y,<br />

f ( x,<br />

y)<br />

) dx = −∫∫⎜<br />

+ ⎟dxdy<br />

= −∫∫<br />

dxdy −<br />

∂ ∫∫ dxdy . (5)<br />

⎝ ∂y<br />

∂z<br />

∂y<br />

⎠<br />

y ∂z<br />

∂y<br />

L<br />

Paskutiniajam reiškiniui pritaikome (2), <strong>ir</strong> tada<br />

D<br />

D<br />

D<br />

(4)


∂X<br />

∂X<br />

∂z<br />

∫ X ( x,<br />

y,<br />

z)<br />

dx = −∫∫<br />

dxdy −<br />

∂ ∫∫ dxdy . (6)<br />

y ∂z<br />

∂y<br />

λ<br />

D D<br />

Integralai (6) gali būti pertvarkomi į <strong>pav</strong><strong>ir</strong>šinius, nes iš anksčiau žinome, kad<br />

Taip gauname<br />

∫∫<br />

σ<br />

( x,<br />

y,<br />

z)<br />

cos dσ<br />

= F(<br />

x,<br />

y z)<br />

© A. Laurutis, D.Šiaučiūnas Paskaitų konspektas 108<br />

∫∫<br />

F γ , dxdy .<br />

⎧<br />

⎪<br />

⎨<br />

⎪<br />

⎪⎩<br />

∫∫<br />

D<br />

∫∫<br />

D<br />

D<br />

∂X<br />

∂X<br />

dxdy =<br />

∂ ∫∫ cosγdσ<br />

,<br />

y ∂y<br />

σ<br />

∂X<br />

∂z<br />

∂X<br />

∂z<br />

dxdy =<br />

∂ ∂ ∫∫ cosγdσ<br />

.<br />

z y ∂z<br />

∂y<br />

cosβ<br />

∂f<br />

∂z<br />

Dabar jau tenka pas<strong>ir</strong>emti formulėmis (1). Iš jų paimame santykį = − = − , o iš<br />

cosγ<br />

∂y<br />

∂y<br />

∂z<br />

∂X<br />

∂z<br />

∂X<br />

čia cosγ = −cosβ<br />

. Tai įvertinę, kadangi cosγdσ cos βdσ<br />

∂y<br />

∫∫ = −<br />

z ∂y<br />

∫∫<br />

, iš (7) gauname<br />

∂<br />

∂z<br />

Formules iš reiškinio (7a) įrašome į (6), todėl<br />

∫<br />

λ<br />

=<br />

X<br />

∫∫<br />

σ<br />

( x,<br />

y,<br />

z)<br />

⎧<br />

⎪<br />

⎨<br />

⎪<br />

⎪⎩<br />

∫∫<br />

D<br />

∫∫<br />

D<br />

dx = −<br />

σ<br />

σ<br />

∂X<br />

∂X<br />

dxdy =<br />

∂ ∫∫ cosγdσ<br />

,<br />

y ∂y<br />

σ<br />

∂X<br />

∂z<br />

∂X<br />

dxdy = −<br />

∂ ∂ ∫∫ cosβdσ<br />

.<br />

z y<br />

∂z<br />

∫∫<br />

σ<br />

∂X<br />

cosγdσ<br />

+<br />

∂y<br />

⎛ ∂X<br />

∂X<br />

⎞<br />

⎜−<br />

cosγ<br />

+ cosβ<br />

⎟dσ<br />

.<br />

⎝ ∂y<br />

∂z<br />

⎠<br />

σ<br />

∫∫<br />

σ<br />

σ<br />

∂X<br />

cosβdσ<br />

=<br />

∂z<br />

Būtina pabrėžti, kad kontūro λ apėjimo kryptis yra suderinta su teigiama normalės →<br />

n<br />

kryptimi, t.y., jei į <strong>pav</strong><strong>ir</strong>šiųσ žiūrėsime iš normalės →<br />

n galo, tai λ apeinama prieš laikrodžio<br />

rodyklę.<br />

Galime užrašyti <strong>ir</strong> kitas, analogiškas reiškiniui (8), formules<br />

λ<br />

∫<br />

λ<br />

∫<br />

σ<br />

(7)<br />

(7a)<br />

⎛ ∂Y<br />

∂Y<br />

⎞<br />

Y ( x,<br />

y,<br />

z)<br />

dy = ∫∫ ⎜−<br />

cosα<br />

+ cosγ<br />

⎟dσ , (8a)<br />

⎝ ∂z<br />

∂x<br />

⎠<br />

⎛ ∂Z<br />

∂Z<br />

⎞<br />

Z( x,<br />

y,<br />

z)<br />

dz = ∫∫⎜<br />

− cos β + cosα<br />

⎟dσ<br />

. (8b)<br />

⎝ ∂x<br />

∂y<br />

⎠<br />

σ<br />

Atsk<strong>ir</strong>ai sudedame formulių (8, 8a <strong>ir</strong> 8b) ka<strong>ir</strong>iąsias <strong>ir</strong> dešiniąsias puses <strong>ir</strong> gausime taip<br />

vadinamą Stokso formulę<br />

∫<br />

λ<br />

⎛ ∂Y<br />

∂X<br />

⎞ ⎛ ∂Z<br />

∂Y<br />

⎞<br />

⎛ ∂X<br />

∂Z<br />

⎞<br />

Xdx + Ydy + Zdz = ∫∫⎜<br />

− ⎟cos<br />

γdσ + ∫∫⎜<br />

− ⎟cosαdσ<br />

+ ∫∫⎜<br />

− ⎟cos βdσ<br />

(9)<br />

⎝ ∂x<br />

∂y<br />

⎠ ⎝ ∂y<br />

∂z<br />

⎠<br />

⎝ ∂z<br />

∂x<br />

⎠<br />

σ<br />

σ<br />

σ<br />

(8)


Ši formulė nustato ryšį tarp <strong>pav</strong><strong>ir</strong>šinio integralo tam tikru <strong>pav</strong><strong>ir</strong>šiumi σ <strong>ir</strong> <strong>kreiv</strong>inio integralo<br />

šio <strong>pav</strong><strong>ir</strong>šiaus ribine linija λ , kai λ apėjimo kryptis yra prieš laikrodžio rodyklę, o cosγ ≥ 0 .<br />

I išvada. Jei <strong>pav</strong><strong>ir</strong>šius σ yra dalis plokštumos, lygiagrečios xOy plokštumai, tai yra, kai<br />

π<br />

⎛ ∂Y<br />

∂X<br />

⎞<br />

∆z = 0 , α = β = , γ = 0 , gauname Grino formulę<br />

2<br />

∫ Xdx + Ydy = ∫∫⎜<br />

− ⎟dσ<br />

, o tai reiškia,<br />

⎝ ∂x<br />

∂y<br />

λ<br />

σ ⎠<br />

kad Grino formulė yra atsk<strong>ir</strong>as Stokso formulės atvejis.<br />

II išvada. Jei formulėje (9) yra įvykdomos visos trys sąlygos<br />

∂Y<br />

∂X<br />

∂Z<br />

∂Y<br />

− = 0 , − = 0,<br />

∂x<br />

∂y<br />

∂y<br />

∂z<br />

∂X<br />

∂Z<br />

− = 0,<br />

(10) tai <strong>kreiv</strong>inis integralas<br />

∂z<br />

∂x<br />

Xdx + Ydy + Zdz = 0 , (11) o tai reiškia, kad jis<br />

© A. Laurutis, D.Šiaučiūnas Paskaitų konspektas 109<br />

∫ λ<br />

nepriklauso nuo integravimo <strong>kreiv</strong>ės formos. Kaip <strong>ir</strong> plokščios <strong>kreiv</strong>ės atveju, sąlygos (10) <strong>ir</strong> (11)<br />

išpildomos, kai pointegralinis reiškinys Xdx + Ydy + Zdz yra kokios nors funkcijos pilnasis<br />

diferencialas, t. y.<br />

( u(<br />

x,<br />

y z)<br />

)<br />

Xdx + Ydy + Zdz = d , , o todėl<br />

( N )<br />

∫<br />

( M )<br />

( N )<br />

Xdx + Ydy + Zdz = ∫ du = u<br />

M <strong>ir</strong> N <strong>kreiv</strong>ės λ taškai.<br />

( M )<br />

§ 18. GAUSO – OSTROGRADSKIO FORMULĖ<br />

( N ) − u(<br />

M )<br />

Tegul trimatėje erdvėje yra teisinga sritis V , apribota uždaru <strong>pav</strong><strong>ir</strong>šiumi σ , kuri<br />

suprojektuota į xOy plokštumą yra taip pat<br />

teisinga dvimatė sritis D . Tada <strong>pav</strong><strong>ir</strong>šių σ<br />

z<br />

galime išdalinti į tris atsk<strong>ir</strong>as dalis: 1. <strong>pav</strong><strong>ir</strong>šių<br />

σ 1 : z = f1(<br />

x,<br />

y)<br />

, apribojantį sritį V iš apačios, 2.<br />

<strong>pav</strong><strong>ir</strong>šių σ 2 : z = f2<br />

( x,<br />

y)<br />

, apribojantį sritį V iš<br />

v<strong>ir</strong>šaus, 3. cilindrinį <strong>pav</strong><strong>ir</strong>šių σ 3 , apribojantį sritį<br />

V iš<br />

šonų. Cilindrinio <strong>pav</strong><strong>ir</strong>šiaus sudarančioji tiesė yra<br />

lygiagreti Oz ašiai, o normalė – jai statmena.<br />

Pas<strong>ir</strong>enkame teigiamą <strong>pav</strong><strong>ir</strong>šiaus σ 2<br />

orientaciją, tada <strong>pav</strong><strong>ir</strong>šiuje σ 2 ,<br />

⎛ ∧<br />

→ →⎞<br />

⎜ ⎟ ⎡ π ⎤<br />

⎜<br />

n , k<br />

⎟<br />

∈<br />

⎢<br />

0,<br />

⎣ 2 ⎥<br />

,<br />

⎦<br />

⎝ ⎠<br />

∧<br />

⎛ → →⎞<br />

t.y. šiame <strong>pav</strong><strong>ir</strong>šiuje cos γ = cos<br />

⎜<br />

n,<br />

k<br />

⎟<br />

> 0 , tada<br />

⎜ ⎟<br />

⎝ ⎠<br />

∧<br />

⎛ → →⎞<br />

<strong>pav</strong><strong>ir</strong>šiuje 1 : cos cos<br />

⎜<br />

,<br />

⎟<br />

π<br />

σ γ = n k < 0 , o <strong>pav</strong><strong>ir</strong>šiuje σ 0<br />

⎜ ⎟<br />

3 : cosγ<br />

= cos = . Nagrinėkime integralą<br />

2<br />

⎝ ⎠<br />

x<br />

n3 r<br />

n2 r<br />

n1 r<br />

λ<br />

D<br />

σ 2<br />

σ3<br />

σ1<br />

, čia<br />

y


∫∫∫<br />

V<br />

=<br />

∫∫<br />

D<br />

( x,<br />

y,<br />

z)<br />

∂Z<br />

∂z<br />

Z<br />

© A. Laurutis, D.Šiaučiūnas Paskaitų konspektas 110<br />

( x,<br />

y )<br />

( x,<br />

y )<br />

( x,<br />

y,<br />

f ( x,<br />

y)<br />

) dxdy − Z(<br />

x,<br />

y,<br />

f ( x,<br />

y)<br />

) dxdy.<br />

2<br />

dxdydz =<br />

⎛<br />

⎜<br />

⎝<br />

f 2<br />

⎜ ∫∫ ∫<br />

D<br />

f1<br />

∫∫<br />

D<br />

∂Z<br />

⎞<br />

dz⎟dxdy<br />

=<br />

∂z<br />

⎟<br />

⎠<br />

Jau žinome, kad Z(<br />

x y,<br />

f ( x,<br />

y)<br />

) dxdy = Z(<br />

x,<br />

y,<br />

z)<br />

∫∫<br />

D<br />

Z<br />

∫∫<br />

D<br />

, 2<br />

∫∫<br />

σ<br />

2<br />

1<br />

cosγdσ<br />

( x y,<br />

f ( x,<br />

y)<br />

) dxdy = Z(<br />

x,<br />

y,<br />

z)(<br />

− cos ) dσ<br />

= − Z(<br />

x,<br />

y,<br />

z)<br />

, 1<br />

Reiškinius (2) <strong>ir</strong> (3) įrašome į formulę (1), tada<br />

∫∫∫<br />

V<br />

( x,<br />

y,<br />

z)<br />

∂Z<br />

∂z<br />

∫∫<br />

Dar pridedame ( x,<br />

y,<br />

z)<br />

∫∫<br />

σ<br />

1<br />

σ<br />

3<br />

dxdydz =<br />

∫∫<br />

σ<br />

1<br />

Z<br />

∫∫<br />

σ<br />

1<br />

∫∫<br />

( x,<br />

y,<br />

z)<br />

cos dσ<br />

+ Z(<br />

x,<br />

y,<br />

z)<br />

γ cosγdσ<br />

∫∫<br />

γ cosγdσ<br />

.<br />

Z cosγdσ<br />

= 0 , nes šiame <strong>pav</strong><strong>ir</strong>šiujeσ 3 – cos γ = 0,<br />

tada<br />

( x,<br />

y,<br />

z)<br />

cos dσ<br />

+ Z(<br />

x,<br />

y,<br />

z)<br />

cosγdσ<br />

+ Z(<br />

x,<br />

y,<br />

z)<br />

cosγdσ<br />

= Z(<br />

x,<br />

y,<br />

z)<br />

∫∫<br />

σ<br />

∫∫<br />

Z γ cosγdσ<br />

<strong>ir</strong><br />

2<br />

∫∫∫<br />

V<br />

∂Z<br />

dxdydz =<br />

∂z<br />

Analogiškai galime gauti šias išraiškas:<br />

∂Y<br />

∫∫∫ dxdydz =<br />

∂y<br />

V<br />

σ<br />

∫∫<br />

σ<br />

∫∫<br />

σ<br />

3<br />

Z<br />

Y<br />

( x y,<br />

z)<br />

σ<br />

2<br />

σ<br />

∫∫<br />

σ<br />

1<br />

(1)<br />

(2)<br />

(3)<br />

, cosγdσ<br />

. (4)<br />

( x y,<br />

z)<br />

, cos βdσ<br />

, (4a)<br />

∂X<br />

∫∫∫ dxdydz = ∫∫ X ( x,<br />

y,<br />

z)<br />

coαβdσ<br />

. (4b)<br />

∂x<br />

V<br />

σ<br />

Sudedame jas visas (4, 4a, 4b) <strong>ir</strong> gauname Gauso – Ostrogradskio formulę<br />

∫∫∫<br />

V<br />

⎛ ∂X<br />

∂Y<br />

∂Z<br />

⎞<br />

⎜ + + ⎟dxdydz<br />

=<br />

⎝ ∂x<br />

∂y<br />

∂z<br />

⎠<br />

∫∫<br />

σ<br />

( X ( x y,<br />

z)<br />

cos + Y ( x,<br />

y,<br />

z)<br />

cos β + Z(<br />

x,<br />

y,<br />

z)<br />

cosγ<br />

)<br />

, α dσ<br />

,<br />

kuri nustato ryšį tarp trilypio integralo srityje V <strong>ir</strong> <strong>pav</strong><strong>ir</strong>šinio integralo <strong>pav</strong><strong>ir</strong>šiumi σ , ribojančiu<br />

šią sritį.<br />

Pavyzdys. Apskaičiuoti <strong>pav</strong><strong>ir</strong>šinį integralą<br />

3<br />

3<br />

3<br />

I = x dydz + y dxdz + z dxdy išoriniu sferos<br />

∫∫<br />

S<br />

2 2 2 2<br />

+ y + z R <strong>pav</strong><strong>ir</strong>šiumi.<br />

φ<br />

x =<br />

> Gauso – Ostrogradskio formulę perrašome į<br />

r<br />

<strong>pav</strong>idalą<br />

θ<br />

⎛ ∂X<br />

∂Y<br />

∂Z<br />

⎞<br />

x<br />

∫∫∫ ⎜ + + ⎟dxdydz<br />

= ∫∫ X ( x,<br />

y,<br />

z)<br />

dydz + Y ( x,<br />

y,<br />

z)<br />

dxdz + Z(<br />

x,<br />

y,<br />

z)<br />

dxdy<br />

⎝ ∂x<br />

∂y<br />

∂z<br />

V<br />

⎠<br />

σ<br />

, nes<br />

z<br />

P(x,y,z)<br />

y


⎧dσ<br />

cosα<br />

= dσ<br />

yz = dydz,<br />

⎪<br />

⎨dσ<br />

cosβ<br />

= dσ<br />

xz = dxdz,<br />

⎪<br />

⎩dσ<br />

cosγ<br />

= dσ<br />

xy = dxdy.<br />

I<br />

=<br />

=<br />

todėl<br />

3<br />

3<br />

3<br />

x dydz + y dxdz + z dxdy =<br />

2 2 2<br />

∫∫∫ ( 3x<br />

+ 3y<br />

+ 3z<br />

) dxdydz.<br />

V<br />

∫∫<br />

S<br />

Pereiname į sferinę koordinačių sistemą<br />

⎧x<br />

= r sinϕ<br />

cosθ<br />

,<br />

⎪<br />

⎨y<br />

= r sinϕ<br />

sinθ<br />

,<br />

⎪<br />

⎩z<br />

= r cosϕ,<br />

Integravimo sritis yra<br />

2π<br />

4<br />

2<br />

I = 3∫<br />

dθ∫<br />

sinϕdϕ<br />

∫ r dr = 3θ<br />

0<br />

0<br />

π<br />

0<br />

2 2 2 2<br />

2<br />

tada x + y + z = r , dxdydz = r sinϕdrdθdϕ.<br />

R<br />

0<br />

[ 0;<br />

R]<br />

,<br />

[ 0;<br />

2π<br />

] ,<br />

[ 0;<br />

π ] ,<br />

⎧r<br />

∈<br />

⎪<br />

V = ⎨θ<br />

∈ todėl<br />

⎪<br />

⎩ϕ<br />

∈<br />

π<br />

5<br />

r<br />

5<br />

12πR<br />

5<br />

π<br />

( − cosϕ<br />

) = 3⋅<br />

2π<br />

( 1+<br />

1)<br />

= . <<br />

0<br />

R<br />

0<br />

© A. Laurutis, D.Šiaučiūnas Paskaitų konspektas 111<br />

R<br />

5<br />

5<br />

5

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!