Submodeling

volně přeložil: Petr Novák,

z knihy Erdogan Madenci, Ibrahim Guven, The Finite Element Method and Applications in Engineering Using ANSYS, 2006, ISBN 0 387 28289 0; kapitola 11: Advanced Topics in Ansys.

 obsah
 

Úvod

Teorie

Příklad

    Analýza globálního modelu

    Skript pro globální model

    Analýza submodelu

    Skript pro submode  

Úvod:

Submodeling je metoda používaná pro získání přesnějších výsledků ve specifické oblasti. Vyžaduje existující model, který je pojmenován „globální model" a korespondující řešení. Globální model používá hrubé síťování pro minimalizaci výpočetní náročnost. Oproti tomu je na submodel použita mnohem jemnější síť pro zvýšení přesnosti výpočtu. Převod výsledků ze submodelu do globálního modelu je zajištěn okrajovými podmínkami aplikovanými na submodel.

Teorie:

Představte si ozubené kolo jako je na obr. 1. Pokud je naším největším zájmem znát hodnotu intenzity napětí v regionu viz. obr. 1., pak samostatný konečně prvkový model (submodel) tohoto regionu je generován s mnohem jemnějším síťováním. Hranice submodelu, které nejsou externí hranice globálního modelu („external boundaries") jsou nazvány „řezné hranice" („Cut boundaries"). Výsledné stupně volnosti z globálního modelu jsou interpolovány jako okrajové podmínky na řeznou hranici. Okrajové podmínky na externích hranicích submodelu musí být stejné jako ty, které byly použity na Globálním modelu. Řezné hranice se nemusí shodovat s globálními elementovými hranicemi.

Síťování Globálního modelu a Submodelu.

 

Obr. 1. Síťování Globálního modelu a Submodelu.

Typické kroky pro analýzu s využitím submodelu sestává z:

  • Vytvoření globálního modelu a získání výsledků.

  • Uložení globálního modelu a jeho řešení.

  • Vytvoření submodelu a jeho uložení.

  • Výběr uzlů podél řezných hranic a jejich zápis do souboru.

  • Vyvolání globálního modelu.

  • Znovu-načtení výsledků globálního modelu.

  • Interpolace na řeznou hranici (zápis okrajových podmínek do souboru s uzly řezné roviny.

  • Znovu-načtení submodelu.

  • Načtení souboru s podmínkami uzlů řezné roviny.

  • Zatížení (zatížení tělesa a/nebo podmínek podél externích hranic).

  • Získání řešení submodelu.

Následuje příklad submodelové analýzy desky s kruhovým otvorem, na kterém je demonstrováno, jak správně použít submodel do analýzy.

...obsah  

Příklad:

Tenká čtvercová plocha (10×10 in) s kruhovým otvorem (o poloměru 2.5 in) je namáhána stejnoměrným tahovým napětím (100 psi) ve vertikálním směru podél její horní hrany. Uchycena je na spodní hraně (viz. obr. 2 ). Protože je deska tenká a namáhána rovinným zatížením, je možná plošná idealizace napětí. Materiálové vlastnosti jsou: Modul pružnosti E=10×106 psi a Poissonova konstanta ν=0.25. Úkolem je vytvořit analýzu s použitím submodelingu v okolí kruhového otvoru pro zjištění přesnějších výsledků.

Plocha s kruhovým otvorem stejnoměrně tahově namáhána.  

Obr. 2. Plocha s kruhovým otvorem stejnoměrně tahově namáhána.

 ...obsah

Analýza globálního modelu:

Pokud chceme, aby naše submodelová analýza byla úspěšná, je nutné specifikovat název pracovního souboru („jobnames") jak pro globální model, tak i pro submodel (použitím /FILNAM příkazu).

V našem příkladě je globální model specifikován jako GLO. Poté co je získáno řešení globálního modelu je použit příkaz SAVE, který ukládá jak model, tak i řešení.

...obsah

Skript pro globální model:

/FILNAM,GLO

/PREP7

ET,1,PLANE42

MP,EX,1,10E6

MP,NUXY,1,0.25

RECTNG,0,5,0,5

PCIRC,2.5

ASBA,1,2

CSYS,1

LSEL,S,LOC,X,2.5

LESIZE,ALL, , ,8

CSYS

ALLSEL

LESIZE,ALL, , ,4

LCCAT,2,3

MSHKEY,1

AMESH,ALL

ARSYM,X,ALL

ARSYM,Y,ALL

NUMMRG,ALL

FINISH

/SOLU

NSEL,S,LOC,Y,-5

D,ALL,ALL

NSEL,S,LOC,Y,5

SF,ALL,PRES,-1000

ALLSEL

SOLVE

SAVE

FINISH

/POST1

PLNSOL,S,Y

FINISH

!Přejmenování pracovního souboru na GLO.

!Vstup do PREPROCESSORu.

!Použít PLANE42 jako element typ 1.

!Specifikace mat. vlastností, (Modul pruž.)

!a poissonovy konst.

!Vytoření obdélníku.

!Vytvoření kruhové plochy.

!Odečtení ploch.

!Přepnutí do globálního cylindrického SS.

!Výběr úseček podle umístění.

!Počet rozdělení úseček.

!Přepnutí do globálního kartézského SS.

!

!

!Zřetězit úsečky 2 a 3.

!Užít "MAPPED" síťování.

!Síťovat všechny plochy.

!Zrcadlit síťované plochy podle YZ-roviny.

!Zrcadlit síťované plochy podle XZ-roviny.

!Sloučit všechny zdvojené entity.

!Ukončit PREPROCESSOR.

!Vstup do řešiče.

!Vybrat uzly podle umístění.

!Okrajové podmínky na vybrané uzly.

!

!Zatížení na vybrané uzly.

!

!Řešit.

!Uložit model a výsledek.

!Ukončit řešící processor.

!Vstup do POSTPROCESSORu.

!Zobrazit rozložení napětí v y-směru.

!Ukončit POSTPROCESSOR.

...obsah

Analýza submodelu:

Na části modelu je vytvořena jemnější síť. Na obr. 2 je uzavřena čárkovanou čarou, obklopující submodel. Obr. 3 zobrazuje síť globálního modelu a region, kde je vytvořen submodel. Zjemněná síť pro submodel je zobrazena na obr. 4 . Jméno pro submodel je v našem případě „SUB". Všechny úsečky kromě čáry ohraničující kruhovou plochu jsou řezné hranice. Po vytvoření sítě pro submodel jsou vybrány uzly na řezné hranici a zapsány do externího ASCII souboru (SBS.NOD) s použitím „NWRITE" příkazu. Následně je model uložen („SAVE" příkaz). Interpolace řezných hranic může být provedena s využitím hraničních uzlových bodů ze souboru SBS.NOD a výsledku řešení globálního modelu (ze souboru GLO.RST). Pro tento účel je načten globalní model (RESUME příkaz) a jsou načteny výsledky.

Interpolace řezných hranic je provedena příkazem CBDOF, který zapíše interpolované výsledky okrajových podmínek (v našem případě do souboru SBS.CBD). Submodel je následně zatížen a SBS.CBD soubor je načten jako vstup.

Napětí v y-směru z globálního modelu a submodelu jsou zobrazena na obr. 5 .

Síť pro globální model (vlevo) a umístění submodelu (vpravo).

Obr. 3. Síť pro globální model (vlevo) a umístění submodelu (vpravo).

Síť submodelu (překrytý na globálním modelu).  

Obr. 4. Síť submodelu (překrytý na globálním modelu).

 Napětí v y-směru získané z globálního modelu a submodelu vztažené k x-souřadnicím.

Obr. 5. Napětí v y-směru získané z globálního modelu a submodelu vztažené k x-souřadnicím.

...obsah

Skript pro submodel:

/CLEAR,START

/FILNAM,SUB

/PREP7

ET,1,PLANE42

MP,EX,1,10E6

MP,NUXY,1,0.25

RECTNG,0,3,0,3

PCIRC,2.5

ASBA,1,2

LESIZE,9, , ,10

LESIZE,10, , ,10

LESIZE,5, , ,40

LESIZE,2, , ,20

LESIZE,3, , ,20

LCCAT,2,3

MSHKEY,1

AMESH,ALL

NSEL,S,LOC,Y,0

NSEL,A,LOC,X,0

NSEL,A,LOC,Y,3

NSEL,A,LOC,X,3

NWRITE,SBS,NOD

ALLSEL

SAVE

FINISH

RESUME,GLO,db

/POST1

FILE,GLO,rst

CBDOF,SBS,NOD, ,SBS,CBD

FINISH

RESUME,SUB,db

/SOLU

/INPUT,SBS,CBD

SOLVE

FINISH

/POST1

PLNSOL,S,Y

!Vyčistit databázi a začít nový.

!Změnit jméno prac. soub. na SUB.

!Vstup do PREPROCESSORu.

!Použít PLANE42 jako element typ 1.

!Definování materiálové vlast. (EX)

!Definování Poissonovy konst.

!Vytvořit obdélník.

!Vytvořit kruh.

!Odečíst kruh z obdélníku.

!Specifikování počtu rozdělení úsečky.

!

!

!

!

!Zřetězit úsečky 2 a 3.

!Užít "MAPPED" síťování.

!Síťovat všechny plochy.

!Vybrat uzly podle umístění.

!Dovybrat úsečky podle umístění.

!

!

!Uložit uzly do souboru "SBS.NOD".

!

!

!

!Načíst soubor "GLO.DB".

!Vstup do POSTPROCESSORu.

!Deklarace výstupního souboru.

!Provést interpolaci řezných hranic.

!Ukončit POSTPROCESSOR.

!Načíst soubor "SUB.DB".

!Vstup do řešiče.

!Načíst okrajové podmínky řezné hranice.

!Řešit.

!Ukončit řešič.

!Vstup do POSTPROCESSORu.

!Zobrazit rozložení napětí v y-směru.

...obsah