kart_kreiv_ir_pav_integralai.pdf
kart_kreiv_ir_pav_integralai.pdf
kart_kreiv_ir_pav_integralai.pdf
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