ansys計(jì)算線程數(shù)的案例
ANSYS模型的Model Order Reduction（提高瞬態(tài)計(jì)算速度數(shù)百倍）
最近在做一些Fast simulation的工作，用到了Model Order Reduction (MOR),覺得十分有用，在這里推薦給大家，并且附上ANSYS 和Matlab的source code.現(xiàn)在最流行的MOR主要基于Krylov subspace projection method, 基本原理如下：
假設(shè)原來的linear system model的node數(shù)是N (NxN的矩陣）, N ~ 100000, 我們可以創(chuàng)建一個(gè)N x r 的矩陣，把原來NxN的矩陣project到 r x r的矩陣上去 (r<<N)，這樣求解過程可以在這個(gè)小矩陣上完成，然后再把結(jié)果project回到原來的NxN的矩陣空間上去。
以下是源程序介紹：
(1) Example.mac: ANSYS mac file 用來建立一個(gè)三維的熱傳導(dǎo)模型，并加上載荷。用戶必需設(shè)置文件里的SOLVE_OPTION變量。
取0時(shí)，輸出模型的heat capacitance 和 heat conductance 矩陣(用HBMAT commands)，分別保存到C_File.dat和K_file.dat（Thermal load 也保存到K_File.dat）。
取1和2時(shí)分別做steady-state和transient simulation。
(2) MOR_ODE.m: matlab 程序，讀取K_File.dat和C_file.dat進(jìn)行MOR,并且和ANSYS 結(jié)果（ansys.dat,在Example.mac里SOLVE_OPTION＝2）的對(duì)比。用戶必需設(shè)置Solver_Option = 1調(diào)用matlab的lu分解函數(shù)。（Solver_Option = 0 是使用external 的 MOR_MUMPS.exe程序做matrix factorization,速度可提高數(shù)倍。
展開 