请大家帮忙看看我的命令流

forseti

上面的体向下移动,与下面的体接触,产生应力和变形.我的想法是由于初始位置有一段距离,因此先加了一个载荷步使上面的体在不建立接触对的情况下向下移动,到了将要发生接触的位置再建立接触对,重新给位移载荷.

现在的情况是第一个载荷步可以运动到理想位置,加第二个载荷的时候已经移动到位的模块突然又回到上面去了,请大家帮我检查一下命令流,看看问题出在哪里吧,有这方面经验的朋友请不吝赐教,谢谢.命令流如下:

/PREP7
BLOCK,0,10,0,1,0,-10,   
BLOCK,4,6,6,7,-4,-6,
et,1,solid45
et,2,solid95
mptemp,1,293,473,873,1273,1473,1773
mpdata,ex,1,1,2.1e11,1.95e11,1.55e11,0.45e11,0.01e11,0.0005e11  
mpdate,nuxy,1,1,0.3,0.3,0.3,0.3,0.3,0.3
mpdate,alpx,1,1,1.23e-5,1.23e-5,1.23e-5,1.23e-5,1.23e-5,1.23e-5
tb,biso,1,6,2   
tbtemp,293,1
tbdata,1,530e6,2.1e10   
tbtemp,473,2
tbdata,1,258e6,1.95e10  
tbtemp,873,3
tbdata,1,137e6,1.55e10  
tbtemp,1273,4   
tbdata,1,30e6,0.45e10   
tbtemp,1473,5   
tbdata,1,10e6,0.01e10   
tbtemp,1773,6   
tbdata,1,1e6,0.0005e10  
mp,dens,1,7815  
mp,kxx,1,43
mp,c,1,580  
mp,mu,1,0.3
MP,MU,1,0.3
MAT,1   
R,3
REAL,3  
et,3,targe170   
et,4,conta174   
KEYOPT,4,9,0
KEYOPT,4,10,1   

lesize,17,,,4   
lesize,18,,,2   
lesize,19,,,4   
lesize,20,,,2   
lesize,22,,,4   
mat,1   
type,1  
vsweep,2,7,8,1  

esize,5,,,10
lesize,6,,,2
lesize,7,,,10   
lesize,8,,,2
lesize,11,,,10  
mat,1   
type,1  
vsweep,1,1,2,1  

ASEL,S,,,9          !         9
NSLA,s,1            !   
cm,mubiao1,node
asel,s,,,4  
NSLA,s,1
cm,jiechu1,node

/SOL
cmsel,s,mubiao1
d,all,ux,0  
d,all,uy,-4.97

asel,s,,,1       !      area1   
nsla,s,1         !  
d,all,uz,0       !z                  0  
d,all,ux,0  
d,all,uy,0  
asel,s,,,2  
nsla,s,1
d,all,uz,0  
d,all,ux,0  
d,all,uy,0  
asel,s,,,5  
nsla,s,1
d,all,uz,0  
d,all,ux,0  
d,all,uy,0  
asel,s,,,6  
nsla,s,1
d,all,uz,0  
d,all,ux,0  
d,all,uy,0  
allsel,all                 !
time,2                   !               2,      20
nsubst,10   
autots,on                  !                  /
nropt,full,,off            !                                                                       /
nropt,unsym,,off
neqit                      !                                                 25?
lnsrch,auto                !
nlgeom,on                  !
outres,all,all
lswrite,1

/prep7
ALLSEL,ALL  
cmsel,s,mubiao1
type,3
mat,1
real,3             !单元类型3,target170,实常数组3,材料属性1
esurf              !在选择的面上生成节点
ALLSEL,ALL         
cmsel,s,jiechu1    !接触面
type,4             !target173
mat,1
real,3
esurf
R,3,0,0,0.01       !组3的实常数
KEYOPT,4,5,1       !定义单元类型为4的key为9,10,值为1
KEYOPT,4,10,1
allsel,all

/SOL
cmsel,s,mubiao1
d,all,ux,0  
d,all,uy,-0.53

asel,s,,,1       !      area1   
nsla,s,1         !  
d,all,uz,0       !z                  0  
d,all,ux,0  
d,all,uy,0  
asel,s,,,2  
nsla,s,1
d,all,uz,0  
d,all,ux,0  
d,all,uy,0  
asel,s,,,5  
nsla,s,1
d,all,uz,0  
d,all,ux,0  
d,all,uy,0  
asel,s,,,6  
nsla,s,1
d,all,uz,0  
d,all,ux,0  
d,all,uy,0  
allsel,all                 !
time,2                   !               2,      20
nsubst,10   
autots,on                  !                  /
nropt,full,,off            !                                                                       /
nropt,unsym,,off
neqit                      !                                                 25?
lnsrch,auto                !
nlgeom,on                  !
outres,all,all  
lswrite,2

LSSOLVE,1,2,1

sysh320

接触问题不是很懂，希望大家多多讨论
像这种有实例的帖子最好了，鼓励鼓励

ymcheng123

不懂，也在搞桩土接触