然后按下按钮【optional settings…】，在弹出来的接触属性设置中，首先在BASIC页中设置法向惩罚刚度因子（Normal Penalty Stiffness）为0.1。

然后点击friction页，设置刚度矩阵（Stiffness matrix）是非对称阵（Unsymmetric）。然后点击【OK】，关闭该接触属性设置对话框，回到前面的接触向导对话框。在接触向导对话框中，点击【Create】,创建接触对。

此时，如果打开单元类型对话框，会发现里面多了两种单元。即：TARGE170是目标单元，它覆盖在销孔内表面；CONTA174是接触单元，它覆盖在钢销外表面。

如果打开实常数会发现多了一个3号实常数，此实常数即代表了上面所创建的接触。

4．设置边界条件

由于模型仅创建了四分之一，因此要对下图所示的4个面施加对称边界条件。

如下图所示约束钢块一个侧面的所有自由度。

5．求解

在Solution Controls中进行求解设置，打开大变形开关，关闭自动时间步长，并设置该载荷步结束时的时间是100，只有一个载荷子步。然后进行求解。

9．后处理

进入后处理，查看Von Mises应力云图。可见，在接触处应力最大，越往外层，应力越来越小。这与实际情况是一致的。

下面是接触处的渗透图，最大渗透量是0.003229.关于更多接触结果的查看方式可参考文章《在ANSYS mechanical APDL经典界面中查看接触分析结果》。