FEM之求解器(13)---利用C++開(kāi)發(fā)振動(dòng)求解器

1. 支持2階四面體 和 六面體單元

2.求解模型的固有頻率,特征值和特征向量

3. 支持 Abaqus *.inp 關(guān)鍵字文件

之所以選擇 Abaqus *.inp作為關(guān)鍵字文件,是因?yàn)榭梢苑奖愕?使用 Abaqus,Calculix與Solver 做比較。這一點(diǎn)很重要,Abaqus Nastran 等通用FEA之所以成為很多行業(yè)的標(biāo)準(zhǔn),就是因?yàn)榉e累了大量經(jīng)過(guò)驗(yàn)證的benchmark case。

使用與 Abaqus,Calculix benchmark的case 來(lái)介紹3D 動(dòng)力學(xué)求解器:

1. FaModel中導(dǎo)入*.step模型

2. 劃分二階四面體網(wǎng)格

3. 不用設(shè)置邊界條件,求固有頻率

4. 設(shè)置材料,求解屬性(如圖1)

5. 導(dǎo)出為Abaqus *.inp格式文件

6. 分別用Nastran, Calculix, FasVib求解

圖1:

Calculix求解結(jié)果:

FasVib求解結(jié)果:

模態(tài)1:最大位移 2.03e-4, 頻率2.05e-6Hz

模態(tài)3:最大位移 1.69-4,頻率3.21e-6Hz

模態(tài)6:最大位移 2.5e-4,頻率3.21e-6Hz

模態(tài)10:最大位移 5.20e-4,頻率6.2Hz

最大位移結(jié)果完全一樣,頻率求得較大誤差,有可能是FEA參數(shù)的設(shè)置不符合實(shí)際情況造成的。