马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
本帖最后由 houdezhiprince 于 2011-5-9 17:34 编辑
请教各位高手,小弟目前在用ANSYS做有关转子动力学的题目,小弟建立了一个梁的模型(BEAM188),然后加约束,得到CAMPBELL图如附件所示,但执行PRCAMP,1,1,RPM,却没有任何显示,小弟现在想要得到临界转速,恳请各路高手给予小弟指点一二,小弟不胜感激。APDL代码如下
/CLEAR
/BATCH,list
/PREP7
/UNITS,SI ! Define units used: meter
*DIM,SPIN,,7 ! SPIN VELOCITY (RPM)
SPIN(1) = 0.
SPIN(2) = 105000.
SPIN(3) = 70000.
SPIN(4) = 80000.
SPIN(5) = 90000.
SPIN(6) = 100000.
SPIN(7) = 105000.
INNER_RAD = 0.05
DIAM = 2*INNER_RAD
L = 0.9
MP,EX,1,2e+11 ! Young's modulus
MP,DENS,1,7806 ! Density
MP,NUXY,1,0.3 ! Possion ratio
MP,GXY,1,7.6923e+010 ! Shear modulus
!DEFINE PARAMETERS OF DISK
OUTER_RAD = 0.24
DISK_THICK = 0.05
DENS = 7806
PI = 4*ATAN(1.0)
DISK_MASS = DENS*PI*(OUTER_RAD**2-INNER_RAD**2)*DISK_THICK
IX = 0.5*DISK_MASS*(OUTER_RAD**2+INNER_RAD**2)
IY = 0.25*DISK_MASS*(OUTER_RAD**2+INNER_RAD**2)+(DISK_THICK**2)/12
IZ = IY
ET,1,BEAM188,,,2
SECTYPE,1,BEAM,CSOLID
SECDATA,INNER_RAD,30
ET,2,MASS21
R,20,DISK_MASS,DISK_MASS,DISK_MASS,IX,IY,IZ
ET,3,COMBI214
KEYOPT,3,2,1 ! YZ PLANE
KEYOPT,3,3,0 ! ELEMENT IS SYMMETRIC
R,30,3E+11,3E+11,5E+10,5E+10
!CREATE NODES
N,1,0,0
N,31,L,0
FILL,1,31
!NODES FOR BEARING
N,32,0,0.2
N,33,L,0.2
TYPE,1
MAT,1
*DO,I,1,30
E,I,I+1
*END DO
TYPE,2
REAL,20
E,10
TYPE,3
REAL,30
E,1,32
E,31,33
!TYPE,3
!REAL,30
FINI
/SOLU
D,ALL,UX
D,ALL,ROTX
D,32,ALL
D,33,ALL
RATIO = PI/30
ANTYPE,MODAL
CORIOLIS,ON,,,ON ! CORIOLIS ON IN A STATIONARY REFERENCE FRAME
NBF = 20
MODOPT,QRDAMP,NBF,,,ON
/OUT,SCRATCH
OMEGA,SPIN(1)*RATIO
MXPAND,NBF
PSTRES,0
SOLVE
OMEGA,SPIN(2)*RATIO
MXPAND,NBF
PSTRES,0
SOLVE
FINI
/POST1
PLCAMP,1,1,RPM
PRCAMP,1,1,RPM
[img]C:\Users\WangDG\Desktop\4NCQ@PCY@]QL[OI_5TI[QDR.jpg[/img]
|