关于接触刚度的讨论[致远书屋]

1、关于接触刚度的讨论（转载）ANSYS基本概念 2008-09-11 10:11 阅读65评论0 字号： 大大 中中 小小 BBS 锦城驿站我最近在做接触分析，老觉得不合理。接触刚度应该是与接触面等材料属性有关，为什么还要自己定义这个刚度？我仿照使用ANSYS6。1进行结构力学分析里面的接触例子，求解时出现realconstant2hasbeenreferencedbyelementtypeselementtypes1and2oneofwhichiscontactelement.书上说的是通过共享实常数来判别接触对，为什么又出现这样的错误提示呢？请大家帮忙。决定接触刚度所有的接触问题都需要定义接

2、触刚度，两个表面之间渗量的大小取决了接触刚度，过大的接触刚度可能会引起总刚矩阵的病态，而造成收敛困难，一般来谘，应该选取足够大的接触刚度以保证接触渗透小到可以接受，但同时又应该让接触刚度足够小以使不会引起总刚矩阵的病态问题而保证收敛性。程序会根据变形体单元的材料特性来估计一个缺省的接触刚度值，你能够用实常数FKN来为接触刚度指定一个比例因子或指定一个真正的值，比例因子一般在0.01和10之间，当避免过多的迭代次数时，应该尽量使渗透到达极小值。为了取得一个较好的接触刚度值，又可需要一些经验，你可以按下面的步骤过行。1、开始时取一个较低的值，低估些值要比高估些值好因为由一个较低的接触刚度导致的渗透

3、问题要比过高的接触刚度导致的收敛性困难，要容易解决。2、对前几个子步进行计算3、检查渗透量和每一子步中的平衡迭代次数，如果总体收敛困难是由过大的渗透引起的（而不是由不平衡力和位移增量引起的），那么可能低估了FKN的值或者是将FTOLN的值取得大小，如果总体的收敛困难是由于不平衡力和位移增量达到收敛值需要过多的迭代次数，而不是由于过大的渗透量，那么FKN的值可能被高估。4、按需要调查FKN或FTOLN的值，重新分析。(ANSYS公司的资料）我的理解：接触刚度与接触面等材料属性无关，理论上接触刚度越大越好，尽量小的接触渗透。但难收敛。通过共享实常数来判别接触对。要注意使用一个contactelem

4、ent和一个targetelement共享实常数。如：type,1!defined1asacontactelementreal,1mat,1!meshtype,2!defined2asatargetelementreal,1mat,1!mesh我还有一点疑惑：接触刚度对整个过程的应力分布有何影响？会不会出现，用不同的接触刚度都可以给出求解结果，但结果随着接触刚度的不同有很大差异的现象？可以这样理解：现实中两个接触的物体，正常是不能渗透的，两个接触的物体由于受接触力发生变形。仿真的接触物体，除了变形外，还会出现接触渗透，这是因为定义了接触刚度，接触力除以接触刚度就是接触渗透量。所以需要大的接触刚

5、度，但要解决收敛问题。我在计算中确实遇到了这样的问题：给出不同的接触刚度，都能给出具体的计算结果。但针对我的具体情况，解除的计算结果相差很大（10%还多一点），按你的观点，我是不是可以认为，我定义的接触刚度大的那个结果是相对准确的？接触刚度这么说我都糊涂了材料肯定是不允许渗透的，这是接触协调条件决定的刚度按照字面的理解应该是抗击变形的能力，这应该就是与材料属性相关的，能解释一下这个概念吗？谢谢了！xidianstuwrote:接触刚度这么说我都糊涂了材料肯定是不允许渗透的！错，为了计算的方便（避免大的单元奇异及重划分）在ANSYS中允许一定量的渗透，这个渗透的大小，由接触刚度决定，不过好像还有

6、其它的一个什么因素（我忘了）。请问dingdingdang你怎么知道定义大的接触刚度计算的数值相对准确呢？我的计算时也出现像你这样的问题！就是在不改变其他的参数只改变接触刚度其迭代次数和结果都没有大的变化！（我做的练习是实体单元！）依我的个人理解：解的结果会随着接触刚度，穿透容忍度的不同而有所不同。但对于穿透容忍度足够小的情况下，解的结果将随接触刚度影响不会很大。不过，在穿透容忍度小特别小的时候当然不容易收敛。因此在穿透容忍度一定的情况下，当然是接触刚度大穿透小的解更加准确。接触刚度是否就是与分析误差有关的一个量，可能接触分析需要迭代，也就是判断什么时候终止迭代？大家看看这个吧！我想对大家的理

7、解可能有帮助摘自ansys中文网站用户专区。（我记得我曾经贴出过）在有限元分析中，接触单元通常用来描述两物体相互接触或滑动的界面。近年来，ANSYS开发了一系列的接触单元。刚开始有节点对节点单元CONTAC12和CONTAC52，接着有节点对地单元CONTAC26，然后有节点对面单元CONTAC48和CONTAC49。最近几年，我们引入一类面对面接触单元CONTA169和CONTA174，同时还有一种新的节点对节点单元CONTA178。虽然接触单元的参数具有多样性，但我们在使用他们时可谨记重要的一点，他们具有一个共同的特点，即除了CONTA178的KEYOPT（2）=0或1外，所有的接触单元都

8、有接触刚度。在现实中实际上相邻结构之间只是一种空隙，但在有限元分析中，这种空隙是一带有刚度的接触单元，这是因为通过刚度矩阵来实现接触算法的。一些接触单元要求使用者输入刚度值，同时另外的接触单元若没有输入则使用缺省值。分析工程师所面对的问题就是针对给定的条件确定一个合理的刚度值。如果过高，问题将会不收敛，如果过低，可能得到错误的结果。那么我们所面对的问题是怎样才能找到一个正确的刚度值？我认为唯一的方法就是我们必须试用不同的值直到找到正确的值。也就是刚开始我们应该使用一个较小的值，然后稳步的增加直到分析的结果不再有什么变化。那么对于我们这一特定分析的问题，这一点就是我们所想要的合适值。我们可举例说

9、明，如图1所示，平行放置两个悬臂梁，并有少许的交迭，下面的左边固支，上面的右边固支，当在上面梁的自由端施加一个向下位移时，梁变形弯曲并接触下面的梁，然后一起向下运动。用SOLID45单元划分梁，用TARGE170和CONTA174面面接触单元来描述相互作用。在此基础上，把CONTA174单元的刚度从非常低变到非常高，从而来观察它对结果的影响和收敛的迭代次数。图2说明了下梁自由端的偏移随接触单元刚度的变化情况，当刚度增加时，偏移量接近一个常数值（我们可以假定它是一个正确的结果。）图3说明求解所需的迭代次数，当接触单元刚度增加时，求解所需的迭代次数也是增加的，并服从指数关系。如果刚度过高，问题很有

10、可能根本就不收敛。图4说明在上梁自由端接触单元的渗透量，当刚度增加时，渗透量降低。从这些图可知，当接触单元的刚度为10e6时，可获得合理精确的结果。任何大于该值的刚度对下梁的偏移量没有什么影响，而求解所需的迭代数却显著的增加。对于这个题目，10e6的刚度是很适合的。但是，如果改变边界条件、网格密度、两梁之间的相对位置、材料特性或梁的几何形状，能获得满意结果的接触刚度值将是不同的。比如，如果网格密度增加，则接触单元数将增加，每一个单元上的载荷将降低。如果接触单元数增加两倍，一个合适的接触单元刚度值应为原来的一半。由于每个题目都是不一样的，所以在求解之前并没有通用的方法来确定接触单元刚度的最佳值。

11、我们不得不试算一个我们认为合适的值然后查看计算结果。一个有经验的分析工程师可能只查看一个计算结果来判定所取值的合适度，但对于大多数情况而言，最好用一个合理而不过度精确的刚度值进行第一次求解，然后用10倍于该值的刚度进行第二次求解，如果两者结果相差很小，而迭代数增加很多，那么我们则正好取得了曲线上的突变点，从而获得相当好的结果。接触单元刚度问题仅仅是一个例子，即对于分析工程师来说，总是置疑于分析结果的正确与否是非常重要的，并要意识到数值仿真的局限性和潜在的假设及他们怎样影响所分析问题的结果。我在计算爆炸动力问题的时候曾经使用过20的比例因子这是接触问题的计算方法。接触问题的关键在于接触体间的相互

12、关系（废话），此关系又可分为在接触前后的法向关系与切向关系。法向关系：在法向，必须实现两点：1）接触力的传递。2）两接触面间没有穿透。ANSYS通过两种算法来实现此法向接触关系：罚函数法和拉格朗日乘子法。1罚函数法是通过接触刚度在接触力与接触面间的穿透值（接触位移）间建立力与位移的线性关系：接触刚度*接触位移=法向接触力对面面接触单元17*，接触刚度由实常数FKN来定义。穿透值在程序中通过分离的接触体上节点间的距离来计算。接触刚度越大，则穿透就越小，理论上在接触刚度为无穷大时，可以实现完全的接触状态，使穿透值等于零。但是显而易见，在程序计算中，接触刚度不可能为无穷大（否则病态），穿透也就不可能

13、真实达到零，而只能是个接近于零的有限值。以上力与位移的接触关系可以很容易地合并入整个结构的平衡方程组K*X=F中去。并不改变总刚K的大小。这种罚函数法有以下几个问题必须解决：1）接触刚度FKN应该取多大？2）接触刚度FKN取大些可以减少虚假穿透，但是会使刚度矩阵成为病态。3）既然与实际情况不符合的虚假穿透既然是不可避免的，那么可以允许有多大为合适？因此，在ANSYS程序里，通常输入FKN实常数不是直接定义接触刚度的数值，而是接触体下单元刚度的一个因子，这使得用户可以方便地定义接触刚度了，一般FKN取0.1到1中间的值。当然，在需要时，也可以把接触刚度直接定义，FKN输入为负数，则程序将其值理解

14、为直接输入的接触刚度值。对于接近病态的刚度阵，不要使用迭代求解器，例如PCG等。它们会需要更多的迭代次数，并有可能不收敛。可以使用直接法求解器，例如稀疏求解器等。这些求解器可以有效求解病态问题。穿透的大小影响结果的精度。用户可以用PLESOL,CONT,PENE来在后处理中查看穿透的数值大小。如果使用的是罚函数法求解接触问题，用户一般需要试用多个FKN值进行计算，可以先用一个较小的FKN值开始计算，例如0.1。因为较小的FKN有助于收敛，然后再逐步增加FKN值进行一系列计算，最后得到一个满意的穿透值。FKN的收敛性要求和穿透太大产生的计算误差总会是一对矛盾。解决此矛盾的办法是在接触算法中采用扩

15、展拉格朗日乘子法。此方法在接触问题的求解控制中可以有更多更灵活的控制。可以更快的实现一个需要的穿透极限。2拉格朗日乘子法与扩展拉格朗日乘子法拉格朗日乘子法与罚函数法不同，不是采用力与位移的关系来求接触力，而是把接触力作为一个独立自由度。因此这里不需要进行迭代，而是在方程里直接求出接触力（接触压力）来。Kx=F+Fcontact从而，拉格朗日乘子法不需要定义人为的接触刚度去满足接触面间不可穿透的条件，可以直接实现穿透为零的真实接触条件，这是罚函数法所不可能实现的。使用拉格朗日乘子法有下列注意事项：1）刚度矩阵中将有零对角元，使有些求解器不克使用。只能使用直接法求解器，例如波前法或系数求解器。而P

16、CG之类迭代求解器是不能用于有零主元问题的。2）由于增加了额外的自由度，刚度阵变大了。3）一个可能发生的严重问题，就是在接触状态发生变化时，例如从接触到分离，从分离到接触，此时接触力有个突变，产生chattering（接触状态的振动式交替改变）。如何控制这种chattering，是纯粹拉格朗日法所难以解决的。因此，为控制chattering，ANSYS采用的是罚函数法与拉格朗日法混合的扩展拉格朗日乘子法。在扩展拉格朗日法中，可以采用实常数TOLN来控制最大允许穿透值。还有最大允许拉力FTOL。这两个参数只对扩展拉格朗日乘子法有效。在扩展拉格朗日乘子法里，程序按照罚函数法开始，与纯粹拉格朗日法类

17、似，用TOLN来控制最大允许穿透值。如果迭代中发现穿透大于允许的TOLN值，（对178单元是TOLN，而对面面接触单元171-174则是FTOLN）则将各个接触单元的接触刚度加上接触力乘以拉格朗日乘子的数值。因此，这种扩展拉格朗日法是不停更新接触刚度的罚函数法，这种更新不断重复，直到计算的穿透值小于允许值为止。尽管与拉格朗日法相比，扩展拉格朗日法的穿透并不是零，与罚函数法相比，可能迭带次数会更多。扩展拉格朗日法有下列优点：1）较少病态，个接触单元的接触刚度取值可能更合理。2）与罚函数法相比较少病态，与单纯的拉格朗日法相比，没有刚度阵零对角元。因此在选择求解器上没有限制，PCG等迭代求解器都可以

18、应用。3）用户可以自由控制允许的穿透值TOLN。（如果输入了TOLN，而使用罚函数法，则程序忽略它）要是在fkn设置中的所设置的值超过100，那样计算是不是合理的？接触分析中确定接触刚度是较关键的一步。 下面是说明书中的介绍，供参考： 所有的接触问题都需要定义接触刚度，两个表面之间穿透量的大小取决于接 触刚度。过大的接触刚度可能回引起总刚矩阵的病态，从而造成收敛困难。一 般来说，应该选取足够大的接触刚度以保证接触穿透小到可以接受，但同时又 应该让接触刚度足够小以不致于引起总刚矩阵的病态而保证收敛性。 ANSYS程序根据下伏柔体单元的材料特性，来估计一个缺省的接触刚度值。 用户可用实常数FKN来为接触刚度指定一个比例因子或指定一个绝对值。比例 因子一般在0.01-10之间。对于大变形问题，通常为0.01-0.1。用户应当总是 检验以使穿透到极小值，而又避免过多的叠代次数。 为了确定一个较好的接触刚度值，可能需要一些经验。用户可以按照下面 的步骤来进行尝试： 1、开始时取一个较低的值。低估值要比高估值好，因为由一个较低的接触刚度 导致的穿透问题，比过高的接触刚度导致的收敛性困难，要