Oddělení fyzikálních praktik při Kabinetu výuky obecné fyziky MFF UK
PRAKTIKUM I. Úloha č. pXXIV pppppp Název:p p pTeplotní p p p p p p p p p p p p proztažnost p p p p p p p p p p p p p p p p pevných p p p p p p p p p p p p p látek pppppppppppppppppppppppppppppppppppppppppppppppppppppppppppppppppppppppppppppppppppppppppppppppppppppppppp Pracoval: p p p p p p p p p p p p p Lukáš p p p p p p p p p pVejmelka p p p p p p p p p p p p p p p p p p p p p stud. skup. p p p p p p p p p p p p p pFMUZV p p p p p p p p p p p p p(73) p p p p p p p p p p p p p p dne p p p p p p p p p p p p p p 27.3.2013 ppppppppppppppppppppppp Odevzdal dne:p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p
Možný počet bodů Práce při měření
0–5
Teoretická část
0–1
Výsledky měření
0–8
Diskuse výsledků
0–4
Závěr
0–1
Seznam použité literatury
0–1
Celkem
max. 20
Posuzoval: p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p dne p p p p p p p p p p p p p p p p p p p p p p p p
Udělený počet bodů
1
Zadání úlohy 1. Určete koeficient teplotní roztažnosti hliníku, mědi, mosazi a oceli. 2. Výsledky měření zpracujte metodou lineární regrese a graficky znázorněte. 3. V případě hliníku zkoumejte závislost prodloužení na počáteční délce.
2
Teoretický úvod měření
V technické praxi je velmi důležité znát vlastnosti používaných materiálů. Například nevhodná kombinace materiálu budované struktury by při vyšších změnách teploty mohla způsobit její potrhání, zhroucení atp. Je třeba se tedy zajímat o chování látek, jsou-li vystaveny změnám teploty. U většiny látek dochází s rostoucí teplotou ke zvětšení objemu1 . Tato informace by však nestačila, je třeba znát vlastní průběh – tedy „ jak mocÿ reaguje na změnu teploty a zda-li reaguje při různých teplotách „stále stejněÿ. Obecně lze zkoumat, jak se při konstantním tlaku se změnou teploty mění objem látky – objemová teplotní roztažnost. Budeme-li se zabývat pevnými izotropními látkami u kterých převládá jeden význačný rozměr, je šikovné řešit pouze jednodimenzionální problém, tedy zkoumat teplotní délkovou změnu – délková teplotní roztažnost a roztažnost v jednom rozměru pak rozšířit na celý objem. Toto měření bude zkoumat délkovou teplotní závislost hliníku, mědi, mosazi a oceli.
Zavedení potřebných veličin a vztahů Pro porovnávání délkových teplotních roztažností různých materiálů zavádíme součinitel délkové roztažnosti α následujícím vztahem [1] 1 ∂l α= , (1) l0 ∂t p kde l0 značí zkoumaný rozměr při teplotě 0◦ C [m] a závorka představuje délkovou teplotní změnu při konstantním tlaku. Předpokládáme-li, že závislost je lineární, tedy že součinitel α není funkcí teploty, lze po úpravě přímou integrací od t0 do t vztahu (1) získat ∆l = α l0 (t − t0 )
(2)
Přitom z linearity závislosti plyne, že l0 můžeme stanovit při obecné teplotě t0 . Připomeňme jen, že vztah byl odvozen pro situaci p = konst. a předpokladu ∂α ∂t = 0, který je však pro reálné látky splněn jen v určitém teplotním oboru. α – teplotní součinitel délkové roztažnosti [K−1 ], l0 – prvotní rozměr při počáteční teplotě t0 [m], ∆l – absolutní prodloužení l(t) − l(t0 ) [m], t – teplota odpovídající prodloužení ∆l [◦ C], t0 – teplota odpovídající l0 [◦ C]. Pro rozdíl teplot ve vztahu (2) použijme pro jednoduchost jednotek ◦ C, jako jednotku součinitele α užívejme již moderně reciproký kelvin K−1 . Můžeme si to dovolit, neboť |◦ C| = |K| Existují pryže, kde je tomu naopak. Stejně tak u vody se objem v oblasti (0 3,98) ◦ C s rostoucí teplotou zmenšuje [2]. Tuto skutečnost oznažujeme jako anomálie vody a má pro život velký význam. 1
2
2.1
Použité přístroje, měřidla, pomůcky
Vzorky - hlíníková, měďěná, mosazná a ocelová trubice s integrovanou konektoráží, dilatometr s indikátorovými hodinami, nádržka s vodou a teploměrem, termostatický ohřívač s čerpadlem, příslušenství (hadice – oběhového okruhu, vypouštěcí, napouštěcí). Tabulka 1: Použité měřící přístroje a jejich mezní chyby měření. Měřidlo Veličina[jednotka] Mezní chyba Pozn. Dilatometr l0 [m] považujme za přesné Indikátorové hodiny ∆l[m] 10−5 dílek ◦ Teploměr t, t0 [ C] 0,5 dílek
2.2
Důležité hodnoty, konstanty, vlastnosti
Důležité hodnoty pro výpočet nebo látkové konstanty pro porovnání výsledků. 0 = 0,0231 · 10−3 K−1 [3] Teplotní součinitel délkové roztažnosti hliníku při 20 ◦ C: αAl 0 Teplotní součinitel délkové roztažnosti mědi při 20 ◦ C: αCu = 0,017 · 10−3 K−1 [3] 0 Teplotní součinitel délkové roztažnosti mosazi při 20 ◦ C: αmosaz = 0,019 · 10−3 K−1 [3] 0 Teplotní součinitel délkové roztažnosti oceli při 20 ◦ C: αocel = (0,011 0,013) · 10−3 K−1 [3]
2.3
Popis postupu vlastního měření
Příprava, zahájení měření Zkoumaný materiál upevníme do aretačního držáku dilatometru deklarující udávanou počáteční délku tyče. Rychlospojkami trubici zapojíme do okruhu teplonosného média. Nádržku naplníme vodou, kterou necháme zahřát na počáteční teplotu t0 . Namontujeme indikátorové hodiny. Po ustálení teploty a vhodné relaxační době promývání trubice nastavíme ručičku indikátorových hodin na nulu. Pro zahájení měření určíme termostatu, do jaké maximální teploty má ohřívat a průběžně sledujeme protažení ∆l trubice při teplotách t. Měření v závislosti na počáteční délce Cílem je určit prodloužení ∆l při konstantní změně teploty 20 ◦ C pro prvotní délky tyče 200 mm, 400 mm a 600 mm. Upevníme proto trubici do jedné z pozic odpovídající prvotní délce a realizujeme ohřev o 20 ◦ C. Pro efektivnější zpracování zapisujeme prodloužení průběžně pro různé teploty před dosažením daného teplotního rozdílu. Proceduru opakujeme pro všechny tři délky. Vzhledem k dalšímu úkolu můžeme v případě délky 600 mm po dosažení teplotního rozdílu 20 ◦ C pokračovat až k 60 ◦ C. Měření v závislosti na změně teploty Trubici ze zkoumaného materiálu ustavíme do pozice odpovídající prvotní délce 600 mm. Realizujeme ohřev trubice v oboru teplot (20 60) ◦ C a zaznamenáváme odpovídající si dvojice prodloužení ∆l a teplot t. Po dosažení maximální teploty nahradíme teplou vodu v nádržce vodou studenou. V případě, že z důvodu nahřáté spirály topného tělesa nedosáhneme potřebné nejnižší teploty, vyměníme vodu znovu. Měření provedeme pro všechny udané materiály.
3
3
Výsledky měření
3.1
Laboratorní podmínky
Teplota v laboratoři: 23,4 ◦ C Atmosférický tlak: 984,2 hPa Vlhkost vzduchu: 20,9 %
3.2
Způsob zpracování dat
Určení αAl0 ze závislosti na počáteční délce Vyjdeme ze vztahu (2). Roznásobením a zavedením koeficientů A a E získáme
∆l = αl0 t − αl0 t0 ∆l = A · t − E.
(3)
Vykreslením závislostí ∆l = f (t) pro l01|2|3 = (200 | 400 | 600) mm a proložením lineárních regresních přímek získáme regresní koeficienty Ai . S jejich pomocí vypočítáme odpovídající absolutní prodloužení ∆li při konstantním teplotním rozdílu t − t0 = ∆t pro jednotlivé počátení délky l0i jako ∆li = α l0i ∆t = Ai ∆t,
(4)
neboť součin α l0 představuje z rovnice (3) právě regresní koeficient A. Nyní konečně vykreslíme zkoumanou závislost ∆l = f (l0 ). Závislostí proložíme lineární regresní křivku. Bude-li rozdíl teplot t − t0 = ∆t konstantní, lze z rovnice (2) psát ∆l = αl0 ∆t = B · l0 .
(5)
Ze směrnice B proložené přímky vypočítáme součinitel αAl0 jako αAl0 =
B . ∆t
(6)
Určení α ze závislosti na teplotě Výpočet sleduje první část předchozího rozboru. Analogicky s rovnicí (3) můžeme psát teoretickou závislost prodloužení ∆l na teplotě t
∆l = C · t − F.
(7)
Opět tedy vykreslíme závislost ∆l = f (t), proložíme lineární regresní křivku a z regresního koeficientu C vypočítáme příslušný součinitel α vztahem α=
C . l0
(8)
Koeficienty C zjistíme pro data z měření jednotlivých vzorků a na základě tohoto vztahu vypočítáme hledané teplotní součinitele délkové roztažnosti α.
4
Určení nejistot měření Nejistotu jednotlivých součinitelů vypočítáme na základě relativních chyb regresních koeficientů, započítáme i chyby určení veličin vystupujících ve výpočetních vztazích. Prvotní délky l0 budume považovat za přesné, v případě určení součinitele αAl0 hliníku ze závislosti na počáteční délce bude započítána chyba měření teploměrem. Do výsledné nejistoty přispějí chyby, které nelze zahrnout pomocí zákona kvadratického hromadění chyb nepřímého měření, jako chyby určení závislých veličin, které mohou být různé a tedy se regresně nevyruší, dále odlišný odvod tepla z trubice do okolí při daných teplotách a komplexně možnost nekorespodence teploty v místě měření v nádržce a teploty materiálu vzorku při různých okamžicích měření. Příspěvek těchto chyb do výsledné nejistoty odhadujme maximálně na 0,05, P ≈ 1.
3.3
Naměřené hodnoty
Naměřené hodnoty potřebné k výpočtům hledaných veličin jsou zachyceny v tabulkách 2 a 3.
Tabulka 2: Naměřená data pro výpočet Číslo l01 = 200 mm ◦ měření t[ C] ∆l[10−5 m] 1. 22,5 0 2. 23,0 1 3. 25,0 2 4. 27,0 3 5. 29,5 4 6. 30,2 5 7. 34,5 6 8. 36,5 7 9. 38,5 8 10. 41,5 9 11. 42,5 9,5 −5 ◦ −1 A[10 m C ] 0,452 sA [10−5 m ◦ C−1 ] 0,017 −5 ∆l [10 m] 8,9
součinitele αAl0 ze závislosti ∆l = f (l0 ). l02 = 400 mm l03 = 600 mm ◦ −5 ◦ t[ C] ∆l[10 m] t[ C] ∆l[10−5 m] 20,5 0 20,0 0 22,0 1 22,5 3 25,0 4 24,5 6 28,5 6 26,0 8 31,0 9 27,0 9 33,0 11 28,5 11 35,1 13 30,2 15 37,0 15 33,0 17 40,0 17 35,0 20 40,5 17,5 37,0 23 42,0 19 42,0 29 0,894 1,332 0,016 0,017 17,9 26,6
Změřené / odečtené hodnoty / hodnoty zadání Rozdíl teplot pro zkoumání závislosti na počáteční délce: ∆t = 20◦ C, Počáteční délky při zkoumání závislosti na počáteční délce: l01|2|3 = (400 |600 |800) mm.
3.4
Zpracování dat a číselné výsledky
Určení součinitele αAl0 v závislosti na počáteční délce Vykreslíme závislosti ∆l = f (t) pro všechny tři počátení délky l0 z dat tabulky 2. Grafy závislostí jsou proloženy lineární přímky, viz graf 1. Regresní koeficienty, vypočítané programem QtiPlot, jsou zapsány v tabulce 2.
5
Tabulka 3: Naměřená data pro Číslo Hliník měření t[◦ C] ∆l[10−5 m] 1. 20,0 0 2. 22,5 3 3. 24,5 6 4. 26,0 8 5. 27,5 10 6. 29,0 12 7. 30,2 15 8. 35,0 20 9. 37,0 23 10. 39,5 26 11. 42,0 29 12. 44,0 32 13. 47,5 38 14. 51,5 42 15. 54,0 45 16. 57,0 49 17. 60,0 52,5
výpočet součinitelů délkové roztažnosti α Měď Mosaz t[◦ C] ∆l[10−5 m] t[◦ C] ∆l[10−5 m] 20,0 0 20,5 0 21,0 3 21,5 3 23,5 5 23,0 4,5 26,0 8 26,0 8 29,0 11 29,0 11 32,0 14 32,0 15 35,0 17 35,0 18 38,0 20 38,5 22 41,0 23 41,0 24,5 44,0 26,5 43,0 27 47,0 29 46,0 30 50,0 32 49,0 33 53,0 35 52,0 36 55,0 37 55,0 40 57,0 39 57,0 42 59,0 41 59,0 44 60,0 42 60,0 45
závislosti ∆l = f (t). Ocel t[◦ C] ∆l[10−5 m] 20,0 0 22,0 1,5 23,0 2 25,0 3,5 27,0 5 30,0 7 33,0 9 37,0 12,5 40,0 14 43,5 17 46,0 19 48,0 20 50,0 21,5 53,0 23 55,0 25 57,5 27 60,0 28
Graf 1: Závislost ∆l = f (t) hliníku pro jednotlivé počáteční délky l0 .
l 01 = 200 mm l 02 = 400 mm l 03 = 600 mm lineární fit
∆l[10−5 m]
30
20
10
0 20
25
30
35
40
45
t[◦ C]
Ze získaných koeficinetů Ai vypočítáme pomocí vztahu (4) prodloužení ∆li odpovídající udané teplotní změně ∆t = 20 ◦ C pro jednotlivé počáteční délky l0 . Tato prodloužení jsou vypočítána v tabulce 2. Nyní konečně vykreslíme zkoumanou závislost ∆l = f (l0 ) a grafem proložíme přímku, z jejíž směrnice B vypočítáme součinitel αAl0 podle vztahu (6). Graf 2 ukazuje toto proložení a současně jsou pro zajímavost světlejší barvou zaneseny výsledky, když by se prodloužení ∆l určila přímo jako prodloužení ∆l odpovídající konstantní změně teploty ∆t.
6
Velikost směrnice B, fakticky představující relativní prodloužení, z programu QtiPlot: B = (4,43 ± 0,03) · 10−4 Ze vztahu (6) vypočítáme hledaný součinitel αAl0 =
B = 0,022 · 10−3 K−1 . ∆l
Chybu vypočítáme z chyby směrnice B a chyby určení ∆t, viz tabulka 1. q ραAl0 = (3ρB )2 + ρ2∆t = 0,054 Absolutní mezní chyba je pak εαAl0 = 0,0012 · 10−3 K−1 Graf 2: Závislost ∆l = f (l0 ) hliníku při ∆t = 20 ◦ C. 30
Z regresního výpočtu Lineárni fit Z přímého výpočtu Lineární fit
25
∆l[10−5 m]
20
15
10
5
0 0
200
400
600
l0 [mm]
Výpočet součinitelů α ze závislostí ∆l = f (t) Naměřená data pro výpočet součinitelů pro jednotlivé vzorky jsou v tabulce 3. Z těchto dat vykreslíme závislosti ∆l = f (t).
7
Graf 3: Graf závislostí ∆l = f (t) pro hliník, měď, mosaz a ocel. 60
Hliník Měď Mosaz Ocel
55 50 45 ∆l[10−5 m]
40 35 30 25 20 15 10 5 20
30
40
50
60
t[◦ C]
Jednotlivé závislosti proložíme lineárními regresními přímkami, viz graf 3. Směrnice regresních přímek, jejich chyby a výpočet součinitelů délkové teplotní roztažnosti s jejich nejistotami určuje tabulka 4. Výpočet součinitelů probíhá dle vztahu (8). K nejistotě měření kromě odchylky směrnice přispívá odhadovaná chyba ostatních vlivů, viz část způsob zpracování dat. Tabulka 4: Výpočet velikostí a nejistot součinitelů délkové teplotní roztažnosti α ze směrnic C. Vzorek C [10−2 m · ◦ C−1 ] sC [10−5 m · ◦ C−1 ] α[10−6 K−1 ] ρC [1] ρ[1] ε[10−6 K−1 ] Hliník 1,32 7,27 22,1 0,0055 0,043 0,96 Měď 1,02 8,62 16,9 0,0085 0,047 0,80 Mosaz 1,11 9,83 18,4 0,0089 0,048 0,89 Ocel 0,71 5,34 11,9 0,0075 0,046 0,54
V tabulce 4 pro daný vzorek značí C směrnici proložené přímky, sC směrodatnou odchylku určení směrnice, α součinitel délkové teplotní roztažnosti, ρC relativní směrodatnou chybu směrnice, ρ celkovou mezní chybu a ε absolutní mezní nejistotu.
3.5
Číselné výsledky měření
Součinitel délkové teplotní roztažnosti hliníku určený ze závislosti na počáteční délce: αAl0 = (22 ± 1) · 10−6 K−1 ,
P ≈ 1.
Součinitelé délkové teplotní roztažnosti zkoumaných vzorků určené ze závislosti na teplotě: αAl = (22,1 ± 1,0) · 10−6 K−1 ,
P ≈ 1,
αCu = (16,9 ± 0,8) · 10−6 K−1 ,
P ≈ 1,
αmosaz = (18,4 ± 0,9) · 10−6 K−1 , αocel = (11,9 ± 0,5) · 10−6 K−1 ,
8
P ≈ 1, P ≈ 1.
3.6
Grafické výsledky měření
Výsledky měření jsou zaznamenány v grafech v průběhu části protokolu zpracování dat. Graf 2 zachycuje závislosti ∆l = f (l0 ), graf 3 pak ∆l = f (t).
4
Diskuze výsledků
Pro zajímavost uveďme přímé porovnání naměřených součinitelů s tabelovanými hodnotami [3], viz tabulka 5. Tabulka 5: Porovnání naměřených hodnot α s tabelovanými [3] . Materiál Naměřeno α[10−6 K−1 ] Tabelováno α0 [10−6 K−1 ] Hliník 22,1 23,1 Měď 16,9 17 Mosaz 18,4 19 Ocel 11,9 11 13
Z porovnání je zřejmé, že materiály vzorků co se týče délkové teplotní roztažnosti si přibližně odpovídají. Měření považuji vzhledem k řádové shodě výsledků za úspěšné. V měření může být obsažena systematická chyba způsobená následujícím. Dilatometr má uchycení zkoumané trubice přizpůsobené tak, aby při teplotě 20 ◦ C byla její aktivní částí délka 600 mm. I kdyby se podařilo dosáhnout přesné teploty 20 ◦ C teplonosné lázně v nádržce, může být teplota tyče nižší a tedy i počáteční rozměr jiný. Chybu může vnášet i nerovnoměrnost roztažnosti v části trubice, kde je bočně navařena příruba pro napojení okruhové hadice. Grafy ukazují, že ve zkoumaných oblastech jsou závislosti lineární. Graf 2 navíc dokazuje, že přímé určení prodloužení ∆l při ohřevu o 20 ◦ C, by bylo chybné – chyby určení odečtení hodnot by měly velký vliv. Výhodnost regresního nalezení závislostí a následujícího výpočtu příslušných absolutních prodloužení je potvrzena výborným „přimknutíÿ fitované přímky v grafu 2. Součinitel teplotní roztažnosti oceli je velice podobný součiniteli betonu (αbeton = 10 · 10−6 K−1 ). Díky tomu je železobeton z tohoto hlediska poměrně kompaktní. [4]
5
Závěr
Závislosti absolutního prodloužení na teplotě i počáteční délce jsou ve zkoumané oblasti teplot a délek dle očekávání lineární. Součinitel délkové teplotní roztažnosti hliníku určený ze závislosti na počáteční délce: αAl0 = (22 ± 1) · 10−6 K−1 ,
P ≈ 1.
Součinitelé délkové teplotní roztažnosti zkoumaných vzorků určené ze závislosti na teplotě: αAl = (22,1 ± 1,0) · 10−6 K−1 ,
P ≈ 1,
αCu = (16,9 ± 0,8) · 10−6 K−1 ,
P ≈ 1,
−6
αmosaz = (18,4 ± 0,9) · 10
−6
αocel = (11,9 ± 0,5) · 10
−1
K
−1
K
,
,
P ≈ 1, P ≈ 1.
Největší součinitel délkové teplotní roztažnosti z měřených materiálů má hlíník, dále následuje mosaz a měď, poslední je pak ocel. 9
Seznam použité literatury [1] H. Valentová: Fyzikální praktikum, studijní text, MFF UK. (29.3.2013). http://physics.mff.cuni.cz/vyuka/zfp/txt_124.pdf [2] ONLINE: Thermal expansion coefficients, wikipedia.org. (30.3.2013) http://cs.wikipedia.org/wiki/Teplotni_roztaznost [3] ONLINE: Thermal expansion coefficients, wikipedia.org. (30.3.2013) http://en.wikipedia.org/wiki/Thermal_expansion [4] ONLINE : Teplotní součinitel délkové roztažnosti betonu. (31.3.2013). http://www.ebeton.cz/pojmy/teplotni-soucinitel-delkove-roztaznosti
10
