Lab Fiz 2
Lab Fiz 2
Laboratorija fizike 2
dr Djordje Spasojević
Beograd, 2014.
iii
Djordje Spasojević
LABORATORIJA FIZIKE 2
Predgovor
Udžbenik Laboratorija fizike 2 sadrži opis i uputstva za izvodjenje studentskih
eksperimenata u okviru predmeta Laboratorija fizike 2 koji prati nastavu iz
oblasti termodinamike i molekularne fizike. Namenjen je studentima I godine
svih smerova fizike na Univerzitetu u Beogradu - Fizičkom fakultetu, kao i stu-
dentima drugih studijskih grupa koji rade iste eksperimente.
Eksperimenti se dominantno izvode na opremi doniranoj od strane Evropske
unije u okviru IPA projekta ”OPREMANjE STUDENTSKIH LABORATORIJA
NA FIZIČKOM FAKULTETU UNIVERZITETA U BEOGRADU”, realizovanog
u jesen 2013. godine.
Tekst posvećen svakom eksperimentu započinje teorijskim uvodom i nastavlja
se opisom aparature i uputstvima za njeno korišćenje. Sledi opis zadatka (cilja)
eksperimenta i upustvo za njegovu realizaciju. Od studenata se očekuje da u
obradi eksperimentalnih rezultata i njihovoj prezentaciji primene znanja stečena
na predmetu Obrada rezultata merenja. Stoga su detalji procene eksperimental-
nih grešaka često izostavljeni. Uputstva za korišćenje pratećeg softvera su data
u posebnom dokumentu.
Pri koncipiranju kursa nastojalo se da budu pokrivene najznačajnije teme iz
termodinamike i molekularne fizike. Svaki eksperiment je postavljen tako da ga je
moguće uraditi u toku jednog termina. Opcioni zadaci su predvidjeni za posebno
zainteresovane studente i rade se u dodatnim terminima. Poglavlja označena
zvezdicom sadrže teorijske detalje koje studenti I godine ne moraju da znaju, a
dati su radi kompletnosti teksta; pri pripremi eksperimenta predlaže se da se ova
poglavlja izostave.
Teme kursa, tip aparatura i eksperimentalne zadatke je izabrao predmetni
nastavnik. U finalizaciji kursa značajan doprinos su dali studenti saradnici Jelena
Pajović, Nora Trklja, Svetislav Mijatović i Petar Bokan, kao i stručni saradnici
Miloš Marinković i Milorad Mijić. Posebno veliki doprinos je dao Edib Dobardžić,
kao i Zoran Popović. Najveći doprinos prilagodjavanju i modifikaciji aparatura
je dao Duško Grujić, na čemu mu se posebno zahvaljujem.
Djordje Spasojević
predmetni nastavnik
Sadržaj
1 Idealni gasovi 1
1.1 Uvod 1
1.2 Opis aparature 3
1.3 Zadatak eksperimenta 4
2 Realni gasovi 7
2.1 Uvod 7
2.2 Opis aparature 13
2.3 Zadatak - snimanje izotermi gasa SF6 14
5 Termalno širenje 41
5.1 Uvod 41
5.2 Opis aparature 43
5.3 Eksperiment 44
6 Površinski napon 52
vi Sadržaj
6.1 Uvod 52
Zavisnost koeficijenta površinskog napona od temperature 54
Koeficijent površinskog napona binarnih rastvora 54
Merenje koeficijenta površinskog napona metodom otkidanja
prestena (du Noüy metod) 56
6.2 Opis aparature 59
Postavka aparature (torzione vage) pre merenja 59
Procedura merenja koeficijenta površinskog napona 60
6.3 Eksperiment 61
Zadatak 1 - odredjivanje koeficijenta površinskog napona za
smešu vode i etil alkohola na sobnoj temperaturi 61
Zadatak 2 - odredjivanje koeficijenta površinskog napona za vodu
u funkciji temperature 63
7 Viskoznost 66
7.1 Uvod 66
Zavisnost koeficijenta viskoznosti od temperature 67
Zavisnost koeficijenta viskoznosti binarnih rastvora od njihovog
sastava 67
Viskoznost nenjutnovskih fluida 68
Metodi merenja koeficijenta viskoznosti 68
7.2 Opis aparature 71
7.3 Eksperiment 72
Priprema merenja: 72
Zadatak 1 - odredjivanje koeficijenta viskoznosti glicerina u
funkciji temperature 75
Zadatak 2 - odredjivanje koeficijenta viskoznosti rastvora šećera
u vodi u funkciji koncentracije šećera 76
Obrada podataka 76
8 Termodinamički ciklusi 77
8.1 Uvod 77
8.2 Stirling mašina 80
8.3 Opis aparature 81
8.4 Odredjivanje P V dijagrama i efikasnosti Stirling mašine 91
Tok eksperimenta 91
10 Zvuk 122
10.1 Uvod 122
10.2 Rezonancija vazdušnog stuba u cevi 130
Opis aparature 130
Proučavanje rezonancije vazdušnog stuba u cevi - tok eksperimenta 131
10.3 Monokord - slobodne i prinudne oscilacije 134
Opis aparature 134
Povezivanje aparature 134
Slobodne oscilacije monokorda - tok eksperimenta 135
Prinudne oscilacije monokorda (opcioni zadatak) 142
10.4 Dodaci 145
Slobodne oscilacije žice sa prigušenjem 146
Prinudne oscilacije žice sa prigušenjem 151
1 Provera zakona idealnog gasnog
stanja
1.1 Uvod
pV = const , (1.2)
2 Francuski naučnik Šarl (Jacques Alexandre César Charles, 1746-1823) je otkrio ovaj zakon
oko 1780. godine ali svoje otkriće nije objavio. Zakon prvi objavljuje Gej-Lisak 1802.
godine, zbog čega se zakon često (npr. u našoj zemlji) naziva Gej-Lisakov zakon.
3 Francuski fizičar Gej-Lisak zakon(Joseph Louis Gay-Lussac, 1778–1850) je otkrio ovaj
zakon 1809. godine, ne znajući da je u periodu 1700—1702 isti zakon već otkrio francuski
fizičar Amonton (Guillaume Amontons, 1663-1705).
1.2 Opis aparature 3
Vazduh se u špric upušta tako što se slavina (7) otvori - okrene da bude
paralelna cevčici slavine. Upušteni vazduh ima atmosferski pritisak i sobnu tem-
peraturu. Kada je slavina normalna na cevčicu gas je zatvoren u špricu i njegov
pritisak zavisi od temperature i zapremine. Pritisak vazduha u špricu se meri
meračem pritiska (4) koji je povezan na laptop (9). Izmereni pritisak prikazuje
PHYWE Measure program; čim se program pokrene prepoznaće merače koji su
uključeni preko USB-a i početi da prikazuje izmerene vrednosti.
Zapremina šprica se fino podešava i kontroliše pomoću zavrtnja (3). Za pre-
cizno očitavanje (do na 1/4 podeoka) skale na špricu se koristiti lupa (10).
4 Idealni gasovi
Obrada podataka
6 U daljem radu nemojte koristiti vrednosti parametara #A, #B i #C, dobijene za ovu, već
za prethodnu modelnu funkciju.
2 Proučavanje realnih gasova
2.1 Uvod
P V = nRT ,
koliki deo sistema se nalazi u tečnoj, odnosno gasovitoj fazi, a odgovarajući deo
izoterme je paralelan V -osi. Najveću zapreminu VG (T ) na ovom delu izoterme
sistem ima u gasovitoj, a najmanju zapreminu VL u tečnoj fazi. Za sve zaprem-
ine VL < Vx < VG sistem nije homogen - deo sistema se nalazi u tečnoj, a
ostatak u gasovitoj fazi i ove dve faze koegzistiraju. Budući da joj je gustina
veća, tečna faza se (zbog gravitacije) grupiše u donjem delu sistemu iznad kojeg
je gas razdvojen od tečnosti meniskom, slika 2.3.
Ako sistemu pri koegzistenciji tečne i gasovite faze dovodimo energiju, dove-
dena energija se troši na pretvaranje dela tečne faze u gas, ali se pritisak pri
fiksiranoj temperaturi ne menja sve dok se ne istroši sva tečna faza. Obratno,
ako odvodimo energiju, deo gasa se pretvara u tečnost (i pritisak se ne menja sve
dok se ne istroši sva gasna faza). Ako pak sistemu smanjimo zapreminu, sistem
neće povećati pritisak, već će porasti udeo tečne faze a višak energije će biti pre-
dat okolini. Tek kada ceo sistem predje u tečno stanje dalje smanjenje zapremine
rezultira povećanjem pritiska i to veoma naglim jer je stišljivost tečnosti mala.
U toku faznog prelaza gustine tečne i gasovite faze, ρL i ρG , odnosno molarne
zapremine2 faza, VmL i VmG , ostaju konstantne; ono što se menja je udeo sistema
koji se nalazi u tečnoj, odnosno gasovitoj fazi. Pošto udeo sistema u tečnoj,
odnosno gasovitoj, fazi iznosi mL /m, odnosno mG /m, gde su mL i mG mase
tečne i gasovite faze pri konstantnoj ukupnoj zapremini m = mL + mG , to je
zapremina tečne faze VL = VmL mL /m, a gasovite VG = VmG mG /m.
Razlika molarnih zapremina VmL i VmG opada sa porastom temperature i na
kritičnoj temperaturi Tc postaje jednaka nuli. To znači da su na Tc molarne za-
2 Molarna zapremina je količnik zapremine i broja molova; brojno je jednaka zapremini
jednog mola.
2.1 Uvod 9
Slika 2.2 Levo: Izoterme realnog gasa - shematski prikaz. Desno: izoterme izopentana.
Slika 2.3 Koegzistencija gasne i tečne faze u kapilari aparature za proučavanje realnih
gasova - vidi sliku 2.6. Faze su razdvojene meniskom.
Slika 2.4 Johannes Diderik van der Waals, holandski fizičar (1837 - 1923).
a2 8a
Pc = , Tc = , Vmc = 3b , (2.3)
27b2 27Rb
ili obratno, van der Waals konstante:
Ispitivani realni gas (npr SF6 ) se nalazi u staklenoj kapilari (1) i odozdo je
zatvoren živom - slika 2.6. Pritisak gasa u kapilari se reguliše okretanjem točka
(5) a očitava na meraču pritiska (4). Zavrtanjem točka (tj okretanjem točka u
smeru kazaljke na satu) se povećava pritisak, a nivo žive u kapilari se podiže.
Dozvoljeni su pritisci do 50 bar (tj 50 × 105 Pa). Pritisak se ne sme povećati
iznad 50 bar jer može doći do prskanja kapilare.
Kapilara sa gasom se nalazi unutar šire cevi za termostatiranje (2). Ova cev je
crevima povezana sa cirkulacionim termostatom (7) i kroz nju teče providna
tečnost za termostatiranje (voda ili smeša vode i alkohola ako su potrebne
temperature ispod 0◦ C). Tečnost za termostatiranje (hladjenjem ili grejanjem)
obezbedjuje zadatu temperaturu u kapilari sa gasom.
14 Realni gasovi
6 Uredjaj može da zapamti tri zadate temperature. U tu svrhu treba pretisnuti jedno od
dugmadi P1, P2 ili P3.
2.3 Zadatak - snimanje izotermi gasa SF6 15
7 Najmanji podelak na skali merača pritiska iznosi 0,5 bar. Po uobičajenim pravilima
očitavanja podataka sa analogne skale očitavanje bi trebalo vršiti zaokruživanjem na
najbliži ceo podelak. Ovde se medjutim može povećati čitljivost skale (tj rezolucija) jer su
fluktuacije pritiska zanemarive, a mi smo prvenstveno zainteresovani za promene pritiska.
16 Realni gasovi
Obrada podataka:
1. Izmerene izoterme nacrtati na istom grafiku u programu Origin ili Table
Curve 2D.
2. Odrediti kritični tačku na izotermi 45,6◦ C - to je tačka u kojoj je tangenta
horizontalna. Koordinate ove tačke su kritični pritisak Pc i kritična zapremina
Vc .
3. Iz vrednosti za kritični pritisak i kritičnu temperaturu izračunati van der
Waals konstante a i b za SF6 .
4. Na posebnom grafiku nacrtajte iznadkritičnu izotermu (koju ste dobili meren-
jem na najvišoj temperaturi) crtajući pritisak u funkciji 1/V . Tangenta na
početni deo ove krive bi trebalo da odgovara veličini nRT , odakle treba da
odredite broj molova gasa n u kapilari.
5. Izračunajte kritičnu molarnu zapreminu Vmc = Vc /n i utvrdite da li se ona
slaže sa teorijskim predvidjanjem Vmc = 3b i vrednošću van der Waals kon-
stante b.
6. U oblasti faznog prelaza na izabranoj podkritičnoj izotermi nacrtajte za-
preminu sistema V u funkciji zapremine gasa Vg . Eksperimentalne tačke bi
trebale da leže na pravoj
( )
VmL
V = nVmL + 1 − Vg ,
VmG
iz čijeg odsečka nVmL i koeficijenta pravca 1 − VVmG
mL
treba da odredite molarnu
zapreminu tečne faze VmL i molarnu zapreminu gasne faze VmG na tempera-
turi izoterme koju ste koristili.
3 Provera Njutnovog zakona
hladjenja
3.1 Uvod
Telo (sistem) čija je temperatura viša od temperature okoline se hladi sve dok
mu se temperatura ne spusti na temperaturu okoline. Pri hladjenju telo predaje
okolini toplotu. Razmena toplote sa okolinom se vrši na sledeća tri načina:
- zračenjem,
- provodjenjem,
- konvekcijom.
Zračenje
Provodjenje
Konvekcija
su ova dva procesa znatno manje efikasna od prenosa toplote mešanjem.9 Drugim
rečima, konvekcija je najvažniji način hladjenja tela potopljenog u fluid sve dok
je prenos energije zračenjem mali.
odnosno
∫t
dTa kt′ ′
A = A0 − e dt ,
dt′
0
te je
∫t
dTa kt′ ′ −kt
Trel = A0 − e dt e .
dt′
0
U početnom trenutku t = 0 je
(Trel )0 = A0
pa je
∫t
dTa kt′ ′ −kt
Trel = (Trel )0 − e dt e .
dt′
0
Kako je
∫t ∫t ∫t
dTa kt′ ′ ′ ′
e dt = ekt dTa = Ta (t)ekt − Ta0 − k Ta (t′ )ekt dt′ ,
dt′
0 0 0
to je
∫t
kt′
T (t) − Ta (t) = T0 − Ta0 − Ta (t)e kt
+ Ta0 + k ′
Ta (t )e ′ −kt
dt e ,
0
te je konačno
∫t
′
T (t) = T0 + k Ta (t′ )ekt dt′ e−kt . (3.8)
0
Kada je telo hladnije od okolnog fluida ono će se grejati. Zakon grejanja ima isti
oblik kao i zakon hladjenja:
dT
= −kg (T − Ta ) , (3.9)
dt
gde je kg konstanta grejanja. Načini grejanja su isti kao i načini hladjenja (radi-
jacija, provodjenje, konvekcija). Napomenimo da sloj vazduha uz telo postaje
hladniji od ostalog vazduha, te budući da mu je gustina veća struji naniže.
Drugim rečima, konvektivno strujanje vazduha je drugačije nego kod grejanja.
24 Njutnov zakon hladjenja
13 Ćelija tabele se formatira tako što se na nju postavi kursor, desnim klikom se otvori
padajući meni, pa se iz njega odabere Format Cells... Nakon toga se otvara prozor za
formatiranje. Da bi se formatiralo više ćelija na isti način ove ćelije se prvo selektuju
(mišem). Da bi se formatirala cela kolona (ili vrsta) potrebno je selektovati zaglavlje
kolone (ili vrste) pa zatim izvršiti formatiranje.
14 Excel funkcija TEXT daje sadržaj ćelije kao string. Funkcija TIMEVALUE pretvara string
u kojem je upisano vreme u broj izmedju 0 i 1, koji kada se pomnoži sa 86400 (broj
sekundi u jednom danu) daje vreme pretvoreno u sekunde.
26 Njutnov zakon hladjenja
Napomena: ovaj zadatak rade grupe sa manje od tri studenta. Ako grupa broji
više od dva studenta umesto ovog zadatka se radi Zadatak 2.
Važno: aparatura ne sme biti izložena sunčevoj svetlosti niti u blizini grejnih
tela; u prostoriji ne sme biti promaje da ne bi došlo do prinudne konvekcije.
Zatim:
T0f = ea + Ta0 ,
a njena greška
∆T0f = ea ∆a + ∆Ta0 .
Napomena: ovaj zadatak rade grupe sa tri ili više studenta. Ako grupa broji
manje od tri studenta umesto ovog zadatka se ponavlja Zadatak 1, ali za drugi
balon.
Zadatak 2 se radi po istoj proceduri kao i Zadatak 1, samo što jedan student
meri vremena na prvom, drugi na drugom balonu, a treći zapisuje rezultate
19 Logaritamska funkcija ima za svoj argument bezdimenzionalne veličine, a veličina T − Ta0
ima dimenziju temperature. Zapis ln(T − Ta0 ) je skraćeni zapis za puni zapis
ln[(T − Ta0 )/1◦ C]
3.6 Zadatak 3 - provera zakona grejanja 29
22 Malo eksperimentisati sa vrednošću za Lag, tj broj susednih tačaka koje program koristi za
izračunavanje izvoda. Manje vrednosti za Lag daju ”šumovitije”, a veće - glatkije podatke;
prvih Lag vrednosti se gube (npr za Lag=4, prvi izvod je sračunat počevši od 5-te tačke)
4 Odredjivanje specifične i latentne
toplote isparavanja vode i latentne
toplote topljenja leda
4.1 Uvod
3 Može se pokazati da je
T α2p
cp = cv +
ρKT
1. specifična toplota c,
2. specifična latentna toplota isparavanja qi i
3. specifična latentna toplota topljenja qt .
P = UI , (4.4)
∆Q = P ∆t . (4.5)
Slika 4.2 Levo: multimer povezan da meri jačinu struje (tj kao ampermetar).
zeleno-žuta žica desno spaja zajednički (COM) priključak sa jednim izlazom kutije
grejača. Crvena žica levo spaja poklopac termos boce sa priključkom za merenje struje
u opsegu do 20A. Na istu tačku je paralelno priključena plava žica čiji je drugi kraj
priključen na COM ulaz multimera desno. Desno: multimer povezan da meri napon
(tj kao volmetar). Crvena žica desno je priključena da VΩHz priključak multimera, a
njen drugi kraj povezan na izlaz grejača na koji nije povezana žica do ampermetra.
4.4 Opis eksperimenta 35
Postupak merenja
Slika 4.3 Početak lista ’Cp data’ Excel tabele za unos podataka. Praznu tabelu
dobijate od nastavnika.
• izmeriti temperaturu okoline i uneti izmerenu vrednost u ’Cp calc’ listu Excel
tabele;
• digitalnom vagom (8) izmeriti masu mb prazne termos boce bez poklopca i
uneti izmerenu vrednost u ’Cp calc’ listu Excel tabele;
• povezati aparturu uz striktan nadzor nastavnika ili stručnog saradnika; pri-
likom povezivanja grejač mora biti isključen.
1. Nakon što ste završili sa isparavanjem sačekajte par minuta pa otvorite termos
bocu i u nju dolijte vode do nivoa od 7-8 cm ispod vrha, pa zatvorite bocu.
2. Sačekajte 2-3 minuta da se temperatura u boci svuda izjednači. Izmerenu
temperaturu tople vode Tt upišite u ćeliju B1 lista ’Topljenje’ Excel tabele.
3. Zatim otvorite bocu i iz nje prospite oko 10cm vode i što brže izmerite masu
boce sa preostalom toplom vodom (bez poklopca) mbt ; upišite ovu vrednost u
ćeliju B2. Excel sračunava masu tople vode u termos boci mt = mbt − mb , gde
je mb ranije izmerena masa prazne boce bez poklopca. Ova vrednost Excel
prikazuje u ćeliji B3.
4. Zatim brzo dodajte kockice leda (bez vode) iz posude sa ledom do nivoa od
7-8 cm ispod vrha, pa zatvorite bocu.
5. Intenzivno mešajte vodu. Izmerena temperatura će da pada. Kada se
temperatura ustali znači da je postignuta toplotna ravnoteža. Upišite tem-
peraturu ravnoteže Tr u ćeliju B4.
6. Izmerite masu termos boce sa vodom i istopljenim ledom (bez poklopca) mbtl
i upišite je u ćeliju B5; Excel sračunava masu leda kojeg ste dodali ml =
mbtl − mbt .
Obrada podataka
1) Specifična toplota vode: Specifična toplota vode je približno konstantna
u intervalu temperatura od 0◦ C do 100◦ C. Ako pretpostavimo da tokom grejanja
termos boca nije predavala toplotu okolini i da nije bilo isparavanja vode, onda
je količina toplote ∆Q = P t koju je odao grejač tokom vremena t utrošena samo
na zagrevanje vode i termos boce te je
P t = (mv cv + Ck )(T − T0 ) ,
∆Q = ml (qt + cv θr ) .
odakle je
(mt cv + Ck ) · (θt − θr )
qt = − cv θr . (4.12)
ml
Izračunatu specifičnu latentnu toplotu topljenja leda koristeći izmerenu vrednost
za specifičnu toplotu vode cv i proceniti grešku za qt .
gde je dQh količina toplote koju sistem oda za vreme dt usled hladjenja, i gde je
dmv masa vode koja ispari za vreme dt (nju prvo treba zagrejati od tekuće tem-
perature T do temperature ključanja T100 , pa onda ispariti, a za ovo je potebno
utrošiti latentnu toplotu isparavanja qi i izvršiti rad po jedinici mase pri ispara-
vanju RT100 /M ).
Masa vode koja ispari tokom eksperimenta je mala. Jednostavnosti radi ćemo
smatrati da je brzina ispravanja m′ ≡ dmv /dt konstantna. Kako je član cv (T100 −
T ) znatno manji od qi smatraćemo aproksimativno da je
cv (T100 − T ) ≈ 50◦ C × cv ,
tj kao da se sva voda koja isparava prvo greje od ”prosečne” temperature 50◦ C.
Tako dobijamo
gde je
q ∗ ≡ qi + 50◦ C × cv + RT100 /M .
8 U ovom delu tekstu koristimo oznaku θ za temperaturu u Celzijusovim stepenima.
40 Karakteristične toplote za vodu
5.1 Uvod
Čvrsta tela se šire kada njihova temperatura raste, odnosno skupljaju kada im
temperatura opada. Empirijski je utvrdjeno da je za svaku linearnu dimenziju
tela l relativna promena ∆l/l0 u dovoljno uskoj oblasti temperatura propor-
cionalna promeni temperature ∆T :
∆l
= α∆T . (5.1)
l0
Ovde je ∆l = l − l0 apsolutna promena posmatrane linearne dimenzije tela (od
vrednosti l0 na temperaturi T0 do vrednosti l na temperaturi T ), ∆T = T − T0
je odgovarajuća promena temperature, dok je koeficijent proporcionalnosti α
linearni koeficijent termalnog širenja.
Za infinitezimalno male promene temperature zakon termalnog širenja (5.1)
poprima oblik
dl d(ln l)
α= = , (5.2)
ldT dT
što znači da se može smatrati da je linearni koeficijent termalnog širenja jednak
izvodu logaritma linearne dimenzije po temperaturi. Ako ovaj izvod ne zavisi od
temperature, iz (5.2) sledi da je
l = l0 [1 + α(T − T0 )] . (5.4)
• pozitivan, α > 0,
• u znatnom opsegu temperatura ne zavisi od temperature,
42 Termalno širenje
zavisi od temperature T . To se javlja npr kod vode, kod koje je γ čak negativno
u intervalu od 0◦ C do 4◦ C, što znači da se voda u ovom intervalu temperatura
skuplja pri porastu temperature.5 Ovakvo anomalno ponašanje vode ima izuzetne
posledice na zbivanja i život na površini Zemlje.6
5.3 Eksperiment
širenje, četvrti list (Voda) za termalno širenje vode u opsegu od sobne temper-
ature do vrha kapilare i peti list (Voda4C) za registrovanje termalne anomalije
vode na oko 4◦ C.
Slika 5.2 Excel tabela za unos izmerenih vrednosti. U listi Voda4C na slici je dat
korak od 1◦ C da bi se videla promena gustine vode, dok je u TE Template.xls izabran
korak od 0, 1◦ C.
=($C$2-Cn)/$A$2
=($C$2-C2)/$A$2
Faza 1:
Jedan deo studenata koji rade eksperiment snima kako se menja zapremina vode
u balonu pri temperaturama oko 4◦ C, dok drugi deo za to vreme snima termalno
izduženje šipke po ubrzanoj proceduri na dužini od 200mm. Zatim studenti men-
jaju uloge: oni koji su snimali termalno širenje šipke prelaze da snimaju vodu, a
oni koji su snimali vodu prelaze na snimanje termalnog širenja šipke na dužini
od 400mm.
1. Iz suda koji sadrži mešavinu vode i leda napunite balon do vrha vodom (bez
leda) na temperaturi ispod 2◦ C i ubacite u balon (beli) magnetić magnetne
mešalice.
2. Ako u balonu ima mehurića vazduha mešajte vodu tako da isterate mehuriće
iz balona. To radite tako što pomerate magnetni hvatač (beli štapić dužine
oko 30cm) oko balona; hvatač povlači za sobom magnetić u balonu, te se tako
voda u balonu meša.
Slika 5.3 a) Balon za vodu: (1) - kapilara; (2) sonda termometra; (3) - čep.
b) Kapilara: (1) plava nit na kapilari izgleda kao da je prekinuta na nivou
vode što omogućava da se nivo vode precizno proceni.
3. Polako začepite balon čepom (3) zajedno sa uvučenom kapilarom (4) pazeći
da se oko donjeg kraja čepa ne formiraju mehurići vazduha. Usled zatvaranja
balona nešto vode će preći u kapilaru a deo vode isteći iz balona; zato obrišite
balon kada ga zatvorite.
Važno: nivo vode u kapilari kada je balon zatvoren mora biti iznad oznake
8ml na kapilari (tj na oko 0,3ml iznad čepa).
Ovime je izvršena priprema za početak merenja. Balon postaviti na radnom
stolu i stalno mešati vodu u njemu. Pošto je temperatura okoline veća, voda
u balonu će se polako grejati.
5.3 Eksperiment 47
4. Pratiti temperaturu vode u balonu. Svaki put kada se ova temperatura povisi
za 0, 1◦ C očitati zapreminu vode u kapilari. Očitanu temperaturu i zapreminu
vode u kapilari upisati u odgovarajuće ćelije liste Voda4C Excel tabele.
5. Merenje je gotovo kada temperatura vode u balonu predje 7◦ C.9
1. Postavite stopicu, slika 5.4 a), tako da legne u žljeb šipke na poziciji 200mm.
Zašrafite stopicu za postolje aparature. Izmedju gornjeg šrafa na stopici i šipke
postavite gumu, pa pritegnite taj šraf.
2. Napunite kadicu cirkulacionog termostata destilovanom vodom sobne temper-
ature 2cm - 3cm iznad nivoa spirale grejača; uverite se da je šipka povezana
crevima sa termostatom.
3. Postavite temperaturu na termostatu na 80◦ C i uključite ga.10 U nastavku ove
faze eksperimenta termostat bi trebalo da greje vodu konstantnom snagom,
što znači da bi temperatura vode u termostatu, a time i temperatura šipke,
trebalo da raste konstantnom brzinom.
Napomena: prethodne tri tačke eksperimenta radite pod kontrolom nastavnika
i/ili stručnog saradnika.
Slika 5.4 a) Stopica za fiksiranje šipke na izabranom žljebu (ovde na 600mm): (1) -
šraf kojim se stopica šrafi za postolje aparature; donji ispust stopice treba da legne u
žljeb šipke (2) i prišrafi gornjim šrafom (3); izmedju gornjeg šrafa i šipke postaviti
gumu (4) da šraf ne ošteti šipku.
b) Komparator kojim se meri promena dužine šipke. Na velikoj skali se
očitavaju stoti delovi milimetra, a na maloj celi milimetri.
5. Kada na komparatoru, slika 5.4 b), primetite da se dužina izabranog dela šipke
povećala za izabrani korak (koji ovde iznosi ∆l = 0, 01mm) očitajte brzo tem-
peraturu na displeju termostata i zapišite očitanu vrednost u odgovarajuću
ćeliju Excel tabele.
6. Merenje je gotovo kada temperatura na termostatu predje 70◦ C.
Važno: ako primetite da se bilo šta čudno dešava tokom eksperimenta, npr da curi
voda, odmah isključite termostat i pozovite stručnog saradnika ili nastavnika.
Zatim snimite termalno širenje šipke na dužini od 400mm uz malo modi-
fikovanu proceduru:
Jedan deo studenata koji rade eksperiment snima kako se menja zapremina vode
u balonu, dok drugi deo za to vreme snima termalno izduženje šipke na dužini
od 600mm.
Obrada podataka
Slika 5.5 Prozor ”Linear Fit” programskog paketa Origin; služi za setovanje opcija
linearnog fitovanja.
=999.97487*Voda!$A$2/(Voda!$A$2+MAX($B$2:$B$200)-Bn)
12 Pogledajte i ostale opcije fitovanja. Uočite npr da su ispod ”Input Data” izlistani svi setovi
podataka koji će se fitovati. Na korisniku je da izabere koje setove želi da fituje.
5.3 Eksperiment 51
13 Numerički izvod nalazite tako što na listi glavnih menija kliknete na ”Analysis”, pa zatim
na ”Mathematics”, ”Differentiate” i ”Open Dialog” u padajućim menijima koji se redom
otvaraju. Otvoriće se prozor ”Mathematics: differentiate” u kojem treba da ozčite da želite
prvi izvod, tj Derivative order 1. Podaci za numerički izvod će biti upisani u kolonu
”Derivative” tabele u kojoj su i originalni podaci čiji izvod tražite.
6 Odredjivanje koeficijenta
površinskog napona
6.1 Uvod
1 Površinski napon se javlja ne samo kod tečnosti već i kod čvrstih tela.
6.1 Uvod 53
je tečnost u kontaktu. Za datu tečnost ova veličina zavisi od njenog stanja, prven-
stveno od temperature. Važno: naročito mnogo koeficijent površinskog napona
zavisi od prisustva nečistoća. Čak i vrlo male količine nečistoća mnogu znatno
da promene površinski napon.
SI jedinica za koeficijent površinskog napona je N/m, ali se njegove vrednosti
u praksi obično navode u 1000 puta manjim jedinicama - mN/m (mili-njutn po
metru).
γ = γA xA + γB xB . (6.3)
Ovo pravilo može poslužiti samo za grubu procenu jer se u gotovo svim situaci-
jama javljaju znatna odstupanja. Osnovni razlog odstupanja je da je površinski
sloj rastvora znatno bogatiji komponentom rastvora koja ima niži površinski
napon nego ostatak rastvora jer se na taj način ostvaruje niža ukupna energija
sistema, slika 6.3.
γ = γA yA + γB yB , (6.4)
Slika 6.4 Shema metoda otkidanja prstena; θ je kontaktni ugao izmedju tečnosti i
prstena. Predstavljen je prsten sa pravougaonim presekom kakav je na aparturi u
studentskoj laboratoriji. U originalnoj metodi prsten ima kružni poprečni presek.
Sila površinskog napona na svaki segment ove linije je usmerena pod uglom θ u
odnosu na vertikalu, gde je θ kontaktni ugao izmedju prstena i tečnosti. Na taj
način je ukupna sila površinskog napona na prsten jednaka
Slika 6.5 Vremenska zavisnost sile F izmedju prstena i površine tečnosti tokom
različitih faza otkidanja prstena; pozitivna smer za F je vertikalno naviše. (1) - prsten
u vazduhu, F = 0. (2) - prsten uz samu površ tečnosti; sila F > 0 je mala i potiče od
adhezije izmedju prstena i tečnosti. (3) - tokom potiskivanja prstena u tečnost
potrebna je mala sila F < 0 da bi se savladao površinski napon. (4) - prsten je
potpuno potpljen u tečnost; mala sila F > 0 potiče od površinskog napona na
držačima prstena. (5) i (6) - tokom podizanja prstena sila F raste sve dok ne dostigne
maksimum (7). Nakon toga se tokom kidanja dodirne linije prstena i tečnosti
smanjuje te pada na nulu (8) kada je prsten otkinut i visoko iznad tečnosti.
Slika 6.6 Kontaktni ugao izmedju prstena (PHYWE torziona vaga) i vode.
Slika 6.8 a) Kraci torzione vage: (1) - kraci; (2) - točak za podešavanje horizontalnosti
krakova; (3) - markeri horizontalnog položaja kazaljke; (4) - štitnik krakova; (5) -
torziona nit; (6) - veza krakova i torzione niti.
b) - Skala aparature pokazuje silu u mN (najmanji podelak na skali
odgovara sili od 0,1mN): (1) - kazaljka skale; čitljivost skale je povećana na 0,025nN
nonijusnom podelom (2) od 4 podeoka na plastičnom štitniku kazaljke; (3) - točak za
okretanje kazaljke skale;
2. Proveriti libelom da li je štitnik krakova (4), slika 6.8 a), horizontalan; ako
ne, postaviti ga u horizontalni položaj pažljivim pomeranjem aparature u
zavrtnjima kojima je torziona vaga pričvršćena na nosač.
3. Okačiti suv prsten na levi krak torzione vage (1), slika 6.8 a).
4. Postaviti kazaljku skale (1), slika 6.8 b), na nulu.
5. Proveriti da li kraci vage stoje horizontalno; ako ne, podesiti ih okretanjem
točka za podešavanje horizontalnosti krakova (2), slika 6.8 a). Kraci torzione
vage stoje horizontalno kada se nalaze unutar polja ograničenih markerima
(3), slika 6.8 a), na štitniku krakova.
Kraci torzione vage su pričvršćeni za torzionu nit (5), slika 6.8 a), u tački
(6). Okretanjem točka za podešavanje položaja krakova (2), slika 6.8 a), se
uvrće (kraći) deo torzione niti izmedju krakova i točka za nivelisanje krekova
(2). Okretanjem točka skale (1), slika 6.8 b), se uvrće duži deo torzione niti
izmedju krakova i skale aparature.
6.3 Eksperiment
3. Dobro isprati sud za ispitivanu tečnost (2), slika 6.7, prvo običnom pa zatim
destilovanom vodom.
4. Napuniti ovaj sud sa 100ml destilovane vode; zapreminu destilovane vode
odmeriti špricem.
5. Postaviti sud sa destilovanom vodom u (prazno) kupatilo i poklopiti ga.
6. Spustiti prsten aparature kroz levak i otvor poklopca; zatvoriti otvor levkom
pazeći da konac ne dodiruje zidove levka.
7. Proveriti libelom da li postolje kupatila (11), slika 6.7, stoji horizontalno; ako
ne, nivelisati ga šrafovima na koje se postolje oslanja.
Merna faza:
1. Izmeriti koeficijent površinskog napona za destilovanu vodu (bez alkohola).
Silu otkidanja izmeriti tri puta i upisati rezultate u Excel tabelu.10 Ako koris-
tite templejt upišite rezultate u listu Water-Ethanol mixture. Upišite i vred-
nost temperature vode.
Važno: upisujete samo vrednosti za tri merenja sile otkidanja u mN (redom
u kolone mer. 1, mer. 2 i mer. 3) i temperaturu u koloni Trastvor [C]; ostale
vrednosti Excel sam računa:
• Srednju silu Fsr računa po formuli:
=(Fn+Gn+Hn)/3
gde je n=2,3,... redni broj vrste u kojoj se upisuju rezultati.
(čistog) alkohola, a VmW = 18, 0677 cm3 molarna zapremina vode, obe
na 25◦ C. Za 96% alkohol koji se koristi je p = 0, 96.
Važno: u listi BasicData su navedeni pomenuti podaci kao i svi ostali
relavantni podaci (npr sve dimenzije prstena) za eksperiment.
• Koeficijent površinskog napona γ Excel računa po formuli
=0.875*En/(0.12235*COS(RADIANS(13.5)))
gde je n=2,3,... redni broj vrste, a En vrednost srednje sile otkidanja u
mN; ovo u stvari formula (6.10) u kojoj se koristi izmerena vrednost
spoljašnjeg+unutrašnjeg obima prstena (=0,12235m), zatim izmereni
kontaktni ugao izmedju prestena i vode (=13.5◦ ), kao i korekcioni faktor
za naš prsten C = 0, 875.11
2. Dodavati malim špricem alkohol u rastvor tako da zapremina 96% alkohola
u rastvoru iznosi redom 5, 10, 15, 20, 30, 40 i 50 ml.12 Svaki put nakon što
dodate alkohol sačekajte 2-3 minuta, pa zatim izmerite koeficijent površinskog
napona rastvora.
3. Na kraju izmerite površinski napon za 96% alkohol (bez vode).
Obrada podataka:
6. Naliti vruću vodu u kupatilo do 1cm ispod vrha; ako sud sa destilovanom
vodom počne da pliva isprazniti nešto vruće vode iz kupatila špricom.
7. Postaviti termometre: jedan u kupatilo, a drugi uz površ destilovane vode
tako da vrh njegove sonde bude u vodi.
8. Pratiti porast temperature na oba termometra (temperatura destilovane vode
će biti niža od temperature vode u kupatilu). Sačekati dok temperatura des-
tilovane vode ne poraste na oko 70◦ C. Ako temperatura destilovane vode ne
uspeva da se povisi do ∼ 70◦ C, izvući par puta špricom vodu iz kupatila pa
dopuniti kupatilo vrućom vodom iz crvenog balona. Ponavljati ovo sve dok se
ne dostigne željena temperatura destilovane vode.
9. Spustiti prsten aparature kroz levak i otvor poklopca; zatvoriti otvor levkom
pazeći da konac ne dodiruje zidove levka.
10. Merenje koeficijenta površinskog napona započeti kada destilovana voda pri-
bližno dostigne željenu početnu temperaturu.
Merna faza:
Obrada podataka:
7.1 Uvod
Pri kretanju realnih fluida, delići fluida koji su u kontaktu sa nepomičnim zi-
dovima (npr cevi kroz koju fluid struji) su ”slepljeni” za zidove dok se ostali
delići kreću brzinama koje rastu kako se udaljavamo od nepomičnih zidova. Pri
relativno sporom strujanju, formiraju se slojevi fluida pri čemu su brzine delića
iz istog sloja iste, dok se brzine slojeva (u opštem slučaju) razlikuju. Za kretanje
fluida kod kojeg su formirani slojevi kažemo da je laminarno.
Izmedju paralelnih slojeva fluida koji se kreću različitim brzinama deluju sile
viskoznog trenja koje nastoje da uspore brži, a ubrzaju sporiji sloj, slika 7.1. Ove
sile su tangencijalne na dodirnu površinu slojeva i stvaraju tangencijalni napon
τ = dF/dS izmedju slojeva. Ovaj napon je, po Njutnovom zakonu viskoznosti,
oblika
dv
τ =η , (7.1)
dz
gde je dv/dz gradijent brzine u pravcu (z) normalnom na brzinu. Koeficijent
proporcionalnosti η se naziva koeficijent viskoznosti fluida. On zavisi od vrste i
stanja fluida. Fluidi čija je viskoznost opisana Njutnovim zakonom viskoznosti
se nazivaju njutnovski fluidi.
7.1 Uvod 67
ose sistema, brzinom v. Zbog simetrije sistema, ova brzina zavisi samo od r, što
znači da svi delići fluida na cilindru radijusa r ravnomerno rotiraju oko ose sis-
tema ugaonom brzinom ω = v/r, koja takodje zavisi samo od r. Sila viskoznog
trenja dF⃗vis , koja deluje na element površi posmatranog cilindra dS
⃗ = dS⃗er , gde
je ⃗er jedinični vektor normalne na površ, je oblika
dF⃗vis = ⃗τvis dS ,
(1) (1)
delovali momenti Mvis (r + dr) ̸= Mvis (r), što znači da bi rotacija posmatranog
sloja bila ubrzana. Tako nalazimo da je
dω
r3 = const ,
dr
odakle je
A
ω(r) = +B, (7.10)
r2
gde su A i B integracione konstante koje nalazimo iz graničnih uslova. Pošto je
fluid ”slepljen” na svojim granicama, to je ω(R2 ) = 0 (jer je spoljašnji cilindar
nepokretan) i ω(R1 ) = Ω (jer tom ugaonom brzinom rotira unutrašnji cilindar).
Odavde je
( )
R12 R22 1 1
w(r) = Ω 2 − 2 , (7.11)
R2 − R12 r2 R2
zbog čega je brzina delića fluida v = rω jednaka
( )
R12 R22 1 r
v(r) = Ω 2 − 2 , (7.12)
R2 − R12 r R2
a intenzitet momenta sile po jedinici dužine
(1) R12 R22
Mvis = 4πη Ω, (7.13)
R22 − R12
(1)
Odavde nalazimo da je za konačne cilindre viskozni moment Mvis = Mvis H
jednak
R2 R2
Mvis = 4πη 2 1 2 2 HΩ , (7.14)
R2 − R1
pod uslovom da je visina cilindra H znatno veća od R1 i R2 te se efekti krajeva
mogu zanemariti.
Iz izraza (7.14) za moment viskozne sile se vidi da je ovaj moment oblika
Mvis = CηΩ , (7.15)
tj da je proporcionalan proizvodu koeficijenta viskoznosti η i brzine rotacije un-
R 2 R2
utrašnjeg cilindra Ω, sa konstantom proporcionalnosti C = 4π R21−R22 H. Iako
2 1
naš zaključak vredi za veoma duge (teorijski - beskonačne) cilindre, propor-
cionalnost oblika (7.15) je ispunjena i za konačne cilindre, jedino što konstanta
proporcionalnosti C ima neku drugu vrednost koja se za dati uredjaj utvrdjuje
eksperimentalno.
7.2 Opis aparature 71
7.3 Eksperiment
Priprema merenja:
Pre početka merenja potebno je obaviti sledeću pripremu:
13. Ući u meni Measurement, slika 7.5 b); displej pokazuje poruku:
”—Measure Config. —”
a ispod nje podatke o priključenom vretenu u polju ”SP:” (u ovom zadatku to
je L1) i tekuću brzinu u polju ”RPM:”, npr 100. Aktivno polje menija blinka.
U njenu se može menjati uneti podatak; promena se vrši strelicama (3) i (4),
slika 7.5. Prelaz sa jednog na drugo polje se vrši pomoću tastera TAB (2). Na
dnu displeja je prikazana maksimalna vrednost koeficijenta viskoznosti koja
se može meriti pri datom izboru vretena i brzine njegove rotacije.
VAŽNO: viskozimetar daje pouzdane rezultate (sa greškom ±1%) samo
u intervalu izmedju 15% i 90% tekućeg opsega. Izvan njega izmerene vred-
nosti imaju znatno veću grešku, a viskozimetar počinje da se zvučno oglašava.
Naglasimo da uredjaj prikazuje usrednjene podatke merenja; potrebno je
sačekati da se izvrši najmanje 5 rotacija vretena da bi se rezultati prikazani
na displeju mogli smatrati verodostojnim.
14. Merenje se startuje pritiskom na taster ON (6), slika 7.5 b). Kada ga pri-
tisnete vreteno počinje da rotira, a kada ga pritisnete još jednom - vreteno
se zaustavlja, a displej prikazuje ”Motor stop”. Ako ponovo pritisnute ON,
rotacija vretena se nastavlja.
15. Iz merenja se izlazi pritiskom na taster QUIT (5), slika 7.5 b).
Postupak merenja:
=(An/342.3)/(An/342.3+(100-An)/18.015)
gde je n redni broj ćelije (n=2,3,...) u koloni A, a 342.3 i 18.015 su molarne mase
šećera i vode u gramima.
Obrada podataka
1. Nacrtajte koeficijent viskoznosti glicerina u funkciji temperature; koristite
program Table Curve 2D.
2. Nacrtajte prirodni logaritam koeficijenta viskoznosti glicerina η(T ) na apso-
lutnoj temperaturi T u funkciji recipročne vrednosti 1/T apsolutne tempera-
ture T .
3. Nafitujte ovako predstavljene eksperimentalne podatke na Andrade formulu
(7.2) kao modelnu funkciju. Eksperimentalne tačke bi trebalo da leže na
pravoj liniji y = a + bx, gde je x = 1/T , a y = ln(η). Koeficijent ove
prave bi trebalo da odgovara količniku Ea /T , a njen odsečak na y-osi veličini
ln(η∞ ) + Ea /(RT0 ), odakle je potrebno odrediti aktivacionu energiju Ea , kao
i koeficijent viskoznosti η∞ na beskonačnoj temperaturi i proceniti njihove
greške.
4. Nacrtajte koeficijent viskoznosti ηx rastvora šećera u funkciji molarnog udela
šećera x.
5. Nacrtajte prirodni logaritam koeficijenta viskoznosti rastvora šećera ln(ηx )
u funkciji molarnog udela šećera x. Nafitujte eksperimentalne podatke na
modelnu funkciju (7.7), tj na
a + bx
y= ,
1 + cx
gde je y = ln(ηx ), te odredite empirijske konstante a, b i c i procenite njihove
greške.
8 Termodinamički ciklusi - Stirling
mašina
8.1 Uvod
1 Tokom relaksacije svaki lokalni parametar a dostiže postepeneo svoju ravnotežnu vrednost
tako što razlika ∆a izmedju njegove tekuće i ravnotežne vrednosti opada vremenom na
nulu po eksponencijalnom zakonu, ∆a = ∆a0 e−t/τa , gde je τa karakteristično vreme
potrebno da razlika za parametar a opadne e puta; ovo vreme se naziva vremenom
relaksacije parametra a. Za vreme relaksacije celog sistema se uzima supremum svih τa i
on je konačan, osim u slučaju kritičnih ravnotežnih stanja. Smatra da je sistem izrelaksirao
nakon isteka nekoliko njegovih vremena relaksacije.
2 Stanje se smatra skoro ravnotežnim ako bi promene u sistemu bile zanemarive tokom
relaksacije iz tog stanja pri zaustavljenim promenama okoline.
78 Toplotne mašine
Idealni Stirling ciklus je reverzibilni ciklus koji se sastoji od dve izoterme i dve
izohore idealnog gasa,4 slika 8.1; broj molova gasa n je fiksan.
pV = nRT ,
ali ovaj rad ne odlazi na povećanje unutrašnje energije (ona ostaje kon-
stantna) već se u vidu toplote QIII = WIII predaje okolini.
Izohorsko grejanje (proces IV na slici 8.1):
Iz stanja 4 do stanja 1, gas se izohorski greje na zapremini V = V1 . Sistem
ne vrši rad već uzima iz okoline količinu toplote QIV = cV ∆T i povećava
svoju unutrašnju energiju i temperaturu sa T2 na T1 .
tokom izohorskog grejanja (proces IV na slici 8.1). Stoga ukupna količina toplote
koju primi gas tokom jednog ciklusa iznosi:
Q = QI + QIV = nRT1 ln(V2 /V1 ) + ncV (T1 − T2 ) . (8.6)
Kako je ukupni rad koji posmatrani sistem izvrši nad okolinom
W = WI − WIII = nR(T1 − T2 )ln(V2 /V1 ) , (8.7)
to koeficijent korisnog dejstva ηs idealnog Stirling ciklusa idealnog gasa iznosi
T1 − T2 1
ηs = · T1 −T2
, (8.8)
T1 1+ T1 · cV
ln(V2 /V1 )
odnosno
ηc
ηs = , (8.9)
1 + ηc cv /ln(V2 /V1 )
gde je
T1 − T2
ηc = , (8.10)
T1
koeficijent korisnog dejstva idealnog Carnot ciklusa5 koji radi na iste dve tem-
perature T1 i T2 .
5 Carnot ciklus je reverzibilni ciklus koji se sastoji iz dve izoterme (na temperaturama T1 i
T2 ) i dve adijabate; radni sistem je proizvoljan. Po Carnot teoremi, koeficijent korisnog
dejstva bilo kog ciklusa koji razmenjuje toplotu samo na ove dve temperature ne može biti
veći od koeficijenta korisnog dejstva Carnot ciklusa ηc = (T1 − T2 )/T1 .
6 Patentirao je 1816. godine škotski sveštenik i izumitelj Robert Stirling (1790. - 1878.);
prvu praktičnu verziju mašine Stirling je napravio 1818. godine i koristio za pogon pumpe
za vodu u kamenolomu.
8.3 Opis aparature 81
Priprema aparature
jednog termopara je dat na slici 8.6. Sondu (1), slika 8.6, termopara za
merenje temperature T1 pažljivo uvući u Stirling mašinu kroz rupicu
(4), slika 8.3, a konektor (3), slika 8.6, tog termopara povezati na T1
utičnicu p-V-n-T metra, (7) na slici 8.7; na isti način sondu termopara
za temperaturu T2 uvući kroz rupicu (5), a konektor uključiti u T2 .
Slika 8.3 PHYWE Stirling mašina: (1) - zamajac; (2) - vertikalni cilindar i radni klip
(2’); (3) - pVn senzor; (4) - rupica za sondu termopara koji meri temperaturu T1 ; (5)
- rupica za sondu termopara koji meri temperaturu T2 ; (6) - polukružni zamajac; (7) -
crevo za povezivanje sa pVn senzorom; (8) - .
2. Zašrafiti motor/generator, slika 8.5, na postolje Stirling mašine, slika 8.3, i gu-
menim remenom (3), slika 8.5, povezati širi kanal osovine motora sa zamajcem
(1) Stirling mašine.
3. Napuniti gorionik, slika 8.4, gorivom:
i) - ukloniti dimnjak (3);
ii) - odvrnuti glavu gorionika (1);
iii)- isprazniti rezervoar;
iv)- izmeriti masu praznog rezervoara;
v) - napuniti rezervoar gorionika (2) alkoholom skoro do vrha;
vi)- izmeriti masu punog rezervoara i izračunatu masu goriva mg u rez-
ervoaru;
vii)- zašrafiti glavu gorionika na rezervoar i zavrtnjem (5) podesiti da fitilj (5)
viri oko 5 mm.
Slika 8.4 Gorionik Stirling mašine: (1) - glava gorionika; (2) - rezervoar za gorivo
(alkohol) na koji se šrafi glava gorionika; (3) - dimnjak; (4) - fitilj; (5) - zavrtanj za
regulisanje visine fitilja; (6) - poklopac glave gorionika.
Slika 8.5 Motor/generator za PHYWE stirling mašinu: (1) - ulaz za napajanje motora
jednosmernim naponom ±3V; (2) - naponski izlaz iz generatora; (3) - gumeni remen
za povezivanje na zamajac Stirling mašine.
84 Toplotne mašine
Slika 8.7 PHYWE p-V-n-T metar: displej (1) prikazuje broj ciklusa Stirling mašine po
minuti, n; kada je dugme (4) u polju Display isključeno, displej (2) pokazuje
temperaturu T1 u ◦ C i svetli gornja dioda, a kada je dugme (4) uključeno displej (2)
pokazuje razliku temperatura ∆T (= T1 − T2 ) i svetli donja dioda; displej (3) pokazuje
temperaturu T2 u ◦ C; u polju Calibration dugme (5) je za kalibraciju razlike
temperatura ∆T , a dugme (6) za kalibraciju zapremine V ; utičnice za dva NiCr/Ni
termopara: (7) - za temperaturu T1 i (8) - za temperaturu T2 ; (9) - DIN utičnica
(8-polna) na koju se povezuje p-V-n senzor (3) na Striling mašini, slika 8.3; preko ove
utičnice se p-V-n senzor napaja i na nju se iz senzora dovode izmereni signali za
pritisak, zapreminu i broj ciklusa koje procesuira i prikazuje p-V-n-T metar; polje
Oscilloscope sadrži dva analogna BNC izlaza za povezivanje na osciloskop ili
kompjuterski interfejs: izlaz (10) daje napon proporcionalan izmerenom pritisku P , a
izlaz (11) napon proporcionalan izmerenoj zapremine V . U konfiguraciji sa Cobra 3
interfejsom, P izlaz se dužim BNC kablom povezuje na analogni ulaz Analog in 2/S2
Cobra 3 interfejsa, (2) na slici 8.8, a V izlaz kraćim BNC kablom na Analog in 1/S1,
(1) na slici 8.8.
86 Toplotne mašine
Slika 8.8 PHYWE interfejs Cobra 3: (1) - analogni ulaz Analog in 1/S1 (on se preko
BNC-banana adaptera, postavljenim tako da je BNC priklučak iznad banana,
povezuje kraćim BNC kablom sa izlazom V za izmerene vrednosti zapremine, (11) na
slici 8.7 p-V-n-T metra); (2) - analogni ulaz Analog in 2/S2 (2) koji se na isti način
dužim BNC kablom povezuje sa izlazom P za izmerene vrednosti pritiska, (10) na
slici 8.7.
Slika 8.9 pVn senzor za Stirling mašinu: (1) - kabl za povezivanje na p-V-n ulaz na
p-V-n-T metru, (9) na slici 8.7; (2) - crevo kojim se vazduh iz Stirling mašine prenosi
do senzora pritiska radi merenja; (3) - inkrementalni prenosnik koji radi merenja
zapremine V i broja obrtaja n mora biti povezan na osovinu polukružnog zamajca (6)
Stirling mašine, slika 8.3.
8.3 Opis aparature 87
Pre upotrebe potrebno je kalibrisati PHYWE p-V-n-T metar, slika 8.7; to se radi
po sledećoj proceduri:
Slika 8.10 Glavni prozor PHYWE programa ”measure”: (1) - linija menija; (2) - linije
alatki; (3) - padajući meni koji se otvara klikom na meni Gauge.
Slika 8.11 Cobra 3 - Univesal Writer prozor PHYWE programa ”measure”. Aktivan je
tab ”Normal Measurement”. Brzina uzorkovanja (Get value) je 2ms, merenje se
započinje (Start) i zaustavlja (End) klikom (on key press).
Tok eksperimenta
Slika 8.14 Postavljanje parametara merenja u prozoru Universal Writer. Pri tekućem
izboru frekvencije (5 kHz) i broju uzoraka (postavljenom na maksimalnu vrednost
2048), jedno merenje traje oko 0,4 sekunde.
moment sila trenja koja nastoji da izvede merač iz vertikalnog položaja. Inten-
zitet ovog momenta zavisi od sile kojom je merač pritegnut na osovinu. Sa druge
strane, kada merač nije u vertikalnom položaju javlja se i moment gravitacione
sile koji nastoji da merač vrati u vertikalni položaj, a čiji je intenzitet propor-
cionalan sinusu ugla odstupanja merača od vertikale. Ravnotežni ugao otklona
se uspostavlja kada se pomenuta dva momenta izjednače i podešava se zavrtan-
jem (odvrtanjem) zavrtnja (3). Napomenimo i da ravnoteža nije idealna zbog
neizbežnih ekscentriteta osovine i ”truckanja” koje se javlja pri proklizavanju.
Zato merač ne miruje na ravnotežnom uglu otklona, već osciluje oko njega.
P V dijagram Stirling mašine merimo pod opterećenjem zadatim momentom
sile trenja (izmedju merača momenta sile i osovine zamajca). Ravnotežni ugao
8.4 Odredjivanje P V dijagrama i efikasnosti Stirling mašine 93
Slika 8.16 Merač momenta sile (1) montiran na osovinu zamajca Stirling mašine. (2) -
skala za očitavanje momenta sile; (3) - zavrtanj za pritezanje merača.
4. Da bi ste precizno odredili brzinu obrtanja ω, isključite sve kanale osim jednog,
npr U1 ; nafitujte ovaj kanal na kao na slici 8.19; dobijenu vrednost za ω
(parametar fita b) upisati na listu Stirling excel tabele.
Slika 8.19
Proračun P V dijagrama
Slika 8.21
Slika 8.22
4. Kada se otvori novi prozor ”Convert relation to function”, slika 8.23, izabrati
opciju Keep measurement in relation mode, pa kliknuti OK.
Slika 8.23
Slika 8.24
8.4 Odredjivanje P V dijagrama i efikasnosti Stirling mašine 99
9.1 Uvod
P⃗ = −κ∇T , (9.4)
T2 − T1
P =κ , (9.8)
L
što se može iskoristiti za merenje koeficijenta toplotne provodljivosti κ pod
uslovom da možemo meriti gustinu fluksa toplotnog provodjenja P i temperature
T1 i T2 .
dok je granici x = L
T2 − T1 ′ T2 − T1
A#
L = k L + T1 + , (9.19)
2 + k′ L 2 + k′ L
Slika 9.2 Crnom linijom je data kriva tg(λL), a crvenom kriva 2/(λ/k − k/λ) za L = 1
i k = 0.25. Svojstvene vrednosti λn se nalaze na preseku ove dve krive: λ1 ≈ 0.693,
λ2 ≈ 3.293, λ3 ≈ 6.362 ...; vidimo da je već za n = 3 svojstvena vrednost
λ3 ≈ 2π/L = 2π = 6.283.
106 Toplotna provodnost
Slika 9.5 Levo - unutrašnjost kućice: (1) - crna kutija sa sijalicom; (2) -
termoparovi nalepljeni na unutrašnje zidove materijala kojem se meri toplotna
provodljivost; (3) - sonde 4 termometra za merenje temperature vazduha u kućici; (4)
- bočni prozori izolovani stiroporom; (5) - senzor za temperaturu vazduha iznad crne
kutije; (6) - ventilatori. Desno: - sijalica za grejanje kućice montirana u crnoj kutiji.
9.3 Odredjivanje koeficijenta toplotne provodljivosti 109
6. Pre nego što uključite reflektor očitajte sve temperature i zapišite očitane
vrednosti u Excel tabelu. Na raspolaganju vam je excel templejt tabela TP Template.xls
koju pre korišćenja treba da preimenujete tako da se vidi koja je grupa i kada
radila eksperiment. Podatke upisujete na list ”Staklena basta” - vidi sliku 9.7.
Slika 9.7 List ”Staklena basta” excel templejta; u koloni J excel prikazuje srednje
temperature, npr ”Tvb s” je srednja temperatura vazduha bliže staklu.
Slika 9.8 List ”Kucica” excel templejta; u koloni M excel prikazuje srednje
temperature, npr ”Tamb s” je srednja ambijentalna temperatura.
Merenje na kućici:
odnosno C # = (C0# + CL #
)/2.
9.4 Dodaci 113
9.4 Dodaci
Po Fourier zakonu je
P⃗ = −κ∇T , (9.30)
Drugim rečima, raspodela Th (t, x) mora biti rešenje jednačine provodjenja toplote
(9.45), sa homogenim graničnim uslovima (9.46) i (9.47), i da zadovoljava mod-
ifikovani početni uslov (9.48).
Partikularna rešenja jednačine (9.45) nalazimo u obliku proizvoda
funkcije Θ(t) (koja zavisi samo od vremena t) i funkcije X(x) (koja zavisi samo
od koordinate x). Smenom (9.49) u (9.45) nalazimo
( )
1 1 dΘ 1 d2 X
· = · , (9.50)
D Θ dt X dx2
što znači da su obe strane jednačine konstantne. Neka je
1 d2 X
· = p = const ,
X dx2
9.4 Dodaci 117
odnosno
d2 X
= pX . (9.51)
dx2
Pored jednačine (9.51), na granicama ploče x = 0 i x = L funkcija X(x) treba
da zadovolji i granične uslove
dX
= k ′ X|x=0 (9.52)
dx x=0
dX
= −k ′ X|x=L (9.53)
dx x=L
koji slede iz graničnih uslova (9.46) i (9.47) za raspodelu temperature Th (t, x).15
Jednačina (9.51) je obična linearna diferencijalna jednačina II reda sa kon-
stantnim koeficijentima. Za p = −λ2 < 0, njeno opšte rešenje je
X(x) = A cos(λx + ϕ) , (9.54)
gde su A i ϕ proizvoljne konstante. Obzirom da je
dX/dx = −λA sin(λx + ϕ) ,
∫
15 Funkcije f (x) na segmentu [0, L] koje su kvadratno integrabilne ( 0L f (x)2 dx < ∞) čine
realan Hilbertov prostor koji se označava sa L2 [0, L]; skalarni proizvod funkcija f (x) i g(x)
∫
je u ovom prostoru dat sa (f, g) = 0L f (x)g(x)dx. Funkcije sa neprekidnim drugim
izvodom čine linearnu mnogostrukost C 2 [0, L] u prostoru L2 [0, L] koja je gusta, tj
najmanji podprostor C 2 [0, L] koji sadrži C 2 [0, L] je L2 [0, L]. Izraz d2 /dx2 je linearni
operator iz C 2 [0, L] u C 1 [0, L] ⊂ L2 [0, L]. Restrikcija D̂2 ≡ d2 /dx2 |D ovog operatora na
podprostor D funkcija iz C 2 [0, L] koje zadovoljavaju granične uslove (9.52) i (9.53):
df /dx = kf za x = 0, odnosno df /dx = −kf za x = L, je operator Sturm-Liouville tipa i
taj operator je hermitski. Funkcija f (x) za koju je d2 f /dx2 = pf , vidi (9.51), je tzv
svojstvena funkcija operatora d2 /dx2 koja odgovara svojstvenoj vrednosti p. Svojstvene
funkcije koje odgovaraju istoj svojstvenoj vrednosti p čine tzv svojstveni podprostor (za
svojstvenu vrednost p); svaki svojstveni podprostor operatora d2 /dx2 je dvodimenzionalan,
npr eλx i e−λx su linearno nezavisne svojstvene funkcije za p = λ2 , a cos(λx) i sin(λx) za
p = −λ2 , dok su const i ax svojstvene funkcije za p = 0. Pokazuje se da su sve svojstvene
vrednosti za Sturm-Liouville operator D̂ (zbog graničnih uslova) negativne p = −λ2n < 0,
da se mogu prebrojati sa n = 1, 2, ... tako da je λ1 < λ2 < ..., a da je svaki svojstveni
podprostor jednodimenzionalan. Svojstvene funkcije fn (x) i fm (x) su ortogonalne, tj
njihov skalarni proizvod (fn , fm ) je nula jer je
p′ (fp′ , fp′′ ) = (p′ fp′ , fp′′ ) = (D̂fp′ , fp′′ ) = (fp′ , D̂fp′′ ) = (fp′ , p′′ fp′′ ) = p′′ fp′ , fp′′ ) (∗)
pa je (p′ − p′′ )(fp′ , fp′′ ) = 0, odakle je (fp′ , fp′′ ) = 0 zbog p′ ̸= p′′ . Ove funkcije čine jedan
ortogonalni bazis po kojem se može razviti svaka funkcija f iz domena D Sturm-Liouville
operatora D̂:
∞
∑
f (x) = An fn (x) ,
n=1
∫L
f (x)fn (x)dx
(f, fn ) 0
An = = .
(fn , fn ) ∫L
fn2 (x)dx
0
118 Toplotna provodnost
15
videti izraz (*) u fusnoti , dok se normalizacioni faktor
∫L
Nn2 ≡ fn2 (x)dx , (9.63)
0
Θn (t) = e−Dλn t ,
2
(9.65)
te su
Tn (t, x) = An e−Dλn t cos (λn x + ϕn ) ,
2
n = 1, 2, ... (9.66)
2
Zbog faktora eDλ t ovo rešenje daje neograničene varijacije temperature tokom
vremena te fizički nije prihvatljivo, sem za A′ = B ′ = 0 kada je Tp>0 (t, x) ≡ 0.
Obzirom da je sistem jednačina (9.45), (9.46) i (9.47) homogen, zbir njegovih
partikularnih rešenja (9.66) je takodje rešenje istog sistema, te opšte rešenje
sistema jednačina (9.45), (9.46) i (9.47) u proizvoljnom trenutku t > 0 glasi
∞
∑
An e−Dλn t cos (λn x + ϕn ) ,
2
Th (t, x) = (9.67)
n=1
120 Toplotna provodnost
dok je za t = 0 ispunjeno
∞
∑ ∞
∑
Th (0, x) = An cos (λn x + ϕn ) = An fn (x) . (9.68)
n=1 n=1
( )
T2 − T1 T2 − T1
+ kx + T1 + , (9.70)
2 + k′ L 2 + k′ L
koje se na granici x = 0 svodi na
( )
T2 − T1
+ T1 + , (9.71)
2 + k′ L
a na granici x = L na
( )
T2 − T1 ′ T2 − T1
+ k L + T1 + . (9.72)
2 + k′ L 2 + k′ L
17 Ako pomnožimo (9.68) sa fm (x) = cos(λm x + ϕm ) i prointegralimo od 0 do L nalazimo
∫L ∞
∑ ∫L
2
Th (0, x)fm (x)dx = An fn (x)fm (x)dx = Nm Am .
0 n=1 0
9.4 Dodaci 121
10.1 Uvod
(10.1)
U čvrstim telima zvukom nazivamo periodično oscilovanje delića koje se javlja
u vidu talasa. Najjednostavniji opis imaju ravni monohromatski (harmonijski)
putujući talasi; oni se prostiru (putuju) duž jednog pravca i smera duž kojeg
trasiramo x-osu i sve veličine zavise samo od x-koordinate i vremena. Pomeranje
delića je dato izrazom
[ ( ) ]
⃗ cos ω t − x + ϕ ,
⃗u(t, x) = A (10.2)
c
gde je ω = 2πf ugaona frekvencija (pri običnoj frekvenciji f ), c je brzina talasa,
dok je ϕ fazni pomak. Kod ovog talasa se pomeranje ⃗u vrši duž pravca vektora
⃗ čiji intenzitet A = |A|
A, ⃗ ima smisao amplitude talasa. U svakoj tački prostora
delić osciluje prosto-periodično po zakonu
( )
⃗ cos(ωt + ϕ0 ) = A
⃗u = A ⃗ cos 2π t + ϕ0 ,
⃗ cos(2πf t + ϕ0 ) = A
T
gde je
x0
ϕ0 = ϕ − ω ,
c
fazni pomak delića, x0 njegova x-koordinata, a
2π 1
T = = ,
ω f
(vremenski) period oscilacija. Sa druge strane, u bilo kom trenutku t0 pomeranje
delića u prostoru je takodje prosto-periodično:
( ) ( )
⃗ cos ω x + ϕ′ = A
⃗u = A ⃗ cos 2π x + ϕ′ ,
⃗ cos (kx + ϕ′ ) = A
0 0 0
c λ
gde je
ϕ′0 = −(ϕ + ωt0 ) ,
fazni pomak u trenutku t0 . Veličina λ je talasna dužina i ona ima smisao perioda
u prostoru, dok se veličina
ω 2π 2πf
k= = = ,
c λ c
10.1 Uvod 123
naziva (ugaoni) talasni broj i stoji u istom odnosu sa periodom u prostoru kao i
ω sa periodom u vremenu. Vredi i
c = fλ .
Ako se pomeranje vrši duž pravca prostiranja talasa, talas je longitudinalan
(A⃗ = A⃗ex ), a ako se vrši normalno na njega (A ⃗ leži u yz ravni) talas je transverza-
lan. Lako se vidi da pri longitudinalnom talasu dolazi do periodičnog zgušnjavanja
i razredjivanja sredine kroz koju se talas prostire te za longitudinalni talas kažemo
da je kompresioni. Pri transverzalnom talasu dolazi pak do periodičnog smicanja
slojeva normalno na pravac talasa.
Zvuk se javlja i u fluidima (gasovima i tečnostima). Za razliku od zvuka u
čvrstim telima, u fluidima se javljaju samo longitudinalni talasi jer se fluidi ne
opiru smicanju. Napomenimo i da u gasovima nema smisla govoriti o oscilacijama
delića; ono što osciluje je lokalna vrednost pritiska i gustine gasa.
Superpozicijom (sabiranjem) talasa nastaje novi talas. Tako na primer super-
pozicijom dva ravna monohromatska talasa iste amplitude A, ⃗ od kojih prvi
⃗u+ (t, x) = A cos [ω (t − x/c) + ϕ+ ]
⃗
putuje u smeru x-ose, dok drugi talas
⃗u− (t, x) = A ⃗ cos [ω (t + x/c) + ϕ− ]
putuje u suprotnom smeru, nastaje tzv stojeći talas
[ ]
⃗u(t, x) = ⃗u+ (t, x) + ⃗u− (t, x) = 2A ⃗ cos(kx − ∆ϕ) cos(ωt + ϕ) ,
Kada se kroz cev dužine L koja je zatvorena na jednom kraju prostire ka tom
kraju ravan monohromatski putujući talas
Stojeći talas se javlja i pri prostiranju zvuka kroz cev koja je otvorena na oba
kraja. Sličnim postupkom kao i u prethodnom odeljku nalazimo da se rezonancija
javlja ako je dužina cevi jednaka celobrojnom umnošku polovina talasnih dužina
n
λ=L, (10.8)
2
slika 10.2, što znači da se pri fiksiranoj dužini cevi rezonancija javlja za niz
frekvencija
c
fn = n , n = 1, 2, ... (10.9)
2L
dok se pri fiksiranoj frekvenciji f0 , a varijabilnoj dužini cevi, rezonancija javlja
za niz dužina cevi
c
Ln = n , n = 1, 2, ... . (10.10)
2f0
Posmatra se zategnuta žica dužine L čiji su krajevi fiksirani, npr žica na gitari.
Kada se ova žica malo izvede iz ravnotežnog položaja i pusti, doći će do njenih
malih transverzalnih oscilacija koje imaju karakter stojećih talasa.1
Stojeći talasi koji se obrazuju u žici moraju imati čvorove na krajevima žice.
Ovo znači da se mogu javiti samo stojeći talasi za koje je celobrojni umnožak
polovine njihove talasne dužine λ jednak dužini žice L:
λn
n = L , n = 1, 2, ... , (10.11)
2
što znači da će dozvoljne talasne dužine iznositi
2L
λn = , n = 1, 2, ... . (10.12)
n
Najveću talasnu dužinu ima talas za n = 1 i ona iznosi λ1 = 2L. Ovaj talas ima
samo jedan trbuh koji se nalazi na polovini žice, slika 10.3 levo. Za talasnu dužinu
λo = λ1 = 2L kažemo da je osnovna, a kraće talasne dužine (koje odgovaraju n =
2, 3, ...) nazivamo njenim harmonicima. Prvom harmoniku (slika 10.3 sredina)
odgovara λ2 = L, dok je za drugi harmonik (slika 10.3 desno) λ3 = 23 L.
1
Slika 10.3 Prva tri stojeća talasa u žici. Za osnovnu talsnu dužinu je L = λ ,
2 o
dok je
3
za njen prvi harmonik i drugi harmonik L = 1λ2 i L = λ ,
2 3
respektivno.
Stojeći talasi u žici se mogu predstaviti kao superpozicija dva putujuća talasa
koji se p̌rostiru duž žice u suprotnim smerovima.2 Brzina ovih talasa iznosi
√
T
c= , (10.13)
µ
gde je T sila zatezanja žice, a µ podužna masa, tj masa po jedinici dužine. Pri
tome, frekvencije f ovih putujućih talasa, a time i frekvencije stojećih talasa koji
nastaju njihovom superpozicijom, stoje u dobro poznatom odnosu sa talasnom
dužinom:
c = fλ . (10.14)
Odavde se vidi da stojećim talasima odgovaraju tzv prirodne frekvencije žice
c
fn = n , (10.15)
2L
1 Oscilacije žice koje nisu male nastaju pri većim pobudama žice i pokoravaju se složenijim
zakonitostima.
2 Amplitude ovih talasa su jednake i prilikom refleksije na krajevima žice njihova faza trpi
skok za 180◦ .
128 Zvuk
Opis aparature
c = 2bL .
a=0,
Opis aparature
Aparatura za izučavanje slobodnih i prinudnih oscilacija monokorda je prikazana
na slici 10.6.
Povezivanje aparature
Pre početka eksperimenta potrebno je:
Povezivanje aparature ne vrše studenti, već stručni saradnik ili nastavnik. Pre
početka eksperimenta potrebno je da se studenti u njihovom prisustvu uvere da
je aparatura povezana.
1”, ”Graph 2”,... itd redom. Ako želite da izmenite naziv grafika postupite
kao i kod izmene naziva snimka.
Napomena: preporučuje se da snimcima date imena po kojima ih lako pre-
poznajete, a da grafike nazovote po snimcima.
12. Iako DataStudio sadrži opciju za nalaženja spektra signala (”FFT” u polju
”Dispalys”), ova opcija nije dovoljno fina za naše potrebe.9 Za potrebe spek-
tralne analize preporučujemo korišćenje programa Origin. U tu svrhu je potrebno
eksportovati podatke. Kliknite na ”File” meni, pa na ”Export Data...”, nakon
čega se otvara prozor ”Choose a data measurement”, slika 10.11, u kojem treba
da izaberete podatke koji se eksportuju iz liste svih dostupnih podataka.
Označite mišem željene podatke, pa kliknite na ”OK” nakončega se otvara
standardni Windows prozor za izbor imena fajla (u kojem će te sačuvati po-
datke) i lokacije (foldera) u kojem će se taj fajl nalaziti. Izberite ekstenziju
”dat” za fajlove sa sačuvanim podacima (to je za Origin podrazumevana ek-
stenzija za fajlove sa podacima).
13. Nakon što ste sačuvali podatke za rad u Originu, izbrišite sve nepotrebne
stvari iz DataStudio aktivnosti i sačuvajte je pod imenom koje ukazuje na tip
aktivnosti (npr SlobOsc) praćeno nazivom vaše grupe (da se zna ko je radio).
9 Dužine sekvenci signala su relativno male (do 4096 tačaka) i što je posebno važno ne
koriste se prozori što unosi deformacije u spektar.
10.3 Monokord - slobodne i prinudne oscilacije 139
tm = m∆Ts , m = 0, 1, 2, ..., N
10 Umesto obične
√ frekvencije f može se koristiti i ugaona frekvencija ω = 2πf .
11 Ovde je c = T /µ brzina transverzalnih talasa pri sili zatezanja T i podužnoj masi µ. Na
osnovnoj frekvenciji f0 talasna dužina je λ = 2L, a breg talasa je na sredini žice.
140 Zvuk
∑
N/2−1
(s) (s)
um = Al sin(2πfl tm + ϕl ) , (10.23)
l=0
(s) (s)
sa amplitudama Al i faznim pomacima ϕl za diskretan niz ekvidistantnih
frekvencija
fl = lfTs , l = 0, 1, 2, ..., N/2 − 1 ,
fN y = fs /2 .
za
f = fl ,
lar”.13 Pošto nas ovde interesuje samo kako se menjaju amplitude harmonika,
isključite u stavkama ”Plot” sve opcije osim ”Amplitude”, da vas Origin ne bi
zatrpao ostalim podacima. Kada kliknete na ”OK” Origin će sračunati spek-
tar. Podatke za spektar će sačuvati na novom listu (Sheet) knjige iz koje su
podaci, a ubaciće i još jedan list na kojem će prikazati spektar. Preporučuje
se da se taj list ne koristi nega da se podaci za spektar nacrtaju na zasebnom
garfiku. Primer spektra je dat na slici 10.12.
4. Vaš naredni zadataj je da sa grafika spektra slobodnih oscilacija žice identi-
fikujete frekvencije na kojima se javljaju harmonici - to su uzane linije (pikovi)
u spektru. Otvorite novu knjigu i nju upišite frekvencije harmonika u funkciji
njihovog rednog broja. Zatim nacrtajte grafik ovih podataka iz kojeg treba
da se vidi da je frekvencija svakog harmonika celobrojni umnožak osnovne
frekvencije.
Slika 10.16 Logaritmovani podaci za napon (snimljen tokom slobodnih oscilacija žice).
Crvena prava prati envelopu podataka.
Slika 10.18 Osciloskopski prikaz ”Scope” u programu DataStudio. (1) i (2) su kontrole
za dva prikazana signala (”Output Voltage” i ”Voltage, ChB”); za svaki od njih se
nezavisno može nameštati skala na Y-osi, tj broj volti po ”kockici”, (1V/div za prvi i
0.02V/div za drugi signal) i ofset. (3) - služi za skalu na X-osi, tj da se zada vreme
koje odgovara jednoj kockici, ovde je to 5ms/div; skala na X-osi je ista za sve signale.
10.4 Dodaci
fazna brzina transverzalnih brzina talasa u žici pri datom k. Smenom (10.32) i
(10.33) u (10.29) nalazimo
koji glasi
uk (t, x) = Ak e−t/τ cos[k(x ± vk t) + ϕk ] , (10.35)
∑∞
∂u
(t, x) = − An e−βct sin(nπx/L)×
∂t n=1
∂u
(0, x) = v0 (x) , (10.46)
∂t
gde su u0 (x) i v0 (x) poznate funkcije koordinate x. Kako su
∞
∑
u(0, x) = An sin(nπx/L) cos(∆ϕn ) , (10.47)
n=1
∞
∑
∂u
(0, x) = − An sin(nπx/L)[βc cos(∆ϕn ) + ωn sin(∆ϕn )] .
∂t n=1
∫L
sin(nπx/L) sin(mπx/L)dx = 0 ,
0
za n ̸= m, kao i
∫L
sin2 (nπx/L)dx = L/2 .
0
150 Zvuk
Tako nalazimo da je
∫L
2
An cos(∆ϕn ) = u0 (x) sin(nπx/L)dx , (10.48)
L
0
∫L
2
An [βc cos(∆ϕn ) + ωn sin(∆ϕn )] = v0 (x) sin(nπx/L)dx ,
L
0
√
gde su c = T /µ i β = b/(2cµ) kao i kod slobodnih oscilacija. Partikularna
rešenja jednačine (10.55) se mogu naću u kompleksnom obliku
(±)
û(±)
p (t, x) = û0 e
i(kn x∓ωt)
, (10.56)
(±)
gde su û0 kompleksne amplitude ovih rešenja. Smenom (10.56) u (10.55) nalaz-
imo da su amplitude rešenja
(±) F0 (±) i(kn x∓ωt)
û0 = R̂ e , (10.57)
2µ ω,kn
gde je
(±) 1
R̂ω,kn = (10.58)
(c2 kn2 − ω 2 ) ∓ 2iβcω
(kompleksni) odziv žice na (kompleksnu) prinudnu ”jediničnu” silu ei(kn x∓ωt) .
(±)
Obzirom da polarni oblik odziva R̂ω,kn glasi
gde su
( )
1 2βcω
Rω,kn = √ , ϕω,kn = arctg ,
(c2 kn2 − ω 2 )2 + 4β 2 c2 ω 2 c2 k 2 − ω 2
(±)
partikularna rešenja ûp se mogu napisati u obliku
F0
û(±)
p (t, x) = Rω,kn ei[kn x∓(ωt−ϕω,kn )] .
2µ
Iz njih formiramo sledeće partikularno rešenje jednačine prinudnih oscilacija
(10.49) za prinudnu silu F0 sin(ωt) sin(kn x):
[ ]
up (t, x) = Re û(+)
p (t, x) − û(−)
p (t, x)
koje glasi
F0 sin(nπx/L)
up (t, x) = √ sin(ωt − ϕω,kn ) . (10.59)
µ (c2 kn2 − ω 2 )2 + 4β 2 c2 ω 2
Ovo partikularno rešenje opisuje stacionarno stanje oscilovanja - oscilovanje čije
se amplituda i faza ne menjaju tokom vremena.
Opšte rešenje prinudnih oscilacija glasi
ug (t, x) = up (t, x) + uh (t, x) , (10.60)
gde je uh (t, x) rešenje jednačine slobodnih prigušenih oscilacija (10.24) koje bi-
ramo tako da budu zadovoljeni početni uslovi (10.26) i (10.27). Obzirom da su
oscilacije žice prigušene, rešenje se pri bilo kakvim početnim uslovima nakon
prelaznog režima svodi na partikularno rešenje (10.59) koje opisuje stacionarno
oscilovanje žice, zato što nakon prelaznog režima uh (t, x) (približno) pada na
nulu.
Iz oblika stacionarnog rešenja (10.59) vidi se da amplituda prinudnih oscilacija
10.4 Dodaci 153