FEM之求解器(13)---利用C++开发振动求解器

2017-01-15 by:CAE仿真在线来源:互联网

结合以前的例子,本文介绍一下利用C++开发的 3D 振动(Vibration)求解器,主要有如下功能

1. 支持2阶四面体和六面体单元

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

3. 支持 Abaqus *.inp 关键字文件

之所以选择 Abaqus *.inp作为关键字文件,是因为可以方便的使用 Abaqus,Calculix与Solver 做比较。这一点很重要,Abaqus Nastran 等通用FEA之所以成为很多行业的标准,就是因为积累了大量经过验证的benchmark case。

使用与 Abaqus,Calculix benchmark的case 来介绍3D 动力学求解器:

1. FaModel中导入*.step模型

2. 划分二阶四面体网格

3. 不用设置边界条件,求固有频率

4. 设置材料,求解属性(如图1)

5. 导出为Abaqus *.inp格式文件

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

图1:

Nastran求解结果:

Calculix求解结果:

FasVib求解结果:

模态1:最大位移 2.03e-4, 频率2.05e-6Hz

模态3:最大位移 1.69-4,频率3.21e-6Hz

模态6:最大位移 2.5e-4,频率3.21e-6Hz

模态10:最大位移 5.20e-4,频率6.2Hz

最大位移结果完全一样,频率求得较大误差,有可能是FEA参数的设置不符合实际情况造成的。

开放分享：优质有限元技术文章,助你自学成才

编辑