ANSYS接触问题的计算方法及参数设置
接触问题的关键在于接触体间的相互关系，此关系又可分为在接触前后的法向关系与切向关系。法向关系：在法向，必须实现两点：1）接触力的传递。2）两接触面间没有穿透。ANSYS通过两种算法来实现此法向接触关系：罚函数法和拉格朗日乘子法。1．罚函数法是通过接触刚度在接触力与接触面间的穿透值（接触位移）间建立力与位移的线性关系：接触刚度接触位移法向接触力对面面接触单元17，接触刚度由实常数FKN来定义。穿透值在程序中通过分离的接触体上节点间的距离来计算。接触刚度越大，则穿透就越小，理论上在接触刚度为无穷大时，可以实现完全的接触状态，使穿透值等于零。但是显而易见，在程序计算中，接触刚度不可能为无穷大（否则病态），穿透也就不可能真实达到零，而只能是个接近于零的有限值。以上力与位移的接触关系可以很容易地合并入整个结构的平衡方程组KXF中去。并不改变总刚K的大小。这种罚函数法有以下几个问题必须解决：1）接触刚度FKN应该取多大？2）接触刚度FKN取大些可以减少虚假穿透，但是会使刚度矩阵成为病态。3）既然与实际情况不符合的虚假穿透既然是不可避免的，那么可以允许有多大为合适？因此，在ANSYS程序里，通常输入FKN实常数不是直接定义接触刚度的数值，而是接触体下单元刚度的一个因子，这使得用户可以方便地定义接触刚度了，一般FKN取01到1中间的值。当然，在需要时，也可以把接触刚度直接定义，FKN输入为负数，则程序将其值理解为直接输入的接触刚度值。对于接近病态的刚度阵，不要使用迭代求解器，例如PCG等。它们会需要更多的迭代次数，并有可能不收敛。可以使用直接法求解器，例如稀疏求解器等。这些求解器可以有效求解病态问题。穿透的大小影响结果的精度。用户可以用PLESOLCONTPENE来在后处理中查看穿透的数值大小。如果使用的是罚函数法求解接触问题，用户一般需要试用多个FKN值进行计算，可以先用一个较小的FKN值开始计算，例如01。因为较小的FKN有助于收敛，然后再逐步增加FKN值进行一系列计算，最后得到一个满意的穿透值。
fFKN的收敛性要求和穿透太大产生的计算误差总会是一对矛盾。解决此矛盾的办法是在接触算法中采用扩展拉格朗日乘子法。此方法在接触问题的求解控制中可以有更多更灵活的控制。可以更快的实现一个需要的穿透极限。2．拉格朗日乘子法与扩展拉格朗日乘子法拉格朗日乘子法与罚函数法不同，不是采用力与位移的关系来求接触力r