Výpočtové systémy - OptiSLang

Citlivostní analýza

V programu OptiSLang nalezneme nejenom optimalizační algoritmy, ale také spoustu metod, které využijeme před samotnou optimalizací k zjištění charakteru daného problému. Jednou z nich je citlivostní analýza sloužící zjištění závislostí (korelace) mezi vstupními i výstupními parametry.

Pro účely citlivostní analýzy přidáme několik dalších parametrů. Vstupními parametry bude kromě orientace láken v lamině ještě její hustota, moduly pružnosti laminy ve směru vláken \(E_1\) a ve směru kolmém na vlákna \(E_2\) a rozměr nosníku \(b\). Dalším výstupním parametrem pak bude hmotnost nosníku. Soubor zmen_uhel.proc používaný v optimalizaci upravíme na zmen_parametry.proc obsahující:

*open_model ..\..\pasek.mud
*set_orient_angle
40.4
*mater_param structural:youngs_modulus1
87000e6
*mater_param structural:youngs_modulus2
5500e6
*mater_param general:mass_density
1500
*geometry_type mech_planar_pstress
*geometry_param norm_to_plane_thick
0.05
*element_mass
all_existing
*write_marc pasek_uprava.dat yes
*quit yes

Z kódu je patrné, že hmotnost nosníku zjistíme pomocí příkazu *element_mass, který vytiskne danou hmotnost do příkazové řádky. Abychom ale mohli hmotnost načíst do programu OptiSLang potřebujeme textový soubor s hodnotou hmotnosti. Řešením je soubor mentat.log generovaný automaticky při úpravě modelu obsahující výpis příkazové řádky. Jelikož bychom tento soubor přepsali při následném vyhodnocování modelu je nutné jej před vyhodnocením přejmenovat. To provedeme přidáním příkazu copy mentat.log hmotnost.log do souboru krok.bat, který uložíme jako krok_parametry.bat s příkazy:

call mentat -bg -pr zmen_parametry.proc
call run_marc -j pasek_uprava
copy ..\..\vyhodnot.proc;
copy mentat.log hmotnost.log
call mentat -bg -pr vyhodnot.proc

Nastavení parametrů, spuštění a vyhodnocení citlivostní analýzy je zobrazeno v animaci 1.

Uvodni obrazek je kdo ví proč fuč, ale animace by tu měla být, stačí kliknout!

Animace 1: Provedení citlivostní analýzy.


Západočeská univerzita v Plzni | Fakulta aplikovaných věd | Katedra mechaniky