圆环 接触问题 算不通？？？
圆环 接触问题 算不通？？？模型就是一个圆环，圆环里面一个圆，想通过平移里面圆的位移来求解外面圆环 接触后的变形量。
也就说里面的圆 只能给圆环接触是施加压力，而不施加拉力。。
但是设置计算后，始终都不是接触计算的结构，也就是所其实圆环无变形？？？
请问问题出在哪了？？？
命令如下：
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 编辑 ]
页:
[1]