ansys圆盘转动例子

2017-03-02 by:CAE仿真在线来源:互联网

/clear

RADIUS_INNER = 10
RADIUS_OUTER = 30
LOAD_OMEGA = 100

/auto
/triad,lbot
/pnum,type,1
/num,1

/prep7
et,1,63
r,1,0.1
mp,ex ,1,10e6
mp,nuxy,1,0.3
mp,dens,1,0.1/386.1

et,2,184
keyopt,2,1,1

pcirc,RADIUS_OUTER,RADIUS_INNER,0,90
pcirc,RADIUS_OUTER,RADIUS_INNER,90,180
pcirc,RADIUS_OUTER,RADIUS_INNER,180,270
pcirc,RADIUS_OUTER,RADIUS_INNER,270,360
nummrg,kp

aatt,1,1,1
mshape,0
mshkey,1
esize,5
amesh,all

! create rigid beams inside
csys,1
nrota,all
lsel,s,loc,x,RADIUS_INNER
nsll,s,1
n,,0,0,0
type,2
real,2
mat,2
*get,NCOUNT,node,,count
NEXTNODE=0
*do,ICOUNT,1,NCOUNT-1
NEXTNODE=ndnext(NEXTNODE)
e,node(0,0,0),NEXTNODE
*enddo
csys,0

! create dummy element
nsel,all
n,,RADIUS_OUTER+5,0,0
e,node(RADIUS_OUTER,0,0),node(RADIUS_OUTER+5,0,0)

d,node(0,0,0),all
ddele,node(0,0,0),rotz

allsel,all
finish

/solu
antype,trans
nlgeom,on
rescon,define,none
outres,all,all
time,0.04
deltim,6e-4,1e-5,0.1

nsel,u,node,,node(0,0,0)
csys,1
*get,NCOUNT,node,,count
NEXTNODE=0
*do,ICOUNT,1,NCOUNT
NEXTNODE=ndnext(NEXTNODE)
ic,NEXTNODE,uy,,LOAD_OMEGA*nx(NEXTNODE)
*enddo
csys,0

allsel,all
solve
finish

/post1
set,last
rsys,1
plns,u,sum
ANDATA,0.1,,2,0,0,1,0,1

开放分享：优质有限元技术文章,助你自学成才