ansys中轴承刚度随转速变化定义与求解方法

我用214单元模拟轴承求解转子系统的临界转速，把刚度设置为转速的函数，但是把命令流输入之后ANSYS到了求解部分就自动停止，不进行计算。下面是定义table数组和求解部分的命令流，请前辈指点下错误出在哪里。

另外，出问题后我查过file.err里面有一个错误提示大意是：omegas missing。

/prep7

l0=1.3

omega1=0

omega2=4000

omega3=8000

kxx1=3.2e6

kxx2=3e6

kxx3=3.5e6

kyy1=8e6

kyy2=8.2e6

kyy3=8.6e6

cx=2e-4

cy=1e-3

*dim,kxx,table,3,1,1,omegas

kxx(1,1)=kxx1,kxx2,kxx3

kxx(1,0)=omega1,omega2,omega3

*dim,kyy,table,3,1,1,zhuansu

kxx(1,1)=kyy1,kyy2,kyy3

kxx(1,0)=omega1,omega2,omega3

et,1,185,,2

et,2,214

keyopt,2,3,1

et,3,214

keyopt,3,3,1

et,4,21

r,1

r,2,%kxx%,%kxx%,,,cx,cx

r,3,%kyy%,%kyy%,,,cy,cy

/solu

nmod=10

antype,modal

modopt,qrdamp,nmod,,,on

mxpand,nmod,,,yes

coriolis,on,,,on

*do,i,1,3

omega,,,rotation(i,1)*2*acos(-1)/60

solve

*enddo

finish