求助！同样的模型用usermat 刚体位移，用标准库材料就没问题

elevenhours

usermat用ansys自带的例子。按说是和 tb,biso 一个效果，可为甚么用到接触分析当中同样的设置tb,biso可以计算，但usermat就算不出来接触？如图真相两图中模型一样，接触面在正中y=0面，两部分左边都约束ux，下半部分底部约束uy，上半部分上边受压（uy=-0.04）

              图1 使用usermat计算结果的位移云图
图中红色位移为0，蓝色位移为-0.04，可以看到两部分根本没有接触，上半部分整体位移了-0.04与下半部分重叠了一些。

              图2 使用tb,biso计算结果的位移云图
可以看到，由于上下两部分在接触面上相接，整个模型呈典型的受压状态。底边为0，上边-0.04。
怎么会这样？难道接触不能用usermat定义基体？
附上命令吧

---

finish
/clear
/config,NRES,5000
/PREP7

/PNUM,KP,1  
/PNUM,LINE,1
/PNUM,AREA,1
/pnum,node,0
/pnum,elem,0
/DSCALE,1,1.0
et,1,183
KEYOPT,1,3,2
et,2,172
KEYOPT,2,2,1
keyopt,2,12,4
et,3,169

*set,INnum,50

*set,Smax,600
*set,DeltaN,1e-3
*set,DeltaT,1e-3

mp,ex,1,2.1e5
mp,nuxy,1,0.3
tb,biso,1,1,4
tbtemp,1.0
tbdata,1,1e3,100
!mp,mu,1,0
!tb,user,1,,4
!tbdata,1,2.1e5, 0.3, 1000,100,       ! E, posn, sigy, H

TB,CZM,2,1,,CBDD       
TBDATA,1,Smax,DeltaN,,,1/INnum/20
mp,mu,2,0.1

*set,a,10
*set,b,5

k,,0,0
k,,a,0
k,,a,b
k,,0,b

l,1,2
l,2,3
l,3,4
l,4,1

al,all
lesize,all,1

mshape,0,2d
mshkey,1
type,1
mat,1

allsel
amap,1,1,2,3,4

lsel,s,,,1
nsll,s,1
esln,all
EREFINE, all,,,1,Clean, off
numcmp,all

allsel
numcmp,all
arsym,y,all
!nummrg,all,,,,low

!!!!!!!!!!!!!!!!!!!!!!
lsel,s,,,1
nsll,s,1
esln,s
mat,2
type,2
esurf,all

lsel,s,,,5
nsll,s,1
esln,s
mat,2
type,3
esurf,all

/SOL
antype,0
NLGEOM,1
NSUBST,INnum
outres,erase
outres,all,1
autots,0
time,1
nsel,s,loc,y,-b
d,all,uy
allsel
nsel,r,loc,x,0
d,all,ux
allsel

nsel,s,loc,y,b
!nsel,r,loc,x,0
d,all,uy,-0.04
allsel
solve

/POST1  
SET,LAST
/EFACET,1   
PLNSOL, u,Y, 0,1.0

dragonsea

楼主解决了吗，我也遇到类似的问题，后改为NLGEOM,1，就没有了，但出现单元One or more elements have become highly distorted.  Excessive         
distortion of elements is usually a symptom indicating the need for     
corrective action elsewhere. 的问题