有限元分析几个重要原则(一)

有限元分析重要原则(1)：尽量把所有不会发生位移的节点都固定住，不要让求解器再去通过迭代计算来确定这些节点的位移。

举个简单例子。一个二维平面应变问题，包含两个弹性体，即圆筒和平板，如图1所示。在圆筒中心的圆孔内壁上定义了固支边界条件，在平板顶部中央的A点给定义了位移U2=-2，希望使平板向正下方移动，和圆筒发生接触。提交分析后，计算可以完成，但在分析结果中看到平板发生了异常的位移，如图2所示。这是什么原因引起的？

图1 定义了位移边界的模型

图2 后处理时看到平板发生了异常的位移

对于三维模型，每个部件都有3个平动自由度和3个转动自由度；对于二维模型，每个部件都有2个平动自由度和1个转动自由度。在建立静力分析模型时，必须在模型每个实体的所有平动和转动自由度上定义足够的边界条件，以避免它们出现不确定的刚体位移，否则将导致分析往往无法收敛，即使能够收敛，结果也往往是错误的。

本例中，圆筒上定义了固支边界条件，不会出现刚体位移。但是平板在x方向上没有定义任何边界条件，因此在x方向上的刚体位移是不确定的；在y方向上，只在一个节点（A点）上给定了位移U2，这时整个平板仍然可以绕A点做刚体转动，即除了A点之外，平板上的其他节点的U2都是不确定的。

尽管整个模型并没有使平板发生转动或x方向平动的载荷，直观感觉上此模型似乎是没问题的，但这样的模型符合有限元分析的要求。这种“因为没有受力，所以不会移动”的因果关系，只是我们根据生活经验在头脑中进行逻辑分析时的思路，而Abaqus/Standard的求解过程恰恰与此相反，其过程是：迭代尝试各种可能的位移状态，检验它们是否能够满足静力平衡方程。

在本实例中，无论平板发生多大的转动或x方向的平动，都可以满足静力平衡方程，即符合静力平衡条件的位移解有无限个，因此会出现“数值奇异”。有限元是一种数值计算方法，计算过程中的微小数值误差会导致平板在缺乏约束的自由度上发生刚体运动，因此会看到如图2所示的异常结果。

解决方法(1)：本实例中的模型是左右对称的，因此圆筒和平板都应该只取一半建模，在整个对称面上定义对称边界条件，即U1=0，这样平板就不会再发生转动或在x方向上产生平动。需要注意的是，一个模型是否具有对称性，不仅取决于它的几何形状，还要看材料、载荷、边界条件和接触等是否都是对称的，即变形后的模型是否是对称的。

解决方法(2)：如果模型不具有对称性，就需要根据具体情况添加适当的边界条件，以消除不确定的刚体位移。本实例中，可以在平板中央对称线上定义边界条件U1=0。需要注意的是，不能只定义A点的U1=0，因为这样整个平板仍然可以绕A点做刚体转动。

 

有限元分析重要原则(2)：模型中仅仅靠两个外力达到静力平衡是不够的，必须要借助于边界条件处的支反力达到平衡。

同样举个简单的例子。二维平板两端受到均布拉力，如图3a所示，如果直接对整个平板建模，没有边界条件，提交分析后，往往会出现“数值奇异”的警告信息。因为这时尽管整个平板处于静力平衡状态，但仍然会出现不确定的刚体位移，因为整个平板是悬浮在空中的，有无数种可能的位移状态。

一个更为合理的建模方法如图3b所示，根据对称性，只取1/4建模，在两个对称面上分别定义对称边界条件。这样就能够保证静力平衡方程的位移解是唯一的，静力分析才能够收敛。

图3a 二维平板两端受拉

图3b 根据对称性取1/4建模

 

相关知识(1)

建立一个模型之前，应该考虑的第一件事就是，这个模型是否具备对称性，是否可以只取1/2、1/4甚至是1/8进行建模。这样做有多方面的重要意义。

1）在对称面上定义对称边界条件，有助于避免刚体位移问题；

2）可以大大减小模型的规模，缩短计算时间；

3）接触面上的节点减少一半，接触分析就更容易收敛；

4）施加了对称边界条件之后，整个模型的支承状况变得更加稳固，可能出现的位移状态大大减少，Abaqus/Standard不用再去反复尝试那些不具备对称性的位移解，这样就更容易找到正确的位移解，会使复杂的非线性分析更容易收敛。

 

相关知识(2)

对于动力分析，不需要在所有自由度上定义足够的边界条件，因为动力分析会考虑惯性力，可以避免产生无限大的瞬时运动。如果在动力分析时看到“数值奇异”的警告信息，往往是由于模型中存在其他问题，例如“过度塑性”等。

来源：DeepFEA