3维实体单元求解临界转速问题

514981310

各位大侠好，我前阵试着用梁单元188模拟转子和214单元模拟轴承计算转子临界转速，现在想用3维实体单元的方法计算临界转速，这次用的185单元，试着按思路编写了命令流，但是求解时候直接弹出，不能计算，若有高手可否指点一二，先谢过各位了，下面附上命令流。
/FILNAME,Critical Speed,0
/PREP7
L1=240e-3   !转轴的长度
D1=130e-3  !转轴的直径
RD1=130E-3   !转子内径
RD2=144E-3   !转子外径
KYY1=1.7E9
KZZ1=1.7E9
ET,1,SOLID185   !定义1号单元
ET,2,COMBI214 !定义2号单元
KEYOPT,2,2,1  !设置3号单元关键字，轴承平行于Y-Z平面
R,1,KYY1,KZZ1,!提取节点号
MPTEMP,,,,,,,,    !赋材料
MPTEMP,1,0  
MPDATA,EX,1,,2.1E11
MPDATA,PRXY,1,,0.30
MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,DENS,1,,7800
MPDATA,DAMP,1,,2E-5
*afun,deg
WPROTA,,,90
PCIRC,0,130E-3,0,360     !建立转轴模型
PCIRC,0,144E-3,0,360  
AOVLAP,ALL               !建立转子轭模型
NUMCMP,ALL                !压缩编号  
WPROTA,,,-90
VEXT,1,,,280E-3
AGEN,2,2,,,20e-3,,,,
ADELE,2,,,1
VEXT,8,,,140E-3
LSEL,S,,,9,12,1
ESIZE,,30
ALLSEL,ALL
LSEL,S,,,15
ESIZE,,60
VSEL,S,,,1,2
VATT,1,,1
MSHAPE,0,3D          ！剖分
VSWEEP,1
VSWEEP,2
LOCAL,11,1,0,0,0,0,0,90
CSYS,11
*DIM,ND,ARRAY,120,1     ！转子上面连接轴承节点数组
*DIM,ND1,ARRAY,120,1    ！转子上面连接轴承节点数组
*DIM,ANG,ARRAY,120,1    ！角度数组
*DIM,DY10,ARRAY,120,1   ！轴承数组，一共两个轴承
*DIM,DY170,ARRAY,120,1
*DO,I,0,119,1
ANG(I+1)=I*360/120           ！角度数组赋值
*ENDDO
*DO,I,1,120,1
ND(I)=NODE(130E-3,ANG(I),(10E-3))        ！转子上面连接轴承节点数组赋值，共两个
ND1(I)=NODE(130E-3,ANG(I),(170E-3))
*ENDDO
NSEL,S,,,ND(1)
*DO,I,2,120,1
NSEL,A,,,ND(I)
*ENDDO
NPLOT
NSEL,S,,,ND1(1)
*DO,I,2,120,1
NSEL,A,,,ND1(I)
*ENDDO
NPLOT
ALLSEL,ALL
*DO,I,1,120,1
N,,D1*1.1,ANG(I),10E-3                   ！在转子外围建立一圈214单元用来模拟轴承，共两个
N,,D1*1.1,ANG(I),170E-3
*ENDDO
*DO,I,1,120,1
DY10(I)=NODE(D1*1.1,ANG(I),10E-3)         ！轴承数组赋值
DY170(I)=NODE(D1*1.1,ANG(I),170E-3)
*ENDDO
*DO,I,1,120,1
TYPE,2   
REAL,1
E,DY10(I),ND(I)                          ！将214轴承单元与转子表面连接，建立单元，
E,DY170(I),ND1(I)
*ENDDO
ALLSEL,ALL
/SOLU               
antype,modal !设置分析类型为模态分析
MODOPT,QRDAMP,20,,,ON !提取20阶
MXPAND,20,            !扩展20阶
!定义位移约束
NSEL,S,,,DY10(1)
NSEL,A,,,DY170(1)
*DO,I,2,120,1
NSEL,A,,,DY10(I)
NSEL,A,,,DY170(I)
*ENDDO
D,ALL,ALL
!考虑陀罗效应
CORIOLIS,1, , ,1,0
!施加转速开始求解
CSYS,0
ALLSEL,ALL
VSEL,S,,,1,2,1
ESLV,S
CM,ROTOR,ELEM
ALLSEL,ALL
CMOMEGA,ROTOR,0
solve

CMOMEGA,ROTOR,1000*2*3.1415/60
solve
CMOMEGA,ROTOR,2000*2*3.1415/60
solve
CMOMEGA,ROTOR,3000*2*3.1415/60
solve
CMOMEGA,ROTOR,4000*2*3.1415/60
solve
CMOMEGA,ROTOR,5000*2*3.1415/60
solve
CMOMEGA,ROTOR,6000*2*3.1415/60
solve
CMOMEGA,ROTOR,7000*2*3.1415/60
solve
CMOMEGA,ROTOR,8000*2*3.1415/60
solve
CMOMEGA,ROTOR,9000*2*3.1415/60
solve
CMOMEGA,ROTOR,10000*2*3.1415/60
solve
CMOMEGA,ROTOR,11000*2*3.1415/60
solve
CMOMEGA,ROTOR,12000*2*3.1415/60
solve
CMOMEGA,ROTOR,13000*2*3.1415/60
solve

/POST1               
PRCAMP,,1,RPM !画坎贝尔图
PLCAMP,,1,RPM
!提取临界转速
*GET,CRIC1,CAMP,1,VCRI,,,
*GET,CRIC2,CAMP,2,VCRI,,,
*GET,CRIC3,CAMP,3,VCRI,,,
*GET,CRIC4,CAMP,4,VCRI,,,
*GET,CRIC5,CAMP,5,VCRI,,,

liuxiaotian

直接弹出？ 毫无征兆 毫无报错 毫无提示？   仔细看看  列出来