关于HYPERMESH中做接触需要注意的一些问题

研究了两天，终于大概搞出个所以然来。 在Hypermesh中做接触分析，要区分两种状态：

1 接触体（GA）与目标体(GB)的距离如果小于1e-14，那就需要定义一个局

部坐标系，坐标系的X轴是GA->GB的方向。如果是圆周接触，那就需要定义一个柱坐标系，X轴是径向方向。

1.1 在

gap单元的property（PGAP）中，U0=AUTO，KA=auto，KB=blank，MU1=blank,MU2=blank。

control card的GAPPRM中，CKGPDIR=REV(因为距离非常近，有可能GA已经渗透入GB，此时X方向就不能代表GA->GB。REV选项可以调整这种反向的状态)，GAPGPRJ=NORM，ERRMSG=FULL（输入所有错误信息，可以看出哪里出了问题）。

1.2 在

2 GA与GB距离大于1e-14，无需定义局部坐标系。

2.1 在

gap单元的property（PGAP）中，U0=AUTO，KA=auto，KB=10（一个小一些的值），MU1=blank,MU2=blank。或U0=0，KA=auto,KB=blank。 control card的GAPPRM中，CKGPDIR=ERR，GAPGPRJ=NORM，ERRMSG=FULL

2.2 在

对于PGAP一些参数的解释： U0 F0 初始间隙。当GA与GB的距离大于1e-14，但又不算太大时，可以将此值设为0，则KB可以设为blank。 预紧力（不太清楚怎么使用） 接触对闭合时，接触刚度大小，一般设为AUTO。KA是为了防止接触相互渗透而设计的，理论上，大一些的KA会使接触条件更加准确，但是过大的KA，会造成收敛性的问题或刚度矩阵变坏（对于KT尤为如此）。如果确实发生了这种情况，那就要适当减小KA的值，比较合理的KA值是（1e3~1e6）*E*h。 E是材料的弹性模量，h是接触单元周围单元的大小。 接触对未闭合时，接触刚度的大小。默认是KA*1e-14。 横向接触刚度（相当于摩擦力的大小）。默认是MU1*KA。一般auto即可。 静摩擦系数。 动摩擦系数。 KA KB KT MU1 MU2 如果要设置非滑移表面，则可以在control card中，设置PARAM,GAPOFFSET = NO。