ABAQUS划分网格的方法(读书总结).doc

划分网格的方
1. 独立实体（independent instance）和非独立实体(dependent instance)
对非独立实体划分网格时，应在窗口顶部的环境栏中把Object选项设为part,即对部件划分网格；对独立实体划分网格时, 应在窗口顶部的环境栏中把Object选项设为assembly, 即对装配件划分网格

2. 网格单元形状
在MESH功能模块中，Mesh—Controls，弹出Mesh Controls对话框，其中可选择单元形状。2D问题，有以下可供选择的单元形状。
1)Quad：网格中完全使用四边形单元；
2)Quad-dominated:网格中主要使用四边形单元，但在过渡区域允许出现三角形单元。选择Quad-dominated类型更容易实现从粗网格到细网格的过渡；
3）Tri:网格中完全使用三角形单元；
对于3D问题，包括以下可供选择的单元形状：
1） Hex:网格中完全使用六面体单元；
2） Hex-dominated:网格中主要使用六面体单元，但在过渡区域允许出现楔形（三棱柱）单元；
3） Tet:网格中完全使用四面体单元；
4） Wedge:网格中完全使用楔形单元；
Quad(2D问题)和Hex（3D问题）可以用较小的计算代价得到较高的精度，应尽可能选择这两种单元。

3. 网格划分技术
Structured(结构化网格)：采用结构化网格的区域显示为绿色；
Sweep(扫掠网格)：采用扫掠网格的区域显示为黄色；
Free(自由网格)：采用自由网格的区域显示为粉红色；
自由网格技术采用Tri和Tet，一般应选择带内部节点的二次单元来保证精度；结构化网格和扫掠网格一般采用Quad和Hex单元，分析精度相对较高。

4. 划分网格的算法
使用Quad和Hex单元划分网格时，有两种可供选择的算法：Medial Axis(中性轴算法)和Advancing Front(进阶算法)。
Medial Axis(中性轴算法)：首先把要划分网格的区域分成一些简单的区域，然后使用结构化网格划分技术来为简单区域划分网格。Medial Axis(中性轴算法)算法有以下特性：
1） 使用Medial Axis(中性轴算法)更容易得到单元形状规则的网格，但网格和种子的位置吻合较差；
2） 在二维模型中使用Medial Axis(中性轴算法)，选择Minimize the mesh transition(最小化网格过渡)，可以提高网格的质量，但使用这种方法生成的网格更容易偏离种子。
3） 如果在模型的一部分边上定义了受完全约束的种子，Medial Axis(中性轴算法)会自动为其他的边选择最佳的种子分布；
4） Medial Axis算法不支持由CAD模型导入的不精确模型和虚拟拓扑。
Advancing Front(进阶算法)：首先在边界上生成四边形网格，然后再向区域内部扩展。具有以下特性：
1） 使用Advancing Front算法得到的网格可以和种子的位置吻合的很好，但在较窄的区域内，精确匹配每粒种子可能使网格歪斜；
2） 使用Advancing Front算法更容易得到单元大小均匀的网格。在Explicite，网格的小单元会限制增量步长。
3） 使用Advancing Front算法更容易得到从粗网格到细网格的过渡；

5. 网格划分失败的解决方法
1）在Mesh模块中，主菜单Tools-Query下的Geometry Diagnostics,检查模型中是否有自由边、短边、小平面、小尖角或微小的缝隙。如果几何部件是由CAD模型导入的，则应注意检查是否模型本身有问题；如果几何部件是ABAQUS/CAE中创建的，应注意拉伸、切割操作，由于几何坐标的误差，出现的问题。
2）在Mesh模块中，主菜单Tools-Virtual Topology(虚拟拓扑)来合并小的边或面，或忽略某些边或顶点；如下图。

3）在Part模块中，主菜单Tools-Repair,可以修复存在的几何问题，例如可以选择Face/Replace faces来合并面
4）无法生成网格的位置加密种子；

7.选择三维实体单元的类型
7.1 节点数目和插值阶数（分以下三种）
1)线性单元：仅在单元角点处布置节点，在各个方向都采用线性插值；
2）二次单元：在每条边上有中间节点，采用二次插值；
3）二次修正单元：只有Tri和Tet使用这种网格。在每条边上有中间节点，采用修正二次插值；

7.2 线性完全积分单元
Meh-element type,保持Linear参数，就可以设置线性完全积分单元。如CPS4单元（四节点四边形双线性平面应力线性单元）；C3D8单元（8节点六面体线性完全积分单元）。所谓“完全积分”指单元具有规则形状时，所用的高斯积分点的数目足以对单元刚度矩阵中多项式进行精确积分。承受弯曲载荷时，线性完全积分单元会出现剪切自锁问题，造成单元过度刚硬，即使划分很细的网格，计算结果也不理想。

7.3二次完全积分单元
CPS8：8节点四边形二次平面应力完全积分单元；
C3D20：20节点六面体二次完全积分单元；
二次完全积分单元优点：
1） 对应力的计算结构很准确，适合模拟应力集中问题；
2） 一般情况下没有剪力自锁问题；
应用二次完全积分单元应注意以下问题：
1） 不能用于接触分析；
2） 对于弹塑性分析，如果材料是不可压缩的（如金属），则容易产生体积自锁；
3） 当单元发生扭曲或弯曲应力有梯度时，有可能出现某种程度的自锁；

7.4 线性减缩积分单元
对于Quad单元和Tex单元，ABAQUS/CAE默认使用线性减缩积分单元。如CPS4R(四节点四边形双线性平面应力线性减缩积分单元)；C3D8R(8节点6面体线性减缩积分单元)。
线性减缩积分单元比普通的完全积分单元在积分方向上少一个积分点。线性减缩积分单元只在单元的中心点有个积分点，故存在“沙漏”问题。采用线性减缩积分单元模拟承受弯曲载荷时，沿厚度方向至少应划分四个单元。
线性减缩积分单元有以下优点：
1） 对位移的求解结果较准确；2）网格存在扭曲变形时，求解的结果不会受到大影响；3）在弯曲载荷下不容易发生剪切自锁。
缺点：1）需要划分较细的网格来克服沙漏现象；2）不适合模拟应力集中问题。因为线性减缩积分只在单元的中心点有一个积分点，相当于常应力单元，它在积分点上的应力结果是准确的，而经过外插值和平均后的节点应力是不准确。

7.2 二次减缩积分单元
对于Quad单元和Tex单元，也可采用二次减缩积分单元。CPS8R(四节点四边形二次平面应力线性减缩积分单元)；C3D20R(20节点6面体二次减缩积分单元)。
优点：1）不会出现严重的沙漏现象；2）复杂应力状态下，对自锁现象不敏感；
同时使用二次减缩积分单元应注意的问题：1）不能在接触分析中使用；2）不适合大应变问题；3）节点应力结果往往低于二次完全积分单元。

7.3 非协调模式单元
对于Quad和Tex单元，也可设置非协调模式单元。仅在ABAQUS/Standard有非协调模式单元，其目的是克服线性完全积分剪切自锁问题。
优点：1）克服剪切自锁问题，在单元扭曲较小情况下，所得的位移和应力结果较准确；2）在模拟弯曲问题，在厚度方向只需很少的单元，就可以得到与二次单元相同的结果；3）使用了增强变形梯度的非协调模式单元，单元交界处不会重叠和开洞，因此容易扩展到非线性、有限应变的位移。