此代码可能有问题，仅供学习之用。

*SET,RADIUS_INNER , 10
*SET,RADIUS_OUTER , 30
*SET,LOAD_OMEGA , 100
/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
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
*SET,NEXTNODE,0
*DO,ICOUNT,1,NCOUNT-1
*SET,NEXTNODE,NDNEXT(NEXTNODE)
E,NODE(0,0,0),NEXTNODE
*ENDDO
CSYS,0
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
*SET,NEXTNODE,0
*DO,ICOUNT,1,NCOUNT
*SET,NEXTNODE,NDNEXT(NEXTNODE)
IC,NEXTNODE,UY,,LOAD_OMEGA*NX(NEXTNODE)
*ENDDO
CSYS,0
ALLSEL,ALL
SOLVE
FINISH
/POST1
SET,LAST
RSYS,1
ANDATA,0.1,,2,0,0,1,0,1