Abaqus有限元分析——有限元网格划分基本原则

个人专栏—ABAQUS专栏

1. Abaqus2023的用法教程——与VS2022、oneAPI 2024子程序的关联方法 Abaqus2023的用法教程——与VS2022、oneAPI 2024子程序的关联方法
2. Abaqus有限元分析——有限元网格划分基本原则 Abaqus有限元分析——有限元网格划分基本原则
3. 各向同性线弹性材料本构模型umat的应用 各向同性线弹性材料本构模型umat的应用
4. ABAQUS用户子程序二次开发——UMAT介绍（一） ABAQUS用户子程序二次开发——UMAT介绍（一）
5. Fortran语法介绍（一） Fortran语法介绍（一）
6. Fortran语法介绍（二） Fortran语法介绍（二）
7. Fortran语法介绍（三） Fortran语法介绍（三）

---

---

$有限元网格划分简介$：

• 有限元的基本思想是将结构离散化，利用简化几何单元近似逼近连续体，根据变形协调条件综合求解。 $网格划分$是进行有限元数值模拟分析至关重要的一步，它直接决定所生成的有限元模型的优劣，进而影响数值分析结果的精确性和经济性。

• 网格的划分没有定式，只能根据经验划分网格。对有限元知识的深刻了解和丰富的经验式保证网格划分优劣的前提。

• 网格划分涉及单元的形状及其拓扑类型、单元类型、网格生成器的选择、网格密度、单元编号及几何体素。

$网格划分遵循的基本原则$：

$网格数量$:

• 网格数量直接影响计算精度和计算耗时，网格数量增加提高计算精度的同时会增加计算耗时。

• 当网格数量较少时增加网格，计算精度可明显提高，计算时耗不会有明显增加；

• 当网格数量增加到一定程度后，再继续增加网格时，精度提高程度减小，计算耗时大幅增加。

• 在确定网格数量时应权衡两者因素的影响。具体应用时采用试算法，比较不同网格数量划分的计算结果，如果两次计算结果相差较大，可以继续增加网格，相反则停止计算。

• 在决定网格数量时应考虑数据的类型：

  1. 在静力分析时，如果仅计算结构的变形，网格数量可以少一些。如果需要计算应力，则在精度要求相同的情况下应取相对较多的网格。

  2. 在响应计算中，计算应力响应所取的网格数应比计算位移响应多。

  3. 在计算结构固有动力特性时，若仅仅是计算少数低阶模态，可以选择较少的网格，如果计算的模态阶次较高，则应选择较多的网格。

  4. 在热分析中，结构内部的温度梯度不大，不需要大量的内部网格，这时可划分较少的网格。

$网格疏密$:

• 在结构不同部位采用大小不同的网格，这是为了适应计算数据的分布特点。

• 在计算数据变化梯度较大的部位（如应力集中处、几何形状、材料、厚度变化的位置），为了较好地反映数据变化规律，需要采用比较密集的网格。而在计算数据变化梯度较小的部位，为减小模型规模，则应划分相对系数的网格，使整个结构表现出疏密不同的网格划分形式。

• 疏密不同的网格主要用于应力分析（包括静应力和动应力），计算固有特性时趋于采用较均匀、规则的网格形式。这是因为固有频率和振型主要取决于结构质量分布和刚度分布，不存在类似应力集中的现象，采用均匀网格可使结构刚度矩阵和质量矩阵的元素不致相差太大，可减少数值计算误差。

• 在结构温度场计算中也趋于采用均匀网格。

• 分析结果完成后，需检查以下各项，当差值与其中最大应力的比值较大时，该位置网格需要细分：

  1. 单元应力的连续性，比较相邻单元应力值的差值；
  2. 应力偏差：节点上的单元节点应力和节点平均应力差值的较大值

$单元阶次$:

• 许多单元都具有线性、二次和三次等形式，其中二次和三次形式的单元称为高阶单元。

• 选用高阶单元可提高计算精度，因为高阶单元的曲线或曲面边界能够更好地逼近结构的曲线和曲面边界，其次高阶插值函数能更高精度地逼近复杂场函数，当结构形状不规则、应力分布或变形很复杂时可以选用高阶单元。但高阶单元的节点数较多，在网格数量相同的情况下由高阶单元组成的模型规模要大很多，在使用时应权衡考虑计算精度和时间。

• $增加网格数量和单元阶次都可以提高计算精度。$在精度一定的情况下，用高阶单元离散结构时应选择适当的网格数量，太多的网格并不能明显提高计算精度，反而会使计算时间大大增加。为兼顾计算精度和计算量，同一结构可以采用不同阶次的单元，即精度要求高的重要部位用高阶单元，精度要求低的次要部位用低阶单元。不同阶次单元之间或采用特殊的过度单元连接，或采用多点约束等式连接。

$网格质量$:

• 网格质量是指网格几何形状的合理性。质量好坏将影响计算精度。质量太差的网格甚至会中止计算。网格各边或各个内角相差不大、网格面不过分扭曲、边节点位于边界等分点附近的网格质量较好。网格质量可用细长比、锥度比、内角、翘曲量、拉伸值、边节点位置偏差等指标度量。

• 形状比（长边与短边距离之比）一般实体单元的形状比越大，分析误差越大。对于板壳单元，评价应力为主时不宜超过1：3，评价位移为主时不宜超过1：5；对于块体单元，评价应力为主时不宜超过1：2，评价位移为主时不宜超过1：3；在应力分布几乎没有变化的区域使用单元，可适当放大形状比。

• 倾角（表示单元偏离直角四边形的程度）四边形的内倾角最好在 $45^\circ C\sim 135^\circ C $，不要超过 $15^\circ C\sim 165^\circ C $.

• 锥度（限于四边形）用几何偏离表示四边形单元的变形程度。

• 翘曲 指四边形单元的四个节点偏离同一平面的程度（限四边形单元）。在两个曲面相连的位置的四边形单元，翘曲比较明显的四边形单元应使用两个三角形单元，精度相对较好。

$网格分界面和分界点$:

• 结构中的一些特殊界面和特殊点应分为网格边界或节点以便定义材料特性、物理特性、载荷和位移约束条件。即应使网格形式满足边界条件特点，而不应让边界条件适应网格。

$位移协调性$:

• 单元上的力和力矩能够通过节点传递相邻单元。为保证位移协调性，一个单元的节点必须同时也是相邻单元的节点，而不应是内点或边界点。相邻单元的共有节点具有相同的自由度性质。否则，单元之间必须用多点约束等式或约束单元进行约束处理。

$网格布局$:

• 当结构形状对称时，其网格也应划分对称网格，以使模型表现处相应的对称性质（如集中质矩阵对称）。不对称布局会引起一定误差。

$节点和单元编号$:

• 节点和单元编号影响结构总刚度矩阵的带宽和波前数，因而影响计算时间和存储容量的大小，合理编号有利于提高计算速度。

$刚性连接$:

• 刚性连接的功能是在不太重要的位置上将结构连接起来，并传递荷载。

• 使用刚性连接时，在连接位置的某一个方向位移不是连续的，应力分布也不圆滑。

• 从属节点不能产生与相邻节点的正常位移–位移不连续。

• 应力的不连续发生在距连接位置单元特性长度（一般为厚度或高度尺寸）的局部范围内，该范围内的应力不可信。

• 在受扭位置最好不要使用刚性连接，因为刚性连接约束了截面的翘曲，夸大了结构的抗扭刚度。

$细节处理$:

• 分析模型中对于构造细节是否应该进行模拟，视分析目的而定。

• 对于分析不重要的细节不应当包含在分析模型中。如结构整体受力分析、结构动力分析时不需要对结构细节进行详细描述。

• 在结构局部应力分析中，如倒角或孔等细节可能时最大应力出现的位置，这些细节对于分析目的十分重要，必须进行模拟。

---

欢迎对Abaqus感兴趣的朋友们查看：Abaqus-UMAT开发精品书籍及umat子程序学习

Abaqus非线性粘弹性模型子程序umat——广义MAXWELL粘弹性模型umat解析（朱-王-唐本构模型）

如果你喜欢以上内容，或者对Abaqus学习有兴趣，欢迎收藏关注，博主将持续更新。你的关注、收藏是我持续创作的动力！