ANSYS 中使用接触向导定义多个接触对详细实例(图文)

ANSYS 中如何使用接触向导定义接触对

在 ANSYS 中定义接触通常有两种方法：

1. 用户自己手工创建接触单元和目标单元。这种方法，在定义接触和目标单元时还比较简单，但是在设置或修改单元属性和定义实常数时却比较复杂。需要用户对接触有较深刻的理解和通过实践积累丰富的经验。

2. 使用接触管理器中的接触向导定义接触对：使用接触管理器 (接触向导) 定义接触对 (即接触单元和目标单元) 时，可以定义除了点-点接触以外的各种接触类型；它可以自动生成接触单元和目标单元，并提供了一组默认的单元属性和实常数值。使用这些默认的设置，加上适当的求解设置，对于多数接触问题都能够获得收敛的结果。而且，如果使用默认设置时，计算不收敛或对结果不太满意，也可以通过接触管理器 (接触向导) 对单元属性和实常数方便的进行修改和调整。

因此，我们推荐，在可能的情况下，尽量使用接触管理器 (接触向导) 来定义接触。本文将通过一个实例介绍接触管理器的基本使用方法。

所使用的例子如下：

两块平板，中间夹一个圆球。上面平板的上表面承受压力，分析模型的变形和应力随压力的变化。

两块平板，尺寸都是 (100*100*20)，相距 100。中间夹一个半径 50 的圆球。两个平板分别与圆球的上下边缘接触。尺寸单位为 mm。几何模型如图 1。

图 1 中，为了能够划分映射网格，分别对体积进行了切割材料属性为： 两块平板： E = 201000 Mpa； μ= 0.3 圆球： E = 70100 Mpa； μ= 0.33

接下来对各个 Volumes 划分网格，单元类型采用 solid186 (20 节点六面体)，单元边长统一取 6 mm。网格划分结果如图 2 所示：

载荷为上平板上表面均布压力，最大值 10 Mpa；约束条件将在后面介绍。

以下介绍如何使用接触管理器 (接触向导) 来分析这一接触问题。 一、 使用接触管理器 (接触向导) 前的准备工作

1. 在 ANSYS 中，接触单元和目标单元必须覆盖在已有的单元表面上。因此，在创建接触之前，应该先尽量完成几何建模和划分网格的工作，特别是在要创建接触的表面 (边界线)，必须以某种方式划分网格，以在相应的边界面 (线) 上生成节点，作为以后创建接触时的基础。划分网格的方式则和几何模型有关。比如，对三维面模型，可以划分壳体单元，后面直接以壳体面作为创建接触对的基础；而三维体模型，可以划分实体单元，后面则以实体单元的表面作为创建接触对的基础。

2. 在使用接触管理器 (接触向导) 创建接触对时，可以选为接触面或目标面的对象有：线、面、节点、节点组 (component) 等。如果模型比较复杂，临时选择不太方便，建议将准备创建接触的实体边界 (面、线) 分别建为单独的实体组 (Component)，或者分别取出其包含的节点 (使用 NSLA、NSLL 命令)，建为单独的节点组件 (Component)。 二、打开接触管理器

在前处理中，点击 GUI 屏幕上命令输入小窗口右边的第三个小图标，就可以打开接触管理器：

使用接触管理器 (接触向导) 定义接触时，接触单元是自动生成的。有关的实常数也在接触向导中定义。在创建接触对时，对话框中有一个 Options 按钮，点击它就会出现各种接触属性的设置，这些属性分别与接触单元和目标单元的有关单元选项以及相应的实常数相对应，可以根据问题的具体情况进行设置。

三、创建接触对

下面来创建接触对。需要创建两个接触对，分别为两个平板与上下两个圆半球之间可能接触的部位。

为了便于创建接触对，先创建 4 个组件，分别包含上下平板与球之间的两个接触对的可能接触面。图 4 中名为 A1、A2、A3 和 A4 的四个不同颜色的 Component 分别为相应的四个组件：

为了创建接触对时容易选择面，只显示这 4 个 Component 中的 Areas：

在其中点击 Create 按钮，创建接触对。如果成功创建了接触对，会出现如下提示，说明该接触对的实常数的 ID 号：

在其中，点击 Finish 按钮，完成接触对的创建，回到接触管理器对话框，并会显示所创建的接触对：

如果接触对创建不成功，也会出现一个提示信息，表示接触对创建不成功，这时则需要坚持接触对的创建过程，找出问题、进行修改后再次创建该接触对。

成功创建接触对之后，首先应该检查一下接触对的情况，特别是接触单元和目标单元的法线方向是否正确。

不要退出接触管理器，选中新创建的接触对，点击 Plot Element 按钮，显示接触单元和目标单元如下：

如果先点击 Show Normals 按钮 ，再点击Plot Element 按钮，在显示接触单

元和目标单元时会同时显示单元的法线：

接触单元和目标单元正确的法线方向应该是相对的。假如发现接触单元和/或目标单元的法线方向不对，则需要进行修改。假设有如下图情况，接触单元 (平板底面处) 的法线方向不对，需要进行修改：

为修改接触单元的法线方向，先单独显示接触单元，然后在其中点击 Flip Normals on Elements：

在弹出的拾取框中点击 Pick All > OK，则所有 (接触单元) 的法线方向被颠倒，与目标单元的法线方向相对：

至此，第一个接触对定义结束。接下来，不要退出接触管理器，定义下半球与下面平板之间的接触。

为了定义第二个接触对，先调出与下半球和下面平板有关的组件 A3 和 A4，并显示 Areas 如下：

然后按照和上面相同的方法定义第二个接触对：

由上图可见，接触单元和目标单元的法线方向互相指向对方，是正确的。

使用接触管理器定义接触对后，如果想修改接触对的设置，可以再次进入接触管理器，选择相应的接触对，然后点击属性 (Propertities) 按钮，就可以对所需的属性重新进行设置和修改。

此外，如果你熟悉实常数定义和单元类型的关键选项设置，在使用接触管理器定义接触对之后，也可以分别到实常数或单元类型中查看和修改相应的参数。

在使用接触管理器导定义接触时，不必预先定义相关的单元类型和实常数，软件会自动生成所需的单元类型和实常数。本例中在定义接触之前仅定义过一种单元类型 – solid186，没有定义过实常数。在定义接触后，分别查看一下单元类型和实常数，可以看到分别新增了一些内容。

对于单元类型，新增了 4 种单元类型，分别为两个接触对所使用的目标单元和接触单元，注意：不同的接触对分别使用不同的接触单元和目标单元：

可以查看有关单元类型的属性。如第一种目标单元类型的关键选项如下：

第一种接触单元类型的关键选项如下：

对于实常数，新增了 2 个实常数，分别对应两个接触对所使用的单元 (同时适用于目标单元和接触单元)：

要查看完整的实常数，应将实常数作为接触对中接触单元的属性来显示：

所显示的实常数如下：

图 32 与单元类型 3 - CONTA714 相关的实常数

这些单元类型的关键选项和实常数，在使用接触向导定义接触对时的选项中，多数都有对应的内容。可以在接触向导中选择，也可以在退出接触向导后在单元类型和实常数中进行修改。即使是手工方式创建的接触对，除了点-点接触外，也同样可以打开接触管理器，然后使用接触向导进行查看和修改。本文就不详细叙述了。

四、 设置分析类型，施加载荷和约束

进入 Solution 模块，设置新的分析类型为 STATIC。

在上平板的上表面施加均布压力 10 Mpa。首先选择要施加压力的 Areas：

图 34 选择要施加压力的 Areas

选择 Areas 后，在上图拾取对话框中点击 OK，弹出施加压力的对话框，在其中输入应力值，然后点击OK：

为了查看压力是否正确，先打开边界条件的显示，并设置压力为箭头显示：

再次显示整个模型如下：

约束条件是：

(1) 下平板的下表面 Uz ＝ 0，先选择要施加约束的 AREAS：

点击 OK 后弹出在 AREAS 上施加约束的对话框，在其中

选择 UZ 作为要约束的自由度，约束位移值 为：常数 0，可以不输入 (默认值)：

点击 OK 后，模型显示如下图：

(2) 对 UX 和 UY的约束

首先设置下平板下表面中间点的 Ux = Uy = 0，后面将根据分析情况进一步进行修改。设置过程如下：

在拾取对话框中点击 OK，弹出 Apply U,ROT on KPs 对话框，在其中选择 Ux 和 Uy，然后点击 OK：

施加 KP 点约束后的模型如下：

五、 设置求解参数

对于接触问题，通常可以设置如下求解参数：

(1) 设置结果输出方式：选择 Every substep，输出所有子步的结果：

(2) 打开求解过程的图形追踪

(3) 激活求解控制

(4) 设置时间步长

(5) 设置收敛准则

还可以设置其他求解选项，如线性搜索等，这里不再详述。

五、 求解

点击 Solve > Current LS，弹出两个表框，一个是 /Status Command，可以查看所有的求解设

置：

如求解设置无误，可以点击该表框右上角的 × 号关闭该边框，然后在另一个表框中点击 OK，开始求解：

开始求解之后，可以显示 Output 窗口，查看求解过程中的一些信息，比如，对于本例陆续出现如下有关存在初始间隙和由于存在摩擦改用非对称求解器的信息：

事实上，我们只对下面平板施加了约束，对小球和上面平板未施加约束，需要靠它们之间及小球与下面平板之间的接触来造成约束。在没有发生接触之前，整个模型确实是缺少约束的。

第一个错误信息说明出错的表现是 Uz 值过大；在第二个错误信息

中对出现这一问题的原因进行了分析 –模型中存在只有靠接触表面闭合 (即发生接触) 才能与其它部分连到一起的部分，并提出使用 CNCHECK 检查模型的建议。

在 Output 窗口中出现错误信息的同时，在 GUI 窗口中会出现 Error 信息框，给出与 Output 窗口中类似的出错信息。

不过，出现出错信息并不一定意味着计算过程就此结束。在两个错误信息框中都有一个 Proceed 按钮，点击该按钮，可以使程序继续执行该分析过程。

这里选择对两个错误信息框都点击 Proceed 使分析继续执行： 以下是两个错误信息框：

但是查看 Output 窗口中的信息，会发现时间值不是一直上升的：

这说明有的时间步计算没有收敛。由于打开了自动时间步，遇到这种情况，程序会自动对上一时间步的时间增量进行二分，然后重新确定计算的时间点 (上一时间点加上新的时间步增量)，继续计算。

本例，最终是收敛的。计算过程中对力收敛情况的追踪曲线如下：

上图中，CRIT 曲线表示力收敛的准则值，L2 表示实际计算得到的准则力值，当 L2 的值小于 CRIT 的值(L2 的曲线低于 CRIT 曲线) 时，计算是收敛的；如果 L2 的值大于 CRIT 的值 (L2 的曲线高于 CRIT 曲线)，程序会继续迭代，直到 L2 的值小于 CRIT 的值 (收敛) 或者迭代次数达到允许的最大值而L2 的值仍大于CRIT 的值 (不收敛)。对后一情况，程序会将对上一子步的时间增量减半，求出新的时间值，然后继续分析，直到时间子步的增量小于设置的时间子步最小值，如仍不收敛，将会发布一个计算不收敛的信息并停止计算。本例的后续计算是收敛的。直到计算结束， 出现如下信息框：

七、 结果后处理

分析结束后，可以进入通用后处理，查看各种计算结果。 以下是最后一个子步的 Uz 位移分布云图：

以下是最后一个子步的 Mises 应力分布云图：

为了更清楚的了解结果，可以只显示部分实体的结果，下面分别显示两块平板和小球的Mises 应力结果：为此，首先显示所有 Volumes：

再选择上下两块平板的所有 Volumes：

选择属于两块平板的所有实体：

本文来源：https://www.bwwdw.com/article/58m3.html