接触疑问

wen19861123tuo

发程序吧，图片不会弄，见谅。
我的疑问是为何求解总是出现刚体位移呢？我尝试将加压力的那个面换成加位移，也就是向下的位移 -0.1，就可以求解，说明接触设定是没有问题，但是换成压力就不好用，哪位大虾帮忙解决一下，万分感激！
finish
/clear
r1=10
r2=13
/prep7
cylind,r1,,0,1,-180
wpstyl,defa
wpoffs,,-23
cylind,r2,,0,1,0,180
wpstyl,defa
et,1,45
MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,EX,1,,2e5
MPDATA,PRXY,1,,0.3  
esize,0.5
vsweep,all
et,2,targe170            !3D接触目标面单元170
et,3,conta173            !3D接触接触面单元174
mp,mu,3,0.3               !设定摩擦系数
real,3
mat,3
keyopt,3,9,0
keyopt,3,12,0
keyopt,3,10,2
asel,s,,,8
nsla,s,1
esln,s,0
type,2
esurf,all
asel,s,,,3
nsla,s,1
esln,s,0
type,3
esurf,all
allsel
!****约束
asel,s,loc,y,-23
nsla,s,1
d,all,all
asel,s,loc,z,0
!asel,a,loc,z,1
nsla,s,1
d,all,uz,0
!da,all,ux,0
asel,s,loc,y,0
nsla,s,1
sf,all,pres,10
allsel
/sol
allsel
nsubst,100,1000,1
outres,all,last
autots,1
nropt,full,,off
lnsrch,on  
pred,on  
nlgeom,off      
kbc,0
antype,0   
time,1
solve

xqm820

本人遇到同样问题，有人建议先施加小位移，再施加面载荷。我不知道对不对，交流啊。我qq529034882

mxjdream8

我也遇到类似的问题。现在的方法是人为的移动模型，是两个接触面靠近，然后再施加压力等载荷。但是发现，只要是距离不足够近就出现刚性位移等报错。还有就是发现两个接触面之间的距离对求解得出的压应力结果影响很大。也很是困惑。希望能多交流。

tyrone

好像一般是通过调整初始接触选项可解决