在有限元分析软件中已划分好网格的实体之间的接触区域设置接触对单元。接触对可以在四面体或六面体三维单元上进行划分。一个接触对包括一个目标单元和一个接触单元，形状可以是三角形或四边形。接触单元和目标单元分别覆盖在两个接触体的接触表面上，目标单元覆盖在刚度较大的实体上，而接触单元覆盖在刚度小的实体上。每个接触单元和目标单元可以覆盖一个或多个接触体的三维单元表面。
      接触单元和目标单元为空间壳单元，其形函数和插值方法与结点无转动自由度的普通空间壳单元安全相同。接触对的结点与被覆盖面的实体单元的表面结点位置相同，但编号不同，编号顺序符合右手法则。目标单元和接触单元取相同的实常数编号。
      接触对单元是一种实际上不存在的假想单元。在未接触区，其拉伸模量和剪切模量都趋向于零；在固定接触区，模量变得足够大，以阻止其侵入。当发生侵入接触单元时，通过调整接触单元的模量，使两者变形协调。接触对单元的刚度矩阵为。
      接触面必须给目标面一定的约束条件才能满足接触相容性，由于这些约束条件随状态的改变而改变，一般采用增量迭代法进行求解。目前求解接触问题常用的算法有罚函数法和拉格朗日乘子法，两者共同点是将约束条件的变分问题转化为无约束条件的变分问题求解。拉格朗日乘子法是通过在泛函式引用拉格朗日乘子项来精确地满足接触约束条件，与罚函数法相比，对接触刚度系数不敏感，但是当引入大量的乘子后，使系统的刚度矩阵不对称，需要很多的迭代次数才能收敛。罚函数法的因子取值对计算结果的精度影响很大，必须根据渗透情况进行多次调整。ANSYS软件推荐使用罚函数法和拉格朗日乘子组合型求解接触问题，其表达式为式中：kn为接触刚度，即罚因子；g为渗透深度。
      两者表面之间穿透量的大小取决于接触刚度，应该有足够大的接触刚度，以保证接触穿透小到可以接受，但接触刚度不能太大，否则会引起总刚度矩阵的病态而不能保证收敛性。
      为了满足相容性，必须在目标面的切向施加一个摩擦力，摩擦力的大小与选择的摩擦类型有关。考虑螺栓连接中的接触和变形，选用库伦摩擦接触。在库伦摩擦接触中，两个接触面在开始相互滑动之前是黏合状态，当剪应力超过定义的等效剪应力时，两个表面之间开始滑动。库伦摩擦定义为式中：τｌｉｍ为极限剪切应力；τ为等效剪切应力；μ为摩擦系数；ｐ为接触面压应力。

                                                                                专业从事机械产品设计│有限元分析│CAE分析│结构优化│技术服务与解决方案
                                                                                                                                                      杭州那泰科技有限公司
                                                                              本文出自杭州那泰科技有限公司www.nataid.com，转载请注明出处和相关链接！