/prep7
*dim,spin,,10 ! spin velocity (RPM)
*do,i,1,10
spin(i)=(i-1)*2000
球墨铸铁管qiumogg*enddo
RO = 7806 ! MATERIAL : STEEL
PEX = 2.078E+11
PGXY = 1.E+20 ! NO SHEAR
MP,EX,1,PEX
MP,PRXY,1,0.3
MP,DENS,1,RO
MP,GXY,1,PGXY
MP,BETD,1,2.0E-4
et,1,188
*do,i,1,21
rs232 ttl
n,i,20*(i-1)/1000
*enddo
n,22,410/1000
sectype,1,beam,csolid
secdata,100/2000
secnum,1
e,1:21,2:22
n,23,430/1000
n,24,450/1000
sectype,2,beam,csolid
secdata,170/2000
secnum,2
e,22:23,23:24
*do,i,1,20
n,i+24,450/1000+20*i/1000
*enddo
n,45,860/1000
secnum,1
e,24:44,25:45
*do,i,1,4
n,45+i,860/1000+20*i/1000
*enddo
secnum,2
e,45:48,46:49
*do,i,1,20
n,i+49,940/1000+20*i/1000
*enddo
n,70,1350/1000
secnum,1
e,49:69,50:70
n,71,1370/1000
n,72,1390/1000
secnum,2
e,70:71,71:72
自动打铃器*do,i,1,20
n,i+72,1390/1000+20*i/1000
*enddo
n,93,1800/1000
secnum,1
e,72:92,73:93
!
et,2,214,,1,0 !YZ plane,not symmetry
!kxx,kyy,kxy,kyx,cxx,cyy,cxy,cyx
!r,2,1e7,3e8,1e6,2e6,1e5,1e5,-1e4,-1e4
!r,2,1e6,1e6
n,1001,0,,200/1000
n,1002,1800/1000,200/1000
n,1003,1800/1000,,200/1000
ET,3,COMBIN14 ! ELEMENT TYPE #3 : BEARINGS
KEYOPT,3,2,2 ! Y DIRECTION
ET,4,COMBIN14
KEYOPT,4,2,3 ! Z DIRECTION
R,30,4E+7
type,3
real,30
e,1,1000甲烷制氢
e,93,1002
type,4
real,30
e,1,1001
e,93,1003
allsel
fini
/solu
d,all,ux
d,all,rotx
d,1000,all
d,1001,all
d,1002,all
d,1003,all
RATIO = 4*ATAN(1)/30
ANTYPE,MODAL
CORIOLIS,ON,,,ON,ON ! CORIOLIS ON IN A STATIONARY REFERENCE FRAME
NBF = 20
电价查询MODOPT,QRDAMP,NBF,,,ON QRDOPT,ON ! REUSE FLAG ON
/OUT,SCRATCH
*DO,I,1,10
OMEGA,SPIN(I)*RATIO
MXPAND,NBF
SOLVE
*ENDDO
FINI
/POST1
PRCAMP,,1.,RPM ! PRINT CAMPBELL VALUES FOR SLOPE=1, UNIT= RPM
PLCAMP,,1.,RPM