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