Piezoelektricita¶
Následující poznámky jsou sestaveny podle [1]. Jde o první kapitoly věnující se
elektrostatice a piezoelektricitě. Ukázka piezoelektrické úlohy z lomové mechaniky řešené
pomocí metody konečných prvků systémem FEniCS je k přečtení zde. Zdrojové
soubory k rozpracovanému článku jsou ke stažení zde
.
Elektrostatika vodičů¶
Základem elektrostatiky jsou Maxwellovy rovnice vakua, které však platí na mikroúrovni,
kde \(\boldsymbol{h}\) je mikroskopické magnetické pole a \(\boldsymbol{e}\) je mikroskopické elektrické pole. Vakuum není kupodivu ve statice předmětem zájmu. Tím je vodič, u kterého však vzniká problém, jelikož už z podstaty jeho existence, není v něm žádné elektrické pole. Nenulovost elektrického pole způsobuje proud a součástí šíření proudu vodičem je nevratná změna energie (její disipace), která se u stacionární úlohy nekoná. Proto ta nulovost elektrického pole ve vodiči.
Odtud také plyne, že náboj, se nachází pouze na povrchu vodiče, jinak by totiž způsoboval elektrické pole uvnitř vodiče. Náboje jsou tedy na povrchu vodiče rozloženy tak, aby byly navzájem v rovnováze. Tuto rovnováhu však jde poznat jen z elektrického pole, které tedy u vodiče není k dispozici, je však ve vakuu, které ho obklopuje a má s vodičem společný povrch. Změna elektrického pole \(\boldsymbol{e}\) ve vakuu je plynulá a jeho střední hodnota \(\boldsymbol{E}\), která odpovídá jeho hodnotě v kontinuu, je s jeho hodnotou téměř shodná.
Předpokládáme, že průměrná hodnota magnetického pole \(\overline{\boldsymbol{h}}\) a jeho změny v čase je nulová, pak pro elektromagnetické pole \(\boldsymbol{E}\) podle (1) a (2) platí
Očividně jde o "klasické" rovnice rovnováhy sil a momentů, pokud by si jeden domyslel místo pole \(\boldsymbol{E}\) třeba napětí \(\boldsymbol{\sigma}\) v klasickém mechanickém smyslu. Elegantní způsob řešení rovnic (3) a (4) nabízí zavedení elektrického potenciálu \(\phi\) tak, že
Pak se dosazením (5) do (3) dostane Laplaceova rovnice
Z rovnice rovnováhy (4) se dostanou okrajové podmínky k rovnici (6). Detaily odvození okrajových podmínek lze nalézt v [1]. Z podmínky (4) zejména plyne, že tečné hodnoty \(\boldsymbol{E}\) k povrchu vodiče musí být nulové
Jinými slovy, směr elektrostatického pole musí být normálový k povrchu vodiče a protože pro potenciál \(\phi\) platí (5), musí být elektrický potenciál \(\phi\) na povrchu vodiče konstantní, nebo-li musí tvořit eqipotencálu elektrostatického pole. To také plyne z vlastnosti harmonických funkcí, pro které platí, že svých extrémních hodnot mohou dosáhnout pouze na hranici oblasti. Důležitý, pro stanovení okrajových podmínek, je vztah mezi průměrnou hodnotou hustoty náboje \(\rho\) a elektrickým polem \(\boldsymbol{E}\)
Aplikací posledního výrazu na nekonečně malý objem mezi jednotkovým povrchem vodiče a vakua a z podmínky, že \(\boldsymbol{E}\) musí mít normálový směr k povrchu vodiče plyne
Integrací předchozího vztahu přes plochu vodiče se dostane celkový náboj na povrchu vodiče
Totální energie \(\mathcal{U}\) elektrostatického pole nabitého vodiče se může psát
kde se integruje přes oblast vně vodiče. Druhý integrál je nulový a první se může převést na plošný integrál za předpokladu, že elektrické pole v nekonečnu mizí. Protože je elektrický potenciál podél hranice konstantní, musí být v nekonečnu nulový. Pak tedy platí
kde se integruje se přes plochu \(i\) vodičů v elektrickém poli \(\boldsymbol{E}\) a kde se použil vztah (10). Mezi náboji a odpovídajícími elektrickými potenciály je lineární vztah
kde \(C_{ij}\) mají dimenzi délky a závisí na tvaru a relativní poloze vodičů. Parametry \(C_{ii}\) jsou koeficienty kapacity, parametry \(C_{ij}\) jsou koeficienty elektrostatické indukce. Dá se ukázat, viz [1], že
Lineární závislost (13) mezi \(e_i\) a \(\phi_i\) umožňuje vyjádřit potenciál a náboj následovně
nebo naopak energii \(\mathcal{U}\) jako
Hodně zajímavý je případ vodiče bez náboje v rovnoměrném vnějším elektrickém poli \(\boldsymbol{\mathcal{E}}\), které si lze představit jako pole generované od nábojů v nekonečnu. Přestože je na povrchu vodiče nulový náboj, vnější elektrické pole v něm vytváří moment elektrického dipólu \(\boldsymbol{\mathcal{P}}\) (viz analogie s bodovým dipólem, kdy dva bodové náboje opačných znamének jsou k sobě limitně blízko, tj. jejich náboje se eliminují, ale zůstává elektrické pole dvou nábojů v jediném bodě). Elektrický potenciál \(\phi\) takového dipólu \(\boldsymbol{\mathcal{P}}\) je při vzdálenosti \(\boldsymbol{r}\) od něj roven
Pokud se dále vezme fakt, že
pak se pro potenciální energii nenabitého vodiče v rovnoměrném elektrickém poli může psát
Podobně jako v případě elektrického potenciálu \(\phi\) a náboje \(e\) je vztah mezi momentem dipólu \(\boldsymbol{\mathcal{P}}\) a elektrickým polem \(\boldsymbol{\mathcal{E}}\) lineární, tj.
kde \(V\alpha_{ik}\) jsou komponenty tzv. tenzoru polarizace oblasti. Tento tenzor je symetrický, tj. platí
Pro energii (19) pak platí
V elektrickém poli působí na vodič síly a momenty. Tyto síly opět vychází z Maxwellových rovnic, které vyjadřují napětí tzv. Maxwellův tenzor napětí
Stejně jako v pružnosti, výraz \(\sigma_{ik}n_k\) je síla \(\boldsymbol{F}_s\) na jednotku plochu povrchu vodiče. Protože však elektrické pole na povrchu vodiče má pouze směr normály \(\boldsymbol{n}\) k povrchu, platí pro sílu \(\boldsymbol{F}_s\) následující
nebo také zavedením hustoty povrchového náboje \(\sigma\)
Elegantní způsob vyjádření síly působící v nějakém směru \(q\) vychází z mechaniky (viz Castiglianova věta) jako derivace potenciální energie podle tohoto směru \(q\), tj.
kde index \(e\) značí, že se derivace počítá při konstantním náboji \(e\). Podobně se síla \(F_q\) může vyjádřit jako derivace podle \(q\) při konstantním potenciálu \(\phi\)
kde \(\tilde{\mathcal{U}}\) je vlastně obdoba doplňkové energie v mechanice, tj.
K vyjádření vztahů (26) a (27) za zavádí také diferenciální vztahy, přičemž první je funkcí náboje \(e\) a souřadnice \(q\) a druhý je funkcí elektrického potenciálu \(\phi\) a souřadnice \(q\)
Pro vodiče bez náboje v elektrickém poli, které se v okolí vodiče mění jen zanedbatelně, se síla působící na vodič může vyjádřit jako gradient energie (22) následovně
Kromě síly působí na nenabitý vodič i moment \(\boldsymbol{K}\), pro který platí
Odvození tohoto vztahu lze nalézt v [1]. Síly (26) a (27) a momenty (31) působící na vodič způsobují jeho deformaci, tzv. elektrostrikci, která může hrát významnou roli, ale musí být svázána s rovnicemi pružnosti.
Elektrostatika dielektrika¶
Doplňkem k vodičům jsou dielektrika, které však naopak žádný proud nevedou. Z tohoto důvodu v nich není elektrické pole nulové
kde \(\overline{\rho}\) je výše zmíněná průměrná hodnota hustoty náboje. Nicméně, rovnici
mají dielektrika a vodiče společnou. Nejdříve budeme předpokládat, že dielektrikum není přípojeno na žádný vnější zdroj elektrického náboje. Pak musí platit
kde se integruje přes objem dielektrika. Jedinou možností řešení této integrální rovnice je zavedení vektoru \(\boldsymbol{P}\) pro který platí
zatímco vně dielektrika je \(\boldsymbol{P}=0\). Fyzikální význam naznačí dosazení (35) do (34)
Vektor \(\boldsymbol{P}\) se nazývá vektorem polarizace dielektrika a necharakterizuje jen rozložení hustoty náboje \(\rho\) uvnitř dielektrika, ale také hustotu náboje \(\sigma\) na jeho povrchu
kde \(P_n\) je složka vektoru \(\boldsymbol{P}\) rovnoběžná s vnější normálou k povrchu dielektrika. Dosazením (35) do (32) se dostane
kde
je tzv. elektrická indukce. Rovnice (38) jde zobecnit na případ, kdy náboj \(\rho_{ex}\) není součásti dielektrika, ale do něj dodáván z nějakého vnějšího zdroje, pak
Okrajové podmínky plynou jak z rovnice (38), tak z rovnice (33). Nejdříve tedy posledně jmenovaná rovnice. Z té plyne, že na rozhraní dvou dielektrik musí být spojitost tangenciálních směrů \(\boldsymbol{E}_t\) podél jejich rozhraní
Okrajová podmínka vyplývající z (38) vyžaduje spojitost normálové složky indukce \(D_n\) podél rozhraní mezi dvěma dielektriky
V případě rozhraní dielektrika s vodičem platí
kde \(\sigma_{ex}\) je hustota náboje na povrchu vodiče. Aby soustava rovnic (33) a (38) byla úplná, musí mít elektrické pole \(\boldsymbol{E}\) a indukce \(\boldsymbol{D}\) mezi sebou vztah. Tento vztah je naštěstí lineární a v případě izotropního dielektrika má velice sympatický tvar
kde \(\varepsilon\) je tzv. permitivita dielektrika, která je stavovou veličinou (závisí na stavu dielektrika z pohledu termodynamiky). Podobně jako indukce \(\boldsymbol{D}\), také polarizace \(\boldsymbol{P}\) má svou závislost, lineární, samozřejmě
kde \(\kappa\) je tzv. koeficient polarizace nebo také faktor citlivosti dielektrika. Vzhledem k okolnostem těchto lineárních závislostí mezi \(\boldsymbol{D}\),:math:boldsymbol{P} a \(\boldsymbol{E}\) pro okrajové podmínky (41) a (42) platí
Jestliže se dále vezme v úvahu vztah (5), potom je rovnice rovnováhy (33) automaticky splněna a rovnice rovnováhy (38) má tvar
s okrajovými podmínkami (46) ve tvaru
Další zajímavý případ, který ovlivňuje zavedení homogenní permitivity \(\varepsilon=konst.\), je případ vodiče v tomto homogenním dielektriku. Jak u vodiče, tak i dielektrika je splněna rovnice rovnováhy
s okrajovou podmínkou, že \(\phi\) je konstantní na povrchu vodiče, ale s tím rozdílem, že místo normálové složky elektrického pole napětí \(E_n\) je nulové normálová složka indukce
kde je opět \(\sigma\) hustota náboje na povrchu vodiče. Protože ve vodiči je nulové elektrické pole, ovlivňuje jeho celkovou energii jen okolní prostředí, které je ale nezávislé na jeho termodynamickém stavu, takže nezvyšuje jeho entropii. Naopak, z důvodu průchodu elektrického pole dielektrikem, má elektrické pole významný vliv na jeho termodynamické vlastnosti. Mechanickou práci, kterou musí vykonat elektrické pole přenesením nekonečně malého náboje \(\delta e\) z nekonečna na povrch vodiče z důvodu rozdílů potenciálu \(\phi\) mezi nimi, lze vyjádřit následovně
Jestliže dielektrikum obklopuje vodič, normálová komponenta jeho indukce \(D_n\) ve směru normály \(\boldsymbol{n}\) k povrchu vodiče (tj. směrem vně z dielektrika a tedy dovnitř vodiče) udává podle (50) hustotu povrchového náboje vodiče a tedy také celkovou hodnotu náboje na povrchu vodiče
Protože je elektrický potenciál na povrchu vodiče konstantní, platí
kde poslední integrál je počítán přes celý objem vně vodiče. Pro variaci \(\delta\boldsymbol{D}\) musí stejně jako pro \(\boldsymbol{D}\) platit rovnice rovnováhy (38), tj.
Takže
a (53) bude mít tvar
Jen pro připomenutí, objemový integrál je počítaný přes celý objem vně vodiče, tedy i přes vakuum, které nevyplňuje dielektrikum. Práce vykonaná na teplotně izolovaném tělese je vlastně změna energie tělesa za konstantní entropie. Proto stačí ke vztahu změny energetické rovnováhy \(\delta\mathcal{U}\) přidat předchozí vztah
kde \(T\) je teplota a \(\mathcal{S}\) je entropie. Podobně to platí i pro celkovou volnou energii
kde po dosazení (56) do \(\delta\mathcal{F}\) se dostane
Podobně lze dostat termodynamické veličiny vztažené na jednotku objemu. Jestliže \(U\), \(S\) a \(\rho\) jsou vnitřní energie, entropie a hmotnost jednotkového objemu, tak pro změnu energie \(\mathrm{d}U\) po přidání elektrického pole a dielektrika platí
kde \(\zeta\) je chemický potenciál. Pro doplňkovou energii \(F=U-TS\) platí podobně
Z těchto vztahů se dá vydedukovat další alternativní výrazy pro elektrické pole při konstantní entropii a hmotnosti, příp. konstantní teplotě a hmotnosti (přičemž druhý případ je praktičtější)
Podobně jako v pružnosti i zde je vhodné kromě energií \(U\) a \(F\) zavést také jejich doplňky
Jejich diferenciály pak jsou
a
Odtud plyne
Elektrické pole \(\boldsymbol{E}\) bez změny teploty a hustotu způsobuje deformaci dielektrika. Obecně deformované těleso není izotropní jak z pohledu mechanických, tak i dielektrických vlastností. Deformaci v bodě tělesa popisuje tenzor deformace \(u_{ij}\), který se, jak je známo z pružnosti, může napsat ve tvaru
kde \(\boldsymbol{u}\) je vektor posunutí. Pro další úvahy je vhodné si tvar (67) upravit. Podobně jako v pružnosti prutů, se bude v dalším uvažovat deformace velice tenké vrstvy tloušťky \(h\), jejíž plochy mají orientaci danou normálou \(\boldsymbol{n}\). Jedna z těchto ploch je podrobena virtuálnímu nekonečně malému posunutí \(\boldsymbol{\xi}\) nezávislému na normále \(\boldsymbol{n}\). Takže zatímco objemová změna vrstvy \(h\) odpovídající diagonálním prvkům tenzoru deformace \(u_{ii}\) je vlastně příspěvek dané složky posuvu \(\xi_i\) do odpovídající složky normály \(n_i\), tj. \(\xi_in_i\), kde se však nesčítá přes index \(i\). Podobně, pro vzájemné skosení (natočení) obou ploch vrstvy \(h\) je příspěvek složek vektoru \(\xi_i\) a \(\xi_k\) do odpovídajících, ale kolmých složek normály \(n_k\) a \(n_i\), tj. \(\xi_in_k\). Ve výsledku se tyto příspěvky deformace vrstvy \(h\) mohou zapsat ve tvaru symetrického tenzoru deformace
Práce síly \(\sigma_{ik}n_k\) na posuvu \(\boldsymbol{\xi}\) vrstvy tloušťky \(h\) je
kde pro \(\delta\tilde{F}\) podle (65) při konstantní teplotě, bez změny hustoty dielektrika a zahrnutí vlivu deformace platí
Elektrické pole se mění nejen vlivem změny tloušťky \(h\) deformované vrstvy, ale také jejím natočením o úhel \(\delta\boldsymbol{\phi}\)
Natočení je samozřejmě zakuklený vektorový součin, takže pro \(\delta\boldsymbol{\phi}\) platí
Dosazením tohoto vztahu do (71) se dostane
a pro výraz \(-\boldsymbol{D}\cdot\delta\boldsymbol{E}/4\pi\) platí
Literatura