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 00:22 | Bezpečnostní upozornění

Společnost ZONER informuje o bezpečnostním incidentu, při kterém došlo ke zcizení a zveřejnění části přihlašovacích údajů zákazníků k elektronické poště a webhostingu CZECHIA.COM.

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

Byla vydána nová stabilní verze 1.13 (1.13.1008.32) webového prohlížeče Vivaldi (Wikipedie). Z novinek vývojáři zdůrazňují možnost zobrazení otevřených i uzavřených listů pomocí ikonky Okno na postranní liště a vylepšené stahování (YouTube). Nejnovější Vivaldi je postaveno na Chromiu 62.0.3202.97.

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

Byla vydána verze 2017.3 dnes již průběžně aktualizované linuxové distribuce navržené pro digitální forenzní analýzu a penetrační testování Kali Linux. Aktualizovat systém lze pomocí příkazů "apt update; apt dist-upgrade; reboot". Z novinek lze zmínit 4 nové nástroje: InSpy, CherryTree, Sublist3r a OSRFramework.

Ladislav Hagara | Komentářů: 0
včera 01:55 | Bezpečnostní upozornění

Společnost Uber potvrdila bezpečnostní incident a únik dat v roce 2016. Unikly údaje o 57 milionech cestujících (jména, emailové adresy a čísla mobilních telefonů) a 600 tisících řidičích (navíc čísla řidičských průkazů).

Ladislav Hagara | Komentářů: 1
21.11. 23:44 | Humor

Co vypíše příkaz man půl hodiny po půlnoci? Text "gimme gimme gimme". Jedná se o virtuální velikonoční vajíčko připomínající skupinu ABBA a její hit Gimme! Gimme! Gimme! (A Man After Midnight). Problém nastane, pokud gimme gimme gimme nabourá automatizované testování softwaru. To se pak příkaz man musí opravit [Bug 1515352] [reddit].

Ladislav Hagara | Komentářů: 10
21.11. 18:11 | Zajímavý článek

Mozilla.cz informuje, že Firefox na Fedoře podporuje Client Side Decorations. Firefox na Linuxu se vykresluje včetně standardního záhlaví okna, které je v případě webového prohlížeče většinou nadbytečné a ubírá drahocenné vertikální místo na obrazovce. Verze distribuovaná uživatelům Fedory však nyní obsahuje experimentální podporu pro takzvané Client Side Decorations, které umožňují vykreslování „oušek“ panelů do záhlaví okna.

Ladislav Hagara | Komentářů: 12
21.11. 05:00 | Bezpečnostní upozornění

Maxim Goryachy a Mark Ermolov ze společnosti Positive Technologies budou mít v prosinci na konferenci Black Hat Europe 2017 přednášku s názvem "Jak se nabourat do vypnutého počítače, a nebo jak v Intel Management Engine spustit vlastní nepodepsaný kód". O nalezeném bezpečnostním problému informovali společnost Intel. Ta bezpečnostní problém INTEL-SA-00086 v Intel Management Engine (ME), Intel Server Platform Services (SPS) a Intel

… více »
Ladislav Hagara | Komentářů: 38
21.11. 01:33 | Zajímavý projekt

Na Humble Bundle byla spuštěna akce Humble Book Bundle: Java. Za 1 dolar a více lze koupit 5 elektronických knih, za 8 dolarů a více 10 elektronických knih a za 15 dolarů a více 15 elektronických knih věnovaných programovacímu jazyku Java od nakladatelství O'Reilly. Peníze lze libovolně rozdělit mezi nakladatelství O'Reilly, neziskovou organizaci Code for America a Humble Bundle.

Ladislav Hagara | Komentářů: 0
21.11. 00:11 | Zajímavý projekt

Článek na OMG! Ubuntu! představuje rodinu písma IBM Plex. Jedná se o open source písmo (GitHub) navržené a uvolněné společností IBM (YouTube, Carbon Design System). Ukázka na Font Squirrel.

Ladislav Hagara | Komentářů: 12
20.11. 23:22 | Komunita

Na Humble Bundle lze získat počítačovou hru Brütal Legend (Wikipedie, YouTube) běžící také v Linuxu zdarma. Speciální akce končí ve středu v 19:00.

Ladislav Hagara | Komentářů: 0
Jak se vás potenciálně dotkne trend odstraňování analogového audio konektoru typu 3,5mm jack z „chytrých telefonů“?
 (9%)
 (1%)
 (1%)
 (1%)
 (75%)
 (14%)
Celkem 749 hlasů
 Komentářů: 37, poslední 21.11. 15:21
    Rozcestník

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

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

    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.