圆环 接触问题 算不通？？？

zhangload

圆环 接触问题 算不通？？？

模型就是一个圆环，圆环里面一个圆，想通过平移里面圆的位移来求解外面圆环 接触后的变形量。

也就说里面的圆 只能给圆环接触是施加压力，而不施加拉力。。


但是设置计算后，始终都不是接触计算的结构，也就是所其实圆环无变形？？？

请问问题出在哪了？？？

命令如下：
fini
/clear
/filname,61Con

hs=0.3e-3
di=20e-3
dt=69e-3
le=4e-3
ex=2e11
dens=7800
nuxy=0.3
force=4
fp=1e8

/prep7

cyl4,,,le/2
cyl4,,,le/2,,di/2
cyl4,,,di/2,,dt/2
aglue,1,2
numcmp,area


et,1,63
et,2,shell63
ET,3,TARGE169          ! 2-D TARGET ELEMENTS
ET,4,CONTA171  
mp,ex,1,ex
mp,nuxy,1,nuxy
mp,dens,1,dens
r,1,hs
mp,ex,2,ex*1e5
mp,nuxy,2,nuxy
mp,dens,2,dens
r,2,hs
!asel,s,loc,z,0
aatt,2,2,2
esize,le
!eshape,0,2d

TSHAPE,LINE  !-------------------------

csys,1
lsel,s,loc,x,le/2
lesize,all,,,2
MSHAPE,0,2d
MSHKEY,0
amesh,1
lsel,s,loc,x,di/2
lesize,all,,,6
MSHAPE,0,2d
MSHKEY,2
amesh,3
esize,le
aatt,1,1,1
amesh,2

asel,s,,,2
lsla,s
LSEL,r,loc,x,di/2   
!NSLL,,1                ! SELECT TARGET NODES ON BIGGER CYLINDER
!cm,con,node
TYPE,4
real,3
!lmesh,all
nsll,s
ESURF  !,,REVERSE


csys,1
asel,s,,,3
lsla,s
LSEL,r,loc,x,di/2  
!lsel,s,,,31,34,1
!NSLL,s,1                ! SELECT TARGET NODES ON BIGGER CYLINDER
!nsel,r,loc,x,di/2
!cm,tag,node
TYPE,3
real,3
lmesh,all
nsll,s
ESURF,,REVERSE


allsel

fini
/solu
antype,STATIC
!outres,basic,all
!outpr,basic,all
pred,on
NLGEOM,ON                     ! INCLUDE LARGE DEFORMATION EFFECTS
NSUBST,6                   ! SPECIFY NUMBER OF SUBSTEPS IN LOAD STEP
!lnsr,on
allsel
csys,1
lsel,s,loc,x,dt/2
nsll,s,1
d,all,all,0
allsel

csys,1
nsel,s,loc,x,0,le/6
cm,nc,node

csys,1
asel,s,,,1,3,2
!lsla,s
!LSEL,r,loc,x,di/2  
nsla,s
!cp,1,ux,all
!cp,1,uy,all

!d,nc,z,0,,,,rotz,rotx,roty
allsel
!*CREATE,SOLV2D,MAC     ! CREATE SOLUTION MACRO FOR 2-D CASE
/SOLU
csys,0
D,nc,ux,100e-6
solve



fini
/post1
pldisp
allsel
csys,1
nsel,s,loc,x,0,le/6
prnsol,u

allsel
plNsol,s,EQV
prrsol

[ 本帖最后由 zhangload 于 2007-1-18 11:37 编辑 ]