异形轴频率分析，大家帮忙看看，有什么不对的地方

sy2690

/clear
/filname,zhou
/prep7
et,1,solid45
mp,ex,1,2e11
mp,prxy,1,0.3
mp,dens,1,7850
mp,ex,2,2e11
mp,prxy,2,0.3
mp,dens,2,12560
                     !第一轴段
!建立关键点，连线
k,1,0,0,0
k,2,0,0.14,0
k,3,0.47,0.14,0
k,4,0.47,0.1575,0
k,5,0.64,0.1575,0
k,6,0.64,0.14,0
k,7,0.955,0.14,0
k,8,0.955,0.1575,0
k,9,1.086,0.1575,0
k,10,1.086,0.1775,0
k,11,1.203,0.1775,0
k,12,1.203,0.18,0
k,13,1.233,0.18,0
k,14,1.233,0.1825,0
k,15,1.35,0.1825,0
k,16,1.35,0.185,0
k,17,1.51,0.185,0
k,18,1.51,0.1875,0
k,19,1.84,0.1875,0
k,20,1.84,0,0
*do,i,1,19
lstr,i,i+1
*enddo
lstr,20,1

                      !第二轴段
k,21,3.4,0,0
k,22,3.4,0.18,0
k,23,3.73,0.18,0
k,24,3.73,0.185,0
k,25,3.89,0.185,0
k,26,3.89,0.1825,0
k,27,4.007,0.1825,0
k,28,4.007,0.18,0
k,29,4.037,0.18,0
k,30,4.037,0.1775,0
k,31,4.154,0.1775,0
k,32,4.154,0.1575,0
k,33,4.285,0.1575,0
k,34,4.285,0.14,0
k,35,4.595,0.14,0
k,36,4.595,0.1575,0
k,37,4.762,0.1575,0
k,38,4.762,0.105,0
k,39,4.91,0.105,0
k,40,4.91,0.1,0
k,41,4.98,0.1,0
k,42,4.98,0.0995,0
k,43,5.14,0.0995,0
k,44,5.14,0.0945,0
k,45,5.15,0.0945,0
k,46,5.15,0.099,0
k,47,5.16,0.099,0
k,48,5.16,0,0
*do,i,21,47
lstr,i,i+1
*enddo
lstr,48,21
!生成面
lsel,s,,,1,20,1
al,all
lsel,s,,,21,48,1
al,all
!旋转成体
vrotat,1,2,,,,,1,48,360
vadd,all
                     !创建块体
block,1.84,3.4,-0.241,0.241,0.241,-0.241
block,1.955,3.285,-0.446,0.446,0.16,-0.16
block,1.955,3.285,-0.16,0.16,0.446,-0.446
vadd,1,2,3
numcmp,all
!粘接体
vglue,1,2,3
vatt,1,1,1
mshape,1,3d
mshkey,0
esize,0.05
vmesh,2
vatt,1,1,1
vmesh,3
vatt,2,1,1
vmesh,4

!旋转坐标系
wprota,,,90
!创建柱坐标系，并激活
cswpla,11,1,1,1
!将所有节点旋转到激活坐标系
nrotat,all
fini
/solu
!指定分析类型为静态分析
antype,0
pstres,on!打开预应力选项
!添加约束
dk,1,uy,,,,uz
dk,48,uy,,,,uz
allsel,all!选择所有元素
omega,200*atan(1),0,0,0!定义绕轴旋转角速度
save
solve
fini
/solu
antype,2!模态分析
modopt,lanb,5
mxpand,5
pstres,1
save
solve
fini
/post1
set,list

沙隆巴斯

起码得给大家看看结果，然后再说有什么不对的地方。咱都不是计算机，直接看程序看不出来结果。