考虑轴力的弹性支撑变截面梁的模态分析问题

trueghost

最近用ansys做了一个变截面梁的模态分析。
在加轴力之前分析结果和理论结果符合，
但加了轴力之后，总感觉不大对劲。轴力在0到2000N之间时，频率与未加轴力一样， 但轴力位于2000到2650N之间时，频率明显变小，直至屈曲。下面是我的命令流，请各位指教。
finish
/CLEAR,START
pd=0
/prep7
ANTYPE, STATIC
PSTRES, ON
*do,i,1,20,1
et,i,beam54okok.org !定义变截面桥墩单元
mp,dens,i,2700
mp,ex,i,7e10
mp,prxy,i,0.33
mp,gxy,i,2.63e10
R,i,2.755e-4*(1-(i-1)*pd/20),2.072e-9*(1-(i-1)*pd/20)**3,0.00475*(1-(i-1)*pd/20),0.00475*(1-(i-1)*pd/20),
RMODIF,i,9, , ,okok.org
RMODIF,i,14, ,
RMODIF,i,5,2.755e-4*(1-i*pd/20),2.072e-9*(1-i*pd/20)**3,0.00475*(1-i*pd/20),0.00475*(1-i*pd/20),
RMODIF,i,11, , ,okok.org
RMODIF,i,15, ,
RMODIF,i,13, ,
RMODIF,i,16, , , ,
*enddo
et,21,combin14
keyopt,21,2,1okok.org !纵向弹簧单元
et,22,combin14
keyopt,22,2,6okok.org !扭转弹簧单元
okok.org
r,21,8.2e3
r,22,8.2e3
r,23,5.63e2
okok.org
n,1,0,0,0
n,21,0,0.262,0
fill
n,22,0,0,0
n,23,0,0,0
n,24,0,0.262,0
okok.org
*do,i,1,20,1
type,i
real,i
e,i,i+1
*enddo
okok.org
type,21
real,21
e,1,22
okok.org
type,21
real,22
e,21,24
okok.org
type,22
real,23
e,1,23
okok.org
f,21,fy,-2600
d,22,allokok.org
d,23,all,all
d,24,allokok.org
d,1,uyokok.org
finish
okok.org
!（4）进入求解模块, 指定求解选项并求解
/SOLU
OUTPR, BASIC, 1
SOLVE
FINISH
! （5）进入一般后处理模块, 提取频率计算结果
/POST1
ETABLE, STRS, LS, 1
*GET, STRSS, ELEM, 20, ETAB, STRS
*GET, DEF, NODE, 21, U, y
! （6）重新进入求解模块, 用模态求解方法求解
/SOLU
ANTYPE, MODAL
PSTRES, ON
MXPAND, 3
MODOPT, LANB, 3
SOLVE
! （7）提取频率计算结果
*GET, FREQ1, MODE, 1, FREQ
*GET, FREQ2, MODE, 2, FREQ
*GET, FREQ3, MODE, 3, FREQ
*STATUS, PARM
FINISH