abclinuxu.cz AbcLinuxu.cz itbiz.cz ITBiz.cz HDmag.cz HDmag.cz abcprace.cz AbcPráce.cz
AbcLinuxu hledá autory!
Inzerujte na AbcPráce.cz od 950 Kč
Rozšířené hledání
×

dnes 15:44 | Komunita

Nadace The Document Foundation (TDF) zastřešující vývoj svobodného kancelářského balíku LibreOffice dnes slaví 6 let od svého oficiálního vzniku. Nadace byla představena 28. září 2010. Formálně ale byla založena až 17. února 2012. Poslední lednový den byl vydán LibreOffice 6.0. Dle zveřejněných statistik byl za dva týdny stažen již cca milionkrát.

Ladislav Hagara | Komentářů: 0
dnes 04:44 | Bezpečnostní upozornění

CSIRT.CZ upozorňuje, že byla vydána nová verze 1.2.3 svobodného routovacího démona Quagga (Wikipedie) přinášející několik bezpečnostních záplat. Při nejhorší variantě může dojít až k ovládnutí běžícího procesu, mezi dalšími možnostmi je únik informací z běžícího procesu nebo odepření služby DoS. Konkrétní zranitelnosti mají následující ID CVE-2018-5378, CVE-2018-5379, CVE-2018-5380 a CVE-2018-5381.

Ladislav Hagara | Komentářů: 0
dnes 00:22 | Pozvánky

V sobota 7. dubna proběhne v Brně na FIT VUT nekomerční konference Security Session '18 věnovaná novinkám a aktuálním problémům v oblasti bezpečnosti mobilních a informačních technologií. Organizátoři vyhlásili CFP. Návrhy přednášek a workshopů lze zaslat do 4. března. Spuštěna byla registrace.

Ladislav Hagara | Komentářů: 0
včera 23:55 | Nová verze

Byla vydána verze 1.10 programovacího jazyka Go (Wikipedie). Přehled novinek v poznámkách k vydání. Zdůraznit lze vylepšené cachování buildů.

Ladislav Hagara | Komentářů: 0
15.2. 22:55 | Komunita

V sobotu 21. dubna proběhne v Bratislavě na Fakultě informatiky a informačních technologií STU konference Bratislava OpenCamp 2018. Organizátoři vyhlásili CFP. Návrhy přednášek a workshopů lze zaslat do 10. března. Spuštěna byla registrace.

Ladislav Hagara | Komentářů: 0
15.2. 21:22 | Nová verze

Byla vydána verze 1.24 programovacího jazyka Rust (Wikipedie). Z novinek vývojáři zdůrazňují nový nástroj rustfmt pro reformátování zdrojových kódů dle požadovaného programovacího stylu a inkrementální překlad. Více v poznámkách k vydání. Vyzkoušet Rust lze například na stránce Rust By Example.

Ladislav Hagara | Komentářů: 5
15.2. 13:00 | Komunita

Desktopovému týmu Ubuntu by se hodilo více informací o desktopech uživatelů Ubuntu. V diskusním listu ubuntu-devel byl k diskusi předložen návrh řešení a seznam odesílaných informací. Ve výchozím stavu by mělo být odesílaní informací povoleno. Uživatel by měl mít možnost odesílaní kdykoli jednoduše zakázat [reddit].

Ladislav Hagara | Komentářů: 20
15.2. 05:55 | Komunita

V úterý 13. února vypršela platnost posledního patentu US 7,334,248 souvisejícího s MPEG-2 [Hacker News].

Ladislav Hagara | Komentářů: 22
15.2. 04:44 | Komunita

Organizace Electronic Frontier Foundation (EFF) na svém blogu oznámila, že otevřená certifikační autorita Let's Encrypt dosáhla dalšího milníku. Počet platných certifikátů překročil 50 milionů, viz statistiky. Cekově bylo vydáno již přes 217 milionů certifikátů.

Ladislav Hagara | Komentářů: 0
15.2. 00:22 | Pozvánky

Konference SUSE Expert Days 2018 proběhnou mimo jiné také 1. března v Bratislavě a 15. března v Praze.

Ladislav Hagara | Komentářů: 0
Který webový vyhledávač používáte nejčastěji?
 (3%)
 (28%)
 (62%)
 (2%)
 (3%)
 (1%)
 (1%)
 (1%)
Celkem 338 hlasů
 Komentářů: 34, poslední 14.2. 18:44
    Rozcestník

    Agros2D – výpočet kapacity jiskřiště a indukčnosti dvojlinky

    19. 5. 2010 | František Mach | Návody | 5644×

    V tomto díle se více zaměříme na samotné řešení (processing) problémů a ukážeme si základy práce s příkazy. Budeme řešit dva příklady, a to výpočet kapacity jiskřiště a výpočet indukčnosti dvojvodičového kabelu.

    Obsah

    Úvod

    link

    Jak již bylo řečeno v minulém díle seriálu, Agros2D využívá pro řešení parciálních diferenciálních rovnic popisujících rozložení řešeného pole knihovnu Hermes2D (http://hpfem.org/hermes2d), která je vyvíjena v rámci hp-FEM Group (http://hpfem.org), do které se řadí Nevadská univerzita v Renu, Akademie věd České republiky a také Západočeská univerzita v Plzni.

    Knihovna Hermes2D umožňuje řešit parciální diferenciální rovnice pomocí hp-FEM, která je obdobou klasické FEM (Finite Element Method). hp-FEM kombinuje změnu velikosti elementů (zjemnění řešené sítě) a změnu řádu polynomu, kterým je každý element proložen. Těmito změnami je zajištěna exponenciální konvergence výsledků. Pro provádění těchto změn využívá Hermes2D automatickou hp-adaptivitu (h-adaptivitu, p-adaptivitu), která se provádí na základě největších gradientů řešené veličiny.

    Agros2D rozlišuje mezi počáteční a řešenou sítí. Počáteční síť může být velmi hrubá. Při řešení je pak upravena díky automatické adaptivitě změnou velikosti jednotlivých elementů. Toho lze dosáhnout jak ručním nastavením určité oblasti, celého problému nebo využitím již zmiňované automatické hp-adaptivity.

    Výpočet kapacity jiskřiště

    link

    Jiskřiště je elektrické zařízení, které se využívá při měření a spínání vysokých a velmi vysokých napětí. Nejčastěji se skládá ze dvou elektrod, které mohou být provedeny v mnoha tvarech. Je možné se setkat s kulovými elektrodami, hroty, elektrodami deskovými a mnoha dalšími.

    Princip činnosti jiskřiště je založen na přeskoku v plynném nebo kapalném prostředí. Tohoto přeskoku lze však dosáhnout až za hranicí přeskokového napětí a je tedy možné toto napětí využít pro spínání nebo právě měření. V tomto příkladu budeme řešit elektrostatické pole mezi hrotem a deskovou elektrodou jiskřiště a následně vypočteme celkovou kapacitu daného uspořádání, která je potřebná například k určení nabíjecích charakteristik.

    Agros2D: Matematický model elektrostatického pole jiskřiště
    Matematický model elektrostatického pole jiskřiště

    Příprava problému

    link

    Založíme nový problém. Jako typ řešeného pole zvolíme Elektrostatické pole a jako typ problému Osově symetrický. Vzhledem k povaze úlohy (velký gradient elektrického potenciálu na ostrých hranách hrotu) využijeme při řešení úlohy automatickou hp-adaptivitu, kterou vybereme v příslušném seznamu. Počet adaptivních kroků nastavíme na 10 a relativní chybu (relativní rozdíl výpočtu energetické normy v aktuálním a předchozím adaptivním kroku) na 5 %. Vzhledem k tomu, že využíváme automatickou adaptivitu, nastavíme počet zjemnění (parametr udávající počet zjemnění celé sítě) na nulu a řád polynomu (nejnižší použitý řád polynomu pro celý problém) na jedničku.

    Agros2D: Nastavení problému
    Nastavení problému

    Podle nákresu geometrie a níže uvedené tabulky přidáme postupně všechny potřebné body a následně je propojíme pomocí hran. Vznikne nám tak geometrie o jedné jediné ploše, která reprezentuje okolní prostředí jiskřiště. Do této plochy je nutné přidat značku plochy.

    r (m) 0 0.15 0.15 0 0 0.015 0 0 0 0.35 0.35
    z (m) 0 0 0.01 0.01 0.06 0.11 0.11 0.31 -0.2 0.31 -0.2

    Dále nadefinujeme okrajové podmínky. Řešený problém popisují čtyři okrajové podmínky, které respektují zdrojovou elektrodu, zemnící elektrodu, hranici symetrie a fiktivní hranici ohraničující řešenou oblast. Hranice zdrojové a zemnící elektrody (Hrot, Deska) popisuje Dirichletova okrajová podmínka a je na nich zadán elektrický potenciál. Pro hranu hrotu použijeme 3 kV a pro hranu desky 0 V. Hranice symetrie i fiktivní hranice jsou popsány Neumannovou okrajovou podmínkou a je na nich zadána derivace elektrického potenciálu ve směru vnější normály. Vzhledem k osové symetrii a nulové hodnotě normálové složky intenzity elektrického pole musí být tato derivace na obou hranicích (Symetrie, Fiktivní hranice) rovna nule.

    Dále nadefinujeme potřebné materiály. V tomto příkladě to bude poměrně jednoduché, protože jediný materiál (Vzduch), který se v problému vyskytuje, je okolní prostředí elektrod, kde budeme uvažovat vzduch. Ten má relativní permitivitu rovnou jedné. Nyní již stačí definované okrajové podmínky a materiály přiřadit. Na druhém obrázku níže je 3D náhled na geometrii řešeného problému.

    Agros2D: Geometrie problému
    Geometrie problému

    Agros2D: 3D model problému
    3D model problému

    Řešení

    link

    Nejprve problém diskretizujeme. Pokud vše proběhne v pořádku, můžeme přistoupit k samotnému řešení. Z výstupu řešiče je patrné, že relativní chyba výpočtu nedosáhla požadovaných 5 % a je tedy nutné zvýšit počet adaptivních kroků v dialogu Nastavení problému. Zkusíme tedy nastavit 20 adaptivních kroků a řešení zopakujeme. Tentokráte již relativní chyba dosáhla méně než 5 % a řešení můžeme považovat za dostatečně přesné. Ze zobrazení počáteční sítě, sítě řešené a řádu polynomu, které si můžete navolit v dialogu Nastavení scény, je patrné, že v okolí ostrých rohů elektrod a v prostoru mezi nimi došlo k automatickému zjemnění a využití různých řádů polynomů. U zobrazeného průběhu potenciálu mezi deskou a hrotem je patrný jeho strmý nárůst a tedy i vysoká hodnota intenzity elektrického pole v těsné blízkosti hrotu.

    Agros2D: Počáteční síť, řešená síť a řád polynomu
    Počáteční síť, řešená síť a řád polynomu

    Agros2D: Detail počáteční sítě, řešená sítě a řádu polynomu
    Detail počáteční sítě, řešené sítě a řádu polynomu

    Agros2D: Průběh elektrického potenciálu mezi deskou a hrotem
    Průběh elektrického potenciálu mezi deskou a hrotem

    Pro kontrolu správnosti řešení u elektrostatických problémů je výhodné využít zobrazení ekvipotenciál (čar spojující místa se stejným potenciálem) a vektorů intenzity elektrického pole, které si můžete opět nastavit v dialogu Nastavení scény. Je patrné, že vektory intenzity elektrického pole směřují vždy od hrotu k desce a ekvipotenciály respektují fyzikální realitu rozložení elektrického potenciálu mezi elektrodami.

    Kapacitu jiskřiště lze spočítat dvěma způsoby. První z nich vychází z definice energie elektrostatického pole v lineárním prostředí a druhý z celkového náboje na jedné z elektrod, který je pro obě elektrody kvantitativně stejný a liší se pouze v polaritě. Využijeme první způsob a spočteme si tedy nejprve celkovou energii v oblasti vzduchu. Tato energie je rovna 171,2 mJ. Na základě vztahu We = ½ CU2 lze tedy již vypočítat kapacitu mezi elektrodami. Tento výpočet můžeme provést přímo v terminálu Agros2D, který využívá skriptovací jazyk Python. Pomocí položky Spustit příkaz… v menu Nástroje nebo klávesové zkratky Alt+C si otevřeme terminál a využijeme pro výpočet příkaz uvedený níže. Kapacita jiskřiště je rovna 11,413 pF. Na druhém obrázku níže je vykresleno rozložení elektrického potenciálu v prostoru jiskřiště.

    print(2*1.712e-05/3e6)
    

    Agros2D: Výpočet kapacity
    Výpočet kapacity

    Agros2D: Rozložení elektrického potenciálu
    Rozložení elektrického potenciálu

    Výpočet indukčnosti dvojlinky

    link

    Při projektování elektrických nebo sdělovacích sítí je velmi důležité znát kromě jiných parametrů také indukčnost použitých kabelů. V tomto příkladě si ukážeme, jak ji vypočítat pro dvojvodičový kabel protékaný stejnosměrným proudem. Budeme tedy řešit magnetostatické pole.

    Agros2D: Matematický model magnetostatického pole dvojlinky
    Matematický model magnetostatického pole dvojlinky

    Příprava problému

    link

    Založíme opět nový problém a jako typ řešeného pole zvolíme Magnetické pole, typ problému Kartézský a druh analýzy Ustálený stav. Ostatní nastavení můžeme nechat bez povšimnutí. Agros2D umožňuje kromě standardního ovládání využívat skriptování a také práci pomocí příkazů. Právě práci pomocí příkazů si ukážeme na vytvoření geometrie tohoto příkladu. Všechny funkce a klíčová slova, která můžete při práci využívat, naleznete v nápovědě v kapitole Scripting.

    Mezi základní funkce pro vytváření geometrie patří addnode() (přidání uzlu), addedge() (přidání hrany) a addlabel() (přidání značky oblastí). Vzhledem k podobě geometrie využijeme ale funkci addcircle(), která patří mezi uživatelské a slouží k přidání kružnice, přičemž umožňuje zároveň přiřadit okrajovou podmínku a značku oblasti (materiál).

    Funkce addcircle() má tři povinné a dva nepovinné parametry. První dva parametry udávají souřadnice středu kružnice a třetí parametr její poloměr. Čtvrtý nepovinný parametr udává, jaká okrajová podmínka se použije na hranici kružnice (pokud není zadána, použije se none a hranice bude při řešení uvažována jako rozhraní). Pátý, také nepovinný parametr, značí materiál, který se přiřadí oblasti uvnitř kružnice (pokud se nepoužije, nebude značka oblasti do kružnice přiřazena). Je důležité upozornit, že pokud je použit pátý parametr, bude čtvrtý parametr povinný.

    Nejprve nadefinujeme okrajové podmínky a materiály, abychom následně mohli plně využít funkci addcircle(). Přidáme jednu okrajovou podmínku, která bude respektovat fiktivní hranici problému, a bude jí popisovat Dirrichletova okrajová podmínka. Na této hranici bude vektorový magnetický potenciál roven nule.

    Při definici materiálů musíme respektovat, že jednotlivými žílami kabelu teče proud vždy opačným směrem. Proto musíme pro jejich definici přidat dva materiály. První z nich (L+) bude mít proudovou hustotu rovnou 2 MA/m2 (tato proudová hustota odpovídá proudu 2 A tekoucího vodičem o průřezu 1 mm2) a druhý (L-) s proudovou hustotou -2 MA/m2. Další materiál, který musíme nadefinovat, je izolace kabelu z PVC (Izolace) a jeho okolní prostředí (Prostredi). Izolace kabelu má stejnou permeabilitu jako okolní vzduch a není ji nutné definovat. Pro přehlednost řešeného modelu to však uděláme. Parametry všech materiálů jsou uvedeny v tabulce níže.

    L+ L- Izolace Prostredi
    Relativní permeabilita 1 1 1 1
    Proudová hustota 2 MA/m2 -2 MA/m2 0 0

    Nyní již s využitím příkazu addcircle() postupně vytvoříme celou geometrii.

    addcircle(-0.000965, 0, 0.000565, "none", "L+")
    addcircle(0.000965, 0, 0.000565, "none", "L-")
    addcircle(-0.000965, 0, 0.000965)
    addcircle(0.000965, 0, 0.000965)
    addcircle(0, 0, 0.025, "Fiktivni_hranice")
    

    Při přidávání geometrie respektující izolaci a okolní prostředí nemůžeme využít automatické přidání značky oblastí, protože ta se vždy přidává do středu kružnice, a tak by se značka přidala do nesprávné oblasti. Tyto značky oblasti tedy přidáme ručně nebo můžeme využít příkaz addlabel().

    addlabel(-0.0002, 0, 0, 0, "Izolace")
    addlabel(0.0002, 0, 0, 0, "Izolace")
    addlabel(0, 0.000965, 0, 0, "Prostredi")
    

    Agros2D: Geometrie problému
    Geometrie problému

    Řešení

    link

    Při výpočtu budeme postupovat obdobně jako u předchozího příkladu. Pokud diskretizace problému a jeho řešení proběhly v pořádku, můžeme přistoupit k výpočtu indukčnosti. Opět využijeme energetickou definici indukčnosti, kterou lze zapsat do vztahu Wm = ½ LI2. Nejprve určíme energii magnetického pole ve všech oblastech. Tato energie je rovna 1,185 μJ. Pro výpočet využijeme níže uvedený příkaz. Indukčnost dvojlinky je rovna 59,25 μH na 1 m délky (standardní hloubka u kartézských problémů). Pro ilustraci jsou na dalších obrázcích vykresleny siločáry magnetického pole a vektory magnetické indukce.

    print(2*1.185e-06/2**2)
    

    Agros2D: Výpočet indukčnosti
    Výpočet indukčnosti

    Agros2D: Siločáry magnetického pole
    Siločáry magnetického pole

    Agros2D: Vektory magnetické indukce
    Vektory magnetické indukce

    Oba dnes řešené příklady si opět můžete stáhnout (agros2d-jiskriste.a2d, agros2d-dvojlinka.a2d). Příští díl seriálu bude zaměřen především na skriptování a řešení sdružených problémů.

           

    Hodnocení: 100 %

            špatnédobré        

    Nástroje: Tisk bez diskuse

    Tiskni Sdílej: Linkuj Jaggni to Vybrali.sme.sk Google Del.icio.us Facebook

    Komentáře

    Vložit další komentář

    19.5.2010 13:25 Petr Tobiška
    Rozbalit Rozbalit vše Re: Agros2D – výpočet kapacity jiskřiště a indukčnosti dvojlinky
    Díky za zajímavý článek. Je možné v Agrosu modelovat i jiné fyzikální pole? Zajímalo by mě vektorové řešení EM pole v optických vláknech, případně se ztrátovými vrstvami.
    19.5.2010 16:57 Pavel Karban | skóre: 7
    Rozbalit Rozbalit vše Re: Agros2D – výpočet kapacity jiskřiště a indukčnosti dvojlinky
    Vysokofrekvenční pole není možné v současné době v Agrosu řešit. Jeden ze studentů Elektrotechnické fakulty v Plzni plánuje začlenit tento typ úlohy v rámci své diplomové práce.
    19.5.2010 19:12 peter
    Rozbalit Rozbalit vše Re: Agros2D – výpočet kapacity jiskřiště a indukčnosti dvojlinky
    Takze si nemozem nasimulovat plosak, ktory ma fungovat povedzme na 500MHz? (alebo niekolko GHz) Neexistuje taky (free) soft vhodny pre amatera? Nieco, na co netreba mat vysku z numeriky a eing. Alebo simulovat presluchy pri audio aplikaciach atd.
    19.5.2010 19:54 Franta Mach | skóre: 10
    Rozbalit Rozbalit vše Re: Agros2D – výpočet kapacity jiskřiště a indukčnosti dvojlinky
    To prozatím možné není. Plošné spoje je navíc takřka ve všech případech nutné modelovat jako 3D problém. Aplikací, které by umožňovali uvedený problém řešit je více, ale nejsem si jistý, zda by se dal najít nějaký volně šiřitelný. Napadá mě jen ANSYS Multiphysics, který je možné za určitých podmínek využít zdarma. Ten však není vhodný pro začátečníky a práce v něm vyžaduje poměrně značné znalosti z oblasti technických výpočtů. Další aplikací, která by byla pro daný účel využitelná je COMSOL Multiphysics. Ta je z pohledu uživatelského značně přívětivější, ale není volně šiřitelná. Možná by nejlepší volbou byl nějaký komplexní balík, který se zabývá návrhem plošných spojů a zároveň umožňuje i částečné řešení jejich modelů.
    Agros2D - Počítejte s námi!
    22.5.2010 00:26 host
    Rozbalit Rozbalit vše Re: Agros2D – výpočet kapacity jiskřiště a indukčnosti dvojlinky
    Mám radost z takovýho projektu, že škola vyvíjí aplikaci, která je kvalitní a je dostupná pro veřejnost (GPL). Takhle nějak jsem si vždycky představoval vysoký školství, že bude dostupným zdrojem informací, které bude společnost využívat ku svému prospěchu a rozvoji :-) Doufám, že takových projektů bude stále přibývat, nejen na ZČU. Jakožto stavař jsem především ocenil článek řešení přestupu tepla ve zdivu, díky.

    Založit nové vláknoNahoru

    ISSN 1214-1267   www.czech-server.cz
    © 1999-2015 Nitemedia s. r. o. Všechna práva vyhrazena.