matlab接触刚度编程,关于接触刚度的定义

所有的接触问题都需要定义接触刚度，两个表面之间渗 量的大小取决了接触刚度，过大的接触刚度可能会引起总刚矩阵的病态，而造成收敛困难，一般来谘，应该选取足够大的接触刚度以保证接触渗透小到可以接受，但同时又应该让接触刚度足够小以使不会引起总刚矩阵的病态问题而保证收敛性。

程序会根据变形体单元的材料特性来估计一个缺省的接触刚度值，你能够用实常数FKN来为接触刚度指定一个比例因子或指定一个真正的值，比例因子一般在0.01和10之间，当避免过多的迭代次数时，应该尽量使渗透到达极小值。

为了取得一个较好的接触刚度值，又可需要一些经验，你可以按下面的步骤过行。

1、 开始时取一个较低的值，低估些值要比高估些值好因为由一个较低的接触刚度导致的渗透问题要比 过高的接触刚度导致的收敛性困难，要容易解决。

2、 对前几个子步进行计算

3、 检查渗透量和每一子步中的平衡迭代次数，如果总体收敛困难是由过大的渗透引起的(而不是由不平衡力和位移增量引起的)，那么可能低估了FKN的值或者是将FTOLN的值取得大小，如果总体的收敛困难是由于不平衡力和位移增量达到收敛值需要过多的迭代次数，而不是由于过大的渗透量，那么FKN的值可能被高估。

4、 按需要调查FKN或FTOLN的值，重新分析。

在有限元分析中，接触单元通常用来描述两物体相互接触或滑动的界面。近年来，ANSYS开发了一系列的接触单元。刚开始有节点对节点单元CONTAC12和CONTAC52，接着有节点对地单元CONTAC26，然后有节点对面单元CONTAC48和CONTAC49。最近几年，我们引入一类面对面接触单元CONTA169和CONTA174，同时还有一种新的节点对节点单元CONTA17