【讨论】面面接触实例对比分析（分析FKT对计算结果的影响及ANSYS对接触分析的过程）

lzc0

下面命令流是面面接触实例
/PREP7   
ET,1,185
MP,EX,1,2.1E5   
MP,PRXY,1,0.3   
CYL4,0,0,34,0,100,90,25  
CYL4,0,0,25,0,35,90,150  
VGEN,,2,,,,,-10,,,1  
LESIZE,17,,,15,,,,,1
LESIZE,19,,,15,,,,,1
LESIZE,18,,,2,,,,,1  
LESIZE,20,,,2,,,,,1  
LESIZE,22,,,20,,,,,1
LESIZE,5,,,10,,,,,1  
LESIZE,7,,,10,,,,,1  
LESIZE,6,,,8,,,,,1   
LESIZE,8,,,8,,,,,1   
LESIZE,10,,,3,,,,,1  
VSWEEP,ALL   
/COM,CONTACT PAIR CREATION-START
MP,MU,1,0.2  
MAT,1   
R,3  
REAL,3   
Et,2,170
ET,3,174
R,3,,,0.1,0.1,,  
NROPT,UNSYM  
ASEL,S,,,4   
CM,_TARGET,AREA  
TYPE,2   
NSLA,S,1
ESLN,S,0
ESURF,ALL   
ASEL,S,,,9   
CM,_CONTACT,AREA
TYPE,3   
NSLA,S,1
ESLN,S,0
ESURF,ALL   
CMDEL,_TARGET   
CMDEL,_CONTACT   
ALLSEL,ALL   
EPLOT   
FINISH   
/SOLU   
DA,5,SYMM   
DA,6,SYMM   
DA,11,SYMM   
DA,12,SYMM   
DA,3,ALL
ANTYPE,0
NLGEOM,1
AUTOTS,0
TIME,100
SOLVE   
esel,all
esel,s,ename,,174   
eplot   
/VIEW, 1 ,1,1,1  
/ANG, 1  
/REP,FAST   
FINISH   
/POST1   
AVPRIN,0,0,  
ETABLE, ,CONT,PENE   
!*   
AVPRIN,0,0,  
ETABLE, ,CONT,PRES   
!*   
AVPRIN,0,0,  
ETABLE, ,CONT,SFRI   
!*   
AVPRIN,0,0,  
ETABLE, ,CONT,SFRI   
!*   
PRETAB,CONTPENE,CONTPRES,CONTSFRI
  
小弟通过改变接触单元的材料,并把法向接触刚度设为上例默认的normal stiffness,发现结果CONTPENE,CONTPRES,CONTSFRI的值不发生改变。
MP,EX,2,1E3  
MP,PRXY,2,0.4   
MP,MU,2,0.2  
MAT,2   
R,3,,,-90586,0.1,,  
在此基础上，小弟又将切向接触刚度的系数分别调为2,3,5,10,发现CONTPENE,CONTPRES值没什么变化，
但CONTSFRI的值变化较大。
小弟验算了一下，CONTPRES值等于CONTPENE乘以法向接触刚度，但只有当切向接触刚度的系数为10时，CONTSFRI＝CONTPRES*摩擦系数。
最后，我在切向接触刚度为5的情况下，把材料2的EX改为1E5，发现结果跟EX为1E5时一样。
请问通过以上实例是否可说明，ANSYS中，接触单元刚度矩阵的形成跟其本身材料无关，只受法向接触刚度，切向接触刚度的系数影响。
  
但我实在不明白CONTSFRI是如何计算的，为何出现上述的结果呢！请各位接触大虾们多多指导。
  

不够再添

七年过去了，无人回答。

taogreat

八年无人管

五磅兔

岁月苍桑。留念。

sunxueshuai

话题很好，竟无人问津！缅怀...

低调的学习

九年神贴

yf06130923

坐等十年神贴！

valen_zt

九年一帖。。。

后天的事情

此贴马上十年

shuiruohan9

这是2014年的第一帖，整整十年了，这帖子都沉不了，命够硬

biyao222

默然发现时04年的贴啊

870813752

10年神贴无人解决。。。

koostreet

楼主自己都忘了吧

包大兴

神贴不能沉

dragon_w

面面接触时算的是应力，所以法向应力CONTPRES等于法向穿透值CONTPENE乘以法向刚度FKN

dragon_w

ANSYS通过两种算法来实现此法向接触关系：罚函数法和拉格朗日乘子法。
1．罚函数法
是通过接触刚度在接触力与接触面间的穿透值（接触位移）间建立力与位移的线性关系：
    接触刚度*接触位移=法向接触力
对面面接触单元17*，接触刚度由实常数FKN来定义。 穿透值在程序中通过分离的接触体上节点间的距离来计算。接触刚度越大，则穿透就越小，理论上在接触刚度为无穷大时，可以实现完全的接触状态，使穿透值等于零。但是显而易见，在程序计算中，接触刚度不可能为无穷大（否则病态），穿透也就不可能真实达到零，而只能是个接近于零的有限值。
以上力与位移的接触关系可以很容易地合并入整个结构的平衡方程组K*X=F中去。并不改变总刚K的大小。

dragon_w

http://blog.sina.com.cn/s/blog_7d9d891c0102uw2v.html这个博客上看的

qmchhz

十年磨一剑  继续修炼:@:@;P;P

闻聪聪

十三年的神贴，厉害了

meiqiuchixiaoji

十五年陈酿